結果
問題 | No.1683 Robot Guidance |
ユーザー | tktk_snsn |
提出日時 | 2021-09-17 22:57:53 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 254 ms / 2,000 ms |
コード長 | 2,206 bytes |
コンパイル時間 | 431 ms |
コンパイル使用メモリ | 86,692 KB |
実行使用メモリ | 108,260 KB |
最終ジャッジ日時 | 2023-09-12 08:59:34 |
合計ジャッジ時間 | 8,997 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 170 ms
106,852 KB |
testcase_01 | AC | 170 ms
106,916 KB |
testcase_02 | AC | 171 ms
107,112 KB |
testcase_03 | AC | 186 ms
108,116 KB |
testcase_04 | AC | 183 ms
107,952 KB |
testcase_05 | AC | 190 ms
108,140 KB |
testcase_06 | AC | 199 ms
108,016 KB |
testcase_07 | AC | 190 ms
107,884 KB |
testcase_08 | AC | 243 ms
108,136 KB |
testcase_09 | AC | 188 ms
107,976 KB |
testcase_10 | AC | 215 ms
107,848 KB |
testcase_11 | AC | 214 ms
108,260 KB |
testcase_12 | AC | 244 ms
107,844 KB |
testcase_13 | AC | 72 ms
71,108 KB |
testcase_14 | AC | 224 ms
107,952 KB |
testcase_15 | AC | 189 ms
108,144 KB |
testcase_16 | AC | 72 ms
71,336 KB |
testcase_17 | AC | 74 ms
71,124 KB |
testcase_18 | AC | 73 ms
71,420 KB |
testcase_19 | AC | 170 ms
106,912 KB |
testcase_20 | AC | 171 ms
107,100 KB |
testcase_21 | AC | 178 ms
107,108 KB |
testcase_22 | AC | 171 ms
107,064 KB |
testcase_23 | AC | 170 ms
106,956 KB |
testcase_24 | AC | 254 ms
108,064 KB |
testcase_25 | AC | 171 ms
106,928 KB |
testcase_26 | AC | 74 ms
71,244 KB |
testcase_27 | AC | 171 ms
106,912 KB |
testcase_28 | AC | 170 ms
106,756 KB |
testcase_29 | AC | 72 ms
71,444 KB |
testcase_30 | AC | 167 ms
106,872 KB |
testcase_31 | AC | 169 ms
106,852 KB |
testcase_32 | AC | 169 ms
106,812 KB |
testcase_33 | AC | 190 ms
107,736 KB |
testcase_34 | AC | 195 ms
107,572 KB |
testcase_35 | AC | 187 ms
107,472 KB |
testcase_36 | AC | 193 ms
107,500 KB |
testcase_37 | AC | 182 ms
107,584 KB |
testcase_38 | AC | 183 ms
108,020 KB |
testcase_39 | AC | 170 ms
106,988 KB |
testcase_40 | AC | 72 ms
71,200 KB |
ソースコード
mod = 10 ** 9 + 7 U = 2 * 10 ** 6 class Combination: """ SIZEが10**6程度以下の二項係数を何回も呼び出したいときに使う 使い方: comb = Combination(SIZE, MOD) comb(10, 3) => 120 """ def __init__(self, N, MOD=10 ** 9 + 7): self.MOD = MOD self.__make_factorial_list(N) def __call__(self, n, k): if k < 0 or k > n: return 0 res = self.fact[n] * self.inv[k] % self.MOD res = res * self.inv[n - k] % self.MOD return res def nPk(self, n, k): if k < 0 or k > n: return 0 return self.fact[n] * self.inv[n - k] % self.MOD def nHk(self, n, k): if k == 0: return 1 return self.__call__(n + k - 1, k) def __make_factorial_list(self, N): self.fact = [1] * (N + 1) self.inv = [1] * (N + 1) MOD = self.MOD for i in range(1, N + 1): self.fact[i] = (self.fact[i - 1] * i) % MOD self.inv[N] = pow(self.fact[N], MOD - 2, MOD) for i in range(N, 0, -1): self.inv[i - 1] = (self.inv[i] * i) % MOD return def main(): A, B, X, Y = map(int, input().split()) if (X + Y + A) % 2 == 1 or A < abs(X) + abs(Y): return 0 tmp = [B // 4] * 4 for i in range(B % 4): tmp[i] += 1 move = [1, 0, 0, 0] for i in range(4): move[(i + 1) % 4] += tmp[i] comb = Combination(U + 5, mod) ans = 0 absX = abs(X) absY = abs(Y) for i in range(absX, A + 1, 2): j = A - i if j < absY: break tmp = 1 ok = (i + absX) // 2 ng = (i - absX) // 2 if X < 0: ok, ng = ng, ok tmp = tmp * comb(ok + move[0] - 1, move[0] - 1) % mod if ng: tmp = tmp * comb(ng + move[2] - 1, move[2] - 1) % mod ok = (j + absY) // 2 ng = (j - absY) // 2 if Y < 0: ok, ng = ng, ok if ok: tmp = tmp * comb(ok + move[1] - 1, move[1] - 1) % mod if ng: tmp = tmp * comb(ng + move[3] - 1, move[3] - 1) % mod ans = (ans + tmp) % mod return ans print(main())