bzoj1005: [HNOI2008]明明的烦恼（prufer+高精度）

1005: [HNOI2008]明明的烦恼

题目：传送门 

题解：

　　 毒瘤题啊天~

　　 其实思考的过程还是比较简单的。。。

　　 首先当然还是要了解好prufer序列的基本性质啦

　　 那么和1211大体一致，主要还是利用组合数学：

　　 首先我们把度数和-n记录为sum，那么根据prufer序列，序列的元素个数就是n-2

　　 那就是要在n-2个位置中选sum个，然后就是分别根据度数要求算每个元素在sum个位置中的方案，然后乘起来。最后还要乘上没有度数要求的元素的方案数就...ok啦

　　 思考两分钟...代码两小时...太菜啦！！！！

代码：

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#define qread(x) x=read()
using namespace std;
inline int read()
{
    int f=1,x=0;char ch;
    while(ch<'0' || ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0' && ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return f*x;
}
struct node
{
    int len,a[11000];
    node(){memset(a,0,sizeof(a));}
}no,n1;
void chengfa(int x)
{
    int i;
    for(i=1;i<=no.len;i++)no.a[i]=no.a[i]*x;
    for(i=1;i<=no.len;i++)
    {
        no.a[i+1]+=no.a[i]/10;
        no.a[i]%=10;
    }
    i=no.len;
    while(no.a[i+1]>0)
    {
        i++;
        no.a[i+1]+=no.a[i]/10;
        no.a[i]%=10;
    }
    no.len=i;
    while(no.a[no.len]==0 && no.len>1)no.len--;
}
bool pd(int x)
{
    if(x<2)return false;
    double t=sqrt(double(x+1));
    for(int i=2;i<=t;i++)
        if(x%i==0)
            return false;
    return true;
}
int n,d[1100],pr[1100],s[1100];
int main()
{
    scanf("%d",&n);int cnt=0,sum=0;
    if(n==1){qread(d[1]);if(d[1]){printf("0
");return 0;}else {printf("1
");return 0;}}
    for(int i=1;i<=n;i++)
    {
        qread(d[i]);
        if(d[i]==0){printf("0
");}
        if(d[i]==-1)cnt++;
        else d[i]-=1,sum+=d[i];
    }
    if(sum>n-2){printf("0
");return 0;}
    if(sum<n-2 && cnt==0){printf("0
");return 0;}
    int len=0;
    for(int i=2;i<=n;i++)if(pd(i))pr[++len]=i;
    for(int i=2;i<=n-2;i++)
    {
        int x=i;
        for(int j=1;j<=len;j++)
            while(x%pr[j]==0 && x!=0)
                s[j]++,x/=pr[j];
    }
    for(int i=2;i<=n-2-sum;i++)
    {
        int x=i;
        for(int j=1;j<=len;j++)
            while(x%pr[j]==0 && x!=0)
                s[j]--,x/=pr[j];
    }
    for(int i=1;i<=n;i++)
        if(d[i]>0)
        {
            for(int k=2;k<=d[i];k++)
            {
                int x=k;
                for(int j=1;j<=len;j++)
                    while(x%pr[j]==0 && x!=0)
                        s[j]--,x/=pr[j];
            }
        }
    no.a[1]=1;no.len=1;
    for(int i=1;i<=len;i++)
        while(s[i]--)
            chengfa(pr[i]);
    for(int i=1;i<=n-2-sum;i++)chengfa(cnt);
    for(int i=no.len;i>=1;i--)
        printf("%d",no.a[i]);
    printf("
");
    return 0;
}