結果
問題 | No.2951 Similar to Mex |
ユーザー |
|
提出日時 | 2024-10-21 23:11:37 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 364 ms / 2,000 ms |
コード長 | 777 bytes |
コンパイル時間 | 268 ms |
コンパイル使用メモリ | 82,408 KB |
実行使用メモリ | 97,908 KB |
最終ジャッジ日時 | 2024-10-25 20:51:34 |
合計ジャッジ時間 | 9,896 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 42 |
ソースコード
N,M,K=map(int,input().split())T=max(M,K)+1mod=998244353dp=[[0]*(T+1) for i in range(N+1)]dp[0][0]=1for i in range(1,T+1):ndp=[[0]*(T+1)for i in range(N+1)]for k in range(i):X=pow(i,max(0,min(i,K)-k),mod)for j in range(min(i,N+1)):if i<=M and j<N:ndp[j+1][k]=(ndp[j+1][k]+dp[j][k])%modndp[j][i]=(ndp[j][i]+dp[j][k]*X%mod)%moddp,ndp=ndp,dpS=[[0]*(N+1)for i in range(N+1)]S[0][0]=1for i in range(N):for j in range(i+1):S[i+1][j]=(S[i+1][j]+S[i][j]*j)%modS[i+1][j+1]=(S[i+1][j+1]+S[i][j])%modfac=[0]*(N+1)fac[0]=1for i in range(1,N+1):fac[i]=(i*fac[i-1])%modans=0for i in range(1,N+1):ans=(ans+dp[i][T]*S[N][i]%mod*fac[i]%mod)%modprint(ans)