結果
問題 | No.1141 田グリッド |
ユーザー | sbite |
提出日時 | 2020-07-31 23:07:13 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 220 ms / 2,000 ms |
コード長 | 2,550 bytes |
コンパイル時間 | 2,272 ms |
コンパイル使用メモリ | 211,348 KB |
実行使用メモリ | 11,592 KB |
最終ジャッジ日時 | 2023-09-21 02:11:04 |
合計ジャッジ時間 | 7,831 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge11 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 1 ms
4,380 KB |
testcase_01 | AC | 2 ms
4,380 KB |
testcase_02 | AC | 1 ms
4,376 KB |
testcase_03 | AC | 2 ms
4,380 KB |
testcase_04 | AC | 2 ms
4,376 KB |
testcase_05 | AC | 2 ms
4,376 KB |
testcase_06 | AC | 2 ms
4,376 KB |
testcase_07 | AC | 2 ms
4,380 KB |
testcase_08 | AC | 5 ms
4,380 KB |
testcase_09 | AC | 3 ms
4,380 KB |
testcase_10 | AC | 5 ms
4,380 KB |
testcase_11 | AC | 7 ms
4,380 KB |
testcase_12 | AC | 5 ms
4,376 KB |
testcase_13 | AC | 196 ms
6,168 KB |
testcase_14 | AC | 220 ms
11,592 KB |
testcase_15 | AC | 200 ms
5,376 KB |
testcase_16 | AC | 201 ms
5,424 KB |
testcase_17 | AC | 204 ms
5,368 KB |
testcase_18 | AC | 151 ms
5,348 KB |
testcase_19 | AC | 148 ms
5,400 KB |
testcase_20 | AC | 149 ms
5,396 KB |
testcase_21 | AC | 200 ms
5,300 KB |
testcase_22 | AC | 204 ms
5,356 KB |
testcase_23 | AC | 203 ms
5,612 KB |
testcase_24 | AC | 149 ms
5,300 KB |
testcase_25 | AC | 152 ms
5,168 KB |
testcase_26 | AC | 148 ms
5,424 KB |
testcase_27 | AC | 150 ms
5,416 KB |
testcase_28 | AC | 152 ms
5,272 KB |
testcase_29 | AC | 147 ms
5,388 KB |
testcase_30 | AC | 149 ms
5,400 KB |
testcase_31 | AC | 149 ms
4,948 KB |
testcase_32 | AC | 147 ms
5,340 KB |
testcase_33 | AC | 150 ms
5,356 KB |
ソースコード
#include <bits/stdc++.h> #define _overload3(_1, _2, _3, name, ...) name #define _rep(i, n) repi(i, 0, n) #define repi(i, a, b) for (int i = (a); i < (b); ++i) #define rep(...) _overload3(__VA_ARGS__, repi, _rep, )(__VA_ARGS__) #define ALL(x) x.begin(), x.end() #define chmax(x, y) x = max(x, y) #define chmin(x, y) x = min(x, y) using namespace std; random_device rnd; mt19937 mt(rnd()); using ll = long long; using lld = long double; using VI = vector<int>; using VVI = vector<VI>; using VL = vector<ll>; using VVL = vector<VL>; using PII = pair<int, int>; const int IINF = 1 << 30; const ll INF = 1ll << 60; const ll MOD = 1000000007; ll mpow(ll base, ll num) { if (num <= 0) return 1; ll prev = mpow(base, num / 2); if (num % 2 == 0) { return (prev * prev) % MOD; } else { return (((prev * prev) % MOD) * base) % MOD; } } ll inv(ll num) { return mpow(num, MOD - 2); } int main() { ll h, w; cin >> h >> w; VVL g(h, VL(w, 0)); VVI zero(h, VI(w, 0)); int zeronum = 0; rep(i, h) rep(j, w) { cin >> g[i][j]; if (g[i][j] == 0) { zero[i][j] = 1; g[i][j] = 1; zeronum++; } } VI zeroy(h, 0); VI zerox(w, 0); rep(i, h) { rep(j, w) { zeroy[i] += zero[i][j]; zerox[j] += zero[i][j]; } } VVL multi(h + 1, VL(w + 1, 1)); rep(i, h) multi[i + 1][1] = (multi[i][1] * g[i][0]) % MOD; rep(i, w) multi[1][i + 1] = (multi[1][i] * g[0][i]) % MOD; rep(i, 1, h) { rep(j, 1, w) { //multi[i + 1][j + 1] = mpow(multi[i][j + 1] * multi[i + 1][j], multi[i][j] - 2) * g[i][j]; multi[i + 1][j + 1] = (((((multi[i][j + 1] * multi[i + 1][j]) % MOD) * inv(multi[i][j])) % MOD) * g[i][j]) % MOD; multi[i + 1][j + 1] %= MOD; } } ll total = multi[h][w]; int q; cin >> q; rep(i, q) { int x, y; cin >> y >> x; int z = zeronum; z -= zeroy[y - 1]; z -= zerox[x - 1]; z += zero[y - 1][x - 1]; //cerr << z << endl; if (z != 0) { cout << 0 << endl; continue; } ll tmp = multi[y][w] * inv(multi[y - 1][w]); tmp %= MOD; tmp *= (multi[h][x] * inv(multi[h][x - 1])) % MOD; tmp %= MOD; tmp *= inv(g[y - 1][x - 1]); tmp %= MOD; cout << (total * inv(tmp)) % MOD << endl; } return 0; }