题目描述
小 C 在自己家的花园里种了一棵苹果树, 树上每个结点都有恰好两个分支. 经过细心的观察, 小 C 发现每一天这棵树都会生长出一个新的结点.
第一天的时候, 果树会长出一个根结点, 以后每一天, 果树会随机选择一个当前树中没有长出过结点 的分支, 然后在这个分支上长出一个新结点, 新结点与分支所属的结点之间连接上一条边.
小 C 定义一棵果树的不便度为树上两两结点之间的距离之和, 两个结点之间 的距离定义为从一个点走到另一个点的路径经过的边数.
现在他非常好奇, 如果 NNN 天之后小 G 来他家摘苹果, 这个不便度的期望 EEE 是多少. 但是小 C 讨厌分数, 所以他只想知道E×N!E imes N !E×N! 对 PPP 取模的结果, 可以证明这是一个整数.
输入输出格式
输入格式：
从标准输入中读入数据. 一行两个整数 NNN (N<=2000), PPP .
输出格式：
输出到标准输出中. 输出一个整数表示答案．
输入输出样例
说明
以上是所有 N=3N = 3N=3 时可能的苹果树形态, 其中编号表示这个结点是第几天生 长出来的, 显然每种情况两两结点的距离均为 444 .
考虑n个节点二叉树方案为n!
考虑边的贡献。
设子树大小为sz,则在不考虑编号的情况下经过此边sz*(n-sz)次
方案数为子树内方案*子树外方案。
由于i节点子树内节点编号>i，所以方案为sz!*C(n-i,sz-1)。
子树外方案为i!*(i+1-2)*(i+2-2)*(i+3-2)*......*(n-sz+1-2)。
化简为
i*(i-1)*(n-sz-1)!
1 #include<iostream> 2 #include<cstring> 3 #include<cstdlib> 4 #include<cstdio> 5 #include<cmath> 6 #include<algorithm> 7 #define ll long long 8 using namespace std; 9 inline ll read() { 10 ll x=0,f=1;char ch=getchar(); 11 for(;!isdigit(ch);ch=getchar()) if(ch=='-') f=-1; 12 for(;isdigit(ch);ch=getchar()) x=x*10+ch-'0'; 13 return x*f; 14 } 15 ll n,p; 16 ll jc[2002],c[2002][2002]; 17 int main() { 18 n=read(),p=read(); 19 c[0][0]=1; 20 jc[0]=1; 21 for(int i=1;i<=n;i++) { 22 c[i][0]=1; 23 jc[i]=jc[i-1]*i;jc[i]%=p; 24 for(int j=1;j<=i;j++) c[i][j]=(c[i-1][j]+c[i-1][j-1])%p; 25 } 26 ll ans=0; 27 for(int i=2;i<=n;i++) { 28 for(int sz=1;sz<=n-i+1;sz++) { 29 ans+=sz*(n-sz)*i*(i-1)%p*c[n-i][sz-1]%p*jc[sz]%p*jc[n-sz-1]%p; 30 ans%=p; 31 } 32 } 33 printf("%lld ",ans); 34 }