洛谷P2290 [HNOI2004]树的计数

bzoj1211 [HNOI2004]树的计数

Description

一个有\(n\)个结点的树,设它的结点分别为\(v_1,v_2,\cdots, v_n\),已知第\(i\)个结点\(v_i\)的度数为\(d_i\)

问满足这样的条件的不同的树有多少棵。

Input

第一行是一个正整数\(n\),表示树有\(n\)个结点。第二行有\(n\)个数,第\(i\)个数表示\(d_i\),即树的第\(i\)个结点的度数。其中\(1\le n\le 150\),输入数据保证满足条件的树不超过\(10^{17}\)个。

Output

输出满足条件的树有多少棵。


可以说是个模板题

prufer 编码,这里只给结论,证明去这里看

我不想复制一遍了

注意一下以下说的“种”和“个”是不同的

由于 prufer 编码的一些性质,其实问题求的就是,总共\(n-2\)个元素,其中有\(n\)种不同元素,每种元素有\(d_i-1\)个,的排列数

所以,如果\(\sum_{i+1}^nd_i-1\neq n-2\)或者\(d_i=0\)则无解,不过要特判\(n=1\)的情况

如何证明在那个链接里都有

首先如果没有这个每种元素多少个的条件,那么\(n-2\)个元素的排列数就是\((n-2)!\)

然后要去重,因为对于每一种元素,这\(d_i-1\)个元素无论怎么变换顺序都是算一种,所以当然要除以\((d_i-1)!\)

那么答案就是:

\[\frac{(n-2)!}{\prod_{i=1}^n d_i-1}
\]

但是题目保证的是最后结果不大于\(10^{17}\),所以中间值可能会爆\(long long\)

那么我们采取分解质因数的方法,最后再乘起来

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<iomanip>
#include<cstring>
#define reg register
#define EN std::puts("")
#define LL long long
inline int read(){
register int x=0;register int y=1;
register char c=std::getchar();
while(c<'0'||c>'9'){if(c=='-') y=0;c=std::getchar();}
while(c>='0'&&c<='9'){x=x*10+(c^48);c=std::getchar();}
return y?x:-x;
}
int n;
int prime[155],notprime[155];
inline void get_prime(){
for(reg int i=2;i<=n;i++){
if(notprime[i]) continue;
prime[++prime[0]]=i;
for(reg int j=i+i;j<=n;j+=i) notprime[j]=1;
}
}
int num[155],d[155];
inline void mul(int x,int k){
for(reg int i=1;i<=prime[0];i++)
while(!(x%prime[i])) num[i]+=k,x/=prime[i];
}
int main(){
n=read();
if(n==1){
d[1]=read();
return std::puts(d[1]?"0":"1"),0;
}
get_prime();
int sum=0;
for(reg int i=1;i<=n;i++){
sum+=d[i]=read()-1;
if(d[i]==-1) return std::puts("0"),0;
}
if(sum!=n-2) return std::puts("0"),0;
for(reg int i=1;i<=n-2;i++) mul(i,1);
for(reg int i=1;i<=n;i++)
for(reg int j=1;j<=d[i];j++) mul(j,-1);
LL ans=1;
for(reg int i=1;i<=prime[0];i++)
for(reg int j=1;j<=num[i];j++) ans*=prime[i];
std::printf("%lld",ans);
return 0;
}

