結果

問題 No.2381 Gift Exchange Party
ユーザー ntuda
提出日時 2023-07-16 15:05:36
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 63 ms / 2,000 ms
コード長 1,068 bytes
コンパイル時間 159 ms
コンパイル使用メモリ 82,424 KB
実行使用メモリ 74,520 KB
最終ジャッジ日時 2024-09-17 12:43:47
合計ジャッジ時間 2,180 ms
ジャッジサーバーID
(参考情報)
judge4 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 22
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

N, P = map(int, input().split())
MOD = 998244353
fact = [1] * (N + 1)
rfact = [1] * (N + 1)
r = 1
for i in range(1, N + 1):
fact[i] = r = r * i % MOD
rfact[N] = r = pow(fact[N], MOD - 2, MOD)
for i in range(N, 0, -1):
rfact[i - 1] = r = r * i % MOD
def comb(n, k):
return fact[n] * rfact[k] * rfact[n - k] % MOD
if N < P:
print(fact[N] - 1)
exit()
'''
PN!
1. A[i] == i (i = 0..N-1)
2. PP
N//P1N//P
A[i] == i
P (P-1)!
2
'''
ans = 0
tmp = 0
tmp1 = 1
npp = N // P
for i in range(N // P):
tmp1 *= comb(N - i * P, P) * fact[P-1]
tmp1 %= MOD
tmp += tmp1 * rfact[i+1]
tmp %= MOD
ans = fact[N] - tmp - 1
ans %= MOD
print(ans)
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0