zoukankan      html  css  js  c++  java
  • BZOJ1211

    Portal

    Description

    一个有(n(nleq150))个结点的树,给出每个节点的度数(d_i),求满足度数要求的无根树的个数。

    Solution

    prufer数列_百度百科
    简单来说,一棵(n)个节点的无根树与一个数值在([1,n])的长度为(n-2)的序列一一对应,这个序列便称作prufer数列。其中若点(i)的度数为(d_i),那么(i)在prufer数列中就出现(d_i-1)次。
    那么这道题就很简单了:向一个长度为(n-2)的序列中填入(d_1-1)(1)(d_2-1)(2),...,(d_n-1)(n)

    [egin{align*} ans &= inom{n-2}{d_1-1} inom{n-2-(d_1-1)}{d_2-1} ...inom{n-2-sum_{i=1}^{n-1}(d_i-1)}{d_n-1} \ &= frac{(n-2)!}{(d_1-1)!(n-2-(d_1-1))!} cdot frac{(n-2-(d_1-1))!}{(d_2-1)!(n-2-sum_{i=1}^2(d_i-1))} ... frac{(n-2-sum_{i=1}^{n-1}(d_i-1))!}{(d_n-1)!0!} \ &= frac{(n-2)!}{prod_{i=1}^n (d_i-1)!} end{align*}$$由于阶乘比较大而且坑爹的不取模,所以通过分解质因数的方法来计算。 > 时间复杂度$O(n^2)$。 ##Code ```cpp //[HNOI2004]树的计数 #include <cstdio> int const N=200; int n,d[N]; int fac[N][N]; void init() { for(int i=1;i<=n;i++) { int x=i; for(int j=2;j<=n;j++) while(x%j==0) fac[i][j]++,x/=j; } for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) fac[i][j]+=fac[i-1][j]; } int ansP[N]; int main() { scanf("%d",&n); init(); int sum=0; for(int i=1;i<=n;i++) { scanf("%d",&d[i]); sum+=d[i]; if(d[i]==0&&n!=1) {puts("0"); return 0;} } if(sum!=n*2-2) {puts("0"); return 0;} if(n==1) {puts("1"); return 0;} for(int i=1;i<=n;i++) ansP[i]=fac[n-2][i]; for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) ansP[j]-=fac[d[i]-1][j]; long long ans=1; for(int i=1;i<=n;i++) for(int j=1;j<=ansP[i];j++) ans*=i; printf("%lld ",ans); return 0; } ``` ##P.S. 注意要特判掉无解和$n=1$的情况。]

  • 相关阅读:
    天下无贼的IT版经典台词(zz)
    matlab混合编程向导(vc,vb,.net...)
    mcc生成的代码加入VC,VC的一些设置
    请问:Win2k下面,怎样取消mediaplayer的预览?
    win2k密码恢复
    **上海铁路局2004年最新时刻发布!**
    VC与Matlab混合编程的快速实现(31)
    校车时刻表
    悼Java师兄(转载)
    制作弹出窗口常用技巧九则
  • 原文地址:https://www.cnblogs.com/VisJiao/p/BZOJ1211.html
Copyright © 2011-2022 走看看