def main0(h,w,a,q,rc): mod=10**9+7 ans=[] for r,c in rc: r,c=r-1,c-1 tmp=1 for i in range(h): if i==r:continue for j in range(w): if j==c:continue tmp*=a[i][j] tmp%=mod ans.append(tmp) return ans def main1(h,w,a,q,rc): mod=10**9+7 mat=[] mat.append([1]*(w+1)) ary0=[] ary0.append([0]*(w+1)) for i in range(h): tmp=[1] tmp0=[0] for j in range(w): if a[i][j]: tmp.append(tmp[-1]*a[i][j]%mod) tmp0.append(tmp0[-1]) else: tmp.append(tmp[-1]) tmp0.append(tmp0[-1]+1) mat.append(tmp) ary0.append(tmp0) for j in range(w+1): for i in range(h): mat[i+1][j]*=mat[i][j] mat[i+1][j]%=mod ary0[i+1][j]+=ary0[i][j] ans=[] #for x in mat:print(x) for r,c in rc: r,c=r-1,c-1 # 0含むかどうかチェック t=ary0[r][c] t+=ary0[-1][c]-ary0[r+1][c] t+=ary0[r][-1]-ary0[r][c+1] t+=ary0[-1][-1]+ary0[r+1][c+1]-ary0[r+1][-1]-ary0[-1][c+1] if t: ans.append(0) continue tmp0=mat[r][c] # 左上 tmp1=mat[-1][c]*pow(mat[r+1][c],mod-2,mod)%mod # 左下 tmp2=mat[r][-1]*pow(mat[r][c+1],mod-2,mod)%mod # 右上 tmp3=mat[-1][-1]*mat[r+1][c+1]%mod # 右下 tmp3=tmp3*pow(mat[r+1][-1],mod-2,mod)*pow(mat[-1][c+1],mod-2,mod)%mod # 右下 ans.append(tmp0*tmp1*tmp2*tmp3%mod) return ans if __name__=='__main__': h,w=map(int,input().split()) a=[list(map(int,input().split())) for _ in range(h)] q=int(input()) rc=[list(map(int,input().split())) for _ in range(q)] #ret0=main0(h,w,a,q,rc) #print(ret0) ret1=main1(h,w,a,q,rc) print(ret1) from random import randint if __name__=='__main__1': for _ in range(100): h=randint(2,6) w=randint(2,6) a=[[randint(0,6) for _ in range(w)] for _ in range(h)] q=randint(2,6) rc=[[randint(1,h),randint(1,w)]for _ in range(q)] ret0=main0(h,w,a,q,rc) ret1=main1(h,w,a,q,rc) if ret0!=ret1: print(h,w) for x in a:print(*x) print(q) for x in rc:print(*x) print(ret0) print(ret1) break