#include<bits/stdc++.h>
using namespace std;

// #define int long long

const int P=1e9+7;
const int mod = 1e9 + 7;

int qpow(int a,int b)
{
    int ret=1;
    while(b)
    {
        if(b&1)
            ret=1ll*ret*a%P;
        b>>=1;
        a=1ll*a*a%P;
    }
    return ret;
}

int s;

struct Matrix {
    int a[71][71];

    Matrix(int _t=0){
        for(int i=0;i<s;i++)
            for(int j=0;j<s;j++)
                a[i][j]=(i==j)*_t;
    }
    void Set(int _t = 0) {
        
        for(int i=0;i<s;i++)
            for(int j=0;j<s;j++)
                a[i][j]=(i==j)*_t;
    }
};

Matrix operator *(const Matrix &a,const Matrix &b)
{
    Matrix c;
    for(int i=0;i<s;i++)
        for(int k=0;k<s;k++)
            for(int j=0;j<s;j++)
                c.a[i][j]=(c.a[i][j]+1ll*a.a[i][k]*b.a[k][j])%P;
    return c;
}

Matrix qpow(Matrix a,int b)
{
    Matrix ret(1);
    while(b)
    {
        if(b&1)
            ret=ret*a;
        b>>=1;
        a=a*a;
    }
    return ret;
}

int n,m,T;

int ans[2000005 + 100];
int f[2000005 + 100];
int sol[2000005];

int val[75][75];
vector<int> bm(vector<int> s) {
    int n = s.size() , L = 0 , m = 0;
    vector<int> C(n) , B(n) , T ;
    C[0] = B[0] = 1;

    int b = 1;
    for(int i = 0;i < n;i++) {
        ++m ;
        int d = s[i] % mod;
        for(int j = 1;j <= L;j++) d = (d + 1LL * C[j] * s[i - j]) % mod;
        if(!d) continue ;
        T = C ;
        int coef = 1LL * d * qpow(b , mod - 2) % mod;
        for(int j = m ; j < n;j++) C[j] = (C[j] - 1LL * coef * B[j - m]%mod + mod) % mod;
        if(2 * L > i) continue ;
        L = i + 1 - L;
        B = T ;
        b = d ;
        m = 0;
    }
    C.resize(L + 1) ;
    C.erase(C.begin()) ;
    for(auto &x : C) x = (mod - x) % mod;
    return C ;
}
Matrix ta[150];
vector<int> cal[2000005 + 100];
signed main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cin>>n>>m>>T;
    s=n;
    Matrix A;
    for(int i=0;i<n;i++)
    {
        string s;
        cin>>s;
        // s=string(n,'1');
        vector<int> v;
        for(int j=0;j<n;j++)
        {
            if(s[j]=='0')
                continue;
            v.push_back(j);
        }
        int p=qpow(v.size(),P-2);
        for(auto j:v)
            A.a[j][i]=p;
    }
    ta[0].Set(1);
    ta[1] = A;
    for(int i = 2;i <= n * 2 +2;i++) {
        ta[i] = ta[i - 1] * A;
        // printf("on %d : %d\n",i,ta[i].a[n-1][0]) ;
    }
    for(int i = 1;i <= n;i++) {
        for(int j = 0 ; j <= n+2;j++) {
            val[i][j] = ta[j].a[n - 1][i - 1] ;
            // printf("%d %d : %d\n",i,j,val[i][j]) ;
        }
    }
    vector<int> v;
    vector<int> h( n + 1);
    mt19937 rnd(0); 
    for(int j = 1;j <= n;j++) h[j] = rnd()%mod ;
    for(int i = 0;i < n* 2 ;i++) {
        int d = 0;
        for(int j = 1;j <= n;j++) {
            d = (d + 1LL * ta[i].a[n - 1][j - 1] * h[j]) % mod;
        }
        v.push_back(d) ;
    }
    vector<int> coe = bm(v) ;
    // for(int i = 0 ;i < coe.size();i++) printf("%d ",coe[i]) ; printf("\n") ;
    // v.clear() ;
    // for(int i = 0;i < n* 2 ;i++) {
    //     v.push_back(ta[i].a[n - 1][2 - 1]) ;
    // }
    // coe = bm(v) ;
    // for(int i = 0 ;i < coe.size();i++) printf("%d ",coe[i]) ; printf("\n") ;

    // for(int i = coe.size() ; i < v.size();i++) {
    //     int d = 0;
    //     for(int j = 0;j < coe.size();j++) {
    //         d = (d + 1LL * coe[j] * v[i - j - 1]) % mod;
    //     }
    //     // printf("%d : %d %d\n",i,v[i],d);
    // }
    for(int i = 1;i <= m;i++) {
        int t , v;
        cin >> t >> v;
        // t=rand()%T+1;v=rand()%n+1;
        for(int j = t ; j < t + coe.size() ; j++) {
            f[j] = (f[j] + val[v][j - t]) % mod;
        }
        // ans[t + coe.size()] = (ans[t + coe.size()] + val[v][coe.size() ]) % mod;
        cal[t + coe.size()].push_back(v);
    }
    int E = 0;
    for(int i = 1;i <= T;i++ ) {
        if(i >= coe.size()) {
            for(auto v : cal[i]) {
                for(int j = 0;j < coe.size();j++) {
                    ans[i - j - 1] = (ans[i - j - 1] + val[v][coe.size() - j - 1]) % mod;
                }
            }
            for(int j = 0;j < coe.size();j++) ans[i] = (ans[i] + 1LL * coe[j] * ans[i - j - 1]) %mod;
        }
        sol[i] = (ans[i] + f[i]) % mod;
        // cout << i <<' ' << sol[i] << '\n';
        E ^= sol[i];
    }
    cout << E << '\n';
    return 0;
}