CF742F，奇妙的计数题

感觉通过重心来限制同构很巧妙啊（可能是套路?或许我做题太少了）。
我一开始的想法是把叶子去掉，这样就变成了所有点度数都小于等于 d <script type="math/tex" id="MathJax-Element-70">d</script>，然后通过度数最大的点来限制同构，可能不太行？于是最后还是抄题解。不过这样可以减少点数，增加运行速度。
中途爆int，WA了一发，感觉自己好菜啊。

#include<cstdio> 
const int N=1005;
int f[N][15][N],n,d,i,j,mo,k,l,inv[N];
inline int C(int x,int y){
    if(x<y)return 0;
    register int ans=1,i=1;
    for(;i<=y;++i)ans=1ll*ans*inv[i]%mo*(x-i+1)%mo;return ans;
}
int main(){
    scanf("%d%d%d",&n,&d,&mo);
    if(n<=2)return puts("1"),0;
    for(i=2,inv[1]=1;i<=d;++i)inv[i]=1ll*(mo-mo/i)*inv[mo%i]%mo;
    for(i=0;i<=n;++i)f[1][0][i]=1;
    for(i=2;i<=n;++i)       
        for(j=1;j<i && j<=d;++j)
            for(k=1;k<=n;++k){
                f[i][j][k]=f[i][j][k-1];
                for(l=1;l*k<i && l<=j;++l)
                    f[i][j][k]=(f[i][j][k]+1ll*f[i-k*l][j-l][k-1]*
                        C(f[k][k>1?d-1:0][k-1]+l-1,l))%mo;
                //printf("f[%d][%d][%d]=%d\n",i,j,k,f[i][j][k]);
            }
    j=f[n][d][n>>1];
    if(!(n&1))j=(j+mo-C(f[n>>1][d-1][n>>1],2))%mo;
    return printf("%d\n",j),0;
}