結果

問題 No.1141 田グリッド
コンテスト
ユーザー k
提出日時 2020-09-05 05:15:04
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 115 ms / 2,000 ms
コード長 1,486 bytes
コンパイル時間 2,304 ms
コンパイル使用メモリ 200,760 KB
最終ジャッジ日時 2025-01-14 06:56:29
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>

using namespace std;

#define REP(i,n) for(int i=0; i<(int)(n); i++)

const long long MOD = 1e9 + 7;

long long modpow(long long x, long long p, long long mod) {
  long long ret = 1;
  while (p) {
    if (p & 1)
      ret = ret * x % mod;
    x = x * x % mod;
    p >>= 1;
  }
  return ret;
}

int main() {
  ios_base::sync_with_stdio(0);
  cin.tie(0);

  int h, w;
  cin >> h >> w;
  vector<vector<int> > a(h, vector<int>(w));

  long long sum = 1;
  int zero = 0;
  REP (i, h) REP (j, w) {
    cin >> a[i][j];
    if (a[i][j])
      sum = sum * a[i][j] % MOD;
    else
      ++zero;
  }
  
  vector<long long> hsum(h), wsum(w);
  vector<int> hzero(h), wzero(w);
  REP (i, h) {
    hsum[i] = 1;
    hzero[i] = 0;
    REP (j, w) {
      if (a[i][j])
        hsum[i] = hsum[i] * a[i][j] % MOD;
      else
        hzero[i]++;
    }
  }
  REP (j, w) {
    wsum[j] = 1;
    wzero[j] = 0;
    REP (i, h) {
      if (a[i][j])
        wsum[j] = wsum[j] * a[i][j] % MOD;
      else
        wzero[j]++;
    }
  }
  
  int q;
  cin >> q;
  while (q--) {
    int r, c;
    cin >> r >> c;
    --r, --c;

    int z = hzero[r] + wzero[c] - !a[r][c];
    if (z < zero)
      cout << 0 << endl;
    else {
      long long ret = 1;
      ret = ret * sum % MOD;
      ret = ret * modpow(hsum[r], MOD-2, MOD) % MOD;
      ret = ret * modpow(wsum[c], MOD-2, MOD) % MOD;
      if (a[r][c])
        ret = ret * a[r][c] % MOD;
      cout << ret << endl;
    }
  }
    
  return 0;
}
0