結果

問題 No.1667 Forest
ユーザー kotatsugamekotatsugame
提出日時 2021-09-03 22:27:40
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 321 ms / 3,000 ms
コード長 717 bytes
コンパイル時間 645 ms
コンパイル使用メモリ 65,944 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-12-15 14:42:47
合計ジャッジ時間 3,625 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 15
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:10:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type]
   10 | main()
      | ^~~~

ソースコード

diff # 

#include<iostream>
using namespace std;
int N,mod;
long modpow(long a,long b){return b?modpow(a*a%mod,b/2)*(b%2?a:1)%mod:1;}
long T[301];
long dp[301][301];
long fac[301],invfac[301];
long ans[301];
long comb(int a,int b){return fac[a]*invfac[b]%mod*invfac[a-b]%mod;}
main()
{
	cin>>N>>mod;
	fac[0]=1;
	for(int i=1;i<=N;i++)fac[i]=fac[i-1]*i%mod;
	invfac[N]=modpow(fac[N],mod-2);
	for(int i=N;i--;)invfac[i]=invfac[i+1]*(i+1)%mod;
	T[1]=1;
	for(int k=2;k<=N;k++)T[k]=modpow(k,k-2);
	dp[0][0]=1;
	for(int c=0;c<N;c++)
	{
		for(int i=0;i<=N;i++)if(dp[c][i]>0)
		{
			for(int k=1;i+k<=N;k++)
			{
				(dp[c+k-1][i+k]+=dp[c][i]*comb(N-i-1,k-1)%mod*T[k]%mod)%=mod;
			}
		}
	}
	for(int i=0;i<N;i++)cout<<dp[i][N]<<endl;
}
0