結果
問題 | No.1141 田グリッド |
ユーザー |
![]() |
提出日時 | 2020-07-31 22:46:52 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 96 ms / 2,000 ms |
コード長 | 2,159 bytes |
コンパイル時間 | 2,304 ms |
コンパイル使用メモリ | 182,096 KB |
実行使用メモリ | 7,296 KB |
最終ジャッジ日時 | 2024-07-06 20:07:33 |
合計ジャッジ時間 | 5,247 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 31 |
ソースコード
#pragma GCC optimize("Ofast")#include <bits/stdc++.h>using namespace std;using ll = long long;using pii = pair<int, int>;using pll = pair<long long, long long>;constexpr char ln = '\n';constexpr long long MOD = 1000000007LL;constexpr long long INF = 1001001001LL;constexpr long long LINF = 1001001001001001001;#define all(x) (x).begin(),(x).end()#define rep(i,n) for(int i=0;i<(n);i++)#define rept(i, j, n) for(int i=(j); i<(n); i++)template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; }template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; }long long modinv(long long a, long long m) {long long b = m, u = 1, v = 0;while (b) {long long 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;vector<vector<ll>> grid(h, vector<ll>(w));vector<ll> countZ_h(h), countZ_w(w);ll sum = 1, zero = 0;rep(i, h){rep(j, w){cin >> grid[i][j];if(grid[i][j])sum = (sum * grid[i][j]) % MOD;else{zero++;countZ_h[i]++;countZ_w[j]++;}}}vector<ll> H(h), W(w);rep(i, h){ll num = 1;rep(j, w){if(grid[i][j])num = (num * grid[i][j]) % MOD;}H[i] = modinv(num, MOD);}rep(j, w){ll num = 1;rep(i, h){if(grid[i][j])num = (num * grid[i][j]) % MOD;}W[j] = modinv(num, MOD);}int Q; cin >> Q;vector<ll> ans(Q);rep(q, Q){ll r, c; cin >> r >> c; r--, c--;ll cnt = countZ_h[r] + countZ_w[c];if(grid[r][c]==0)cnt--;if(cnt < zero)ans[q] = 0;else{ll res = sum;res = (res * H[r]) % MOD;res = (res * W[c]) % MOD;if(grid[r][c])res = (res * grid[r][c]) % MOD;ans[q] = res;}}for(auto a: ans)cout << a << ln;}