bzoj1005明明的烦恼(多少棵度数满足要求的树 Purfer+质因子分解高精计算）

题目
学习链接：
https://www.cnblogs.com/zhj5chengfeng/archive/2013/08/23/3278557.html
http://www.matrix67.com/blog/archives/682

一定是个整数。。

//难道数据范围有问题???1e3+5为什么不可以
#include<cstdio>
#include<iostream>
#include<cstring>
#define en '\n'
#define m(a,b) memset(a,b,sizeof a)
using namespace std;
typedef long long ll;
const int N=1e5+5;
int is[N],prime[N],tot;//质因子分解
int d[N],tmp[N],pri[N],cnt;//储存prufer需要的
int ans[N];//高精*低精 ans的空间需要>1e3;
void getprime()
{
    for(int i=2;i<=N-5-2;++i)
    {
        if(!is[i]) prime[++tot]=i;
        for(int j=1;j<=tot&&i*prime[j]<=N-5-2;++j)
        {
            is[i*prime[j]]=1;
            if(i%prime[j]==0) break;
        }
    }
}
inline void cheng(int x)
{
    for(int i=1;i<=ans[0];++i) ans[i]*=x;
    for(int i=1;i<=ans[0];++i) ans[i+1]+=ans[i]/10,ans[i]%=10;
    while(ans[ans[0]+1])
    {
        ++ans[0];
        ans[ans[0]+1]=ans[ans[0]]/10,ans[ans[0]]%=10;
    }
}
int main()
{
    getprime();
    int n;scanf("%d",&n);
    if(n==1)
    {
        scanf("%d",&d[1]);
        puts(d[1]<=0?"1":"0");
        return 0;
    }
    int sum=0;
    for(int i=1;i<=n;++i)
    {
        scanf("%d",&d[i]);
        if(!d[i]||d[i]==n) {puts("0");return 0;}//不可能存在无度的点
        if(d[i]!=-1) tmp[++cnt]=d[i]-1,sum+=tmp[cnt];
    }
    if(n-2<sum) {puts("0");return 0;}//要求过高 不得已实现
    for(int i=1;i<=cnt;++i)
    {
        for(int j=2;j<=tmp[i];++j)
        {
            int jj=j;
            for(int k=1;prime[k]<=jj;++k)
                while(jj%prime[k]==0) --pri[k],jj/=prime[k];
        }
    }
    for(int i=n-2-sum+1;i<=n-2;++i)
    {
        int jj=i;
        for(int k=1;prime[k]<=jj;++k)
            while(jj%prime[k]==0) ++pri[k],jj/=prime[k];
    }
    ans[0]=1,ans[1]=1;
    for(int i=1;i<=tot;++i)
    {
        while(pri[i])// 一定为正值!!!
        {
            --pri[i];
            cheng(prime[i]);
        }
    }
    for(int i=1;i<=n-2-sum;++i) cheng(n-cnt);
    for(int i=ans[0];i>=1;--i) printf("%d",ans[i]);
    puts("");
}