BZOJ 1005 明明的烦恼(prufer序列+高精度)

有一种东西叫树的prufer序列，一个树的与一个prufer序列是一一对应的关系。

设有m个度数确定的点，这些点的度为dee[i],那么每个点在prufer序列中出现了dee[i]-1次。

由排列组合可以推出公式。需要用高精度。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define M 1100
using namespace std;
typedef long long ll;
struct abcd{
    ll xx[400];
    int cnt;
    abcd(int x=0)
    {
        memset(xx,0,sizeof xx);
        xx[1]=x;
        cnt=1;
    }
    ll& operator [] (int x)
    {
        return xx[x];
    }
}ans(1);
abcd operator *= (abcd &x,abcd &y)
{
    int i,j;
    abcd z;
    for(i=1;i<=x.cnt;i++)
        for(j=1;j<=y.cnt;j++)
            z[i+j-1]+=x[i]*y[j],z[i+j]+=z[i+j-1]/100000000,z[i+j-1]%=100000000;
    z.cnt=x.cnt+y.cnt;
    if(!z[z.cnt])
        --z.cnt;
    x=z;
}
ostream& operator << (ostream& os,abcd &x)
{
    int i;
    printf("%lld",x[x.cnt]);
    for(i=x.cnt-1;i;i--)
        printf("%08lld",x[i]);
    return os;
}
int n,m,remain,cnt[M],stack[M],top;
void Decomposition(int x,int y)
{
    int i;
    for(i=2;i*i<=x;i++)
        while(x%i==0)
            cnt[i]+=y,x/=i;
    if(x^1)
        cnt[x]+=y;
}
void Quick_Power(int i,int y)
{
    abcd x(i);
    while(y)
    {
        if(y&1)ans*=x;
        x*=x;
        y>>=1;
    }
}
int main()
{
    int i,x;
    cin>>n;remain=n-2;
    for(i=1;i<=n;i++)
    {
        scanf("%d",&x);
        if(x==-1)
            ++m;
        else if(x>1)
            stack[++top]=x-1,remain-=x-1;
    }
    for(i=2;i<=n-2;i++)
        Decomposition(i,1);
    while(top)
    {
        for(i=2;i<=stack[top];i++)
            Decomposition(i,-1);
        stack[top--]=0;
    }
    for(i=2;i<=remain;i++)
        Decomposition(i,-1);
    Decomposition(m,remain);
    for(i=1;i<=n;i++)
        if(cnt[i])
            Quick_Power(i,cnt[i]);
    cout<<ans<<endl;
}