結果
| 問題 |
No.1141 田グリッド
|
| ユーザー |
umezo
|
| 提出日時 | 2020-08-02 03:42:13 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 179 ms / 2,000 ms |
| コード長 | 1,564 bytes |
| コンパイル時間 | 2,155 ms |
| コンパイル使用メモリ | 198,656 KB |
| 最終ジャッジ日時 | 2025-01-12 13:06:10 |
|
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 31 |
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:57:6: warning: ‘tmp’ may be used uninitialized [-Wmaybe-uninitialized]
57 | ll tmp;
| ^~~
ソースコード
#define rep(i, n) for (int i = 0; i < (int)(n); i++)
#define ALL(v) v.begin(), v.end()
typedef long long ll;
#include <bits/stdc++.h>
using namespace std;
const int MOD=1e9+7;
ll modinv(ll a,ll m){
ll b=m,u=1,v=0;
while(b){
ll 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(){
ll h,w;
cin>>h>>w;
ll A[h+1][w+1],Z[h+1][w+1];
ll x;
for(int i=1;i<=h;i++){
for(int j=1;j<=w;j++){
cin>>x;
if(x==0){
A[i][j]=1;
Z[i][j]=1;
}
else{
A[i][j]=x;
Z[i][j]=0;
}
}
}
ll s[h+2][w+2];
rep(i,h+2){
rep(j,w+2) s[i][j]=0;
}
s[1][1]=0;
for(int i=1;i<h+1;i++){
for(int j=1;j<w+1;j++) s[i+1][j+1]=s[i+1][j]+s[i][j+1]-s[i][j]+Z[i][j];
}
ll all=1;
for(int i=1;i<=h;i++){
for(int j=1;j<=w;j++) all=(all*A[i][j])%MOD;
}
vector<ll> B(h+1),C(w+1);
ll tmp;
for(int i=1;i<=h;i++){
for(int j=1;j<=w;j++){
if(j==1) tmp=A[i][1];
else tmp=(tmp*A[i][j])%MOD;
}
B[i]=modinv(tmp,MOD);
}
for(int j=1;j<=w;j++){
for(int i=1;i<=h;i++){
if(i==1) tmp=A[1][j];
else tmp=(tmp*A[i][j])%MOD;
}
C[j]=modinv(tmp,MOD);
}
ll q;
cin>>q;
ll r,c;
rep(i,q){
cin>>r>>c;
ll a=0;
a+=s[r][c];
a+=s[r][w+1]-s[r][c+1];
a+=s[h+1][c]-s[r+1][c];
a+=s[h+1][w+1]-s[r+1][w+1]-s[h+1][c+1]+s[r+1][c+1];
if(a>0) cout<<0<<endl;
else{
ll ans=(((all*B[r]%MOD)*C[c]%MOD)*A[r][c])%MOD;
cout<<ans<<endl;
}
}
return 0;
}
umezo