牛客（多校7）：Valuable Forests

示例1：
输入

5 1000000007
2
3
4
5
107

输出

2
24
264
3240
736935633

题意：
一个森林的代价为内部每个节点度数的平方和。问所有带标号的n 个点的森林的代价和 。
代码：

#include 
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]);
    }
}