#include <bits/stdc++.h>
using namespace std;
using ll=long long;
using vin=vector<int>;
using vll=vector<long long>;
using vvin=vector<vector<int>>;
using vvll=vector<vector<long long>>;
using vstr=vector<string>;
using vvstr=vector<vector<string>>;
using vch=vector<char>;
using vvch=vector<vector<char>>;
using vbo=vector<bool>;
using vvbo=vector<vector<bool>>;
using vpii=vector<pair<int,int>>;
using pqsin=priority_queue<int,vector<int>,greater<int>>;
#define mp make_pair
#define rep(i,n) for(int i=0;i<(int)(n);i++)
#define rep2(i,s,n) for(int i=(s);i<(int)(n);i++)
#define all(v) v.begin(),v.end()
#define decp(n) cout<<fixed<<setprecision((int)n)
const ll inf=1e9+7;
const ll INF=1e18;

ll mondiv(ll a,ll m=inf){
    ll b=m;ll u=1;ll v=0;
    while(b){
        ll t=a/b;
        a-=t*b;swap(a,b);
        u-=t*v;swap(u,v);
    }
    u%=m;
    if(u<0)u+=m;
    return u;
}

int main(){
    int h,w;cin>>h>>w;
    vvll a(h,vll(w));
    vll row(h,(ll)1),column(w,(ll)1);
    vin zerorow(h),zerocolumn(w);
    ll allzero=(ll)0;
    bool zero=false;
    ll ap=(ll)1;ll tmp;
    rep(i,h)rep(j,w){
        cin>>a[i][j];
        if(a[i][j]==0){
            zerorow[i]++;
            zerocolumn[j]++;
            allzero++;
            zero=true;
            continue;
        }
        ap*=a[i][j];ap%=inf;
        tmp=mondiv(a[i][j]);
        row[i]*=tmp;row[i]%=inf;
        column[j]*=tmp;column[j]%=inf;
    }
    int q;cin>>q;
    vll ans;int r,c;
    rep(i,q){
        cin>>r>>c;r--;c--;
        if(zero){
             tmp=zerorow[r]+zerocolumn[c];
             if(a[r][c]==0)tmp--;
             if(tmp!=allzero){
                 ans.push_back(0);
                 continue;
             }
        }
        tmp=ap*row[r];tmp%=inf;
        tmp*=column[c];tmp%=inf;
        if(a[r][c]!=0)tmp*=a[r][c];tmp%=inf;
        ans.push_back(tmp);
    }
    rep(i,ans.size())cout<<ans[i]<<endl;
}