【51NOD1806】wangyurzee的树(Prufer编码,容斥原理,组合计数)
题意:有n个点和m条限制,每条限制限制了一个点的度数不能为某个数。
求合法的树的个数模10^9+7
n<=10^6
m<=17
思路:WYZ作业
首先m<=17显然是2^m容斥
枚举每一个限制有用或没用,考虑某一个约束情况下的方案数
Caylay定理:n个点的生成树的个数=n^(n-2)
Prufer序列:一个长度为n-2的Prufer序列对应唯一的一棵n个节点的树,度数为a[i]的点在其中出现了(a[i]-1)次
考虑先在序列中填上所有的受约束条件的点,它们的方案数是C(n-2,a[1]-1)*C(n-2-(a[1]-1),a[2]-1)*……
化简得:
(n-2)!/((a[1]-1)!(a[2]-1)!……*(n-2-sigma(a[i]-1))!)
考虑没有约束的点可以随意在剩下的序列中出现,它们的方案数为:
(n-p)^(n-2-sigma(a[i]-1))
其中p为选中的约束条件个数
将两部分相乘,根据p的奇偶性进行容斥计算即可
几个hack点:
1.对于同一个点可能会有多个约束条件,需要判断某一个点是否被多次限制,如果被多次限制就不能计入答案
2.n=1需要特判ans=1(真的吗,这个点可是坑了我5点头盾)
const mo=;
var flag:array[..]of longint;
fac,exf:array[..]of int64;
a,b,c:array[..]of longint;
n,m,i:longint;
sun:int64; function mult(x,y:longint):int64;
var tmp:int64;
begin
tmp:=x; mult:=;
while y> do
begin
if y and = then mult:=mult*tmp mod mo;
tmp:=tmp*tmp mod mo;
y:=y>>;
end;
end; procedure dfs(k:longint);
var i,j,tot,sum:longint;
ans:int64;
begin
if k=m+ then
begin
tot:=; ans:=fac[n-]; sum:=n-;
for i:= to m do
if c[i]> then
begin
inc(tot);
inc(flag[a[i]]);
sum:=sum-(b[i]-);
end;
for i:= to m do
if (c[i]>)and(flag[a[i]]>) then
begin
for j:= to m do
if c[j]> then dec(flag[a[j]]);
exit;
end;
for i:= to m do
if c[i]> then ans:=ans*exf[b[i]-] mod mo;
//writeln(sum);
if sum>= then
begin
ans:=ans*exf[sum] mod mo;
ans:=ans*mult(n-tot,sum) mod mo;
if tot and = then sun:=(sun+ans) mod mo
else sun:=(sun-ans) mod mo;
sun:=(sun+mo) mod mo;
end;
for i:= to m do
if c[i]> then dec(flag[a[i]]);
exit;
end;
if b[k]= then dfs(k+)
else
begin
c[k]:=;
dfs(k+);
c[k]:=;
dfs(k+);
end;
end; begin
assign(input,'51nod1806.in'); reset(input);
assign(output,'51nod1806.out'); rewrite(output);
readln(n,m);
if n= then
begin
writeln();
exit;
end;
for i:= to m do read(a[i],b[i]);
exf[]:=; exf[]:=; fac[]:=;
for i:= to do fac[i]:=fac[i-]*i mod mo;
for i:= to do exf[i]:=exf[mo mod i]*(mo-mo div i) mod mo;
for i:= to do exf[i]:=exf[i-]*exf[i] mod mo;
dfs();
writeln(sun);
close(input);
close(output);
end.
【51NOD1806】wangyurzee的树(Prufer编码,容斥原理,组合计数)的更多相关文章
- 【Foreign】树 [prufer编码][DP]
树 Time Limit: 10 Sec Memory Limit: 256 MB Description Input Output Sample Input 3 2 2 1 Sample Outp ...
- 「模拟赛20180406」膜树 prufer编码+概率
题目描述 给定一个完全图,保证\(w_{u,v}=w_{v,u}\)且\(w_{u,u}=0\),等概率选取一个随机生成树,对于每一对\((u,v)\),求\(dis(u,v)\)的期望值对\(998 ...
- 树的Prufer 编码和最小生成树计数
Prufer数列 Prufer数列是无根树的一种数列.在组合数学中,Prufer数列由有一个对于顶点标过号的树转化来的数列,点数为n的树转化来的Prufer数列长度为n-2.它可以通过简单的迭代方 ...
- BSOJ 5553 wangyurzee的树 prufer序列 容斥
BSOJ我也不知道在哪. 容易想到容斥. 考虑不合法的方案 想到强制某个点的度数为限制即可. 这样就变成了了总方案-一个不合法+两个不合法-3个......的模型了. 坑点 当强制两个相同的点时 方案 ...
- HDU 4390 Number Sequence (容斥原理+组合计数)
HDU 4390 题意: 大概就是这样.不翻译了: Given a number sequence b1,b2-bn. Please count how many number sequences a ...
- 【BZOJ 3027】 3027: [Ceoi2004]Sweet (容斥原理+组合计数)
3027: [Ceoi2004]Sweet Time Limit: 1 Sec Memory Limit: 128 MBSubmit: 71 Solved: 34 Description John ...
- 51nod 1806 wangyurzee的树
基准时间限制:1 秒 空间限制:131072 KB wangyurzee有n个各不相同的节点,编号从1到n.wangyurzee想在它们之间连n-1条边,从而使它们成为一棵树.可是wangyur ...
- Luogu2290 [HNOI2004]树的计数 (组合计数,prufer编码)
这不prufer编码吗,防爆long long就行了啊 #include <iostream> #include <cstdio> #include <cstring&g ...
- BZOJ1005--[HNOI2008]明明的烦恼(树的prufer编码)
1005: [HNOI2008]明明的烦恼 Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 5768 Solved: 2253[Submit][Stat ...
- 树的prufer编码
prufer是无根树的一种编码方式,一棵无根树和一个prufer编码唯一对应,也就是一棵树有唯一的prufer编码,而一个prufer编码对应一棵唯一的树. 第一部分:树编码成prufer序列. 树编 ...
随机推荐
- 转 windows10 U盘系统启动盘怎么制作
转 windows10 U盘系统启动盘怎么制作 现将http://jingyan.baidu.com/article/9f7e7ec05e24d56f29155455.html 将dvd 写入 is ...
- mvc的生命周期
尽管每次都是从msdn里复制粘贴,但是还是有必要的,加深对mvc的理解和官方的表述. 了解 MVC 应用程序执行过程 发送给基于 ASP.NET MVC 的 Web 应用程序的请求首先通过 UrlRo ...
- 使用Kotlin,抛弃findViewById
有没有觉得Android的findViewById挺烦人的.使用Kotlin可以让你彻底抛弃这个烦恼 步骤1.在build.gradle(Module:app)中添加如下一句话 这个在老一点版本的An ...
- Vue.js——打包之后资源路径产生问题
https://blog.csdn.net/qq_30632003/article/details/79353035 https://www.cnblogs.com/diantao/p/7776523 ...
- ES6扩展运算符的使用
<!DOCTYPE html> <html> <head> <meta charset="UTF-8"> <title> ...
- <jsp:include page="${pageContext.request.contextPath/index.jsp" ></jsp:include> 引发的错误
路径引发的错误 如下使用项目路径对jsp页面进行获取,会报javax.servlet.ServletException: File "/web/dbwx/web/public/page_to ...
- SpringBoot 快速开发框架
学习资源:https://ke.qq.com/course/260513(这是Springboot升级版本教程,里面还有一个初级版本的) 1.第一个测试程序,那个覆盖方法加上@Override会报错, ...
- Java基础(十三)--深拷贝和浅拷贝
在上篇文章:Java基础(十二)--clone()方法,我们简单介绍了clone()的使用 clone()对于基本数据类型的拷贝是完全没问题的,但是如果是引用数据类型呢? @Data @NoArgsC ...
- vue组件之this指向问题
[问题描述] 返回顶部组件里,用到数据操作.通过方法里改动this数据,但发现直接使用失效 mounted() { window.onscroll=function(){ ) { this.isAct ...
- JSP的九大对象和四大作用域
1.JSP中九大内置对象为: request 请求对象 类型 javax.servlet.ServletRequest 作用域 Request ...