結果

問題 No.1683 Robot Guidance
ユーザー kai4096_donkai4096_don
提出日時 2021-09-17 22:34:28
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 747 ms / 2,000 ms
コード長 2,210 bytes
コンパイル時間 3,913 ms
コンパイル使用メモリ 227,668 KB
実行使用メモリ 34,468 KB
最終ジャッジ日時 2023-09-12 08:35:13
合計ジャッジ時間 34,003 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge15
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 667 ms
34,180 KB
testcase_01 AC 669 ms
34,180 KB
testcase_02 AC 667 ms
34,108 KB
testcase_03 AC 671 ms
34,180 KB
testcase_04 AC 669 ms
34,360 KB
testcase_05 AC 673 ms
34,184 KB
testcase_06 AC 684 ms
34,212 KB
testcase_07 AC 677 ms
34,124 KB
testcase_08 AC 736 ms
34,468 KB
testcase_09 AC 676 ms
34,228 KB
testcase_10 AC 703 ms
34,284 KB
testcase_11 AC 701 ms
34,264 KB
testcase_12 AC 735 ms
34,176 KB
testcase_13 AC 668 ms
34,296 KB
testcase_14 AC 714 ms
34,264 KB
testcase_15 AC 675 ms
34,172 KB
testcase_16 AC 668 ms
34,368 KB
testcase_17 AC 667 ms
34,368 KB
testcase_18 AC 667 ms
34,364 KB
testcase_19 AC 668 ms
34,328 KB
testcase_20 AC 667 ms
34,224 KB
testcase_21 AC 667 ms
34,120 KB
testcase_22 AC 668 ms
34,296 KB
testcase_23 AC 668 ms
34,288 KB
testcase_24 AC 747 ms
34,288 KB
testcase_25 AC 672 ms
34,180 KB
testcase_26 AC 670 ms
34,352 KB
testcase_27 AC 669 ms
34,344 KB
testcase_28 AC 670 ms
34,296 KB
testcase_29 AC 670 ms
34,328 KB
testcase_30 AC 669 ms
34,260 KB
testcase_31 AC 669 ms
34,180 KB
testcase_32 AC 669 ms
34,232 KB
testcase_33 AC 670 ms
34,376 KB
testcase_34 AC 669 ms
34,372 KB
testcase_35 AC 670 ms
34,292 KB
testcase_36 AC 670 ms
34,296 KB
testcase_37 AC 669 ms
34,180 KB
testcase_38 AC 670 ms
34,132 KB
testcase_39 AC 668 ms
34,112 KB
testcase_40 AC 667 ms
34,180 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
const long nPrime = 1000000007;
//const long nPrime = 998244353;
typedef long long ll;
long Inverse(const long &n){
    if (n % nPrime == 0){
        return 0;
    }
    long a = n, b = nPrime, u = 1, v = 0;
    while (b > 0) {
        long t = a / b;
        a -= t * b;
        swap(a, b);
        u -= t * v;
        swap(u, v);
    }
    u %= nPrime;
    if (u < 0){
        u += nPrime;
    }
    return u;    
}
int main() {
    vector<long> viTable(2000001,1), viInv(2000001,1);
    for(long i = 1; i <= 2000000; i++){
        viTable[i] = (viTable[i-1] * i) % nPrime;
        viInv[i] = Inverse(viTable[i]);
    }
    
    long a,b,x,y;
    cin >> a >> b >> x >> y;
    if(((a+x+y)&1) == 1){
        cout << 0 << endl;
        return 0;
    }
    
    long i1 = (b+4)/4;
    long i2 = (b+2)/4;
    long i3 = (b+3)/4;
    long i4 = (b+1)/4;
    
    long nAns = 0;
    for(long s = 0; s <= a; s++){
        long t = s-x;
        long u = (a-2*s+x+y)/2;
        long v = u-y;
        if(t < 0 || u < 0 || v < 0){
            continue;
        }
        if((i2 == 0 && t > 0) ||(i3 == 0 && u > 0) ||(i4 == 0 && v > 0) ){
            continue;
        }
        
        long iThis = 1;
        iThis *= viTable[s+i1-1];
        iThis %= nPrime;
        iThis *= viInv[s];
        iThis %= nPrime;
        iThis *= viInv[i1-1];
        iThis %= nPrime;
        if(i2 > 0){
            iThis *= viTable[t+i2-1];
            iThis %= nPrime;
            iThis *= viInv[t];
            iThis %= nPrime;
            iThis *= viInv[i2-1];
            iThis %= nPrime;
        }
        if(i3 > 0){
            iThis *= viTable[u+i3-1];
            iThis %= nPrime;
            iThis *= viInv[u];
            iThis %= nPrime;
            iThis *= viInv[i3-1];
            iThis %= nPrime;
        }
        if(i4 > 0){
            iThis *= viTable[v+i4-1];
            iThis %= nPrime;
            iThis *= viInv[v];
            iThis %= nPrime;
            iThis *= viInv[i4-1];
            iThis %= nPrime;
        }
        nAns += iThis;
        nAns %= nPrime;
    }
    
    cout << nAns << endl;
    return 0;
}
0