結果

問題 No.1075 木の上の山
ユーザー kotatsugame
提出日時 2020-06-05 22:56:25
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 48 ms / 2,000 ms
コード長 1,261 bytes
コンパイル時間 789 ms
コンパイル使用メモリ 76,564 KB
実行使用メモリ 42,836 KB
最終ジャッジ日時 2024-12-17 17:17:25
合計ジャッジ時間 2,407 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 30
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:70:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type]
   70 | main()
      | ^~~~

ソースコード

diff # 

#include<iostream>
#include<vector>
using namespace std;
const long mod=1e9+7;
int N,K;
vector<int>G[1010];
long dp[1010][1010];
void dfs1(int u,int p)
{
	for(int i=0;i<K;i++)dp[u][i]=1;
	for(int v:G[u])if(v!=p)
	{
		dfs1(v,u);
		int t=0;
		for(int i=0;i<K;i++)
		{
			t+=dp[v][i];
			if(t>=mod)t-=mod;
			(dp[u][i]*=t)%=mod;
		}
	}
}
long ans;
void dfs2(int u,int p)
{
	vector<int>cs;
	for(int v:G[u])if(v!=p)cs.push_back(v);
	vector<vector<long> >L(cs.size()+1,vector<long>(K,1));
	vector<vector<long> >R(cs.size()+1,vector<long>(K,1));
	long t=p==-1;
	for(int i=0;i<K;i++)
	{
		(ans+=dp[u][i]*t)%=mod;
		if(p!=-1)
		{
			t+=dp[p][i];
			if(t>=mod)t-=mod;
		}
		L[0][i]=t;
	}
	for(int i=0;i<cs.size();i++)
	{
		long t=0;
		for(int j=0;j<K;j++)
		{
			t+=dp[cs[i]][j];
			if(t>=mod)t-=mod;
			L[i+1][j]=L[i][j]*t%mod;
		}
	}
	for(int i=cs.size();i--;)
	{
		long t=0;
		for(int j=0;j<K;j++)
		{
			t+=dp[cs[i]][j];
			if(t>=mod)t-=mod;
			R[i][j]=R[i+1][j]*t%mod;
		}
	}
	for(int i=0;i<cs.size();i++)
	{
		for(int j=0;j<K;j++)
		{
			dp[u][j]=L[i][j]*R[i+1][j]%mod;
		}
		dfs2(cs[i],u);
	}
}
main()
{
	cin>>N>>K;
	for(int i=1;i<N;i++)
	{
		int a,b;cin>>a>>b;
		a--,b--;
		G[a].push_back(b);
		G[b].push_back(a);
	}
	dfs1(0,-1);
	dfs2(0,-1);
	cout<<ans<<endl;
}
0