【bzoj1005】 HNOI2008—明明的烦恼

http://www.lydsy.com/JudgeOnline/problem.php?id=1005 (题目链接)

题意

　　给出标号为1到N的点,以及某些点最终的度数，允许在任意两点间连线，可产生多少棵度数满足要求的树？

Solution 

　　prufer编码。关于prufer，有一篇博客写得很好，还运用组合数学求出了公式，可惜代码是java。这道题还要写高精度，高精度除法太麻烦了，因为组合数一定是整数，所以我们选择将n！分解质因数，最后相加减就可以了。如何分解质因数呢，我当然不会去暴力分解，一次考试中就是暴力分解TLE了= =。 

　　若P为质数，则n!可分解为 一个数×P^x。其中 
　　 

　　并且P^t小于n，所以就很好做了。

代码

// bzoj1005
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<map>
#define inf 2147483640
#define LL long long
#define free(a) freopen(a".in","r",stdin);freopen(a".out","w",stdout);
using namespace std;
inline LL getint() {
    LL x=0,f=1;char ch=getchar();
    while (ch>'9' || ch<'0') {if (ch=='-') f=-1;ch=getchar();}
    while (ch>='0' && ch<='9') {x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

const int maxn=10010;
int ans[maxn*100],np[maxn],a[maxn],d[maxn];
int l=1,n,sum;

void pls(int x,int f) {
    for (int i=2;i<=x;i++)
        if (!np[i]) for (int j=i;j<=x;j*=i) a[i]+=f*x/j;
}
void mul(int x) {
    for (int i=1;i<=l;i++) ans[i]*=x;
    for (int i=1;i<=l;i++) {
        ans[i+1]+=ans[i]/10;
        ans[i]%=10;
    }
    while (ans[l+1]>0) {l++;ans[1+l]+=ans[l]/10;ans[l]%=10;}
}
int main() {
    for (int i=2;i<=1000;i++) if (!np[i]) {
            int x=i+i;
            while (x<=1000) np[x]=1,x+=i;
        }
    scanf("%d",&n);
    if (n==1) {
        int x;scanf("%d",&x);
        if (!x) printf("1");
        else printf("0");
        return 0;
    }
    int cnt=0;ans[1]=1;
    for (int i=1;i<=n;i++) {
        scanf("%d",&d[i]);
        if (!d[i]) {printf("0");return 0;}
        if (d[i]>0) {d[i]--;sum+=d[i];}
        else cnt++;
    }
    if (sum>n-2) {printf("0");return 0;}
    pls(n-2,1);
    pls(n-2-sum,-1);
    for (int i=1;i<=n;i++) if (d[i]) pls(d[i],-1);
    for (int i=2;i<=1000;i++)
        for (int j=1;j<=a[i];j++) mul(i);
    for (int i=1;i<=n-2-sum;i++) mul(cnt);
    for (int i=l;i>=1;i--) printf("%d",ans[i]);
    return 0;
}