結果
問題 | No.2383 Naphthol |
ユーザー | かつっぱ競プロ |
提出日時 | 2023-07-14 22:44:34 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 113 ms / 2,000 ms |
コード長 | 941 bytes |
コンパイル時間 | 85 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 19,456 KB |
最終ジャッジ日時 | 2024-09-16 07:47:55 |
合計ジャッジ時間 | 3,180 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge6 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 18 |
ソースコード
MOD = 998244353 fact = [1] * 220000 def inv(x): return pow(x, MOD-2, MOD) for i in range(1, 220000): fact[i] = fact[i-1] * i % MOD def nCk(n, k): if k < 0 or n - k < 0:return 0 return fact[n] * inv(fact[k] * fact[n-k]) % MOD N, K = map(int, input().split()) if N == 1: print([1,1,3,3,3,1,1][K]) exit(0) abcd = nCk(2 * N + 4, K) abba = nCk(N + 2, K // 2) if K % 2 == 0 else 0 abab = nCk(N + 2, K // 2) if K % 2 == 0 else 0 aabb = 0 if N % 2 == 0 and K % 2 == 0: aabb = nCk(N + 2, K // 2) if N % 2 == 1: if K % 2 == 0: aabb = nCk(N + 1, K // 2) * nCk(2, 0) % MOD aabb += nCk(N + 1, (K // 2) - 1) * nCk(2, 2) % MOD if K % 2 == 1: aabb = nCk(N + 1, K // 2) * nCk(2, 1) % MOD aaaa = 0 if N % 2 == 0 and K % 4 == 0: aaaa = nCk(N // 2 + 1, K // 4) if N % 2 == 1 and K % 2 == 0: aaaa = nCk((N - 1) // 2 + 1, K // 4) print((abcd + aabb + abab + abba) * inv(4) % MOD)