結果

問題 No.1683 Robot Guidance
ユーザー KudeKude
提出日時 2021-09-17 22:09:06
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 116 ms / 2,000 ms
コード長 1,009 bytes
コンパイル時間 295 ms
コンパイル使用メモリ 81,936 KB
実行使用メモリ 96,908 KB
最終ジャッジ日時 2024-06-29 20:43:14
合計ジャッジ時間 5,127 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 85 ms
90,160 KB
testcase_01 AC 81 ms
90,496 KB
testcase_02 AC 88 ms
92,148 KB
testcase_03 AC 99 ms
96,580 KB
testcase_04 AC 94 ms
95,464 KB
testcase_05 AC 97 ms
95,340 KB
testcase_06 AC 100 ms
96,716 KB
testcase_07 AC 96 ms
96,204 KB
testcase_08 AC 115 ms
96,464 KB
testcase_09 AC 98 ms
95,984 KB
testcase_10 AC 104 ms
95,736 KB
testcase_11 AC 103 ms
96,252 KB
testcase_12 AC 116 ms
96,500 KB
testcase_13 AC 86 ms
91,912 KB
testcase_14 AC 106 ms
95,884 KB
testcase_15 AC 102 ms
96,852 KB
testcase_16 AC 86 ms
91,196 KB
testcase_17 AC 85 ms
91,836 KB
testcase_18 AC 84 ms
90,872 KB
testcase_19 AC 87 ms
92,016 KB
testcase_20 AC 87 ms
91,820 KB
testcase_21 AC 89 ms
93,532 KB
testcase_22 AC 84 ms
89,936 KB
testcase_23 AC 88 ms
90,936 KB
testcase_24 AC 113 ms
95,320 KB
testcase_25 AC 88 ms
91,652 KB
testcase_26 AC 86 ms
91,172 KB
testcase_27 AC 90 ms
91,832 KB
testcase_28 AC 88 ms
91,176 KB
testcase_29 AC 87 ms
91,164 KB
testcase_30 AC 87 ms
92,492 KB
testcase_31 AC 90 ms
92,612 KB
testcase_32 AC 89 ms
92,172 KB
testcase_33 AC 102 ms
96,808 KB
testcase_34 AC 103 ms
96,792 KB
testcase_35 AC 99 ms
95,416 KB
testcase_36 AC 99 ms
96,908 KB
testcase_37 AC 97 ms
96,792 KB
testcase_38 AC 97 ms
95,708 KB
testcase_39 AC 89 ms
93,472 KB
testcase_40 AC 83 ms
90,848 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

MOD = 10 ** 9 + 7

fact_range = 2 * 10**6 + 100
fact = [1] * (fact_range + 1)
for i in range(0, fact_range):
    fact[i+1] = fact[i] * (i + 1) % MOD

ifact = [1] * (fact_range + 1)
ifact[fact_range] = pow(fact[fact_range], MOD - 2, MOD)
for i in range(fact_range, 0, -1):
    ifact[i-1] = ifact[i] * i % MOD

def comb(n, k):
    if k == 0:
        return 1
    if k < 0 or n < k:
        return 0
    else:
        return fact[n] * ifact[n-k] % MOD * ifact[k] % MOD


g, t, x, y = map(int, input().split())
ans = 0
t += 1
cnt_r = t // 4 + (t % 4 >= 1)
cnt_u = t // 4 + (t % 4 >= 2)
cnt_l = t // 4 + (t % 4 >= 3)
cnt_d = t // 4
for k1 in range(g + 1):
    k2 = g - k1
    if (x + k1) % 2 or (y + k2) % 2 or k1 < abs(x) or k2 < abs(y):
        continue
    x_r = (x + k1) // 2
    x_l = (k1 - x) // 2
    y_u = (y + k2) // 2
    y_d = (k2 - y) // 2
    ans += comb(x_r + cnt_r - 1, x_r) * comb(x_l + cnt_l - 1, x_l) % MOD * comb(y_u + cnt_u - 1, y_u) % MOD * comb(y_d + cnt_d - 1, y_d)
    ans %= MOD
print(ans)
0