結果
問題 | 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 = 998244353fact = [1] * 220000def inv(x):return pow(x, MOD-2, MOD)for i in range(1, 220000):fact[i] = fact[i-1] * i % MODdef nCk(n, k):if k < 0 or n - k < 0:return 0return fact[n] * inv(fact[k] * fact[n-k]) % MODN, 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 0abab = nCk(N + 2, K // 2) if K % 2 == 0 else 0aabb = 0if 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) % MODaabb += nCk(N + 1, (K // 2) - 1) * nCk(2, 2) % MODif K % 2 == 1:aabb = nCk(N + 1, K // 2) * nCk(2, 1) % MODaaaa = 0if 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)