結果

問題 No.1141 田グリッド
ユーザー どらら
提出日時 2020-08-01 01:26:31
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 148 ms / 2,000 ms
コード長 2,275 bytes
コンパイル時間 1,823 ms
コンパイル使用メモリ 177,312 KB
実行使用メモリ 16,896 KB
最終ジャッジ日時 2024-07-07 00:01:43
合計ジャッジ時間 5,915 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define REP(i,a,n) for(int i=(a); i<(int)(n); i++)
#define rep(i,n) REP(i,0,n)
#define FOR(it,c) for(__typeof((c).begin()) it=(c).begin(); it!=(c).end(); ++it)
#define ALLOF(c) (c).begin(), (c).end()
typedef long long ll;
typedef unsigned long long ull;

static const ll MOD = 1000000007;

int main(){
  int H, W;
  cin >> H >> W;
  vector<vector<int>> v(H, vector<int>(W));
  rep(i,H){
    rep(j,W){
      cin >> v[i][j];
    }
  }

  vector<vector<ll>> dp1(H, vector<ll>(W));
  vector<vector<ll>> dp2(H, vector<ll>(W));
  vector<vector<ll>> dp3(H, vector<ll>(W));
  vector<vector<ll>> dp4(H, vector<ll>(W));

  {
    ll p = 1;
    for(int y=0; y<H; y++){
      p *= v[y][0];
      p %= MOD;
      dp1[y][0] = p;
    }
    for(int x=1; x<W; x++){
      p = 1;
      for(int y=0; y<H; y++){
        p *= v[y][x];
        p %= MOD;
        dp1[y][x] = (dp1[y][x-1] * p)%MOD;
      }
    }
  }
  {
    ll p = 1;
    for(int y=0; y<H; y++){
      p *= v[y][W-1];
      p %= MOD;
      dp2[y][W-1] = p;
    }
    for(int x=W-2; x>=0; x--){
      p = 1;
      for(int y=0; y<H; y++){
        p *= v[y][x];
        p %= MOD;
        dp2[y][x] = (dp2[y][x+1] * p)%MOD;
      }
    }
  }
  {
    ll p = 1;
    for(int y=H-1; y>=0; y--){
      p *= v[y][0];
      p %= MOD;
      dp3[y][0] = p;
    }
    for(int x=1; x<W; x++){
      p = 1;
      for(int y=H-1; y>=0; y--){
        p *= v[y][x];
        p %= MOD;
        dp3[y][x] = (dp3[y][x-1] * p)%MOD;
      }
    }
  }
  {
    ll p = 1;
    for(int y=H-1; y>=0; y--){
      p *= v[y][W-1];
      p %= MOD;
      dp4[y][W-1] = p;
    }
    for(int x=W-2; x>=0; x--){
      p = 1;
      for(int y=H-1; y>=0; y--){
        p *= v[y][x];
        p %= MOD;
        dp4[y][x] = (dp4[y][x+1] * p)%MOD;
      }
    }
  }
  
  int Q;
  cin >> Q;
  rep(i,Q){
    int r, c;
    cin >> r >> c;
    r--; c--;

    ll ret = 1;
    if(r != 0 && c != 0){
      ret *= dp1[r-1][c-1];
      ret %= MOD;
    }
    if(r != 0 && c != W-1){
      ret *= dp2[r-1][c+1];
      ret %= MOD;
    }
    if(r != H-1 && c != 0){
      ret *= dp3[r+1][c-1];
      ret %= MOD;
    }
    if(r != H-1 && c != W-1){
      ret *= dp4[r+1][c+1];
      ret %= MOD;
    }
    cout << ret << endl;
  }
  

  return 0;
}
0