示例1:
输入
5 1000000007
2
3
4
5
107
输出
2
24
264
3240
736935633
题意:
一个森林的代价为内部每个节点度数的平方和。问所有带标号的n 个点的森林的代价和 。
代码:
#include <bits/stdc++.h>
using namespace std;
const int N = 5000+1;
int MOD,n,t;
int C[N][N],F[N][N];
int f[N],g[N],h[N],dp[N];
int main(){
scanf("%d %d",&t,&MOD);
C[0][0]=1;
for(int i=1;i<N;++i) C[0][i]=0,C[i][0]=1;
for(int i=1;i<N;++i)
for(int j=1;j<=i;++j) C[i][j]=(C[i-1][j]+C[i-1][j-1])%MOD;
F[0][0]=1;
for(int i=1;i<N;++i){
F[i][0]=1;
for(int j=1;j<N;++j) F[i][j]=1LL*F[i][j-1]*i%MOD;
}
h[2]=2;
for(int i=3;i<N;++i){
for(int j=0;j<i-1;++j){
h[i]=(h[i]+1LL*C[i-1][j]*j*j%MOD*F[i-1][i-j-2]%MOD*j)%MOD;
}
h[i]=(h[i]+(i-1)*(i-1))%MOD;
h[i]=1LL*h[i]*i%MOD;
// if(i<20) printf("%lld ",h[i]);
}
f[0]=1; g[0]=1; f[1]=1; g[1]=1;
for(int i=2;i<N;++i){
f[i]=F[i][i-2]; g[i]=f[i];
for(int j=1;j<i;++j)
g[i]=(g[i]+1LL*C[i-1][j-1]*f[j]%MOD*g[i-j])%MOD;
// if(i<20) printf("%lld ",[i]);
}
for(int i=2;i<N;++i){
dp[i]=h[i];
for(int j=1;j<i;++j){
dp[i]=(dp[i]+1LL*C[i-1][j-1]*h[j]%MOD*g[i-j])%MOD;
dp[i]=(dp[i]+1LL*C[i-1][j-1]*f[j]%MOD*dp[i-j]%MOD)%MOD;
}
}
while(t--){
scanf("%d",&n);
printf("%d\n",dp[n]);
}
}