結果

問題 No.1141 田グリッド
ユーザー sbitesbite
提出日時 2020-07-31 23:07:13
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 220 ms / 2,000 ms
コード長 2,550 bytes
コンパイル時間 2,272 ms
コンパイル使用メモリ 211,348 KB
実行使用メモリ 11,592 KB
最終ジャッジ日時 2023-09-21 02:11:04
合計ジャッジ時間 7,831 ms
ジャッジサーバーID
(参考情報)
judge12 / judge11
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,380 KB
testcase_01 AC 2 ms
4,380 KB
testcase_02 AC 1 ms
4,376 KB
testcase_03 AC 2 ms
4,380 KB
testcase_04 AC 2 ms
4,376 KB
testcase_05 AC 2 ms
4,376 KB
testcase_06 AC 2 ms
4,376 KB
testcase_07 AC 2 ms
4,380 KB
testcase_08 AC 5 ms
4,380 KB
testcase_09 AC 3 ms
4,380 KB
testcase_10 AC 5 ms
4,380 KB
testcase_11 AC 7 ms
4,380 KB
testcase_12 AC 5 ms
4,376 KB
testcase_13 AC 196 ms
6,168 KB
testcase_14 AC 220 ms
11,592 KB
testcase_15 AC 200 ms
5,376 KB
testcase_16 AC 201 ms
5,424 KB
testcase_17 AC 204 ms
5,368 KB
testcase_18 AC 151 ms
5,348 KB
testcase_19 AC 148 ms
5,400 KB
testcase_20 AC 149 ms
5,396 KB
testcase_21 AC 200 ms
5,300 KB
testcase_22 AC 204 ms
5,356 KB
testcase_23 AC 203 ms
5,612 KB
testcase_24 AC 149 ms
5,300 KB
testcase_25 AC 152 ms
5,168 KB
testcase_26 AC 148 ms
5,424 KB
testcase_27 AC 150 ms
5,416 KB
testcase_28 AC 152 ms
5,272 KB
testcase_29 AC 147 ms
5,388 KB
testcase_30 AC 149 ms
5,400 KB
testcase_31 AC 149 ms
4,948 KB
testcase_32 AC 147 ms
5,340 KB
testcase_33 AC 150 ms
5,356 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
#define _overload3(_1, _2, _3, name, ...) name
#define _rep(i, n) repi(i, 0, n)
#define repi(i, a, b) for (int i = (a); i < (b); ++i)
#define rep(...) _overload3(__VA_ARGS__, repi, _rep, )(__VA_ARGS__)
#define ALL(x) x.begin(), x.end()
#define chmax(x, y) x = max(x, y)
#define chmin(x, y) x = min(x, y)
using namespace std;
random_device rnd;
mt19937 mt(rnd());
using ll = long long;
using lld = long double;
using VI = vector<int>;
using VVI = vector<VI>;
using VL = vector<ll>;
using VVL = vector<VL>;
using PII = pair<int, int>;
const int IINF = 1 << 30;
const ll INF = 1ll << 60;
const ll MOD = 1000000007;

ll mpow(ll base, ll num)
{
    if (num <= 0)
        return 1;
    ll prev = mpow(base, num / 2);
    if (num % 2 == 0)
    {
        return (prev * prev) % MOD;
    }
    else
    {
        return (((prev * prev) % MOD) * base) % MOD;
    }
}

ll inv(ll num)
{
    return mpow(num, MOD - 2);
}

int main()
{
    ll h, w;
    cin >> h >> w;
    VVL g(h, VL(w, 0));
    VVI zero(h, VI(w, 0));

    int zeronum = 0;
    rep(i, h) rep(j, w)
    {
        cin >> g[i][j];
        if (g[i][j] == 0)
        {
            zero[i][j] = 1;
            g[i][j] = 1;
            zeronum++;
        }
    }
    VI zeroy(h, 0);
    VI zerox(w, 0);

    rep(i, h)
    {
        rep(j, w)
        {
            zeroy[i] += zero[i][j];
            zerox[j] += zero[i][j];
        }
    }
    VVL multi(h + 1, VL(w + 1, 1));

    rep(i, h) multi[i + 1][1] = (multi[i][1] * g[i][0]) % MOD;
    rep(i, w) multi[1][i + 1] = (multi[1][i] * g[0][i]) % MOD;
    rep(i, 1, h)
    {
        rep(j, 1, w)
        {
            //multi[i + 1][j + 1] = mpow(multi[i][j + 1] * multi[i + 1][j], multi[i][j] - 2) * g[i][j];
            multi[i + 1][j + 1] = (((((multi[i][j + 1] * multi[i + 1][j]) % MOD) * inv(multi[i][j])) % MOD) * g[i][j]) % MOD;

            multi[i + 1][j + 1] %= MOD;
        }
    }

    ll total = multi[h][w];
    int q;
    cin >> q;
    rep(i, q)
    {
        int x, y;
        cin >> y >> x;
        int z = zeronum;
        z -= zeroy[y - 1];
        z -= zerox[x - 1];

        z += zero[y - 1][x - 1];
        //cerr << z << endl;
        if (z != 0)
        {
            cout << 0 << endl;
            continue;
        }

        ll tmp = multi[y][w] * inv(multi[y - 1][w]);
        tmp %= MOD;

        tmp *= (multi[h][x] * inv(multi[h][x - 1])) % MOD;
        tmp %= MOD;

        tmp *= inv(g[y - 1][x - 1]);
        tmp %= MOD;

        cout << (total * inv(tmp)) % MOD << endl;
    }

    return 0;
}
0