hdu 5917
题意:给你一个无向图,问图中有多少个符合条件的集合?条件为这个集合里面存在一个子集(大小>=3)为团或者都是孤立点。答案mod1e9+7;
根据 Ramsey定理,大于等于6个的集合,肯定存在一个子集的边都是红色或者都是蓝色,即为团还是为孤立点;
所以当n大于等于6的时候,所有的取6个或六个以上的子集的集合都是符合的,所以将这些排列组合的方式全部都计算在内;
即C(n,i) i的取值范围为(6~n) 但是这样子算会超时,我们可以计算C(n,i) i从0开始计算,这样子所有的数加起来,就是2^n
然后再减去C(n,i)i从0到5即可;
n的取值范围到50,所以暴力不超时;
然后还剩下点集为3到5的情况,这个时候,我们就分别枚举点集为3,4,5,判断是否满足情况即可(题目给出的m条边就在这个时候起作用)
只要满足图中有3个为孤立点即可;
代码如下:
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
const int M = 1e2+;
const LL mod = 1e9+;
int ncase,n,m,cas=;
int vis[M][M];
LL fac[M];
void init()
{ ///预处理阶乘
fac[]=;fac[]=;
for (int i=;i<M;i++)
fac[i]=fac[i-]*i%mod;
}
LL quick_pow(LL p,LL k) ///快速幂
{
LL res=,tp=p;
if(k<) return ;
while(k){
if(k&) res=res*tp%mod;
tp=tp*tp%mod;
k>>=;
}
return res;
}
LL C(int n,int m)
{
return fac[n]*quick_pow(fac[m],mod-)%mod*quick_pow(fac[n-m],mod-)%mod;
}///费马小定理求逆元+快速幂
int judge(int i,int j,int k){ ///判断三点之间满不满足不稳定点集
if(vis[i][j]&&vis[i][k]&&vis[j][k]) return ;///三点之间相连
if(!vis[i][j]&&!vis[i][k]&&!vis[j][k]) return ; ///三点之间不互联
return ;
}
int judge1(int i,int j,int k,int l)
{
if(judge(i,j,k)||judge(i,j,l)||judge(j,k,l)||judge(i,k,l)) return ; ///表示4个点中有三个点为不稳定点集就行,为什么呢?
///因为时题目要求的,不稳定点集为3个点。
return ;
}
int judge2(int i,int j,int k,int l,int o){
///表示5个点中有4个点满足就行
if(judge1(i,j,k,l)||judge1(i,j,k,o)||judge1(i,j,l,o)||judge1(j,k,l,o)||judge1(i,k,l,o)) return ;
return ;
} int main(){ init();
scanf("%d",&ncase);
while(ncase--){
scanf("%d%d",&n,&m);
memset(vis,false,sizeof(vis));
for(int i=;i<=m;i++){
int u,v;
scanf("%d%d",&u,&v);
vis[u][v]=vis[v][u]=;
}
LL ans=;
if(n>=){
ans=(ans+quick_pow(,n))%mod; ///多项式定理C(n,0)+C(n,1)+...+C(n,n)=2^n
///C(n,k)(k>=6)表示,从n个中取k个出来,总存在一个不稳定点集的个数(三点之间互联或三点之间不连)
for(int i=;i<;i++) ///减去前5个
ans=(ans-C(n,i)+mod)%mod;
}
if(n>=){ ///暴力取3个
for(int i=;i<=n;i++)
for(int j=i+;j<=n;j++)
for(int k=j+;k<=n;k++)
if(judge(i,j,k)) ans=(ans+)%mod;
}
if(n>=){ ///暴力取4个
for(int i=;i<=n;i++)
for(int j=i+;j<=n;j++)
for(int k=j+;k<=n;k++)
for(int l=k+;l<=n;l++)
if(judge1(i,j,k,l)) ans=(ans+)%mod;
}
if(n>=){ ///暴力取5个
for(int i=;i<=n;i++)
for(int j=i+;j<=n;j++)
for(int k=j+;k<=n;k++)
for(int l=k+;l<=n;l++)
for(int o=l+;o<=n;o++)
if(judge2(i,j,k,l,o)) ans=(ans+)%mod;
}
printf("Case #%d: %lld\n",cas++,ans);
}
return ;}
hdu 5917的更多相关文章
- HDU 5917 Instability ramsey定理
http://acm.hdu.edu.cn/showproblem.php?pid=5917 即世界上任意6个人中,总有3个人相互认识,或互相皆不认识. 所以子集 >= 6的一定是合法的. 然后 ...
- HDU - 5917 水题
题意:n个点m条边,找点集个数,点集满足有任意三个点成环,或者三个点互不相连 题解:暴力复杂度O(n^5/120*O(ok))==O(能过) //#pragma comment(linker, &qu ...
- 近几年ACM/ICPC区域赛铜牌题
2013 changsha zoj 3726 3728 3736 3735 2013 chengdu hud 4786 4788 4790 2013 hangzhou hdu 4770 4771 47 ...
- HDOJ 2111. Saving HDU 贪心 结构体排序
Saving HDU Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total ...
- 【HDU 3037】Saving Beans Lucas定理模板
http://acm.hdu.edu.cn/showproblem.php?pid=3037 Lucas定理模板. 现在才写,noip滚粗前兆QAQ #include<cstdio> #i ...
- hdu 4859 海岸线 Bestcoder Round 1
http://acm.hdu.edu.cn/showproblem.php?pid=4859 题目大意: 在一个矩形周围都是海,这个矩形中有陆地,深海和浅海.浅海是可以填成陆地的. 求最多有多少条方格 ...
- HDU 4569 Special equations(取模)
Special equations Time Limit:1000MS Memory Limit:32768KB 64bit IO Format:%I64d & %I64u S ...
- HDU 4006The kth great number(K大数 +小顶堆)
The kth great number Time Limit:1000MS Memory Limit:65768KB 64bit IO Format:%I64d & %I64 ...
- HDU 1796How many integers can you find(容斥原理)
How many integers can you find Time Limit:5000MS Memory Limit:32768KB 64bit IO Format:%I64d ...
随机推荐
- C#中的异步编程--探索await与async关键字的奥妙之处,原来理解和使用异步编程可以这么简单
前言 await与async是C#5.0推出的新语法,关于await与async有很多文章讲解.但看完后有没有这样一种感觉,感觉这东西像是不错,但好像就是看不太懂,也不清楚该怎么使用.虽然偶有接触,但 ...
- Hadoop学习之路(8)Yarn资源调度系统详解
文章目录 1.Yarn介绍 2.Yarn架构 2.1 .ResourceManager 2.2 .ApplicationMaster 2.3 .NodeManager 2.4 .Container 2 ...
- too many open files异常及解决办法
1.问题描述: 最近项目中由于大批量的操作文件,项目部署在Linux服务器上,突然发现出现了 too many open files异常. 经过一番查阅,(直译过来就是文件打开数过多)但是查阅了相关的 ...
- 浅谈python的第三方库——pandas(终)
作为pandas系列的最终章,本文引出一个数据"复制"问题. 示例如下: 从上图中可以看到:我们对data_pd做了删除一行的操作,但是这并没有改变变量data_pd在内存中的值, ...
- ES6基础与解构赋值(高颜值弹框小案例!)
let只作用在当前块级作用域内使用let或者const声明的变量,不能再被重新声明let不存在`变量提升` console.log(dad); var dad = '我是爸爸!';//预定义undef ...
- C#常规TcpListener
1.Xaml <Window x:Class="Server.MainWindow" xmlns="http://schemas.microsoft.com/win ...
- 获取redis指定实例中所有的key
需求:获取redis指定的实例中所有的key的名字. 千万不要使用keys *,可以使用scan命令的递归方式获取. 以下给出自己写的脚本,经过测试效果还可以. db_ip=5.5.5.101 db_ ...
- 第十周 11.28 psp0
课堂测试总结 学生:马小心 日期:2017.11.28 作业号 日期 过程 估计数据 实际数据 累计数据 时间 ...
- Vue中axios有关请求头的几点小结
在Vue前端中向后端发起http请求会有着两种写法:一种是在vue文件中直接导入axios模板,另外一种是使用Vue的属性$http. 1.在第一种方式中,在同一个工程中所添加的vue文件直接使用ax ...
- Python入门10 —— for循环
1.字符串依次取值 students = ['egon', 'lxx', 'alex'] i = 0 while i < 3: print(students[i]) i += 1 2.针对循环取 ...