結果

問題 No.1141 田グリッド
ユーザー yuji9511
提出日時 2020-07-31 23:20:28
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 59 ms / 2,000 ms
コード長 2,237 bytes
コンパイル時間 2,272 ms
コンパイル使用メモリ 199,616 KB
最終ジャッジ日時 2025-01-12 11:44:46
ジャッジサーバーID
(参考情報)
judge5 / judge5
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

diff #

/*** author: yuji9511 ***/
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using lpair = pair<ll, ll>;
const ll MOD = 1e9+7;
const ll INF = 1e18;
#define rep(i,m,n) for(ll i=(m);i<(n);i++)
#define rrep(i,m,n) for(ll i=(m);i>=(n);i--)
#define printa(x,n) for(ll i=0;i<n;i++){cout<<(x[i])<<" \n"[i==n-1];};
void print() {}
template <class H,class... T>
void print(H&& h, T&&... t){cout<<h<<" \n"[sizeof...(t)==0];print(forward<T>(t)...);}
    ll power(ll x, ll n){
        if(n == 0) return 1LL;
        ll res = power(x * x % MOD, n/2);
        if(n % 2 == 1) res = res * x % MOD;
        return res;
    }

void solve(){
    ll H,W;
    cin >> H >> W;
    vector< vector<ll> > A(H+1, vector<ll>(W+1));
    rep(i,0,H){
        rep(j,0,W){
            cin >> A[i][j];
        }
    }
    ll mul = 1;
    rep(i,0,H){
        rep(j,0,W){
            if(A[i][j] != 0){
                mul *= A[i][j];
                mul %= MOD;
            }
        }
    }
    ll tate[100010], yoko[100010];
    rep(i,0,H){
        ll res = 1;
        rep(j,0,W){
            if(A[i][j] != 0){
                res *= A[i][j];
                res %= MOD;
            }
        }
        tate[i] = res;
    }
    rep(j,0,W){
        ll res = 1;
        rep(i,0,H){
            if(A[i][j] != 0){
                res *= A[i][j];
                res %= MOD;
            }
        }
        yoko[j] = res;
    }
    vector< vector<ll> > sum(H+2, vector<ll>(W+2, 0));
    ll n1[100010] = {}, n2[100010] = {};
    ll zero_cnt = 0;
    rep(i,0,H){
        rep(j,0,W){
            if(A[i][j] == 0){
                n1[i]++;
                n2[j]++;
                zero_cnt++;
            }
        }
    }

    ll Q;
    cin >> Q;
    while(Q--){
        ll r,c;
        cin >> r >> c;
        r--; c--;
        ll num = n1[r] + n2[c];
        if(A[r][c] == 0) num--;
        if(num != zero_cnt){
            print(0);
        }else{
            ll ans = mul * power(tate[r], MOD-2) % MOD * power(yoko[c], MOD-2) % MOD;
            if(A[r][c] != 0){
                ans *= A[r][c];
                ans %= MOD;
            }
            print(ans);
        }

    }


}

int main(){
    cin.tie(0);
    ios::sync_with_stdio(false);
    solve();
}
0