結果
問題 | No.1667 Forest |
ユーザー | とりゐ |
提出日時 | 2022-11-29 14:35:29 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 322 ms / 3,000 ms |
コード長 | 763 bytes |
コンパイル時間 | 388 ms |
コンパイル使用メモリ | 81,792 KB |
実行使用メモリ | 76,132 KB |
最終ジャッジ日時 | 2024-10-07 01:53:23 |
合計ジャッジ時間 | 3,841 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 321 ms
75,776 KB |
testcase_01 | AC | 322 ms
75,904 KB |
testcase_02 | AC | 313 ms
75,904 KB |
testcase_03 | AC | 77 ms
71,040 KB |
testcase_04 | AC | 303 ms
76,132 KB |
testcase_05 | AC | 224 ms
74,624 KB |
testcase_06 | AC | 165 ms
73,344 KB |
testcase_07 | AC | 136 ms
72,576 KB |
testcase_08 | AC | 102 ms
71,808 KB |
testcase_09 | AC | 92 ms
71,296 KB |
testcase_10 | AC | 83 ms
71,168 KB |
testcase_11 | AC | 76 ms
70,656 KB |
testcase_12 | AC | 65 ms
65,664 KB |
testcase_13 | AC | 65 ms
65,664 KB |
testcase_14 | AC | 65 ms
66,048 KB |
testcase_15 | AC | 64 ms
65,792 KB |
testcase_16 | AC | 64 ms
65,536 KB |
testcase_17 | AC | 64 ms
66,048 KB |
ソースコード
n,mod=map(int,input().split()) M=(10**5)*3+1 fac=[1]*M ninv=[1]*M finv=[1]*M for i in range(2,M): fac[i]=fac[i-1]*i%mod ninv[i]=(-(mod//i)*ninv[mod%i])%mod finv[i]=finv[i-1]*ninv[i]%mod def binom(n,k): if n<0 or k<0: return 0 if k>n: return 0 return (fac[n]*finv[k]%mod)*finv[n-k]%mod def fpow(x,k): res=1 while k: if k&1: res=res*x%mod x=x*x%mod k>>=1 return res cnt=[0]*(n+1) for i in range(1,n+1): if i==1: cnt[i]=1 else: cnt[i]=fpow(i,i-2) dp=[[0]*(n+1) for i in range(n+1)] dp[0][0]=1 for i in range(n): for j in range(i+1): for k in range(1,n+1): if i+k<=n: dp[i+k][j+k-1]+=binom(n-i-1,k-1)*dp[i][j]%mod*cnt[k] dp[i+k][j+k-1]%=mod for i in range(n): print(dp[n][i])