結果
問題 | No.2468 Mercurialist |
ユーザー | 👑 rin204 |
提出日時 | 2023-08-13 01:40:27 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 148 ms / 2,000 ms |
コード長 | 1,930 bytes |
コンパイル時間 | 370 ms |
コンパイル使用メモリ | 82,520 KB |
実行使用メモリ | 72,348 KB |
最終ジャッジ日時 | 2024-04-30 21:04:38 |
合計ジャッジ時間 | 3,853 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 41 ms
53,684 KB |
testcase_01 | AC | 39 ms
52,588 KB |
testcase_02 | AC | 39 ms
52,680 KB |
testcase_03 | AC | 131 ms
71,280 KB |
testcase_04 | AC | 114 ms
69,440 KB |
testcase_05 | AC | 148 ms
70,984 KB |
testcase_06 | AC | 63 ms
66,932 KB |
testcase_07 | AC | 63 ms
66,628 KB |
testcase_08 | AC | 109 ms
68,408 KB |
testcase_09 | AC | 50 ms
61,660 KB |
testcase_10 | AC | 95 ms
63,400 KB |
testcase_11 | AC | 58 ms
64,440 KB |
testcase_12 | AC | 38 ms
53,420 KB |
testcase_13 | AC | 39 ms
53,420 KB |
testcase_14 | AC | 109 ms
71,240 KB |
testcase_15 | AC | 63 ms
66,888 KB |
testcase_16 | AC | 60 ms
66,164 KB |
testcase_17 | AC | 59 ms
66,020 KB |
testcase_18 | AC | 84 ms
69,288 KB |
testcase_19 | AC | 84 ms
70,488 KB |
testcase_20 | AC | 101 ms
70,284 KB |
testcase_21 | AC | 115 ms
72,348 KB |
testcase_22 | AC | 52 ms
63,208 KB |
testcase_23 | AC | 73 ms
70,320 KB |
testcase_24 | AC | 81 ms
69,780 KB |
testcase_25 | AC | 78 ms
71,672 KB |
testcase_26 | AC | 73 ms
67,760 KB |
testcase_27 | AC | 74 ms
68,984 KB |
testcase_28 | AC | 114 ms
72,164 KB |
testcase_29 | AC | 71 ms
67,636 KB |
testcase_30 | AC | 68 ms
68,116 KB |
testcase_31 | AC | 101 ms
68,100 KB |
testcase_32 | AC | 78 ms
67,156 KB |
testcase_33 | AC | 93 ms
71,096 KB |
ソースコード
class Combination: def __init__(self, n, MOD=998244353): self.fact = [1] * (n + 1) self.invfact = [1] * (n + 1) self.MOD = MOD for i in range(1, n + 1): self.fact[i] = self.fact[i - 1] * i % MOD self.invfact[n] = pow(self.fact[n], MOD - 2, MOD) for i in range(n - 1, -1, -1): self.invfact[i] = self.invfact[i + 1] * (i + 1) % MOD def extend(self, n): le = len(self.fact) if n < le: return self.fact.extend([1] * (n - le + 1)) self.invfact.extend([1] * (n - le + 1)) for i in range(le, n + 1): self.fact[i] = self.fact[i - 1] * i % self.MOD self.invfact[n] = pow(self.fact[n], self.MOD - 2, self.MOD) for i in range(n - 1, le - 1, -1): self.invfact[i] = self.invfact[i + 1] * (i + 1) % self.MOD def nPk(self, n, k): if k < 0 or n < k: return 0 if n >= len(self.fact): self.extend(n) return self.fact[n] * self.invfact[n - k] % self.MOD def nCk(self, n, k): if k < 0 or n < k: return 0 if n >= len(self.fact): self.extend(n) return (self.fact[n] * self.invfact[n - k] % self.MOD) * self.invfact[k] % self.MOD def nHk(self, n, k): if n == 0 and k == 0: return 1 return self.nCk(n + k - 1, k) def Catalan(self, n): return (self.nCk(2 * n, n) - self.nCk(2 * n, n - 1)) % self.MOD MOD = 998244353 x, y, z, k = map(int, input().split()) Comb = Combination(x + y + z) ans = 0 n = x + y + z for i in range(1, y + z + 2): mi = min(y + z - i + 1 + k, y + z) mi2 = min(y, y + z - mi) tmp = pow(mi, mi2, MOD) * Comb.fact[y + z - mi2] % MOD tmp *= Comb.nCk(n - i, y + z - (i - 1)) ans += tmp % MOD ans %= MOD ans *= pow(Comb.nCk(n, x) * Comb.fact[y + z] % MOD, MOD - 2, MOD) print(ans % MOD)