Prufer序列/排列组合+高精度

　　窝不会告诉你我是先做了BZOJ1211然后才来做这题的>_>（为什么？因为我以前不会高精度呀……）

　　在A了BZOJ 1211和1089之后，蒟蒻终于有信心来写这道神题啦= =

　　嗯还是先说下做法吧~

　　……

　　还是出门左转去看黄学长的博客吧……我懒得写了……其实就是Prufer序列+高精度= =嗯就是之前说的那两道题的加和……

 /**************************************************************

     Problem: 1005

     User: Tunix

     Language: C++

     Result: Accepted

     Time:188 ms

     Memory:1304 kb

 ****************************************************************/

 //BZOJ 1005

 #include<cstdio>

 #include<cstring>

 #include<cstdlib>

 #include<iostream>

 #include<algorithm>

 #define rep(i,n) for(int i=0;i<n;++i)

 #define F(i,j,n) for(int i=j;i<=n;++i)

 #define D(i,j,n) for(int i=j;i>=n;--i)

 #define pb push_back

 using namespace std;

 typedef long long LL;

 inline int getint(){

     int r=,v=; char ch=getchar();

     for(;!isdigit(ch);ch=getchar()) if (ch=='-') r=-;

     for(; isdigit(ch);ch=getchar()) v=v*-''+ch;

     return r*v;

 }

 const int N=;

 /*******************template********************/

 struct bint{

     int l,v[];

     bint(){l=; memset(v,,sizeof v);}

     int& operator [] (int x){return v[x];}

 }ans;

 const int Limit=;

 void print(bint a){

     printf("%d",a[a.l]);

     D(i,a.l-,) printf("%03d",a[i]);

     puts("");

 }

 bint operator * (bint a,int p){

     int tmp=;

     F(i,,a.l){

         a[i]=a[i]*p+tmp;

         tmp=a[i]/Limit;

         a[i]%=Limit;

     }

     if (tmp) a[++a.l]=tmp;

     return a;

 }

 int n,a[N],b[N],prime[N],tot,m,cnt;

 bool vis[N];

 void ready(int n){

     F(i,,n){

         if (!vis[i]) prime[++tot]=i;

         F(j,,tot){

             if (i*prime[j]>n) break;

             vis[i*prime[j]]=;

             if (i%prime[j]==) break;

         }

     }

 }

 void add(int k,int v){

     F(j,,tot){

         int x=k;

         while (x){

             b[j]+=x/prime[j]*v;

             x/=prime[j];

         }

     }

 }

 int main(){

 #ifndef ONLINE_JUDGE

     freopen("1005.in","r",stdin);

     freopen("1005.out","w",stdout);

 #endif

     ready();

     n=getint();

     if (n==){

         a[]=getint();

         if (a[]<=) puts("");

         else puts("");

         return ;

     }

     F(i,,n){

         a[i]=getint();

         if (a[i]== || a[i]>n-) {puts(""); return ;}

     }

     F(i,,n){

         if (a[i]>){

             a[i]-=;

             m+=a[i];

         }

         else cnt++;

     }

     if (m>n-){ puts(""); return ;}

     int tmp=n-;

     F(i,,n){

         if (a[i]>){

             add(tmp,);

             add(a[i],-);

             add(tmp-a[i],-);

             tmp-=a[i];

         }

     }

 //  F(i,1,tot) printf("%d ",b[i]); puts("");

     ans[]=;

     F(i,,tot) F(j,,b[i]) ans=ans*prime[i];

     F(j,,tmp) ans=ans*cnt;

     print(ans);

     return ;

 }

1005: [HNOI2008]明明的烦恼

Time Limit: 1 Sec Memory Limit: 162 MB
Submit: 2964 Solved: 1182
[Submit][Status][Discuss]

Description

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

Input

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

Output

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

Sample Input

3
1
-1
-1

Sample Output

HINT

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

Source

[Submit][Status][Discuss]

明明

【BZOJ】【1005】【HNOI2008】明明的烦恼