[bzoj1005]明明的烦恼

根据purfer序列的原理，每一个purfer序列都一一对应了一棵树，每一个点在purfer序列中出现的次数就是它的度数，那么直接用组合数去计算即可，注意要加高精度

#include<cstdio>
#include<cstring>
using namespace std;
struct ji{
    int len,a[10001];
}a;
int n,m,k,d[1001],vis[1001];
void cheng(int b){
    a.a[1]*=b;
    for(int i=2;i<=a.len;i++){
        a.a[i]=a.a[i]*b+a.a[i-1]/10;
        a.a[i-1]%=10;
    }
    while (a.a[a.len]>9){
        a.a[a.len+1]=a.a[a.len]/10;
        a.a[a.len++]%=10;
    }
}
void calc(int k,int p){
    for(int i=2;i*i<=k;i++)
        while (k%i==0){
            vis[i]+=p;
            k/=i;
        }
    if (k>1)vis[k]+=p;
}
void c(int a,int b){
    for(int i=b+1;i<=a;i++)calc(i,1);
    for(int i=1;i<=a-b;i++)calc(i,-1);
}
int main(){
    scanf("%d",&n);
    m=n-2;
    k=n;
    for(int i=1;i<=n;i++)scanf("%d",&d[i]);
    a.len=a.a[1]=1;
    for(int i=1;i<=n;i++)
        if (d[i]!=-1){
            c(m,d[i]-1);
            m-=d[i]-1;
            k--;
        }
    for(int i=1;i<=m;i++)calc(k,1);
    for(int i=1;i<=1000;i++)
        for(int j=1;j<=vis[i];j++)cheng(i);
    for(int i=a.len;i;i--)printf("%d",a.a[i]);
}