P2290 [HNOI2004]树的计数(bzoj1211)的更多相关文章

  1. Luogu P2290 [HNOI2004]树的计数 Prufer序列+组合数

    最近碰了$prufer$ 序列和组合数..于是老师留了一道题:P2624 [HNOI2008]明明的烦恼 qwq要用高精... 于是我们有了弱化版:P2290 [HNOI2004]树的计数(考一样的可 ...

  2. P2290 [HNOI2004]树的计数

    P2290 [HNOI2004]树的计数prufer序列模板题 #include <iostream> #include <cstdio> #include <queue ...

  3. 洛谷 P2290 [HNOI2004]树的计数

    题目描述 输入输出格式 输入格式: 输入文件第一行是一个正整数n,表示树有n个结点.第二行有n个数,第i个数表示di,即树的第i个结点的度数.其中1<=n<=150,输入数据保证满足条件的 ...

  4. LUOGU P2290 [HNOI2004]树的计数(组合数,prufer序)

    传送门 解题思路 \(prufer\)序,就是所有的不同的无根树,都可以转化为唯一的序列.做法就是每次从度数为\(1\)的点中选出一个字典序最小的,把这个点删掉,并把这个点相连的节点加入序列,直到只剩 ...

  5. bzoj1211: [HNOI2004]树的计数 prufer编码

    题目链接 bzoj1211: [HNOI2004]树的计数 题解 prufer序 可重排列计数 代码 #include<bits/stdc++.h> using namespace std ...

  6. BZOJ1211: [HNOI2004]树的计数

    1211: [HNOI2004]树的计数 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 1245  Solved: 383[Submit][Statu ...

  7. prufer BZOJ1211: [HNOI2004]树的计数

    以前做过几题..好久过去全忘了. 看来是要记一下... [prufer] n个点的无根树(点都是标号的,distinct)对应一个 长度n-2的数列 所以 n个点的无根树有n^(n-2)种 树 转 p ...

  8. bzoj1211: [HNOI2004]树的计数(prufer序列+组合数学)

    1211: [HNOI2004]树的计数 题目:传送门 题解: 今天刚学prufer序列,先打几道简单题 首先我们知道prufer序列和一颗无根树是一一对应的,那么对于任意一个节点,假设这个节点的度数 ...

  9. [bzo1211][HNOI2004]树的计数_prufer序列

    树的计数 bzoj-1211 HNOI-2004 题目大意:题目链接. 注释:略. 想法: prufer序列有一个性质就是一个数在prufer序列中出现的次数等于这个prufer序列生成的树中它的度数 ...

随机推荐

  1. MTK Android 设置下添加一级菜单[ZedielPcbTest]

    功能描述:Android7.1.2 设置下添加一级菜单[ZedielPcbTest],点击ZedielPcbTest,启动ZedielPcbTest.apk应用. 编译:需要在out目录删除Settt ...

  2. javascript中的constructor

    constructor,构造函数,对这个名字,我们都不陌生,constructor始终指向创建当前对象的构造函数. 这里有一点需要注意的是,每个函数都有一个prototype属性,这个prototyp ...

  3. YII基础二

    YII多表关联查询 ->select("{{%relation_detail}}.name") ->from("{{%user}}") ->l ...

  4. Jmeter常用元件

    1.测试计划:测试元件的容器,相当于一个项目名称 线程组:  2.监听器:负责收集测试结果,同时也被告知了结果显示的方式 (1)查看结果树:看具体某个请求——请求响应,结果明细 (2)聚合报告:汇总报 ...

  5. OO第三单元作业分析

    一.JML的理论基础应用工具链 JML是用于对Java程序进行规格化设计的一种表示语言.基于Larch方法构建. (1)注释 JML以javadoc注释的方式来表示规格,每行都以@起头.有两种注释方式 ...

  6. alg-链表中有环

    typedef struct ListNode { int val; ListNode *next; ListNode(int x) : val(x), next(nullptr) {} }ListN ...

  7. Python——flask漏洞探究

    python的用途是真的多,就连网站也能做,这个有点像Java的Servlet flask基础 hello world 我们先从基础的开始,在网页上打出hello world,python代码如下: ...

  8. 利用SSIS进行SharePoint 列表数据的ETL

    好几年前写了一篇<SSIS利用Microsoft Connector for Oracle by Attunity组件进行ETL!>,IT技术真是日新月异,这种方式对于新的SQL SERV ...

  9. HashMap主要方法源码分析(JDK1.8)

    本篇从HashMap的put.get.remove方法入手,分析源码流程 (不涉及红黑树的具体算法) jkd1.8中HashMap的结构为数组.链表.红黑树的形式     (未转化红黑树时)   (转 ...

  10. L - Neko does Maths CodeForces - 1152C 数论(gcd)

    题目大意:输入两个数 a,b,输出一个k使得lcm(a+k,b+k)尽可能的小,如果有多个K,输出最小的. 题解: 假设gcd(a+k,b+k)=z; 那么(a+k)%z=(b+k)%z=0. a%z ...