結果

問題 No.1683 Robot Guidance
ユーザー ktr216ktr216
提出日時 2021-09-17 23:15:06
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 2,494 bytes
コンパイル時間 1,498 ms
コンパイル使用メモリ 167,148 KB
実行使用メモリ 23,436 KB
最終ジャッジ日時 2023-09-12 09:17:22
合計ジャッジ時間 12,105 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 40 ms
18,448 KB
testcase_01 AC 40 ms
18,676 KB
testcase_02 AC 42 ms
18,512 KB
testcase_03 AC 175 ms
18,748 KB
testcase_04 AC 75 ms
18,488 KB
testcase_05 AC 228 ms
18,448 KB
testcase_06 AC 536 ms
18,600 KB
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 AC 1,079 ms
18,592 KB
testcase_11 WA -
testcase_12 AC 1,967 ms
18,604 KB
testcase_13 AC 40 ms
18,508 KB
testcase_14 AC 1,388 ms
18,624 KB
testcase_15 AC 207 ms
18,536 KB
testcase_16 AC 40 ms
18,636 KB
testcase_17 AC 40 ms
18,648 KB
testcase_18 AC 40 ms
18,688 KB
testcase_19 AC 40 ms
18,732 KB
testcase_20 AC 40 ms
18,644 KB
testcase_21 AC 41 ms
18,516 KB
testcase_22 AC 40 ms
18,540 KB
testcase_23 AC 40 ms
18,548 KB
testcase_24 TLE -
testcase_25 AC 39 ms
18,624 KB
testcase_26 AC 40 ms
18,624 KB
testcase_27 AC 39 ms
18,752 KB
testcase_28 AC 40 ms
18,608 KB
testcase_29 AC 39 ms
18,764 KB
testcase_30 AC 40 ms
18,752 KB
testcase_31 AC 40 ms
18,676 KB
testcase_32 AC 41 ms
18,748 KB
testcase_33 AC 40 ms
18,736 KB
testcase_34 AC 40 ms
18,672 KB
testcase_35 AC 41 ms
18,796 KB
testcase_36 AC 40 ms
18,816 KB
testcase_37 AC 40 ms
18,628 KB
testcase_38 WA -
testcase_39 AC 39 ms
18,680 KB
testcase_40 AC 40 ms
18,748 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define rep(i, l, r) for (int i = (l); i < (r); i++)
using namespace std;

typedef long long ll;

vector<ll> f(2000001, 1);
ll MOD = 1000000007;

ll mod_pow(ll x, ll n, ll mod) {
    ll res = 1;
    while (n > 0) {
        if (n & 1) res = res * x % mod;
        x = x * x % mod;
        n >>= 1;
    }
    return res;
}

ll comb(ll x, ll y) {
    ll ret = f[x];
    ret = ret * mod_pow(f[y], MOD - 2, MOD) % MOD;
    ret = ret * mod_pow(f[x - y], MOD - 2, MOD) % MOD;
    return ret;
}

int main() {
    ll A, B, X, Y;
    rep(i, 0, 2000000) f[i + 1] = f[i] * (i + 1) % MOD;
    cin >> A >> B >> X >> Y;
    if (A < abs(X) + abs(Y)) {
        cout << 0 << endl;
        return 0;
    }
    if (A % 2 != (X + Y) % 2) {
        cout << 0 << endl;
        return 0;
    }
    if (B == 0) {
        if (X == A) cout << 1 << endl;
        else cout << 0 << endl;
        return 0;
    }
    if (B == 1) {
        if (X >= 0 && Y >= 0 && X + Y == A) cout << 1 << endl;
        else cout << 0 << endl;
        return 0;
    }
    if (B == 2) {
        if (Y >= 0 && abs(X) + Y <= A) cout << 1 << endl;
        else cout << 0 << endl;
        return 0;
    }
    if (B == 3) {
        cout << A - abs(X) - abs(Y) + 1 << endl;
        return 0;
    }
    ll R = (B + 4) / 4, U = (B + 3) / 4, L = (B + 2) / 4, D = (B + 1) / 4;
    ll ans = 0, C = (A - abs(X) - abs(Y)) / 2;
    rep(i, 0, C + 1) {
        //cout << "i=" << i << endl;
        ll P = 1;
        if (X > 0) {
            P = P * comb(R - 1 + X + i, R - 1) % MOD;
            //cout << "ans=" << ans << endl;
            P = P * comb(L - 1 + i, L - 1) % MOD;
            //cout << "ans=" << ans << endl;
        } else {
            P = P * comb(R - 1 + i, R - 1) % MOD;
            //cout << "ans=" << ans << endl;
            P = P * comb(L - 1 - X + i, L - 1) % MOD;
            //cout << "ans=" << ans << endl;
        }
        if (Y > 0) {
            P = P * comb(U - 1 + Y + C - i, U - 1) % MOD;
            //cout << "ans=" << ans << endl;
            P = P * comb(D - 1 + C - i, D - 1) % MOD;
            //cout << "ans=" << ans << endl;
        } else {
            P = P * comb(U - 1 + C - i, U - 1) % MOD;
            //cout << "ans=" << ans << endl;
            P = P * comb(D - 1 - Y + C - i, D - 1) % MOD;
            //cout << "ans=" << ans << endl;
        }
        ans = (ans + P) % MOD;
    }
    cout << ans << endl;
    //cout << "R=" << R << " U=" << U << " L=" << L << " D=" << D << endl;
}
0