Description
Solution
\\(n\\) 个点的二叉树的方案数是 \\(n!\\)
证明十分显然:新加入的点占掉了 \\(1\\) 个位置,新加了 \\(2\\) 个位置,那么多出来一个位置,所以第 \\(i\\) 个点有 \\(i\\) 种放法
考虑每条边被经过的次数,设子树大小为 \\(size\\),就是 \\(size*(n-size)\\)
以此考虑每个点父边被经过的次数,枚举子树大小
然后贡献就是子树内部形态的方案数乘以外部形态的方案数
内部的显然就是 \\(size!\\) , 但是编号还不确定,所以是 \\(size!*C_{n-i}^{size-1}\\)
外部的我们先确定一个大小为 \\(i\\) 的树,再把多出来的 \\(n-size-i+1\\) 个拼上去,方案数为 \\(\\frac{(n-j-1)!}{(i-2)!}\\)
#include<bits/stdc++.h>
using namespace std;
const int N=2010;
int n,mod,c[N][N],Fac[N],w[N][N];
inline int F(int x,int y){
if(y<=0)return 1;
return w[x][y];
}
int main(){
freopen("pp.in","r",stdin);
freopen("pp.out","w",stdout);
cin>>n>>mod;
for(int i=0;i<=n;i++){
c[i][0]=1;
for(int j=1;j<=i;j++)c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod;
}
Fac[0]=1;
for(int i=1;i<=n;i++)Fac[i]=1ll*Fac[i-1]*i%mod;
for(int i=1;i<=n;i++){
w[i][0]=1;
for(int j=1;j<=n;j++)w[i][j]=1ll*w[i][j-1]*(i+j-2)%mod;
}
int ans=0;
for(int i=2;i<=n;i++)
for(int j=n-i+1;j>=1;j--)
ans=(ans+1ll*Fac[j]*c[n-i][j-1]%mod*Fac[i]%mod*F(i,n-j-i+1)%mod*(n-j)%mod*j)%mod;
cout<<ans<<endl;
return 0;
}