結果
問題 |
No.1141 田グリッド
|
ユーザー |
![]() |
提出日時 | 2020-07-31 22:05:12 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,490 bytes |
コンパイル時間 | 1,919 ms |
コンパイル使用メモリ | 123,884 KB |
最終ジャッジ日時 | 2025-01-12 10:11:57 |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 21 WA * 10 |
ソースコード
#include<iostream> #include<string> #include<vector> #include<algorithm> #include<bitset> #include<set> #include<map> #include<stack> #include<queue> #include<deque> #include<list> #include<iomanip> #include<cmath> #include<cstring> #include<functional> #include<cstdio> #include<cstdlib> using namespace std; #define repr(i, a, b) for (int i = (int)(a); i < (int)(b); i++) #define rep(i, n) repr(i, 0, n) #define INF 2e9 #define MOD 1000000007 //#define MOD 998244353 #define LINF (long long)4e18 #define jck 3.141592 #define PI acos(-1) const double EPS = 1e-10; using ll = long long; using Pi = pair<int,int>; using Pl = pair<ll,ll>; 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(){ int h,w; cin >> h >> w; vector<vector<int>> a(h,vector<int>(w)); rep(i,h)rep(j,w) cin >> a[i][j]; int q; cin >> q; vector<vector<ll>> seki(h+1,vector<ll>(w+1)); rep(i,h+1) seki[i][0] = 1; rep(i,w+1) seki[0][i] = 1; rep(i,h)rep(j,w){ seki[i+1][j+1] = a[i][j]*seki[i+1][j]%MOD*seki[i][j+1]%MOD*modinv(seki[i][j],MOD)%MOD; } while(q--){ int r,c; cin >> r >> c; r--; c--; ll ans = seki[h][w]*modinv(seki[r+1][w]*modinv(seki[r][w],MOD)%MOD,MOD)%MOD*modinv(seki[h][c+1]*modinv(seki[h][c]%MOD,MOD),MOD)%MOD*a[r][c]%MOD; cout << ans << endl; } }