結果

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

ソースコード

diff #
プレゼンテーションモードにする

#include <iostream>
#include <algorithm>
#include <numeric>
#include <vector>
#include <string>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <stack>
#include <iomanip>
#include <functional>
#include <bitset>
#include <limits>
#include <cstdio>
#include <cmath>
#include <cassert>
#include <random>
#ifdef DEBUG
#include "library/Utility/debug.cpp"
#else
#define debug(...)
#endif
#define rep(i,n) for(int i=0;i<(n);++i)
#define EL '\n'
#define print(i) std::cout << (i) << '\n'
#define all(v) (v).begin(), (v).end()
using lnt = long long;
struct FIO{FIO(){std::cin.tie(0);std::cout.tie(0);std::ios_base::sync_with_stdio(0);std::cout<<std::fixed<<std::setprecision(15);}}fIO;
/*-*/
constexpr lnt MOD = 1e9+7;
struct mint
{
    lnt v;
    mint():v(0){}
    mint(lnt v):v((v+MOD)%MOD){}
    mint operator-()const{ return mint(0) - *this; }
    mint& operator+=(const mint& a){ if((v+=a.v)>=MOD) v-=MOD; return *this; }
    mint& operator-=(const mint& a){ if((v+=MOD-a.v)>=MOD) v-=MOD; return *this; }
    mint& operator*=(const mint& a){ (v*=a.v)%=MOD; return *this; }
    mint& operator/=(const mint& a){ (*this) *= a.inv(); return *this; }
    mint operator+(const mint& a)const{ return mint(*this) += a; }
    mint operator-(const mint& a)const{ return mint(*this) -= a; }
    mint operator*(const mint& a)const{ return mint(*this) *= a; }
    mint operator/(const mint& a)const{ return mint(*this) /= a; }
    bool operator<(const mint& a)const{ return v < a.v; }
    bool operator==(const mint& a)const{ return v == a.v; }
    mint pow(lnt k)const{ mint r(1),t(v); while(k){ if(k&1) r*=t; t*=t; k>>=1; } return r; }
    mint inv()const{ return pow(MOD-2); }
    static mint comb(lnt n, lnt k) { if(n-k<k) k=n-k; mint num(1), dom(1); for(int i=0;i<k;i++) { num*=n-i; dom*=i+1; } return num/dom; }
    static std::vector<mint> construct_comb(int n) {
        std::vector<mint> c(n+1); mint a = 1; c[0] = a; for(int i=1;i<=n;i++) { a = a*mint(n+1-i)/i; c[i] = a; } return c;
    }
    static std::vector<mint> construct_fact(int n) { std::vector<mint> f(n+1,1); for(int i=2;i<=n;i++) f[i]=f[i-1]*i; return f; }
};
std::istream& operator>>(std::istream&i,mint&a){ lnt t; i>>t; a=mint(t); return i; }
std::ostream& operator<<(std::ostream&o,const mint&a){ o<<a.v; return o; }
int main() {
    int h,w;
    std::cin >> h >> w;
    std::vector<std::vector<lnt> > g(h,std::vector<lnt>(w));
    rep(i,h) rep(j,w) std::cin >> g[i][j];
    std::vector<std::vector<std::vector<mint> > > A(4,std::vector<std::vector<mint> >(h,std::vector<mint>(w)));
    rep(k,4) rep(i,h) rep(j,w) A[k][i][j]=mint(g[i][j]);
    rep(k,4) {
        rep(x,h) rep(y,w) {
            int i=(k&1)?x:h-x-1;
            int j=(k>=2)?y:w-y-1;
            if(x==0&&y==0) A[k][i][j]=mint(g[i][j]);
            else if(x==0) A[k][i][j]=mint(g[i][j])*A[k][i][(k>1)?j-1:j+1];
            else if(y==0) A[k][i][j]=mint(g[i][j])*A[k][(k&1)?i-1:i+1][j];
            else {
                A[k][i][j]=A[k][(k&1)?i-1:i+1][j]*mint(g[i][j])*A[k][i][(k>=2)?j-1:j+1]/A[k][(k&1)?i-1:i+1][(k>=2)?j-1:j+1];
            }
        }
    }
    debug(A[0]);
    int q;
    std::cin >> q;
    rep(qq,q) {
        int r,c;
        std::cin >> r >> c;
        r--;c--;
        mint ans=1;
        debug(r,c);
        if(r>0&&c>0) ans*=A[3][r-1][c-1];
        if(r>0&&c<w-1) ans*=A[1][r-1][c+1];
        if(r<h-1&&c>0) ans*=A[2][r+1][c-1];
        if(r<h-1&&c<w-1) ans*=A[0][r+1][c+1];
        print(ans);
    }
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0