結果
| 問題 |
No.1667 Forest
|
| コンテスト | |
| ユーザー |
👑 Kazun
|
| 提出日時 | 2021-09-04 02:35:44 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 149 ms / 3,000 ms |
| コード長 | 914 bytes |
| コンパイル時間 | 198 ms |
| コンパイル使用メモリ | 82,392 KB |
| 実行使用メモリ | 69,128 KB |
| 最終ジャッジ日時 | 2024-12-16 00:46:52 |
| 合計ジャッジ時間 | 2,430 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 15 |
ソースコード
def convolution(a,b,length):
c=[0]*length
for i in range(len(a)):
for j in range(len(b)):
if i+j>=length:
break
c[i+j]+=a[i]*b[j]
c[i+j]%=Mod
#==================================================
N,Mod=map(int,input().split())
X=[]
for i in range(N+1):
if i==0:
X.append(0)
elif i==1:
X.append(1)
else:
X.append(pow(i,i-2,Mod))
Binom=[[0]*(N+1) for _ in range(N+1)]
for n in range(N+1):
for r in range(n+1):
if r==0 or r==n:
Binom[n][r]=1
else:
Binom[n][r]=(Binom[n-1][r]+Binom[n-1][r-1])%Mod
DP=[[0]*N for _ in range(N+1)]
DP[0][0]=1
for n in range(N):
for m in range(n+1):
for k in range(1,N-n+1):
alpha=(Binom[N-n-1][k-1]*X[k])%Mod
DP[n+k][m+k-1]+=(DP[n][m]*alpha)%Mod
DP[n+k][m+k-1]%=Mod
print(*DP[N],sep="\n")
Kazun