結果

問題 No.1667 Forest
ユーザー 👑 Kazun
提出日時 2021-09-04 16:28:36
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 101 ms / 3,000 ms
コード長 670 bytes
コンパイル時間 317 ms
コンパイル使用メモリ 82,188 KB
実行使用メモリ 68,824 KB
最終ジャッジ日時 2024-12-21 05:27:31
合計ジャッジ時間 2,544 ms
ジャッジサーバーID
(参考情報)
judge2 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 15
権限があれば一括ダウンロードができます

ソースコード

diff #

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):
    B=Binom[n]
    C=Binom[n-1]
    for r in range(n+1):
        if r==0 or r==n:
            B[r]=1
        else:
            B[r]=(C[r]+C[r-1])%Mod

DP=[[0]*N for _ in range(N+1)]
DP[0][0]=1

for n in range(N):
    B=Binom[N-n-1]
    D=DP[n]
    for k in range(1,N-n+1):
        alpha=(B[k-1]*X[k])%Mod
        E=DP[n+k]
        for m in range(n+1):
            E[m+k-1]+=(D[m]*alpha)%Mod
            E[m+k-1]%=Mod

print(*DP[N],sep="\n")
0