結果
問題 |
No.1141 田グリッド
|
ユーザー |
![]() |
提出日時 | 2021-03-20 14:56:19 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 604 ms / 2,000 ms |
コード長 | 1,070 bytes |
コンパイル時間 | 161 ms |
コンパイル使用メモリ | 82,300 KB |
実行使用メモリ | 95,300 KB |
最終ジャッジ日時 | 2024-11-21 00:43:09 |
合計ジャッジ時間 | 9,642 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 31 |
ソースコード
import sys import io, os input = io.BytesIO(os.read(0,os.fstat(0).st_size)).readline h, w = map(int, input().split()) A = [list(map(int, input().split())) for _ in range(h)] mod = 1000000007 X = [[1]*(w+1) for i in range(h+1)] total0 = 0 row0 = [0]*h col0 = [0]*w for i in range(h): for j in range(w): if A[i][j] != 0: X[i+1][j+1] = X[i+1][j]*X[i][j+1]*pow(X[i][j],mod-2,mod)*A[i][j] else: X[i+1][j+1] = X[i+1][j]*X[i][j+1]*pow(X[i][j],mod-2,mod)*1 total0 += 1 row0[i] += 1 col0[j] += 1 X[i+1][j+1] %= mod def rect(r1, c1, r2, c2): res = X[r2][c2]*pow(X[r2][c1],mod-2, mod)*pow(X[r1][c2],mod-2, mod)*X[r1][c1] return res%mod #print(X) #print(X[h][w]%mod) q = int(input()) for i in range(q): r, c = map(int, input().split()) r, c = r-1, c-1 z = total0-row0[r]-col0[c] if A[r][c] == 0: z += 1 if z > 0: print(0) else: ans = rect(r+1, c+1, h, w)*rect(0, c+1, r, w)*rect(r+1, 0, h, c)*rect(0, 0, r, c) print(ans%mod)