結果
問題 | No.1683 Robot Guidance |
ユーザー |
![]() |
提出日時 | 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 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 38 |
ソースコード
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)