import times, strutils, sequtils, math, algorithm, tables, sets, lists, intsets import critbits, future, strformat, deques template `max=`(x,y) = x = max(x,y) template `min=`(x,y) = x = min(x,y) template `mod=`(x,y) = x = x mod y template scan2 = (scan(), scan()) template scan3 = (scan(), scan()) let read* = iterator: string {.closure.} = while true: (for s in stdin.readLine.split: yield s) proc scan(): int = read().parseInt proc scanf(): float = read().parseFloat proc toInt(c:char): int = return int(c) - int('0') # 拡張ユークリッド互除法 proc xgcd(a,b:int):(int,int,int)= var b = b x0=1 y0=0 x1=0 y1=1 q:int a = a while b!=0: (q,a,b) = (a.div(b),b,a.mod(b)) (x0,x1)=(x1,x0-q*x1) (y0,y1)=(y1,y0-q*y1) return (a,x0,y0) # modの逆元 proc modinv(a,m:int):int= var (g,x,y)=xgcd(a,m) if g!=1: return -1 else: return x mod m proc solve()= var h = scan() w = scan() a = newseqwith(h,newseqwith(w,scan())) board_UL = newseqwith(h+1,newseqwith(w+1,1)) board_UR = newseqwith(h+1,newseqwith(w+1,1)) board_LL = newseqwith(h+1,newseqwith(w+1,1)) board_LR = newseqwith(h+1,newseqwith(w+1,1)) q = scan() yx = newseqwith(q,(scan(),scan())) md = 10^9+7 for y in 1..h: for x in 1..w: boardUL[y][x] = ((((a[y-1][x-1]*boardUL[y-1][x]).mod(md)*(boardUL[y][x-1]*modinv(boardUL[y-1][x-1],md)).mod(md)).mod(md))+md).mod(md) for y in 1..h: for x in countdown(w-1,0): boardUR[y][x] = ((((a[y-1][x]*boardUR[y-1][x]).mod(md)*(boardUR[y][x+1]*modinv(boardUR[y-1][x+1],md)).mod(md)).mod(md))+md).mod(md) for y in countdown(h-1,0): for x in 1..w: boardLL[y][x] = ((((a[y][x-1]*boardLL[y+1][x]).mod(md)*(boardLL[y][x-1]*modinv(boardLL[y+1][x-1],md)).mod(md)).mod(md))+md).mod(md) for y in countdown(h-1,0): for x in countdown(w-1,0): boardLR[y][x] = ((((a[y][x]*boardLR[y+1][x]).mod(md)*(boardLR[y][x+1]*modinv(boardLR[y+1][x+1],md)).mod(md)).mod(md))+md).mod(md) #echo boardUL.join("\n") #echo "" #echo boardUR.join("\n") #echo "" #echo boardLL.join("\n") #echo "" #echo boardLR.join("\n") for (y, x) in yx: var t = 1 t*=board_UL[y-1][x-1] t.mod=md t*=board_UR[y-1][x] t.mod=md t*=board_LL[y][x-1] t.mod=md t*=board_LR[y][x] t.mod=md echo t solve()