結果
問題 | No.1141 田グリッド |
ユーザー |
|
提出日時 | 2020-07-31 22:04:06 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 194 ms / 2,000 ms |
コード長 | 2,269 bytes |
コンパイル時間 | 2,276 ms |
コンパイル使用メモリ | 202,524 KB |
最終ジャッジ日時 | 2025-01-12 10:09:49 |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 31 |
ソースコード
#include <bits/stdc++.h>using namespace std;#define rep(i, n) for(ll i = 0, i##_len = (n); i < i##_len; ++i)#define rep2(i, x, n) for(ll i = x, i##_len = (n); i < i##_len; ++i)#define all(n) begin(n), end(n)using ll = long long;using P = pair<ll, ll>;using vi = vector<int>;using vl = vector<ll>;using vs = vector<string>;using vc = vector<char>;using vb = vector<bool>;using vd = vector<double>;vi dir = {-1, 0, 1, 0, -1, -1, 1, 1, -1};const int mod = 1000000007;struct mint {ll x; // typedef long long ll;mint(ll x = 0) : x((x % mod + mod) % mod) {}mint operator-() const { return mint(-x); }mint& operator+=(const mint a) {if((x += a.x) >= mod) x -= mod;return *this;}mint& operator-=(const mint a) {if((x += mod - a.x) >= mod) x -= mod;return *this;}mint& operator*=(const mint a) {(x *= a.x) %= mod;return *this;}mint operator+(const mint a) const { return mint(*this) += a; }mint operator-(const mint a) const { return mint(*this) -= a; }mint operator*(const mint a) const { return mint(*this) *= a; }mint pow(ll t) const {if(!t) return 1;mint a = pow(t >> 1);a *= a;if(t & 1) a *= *this;return a;}// for prime modmint inv() const { return pow(mod - 2); }mint& operator/=(const mint a) { return *this *= a.inv(); }mint operator/(const mint a) const { return mint(*this) /= a; }};istream& operator>>(istream& is, const mint& a) { return is >> a.x; }ostream& operator<<(ostream& os, const mint& a) { return os << a.x; }mint fact(int n) {if(n == 1 || n == 0) return 1;mint res = n;return res * fact(n - 1);}int main() {ll h, w;cin >> h >> w;vector<mint> r(h, 1), c(w, 1);vector<vl> a(h, vl(w, 1)), za(h, vl(w, 0));vl zr(h, 0), zc(w, 0);ll zero = 0;mint s = 1;rep(i, h) rep(j, w) {ll b;cin >> b;if(b) {a[i][j] = b;r[i] *= b;c[j] *= b;s *= b;} else {za[i][j] = 1;zr[i]++;zc[j]++;zero++;}}ll q;cin >> q;rep(_, q) {ll i, j;cin >> i >> j;i--;j--;ll nz = zr[i] + zc[j] - za[i][j];if(nz == zero)cout << s / r[i] / c[j] * a[i][j] << '\n';elsecout << 0 << '\n';}}