[BZOJ1005] [HNOI2008] 明明的烦恼 (prufer编码)

Description

　　自从明明学了树的结构,就对奇怪的树产生了兴趣......给出标号为1到N的点,以及某些点最终的度数,允许在
任意两点间连线,可产生多少棵度数满足要求的树?

Input

　　第一行为N(0 < N < = 1000),
接下来N行,第i+1行给出第i个节点的度数Di,如果对度数不要求,则输入-1

Output

　　一个整数,表示不同的满足要求的树的个数,无解输出0

Sample Input

3
1
-1
-1

Sample Output

2

HINT

　　两棵树分别为1-2-3;1-3-2

Source

Solution

　　懒，不想写大段的证明（况且我觉得我讲不懂），可以戳这，该大神讲的非常清楚

　　或许$prufer$的魅力就在于每一个地方可以放任意的数吧，这倒是解决了$BZOJ1430$的疑问

 #include <bits/stdc++.h>
 using namespace std;
 int d[];
 struct bigint
 {
     int a[], len;

     bigint()
     {
         memset(a, , ), len = ;
     }

     bigint operator* (const int &rhs) const
     {
         bigint ans;
         ans.len = len + ;
         for(int i = ; i <= len; ++i)
             ans.a[i] += a[i] * rhs;
         for(int i = ; i < ans.len; ++i)
             if(ans.a[i] > )
             {
                 ans.a[i + ] += ans.a[i] / ;
                 ans.a[i] %= ;
             }
         while(!ans.a[--ans.len]);
         return ans;
     }

     bigint operator/ (const int &rhs) const
     {
         bigint ans;
         ans = *this, ++ans.len;
         for(int i = ans.len; i; --i)
         {
             ans.a[i - ] += ans.a[i] % rhs * ;
             ans.a[i] /= rhs;
         }
         while(!ans.a[--ans.len]);
         return ans;
     }
 };

 int main()
 {
     int n, sum = , cnt = ;
     bigint ans;
     scanf("%d", &n);
     for(int i = ; i <= n; ++i)
     {
         scanf("%d", d + i);
         if(!d[i])
         {
             puts("");
             return ;
         }
         if(~d[i]) ++cnt, sum += d[i] - ;
     }
     if(sum >  * n - )
     {
         puts("");
         return ;
     }
     ans.a[] = ;
     for(int i = n -  - sum; i < n - ; ++i)
         ans = ans * i;
     for(int i = ; i <= n -  - sum; ++i)
         ans = ans * (n - cnt);
     for(int i = ; i <= n; ++i)
         for(int j = ; j <= d[i] - ; ++j)
             ans = ans / j;
     for(int i = ans.len; i; --i)
         printf("%d", ans.a[i]);
     puts("");
     return ;
 }