题意:给你一个无向图,问图中有多少个符合条件的集合?条件为这个集合里面存在一个子集(大小>=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的更多相关文章

  1. HDU 5917 Instability ramsey定理

    http://acm.hdu.edu.cn/showproblem.php?pid=5917 即世界上任意6个人中,总有3个人相互认识,或互相皆不认识. 所以子集 >= 6的一定是合法的. 然后 ...

  2. HDU - 5917 水题

    题意:n个点m条边,找点集个数,点集满足有任意三个点成环,或者三个点互不相连 题解:暴力复杂度O(n^5/120*O(ok))==O(能过) //#pragma comment(linker, &qu ...

  3. 近几年ACM/ICPC区域赛铜牌题

    2013 changsha zoj 3726 3728 3736 3735 2013 chengdu hud 4786 4788 4790 2013 hangzhou hdu 4770 4771 47 ...

  4. HDOJ 2111. Saving HDU 贪心 结构体排序

    Saving HDU Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total ...

  5. 【HDU 3037】Saving Beans Lucas定理模板

    http://acm.hdu.edu.cn/showproblem.php?pid=3037 Lucas定理模板. 现在才写,noip滚粗前兆QAQ #include<cstdio> #i ...

  6. hdu 4859 海岸线 Bestcoder Round 1

    http://acm.hdu.edu.cn/showproblem.php?pid=4859 题目大意: 在一个矩形周围都是海,这个矩形中有陆地,深海和浅海.浅海是可以填成陆地的. 求最多有多少条方格 ...

  7. HDU 4569 Special equations(取模)

    Special equations Time Limit:1000MS     Memory Limit:32768KB     64bit IO Format:%I64d & %I64u S ...

  8. HDU 4006The kth great number(K大数 +小顶堆)

    The kth great number Time Limit:1000MS     Memory Limit:65768KB     64bit IO Format:%I64d & %I64 ...

  9. HDU 1796How many integers can you find(容斥原理)

    How many integers can you find Time Limit:5000MS     Memory Limit:32768KB     64bit IO Format:%I64d ...

随机推荐

  1. C#中的异步编程--探索await与async关键字的奥妙之处,原来理解和使用异步编程可以这么简单

    前言 await与async是C#5.0推出的新语法,关于await与async有很多文章讲解.但看完后有没有这样一种感觉,感觉这东西像是不错,但好像就是看不太懂,也不清楚该怎么使用.虽然偶有接触,但 ...

  2. Hadoop学习之路(8)Yarn资源调度系统详解

    文章目录 1.Yarn介绍 2.Yarn架构 2.1 .ResourceManager 2.2 .ApplicationMaster 2.3 .NodeManager 2.4 .Container 2 ...

  3. too many open files异常及解决办法

    1.问题描述: 最近项目中由于大批量的操作文件,项目部署在Linux服务器上,突然发现出现了 too many open files异常. 经过一番查阅,(直译过来就是文件打开数过多)但是查阅了相关的 ...

  4. 浅谈python的第三方库——pandas(终)

    作为pandas系列的最终章,本文引出一个数据"复制"问题. 示例如下: 从上图中可以看到:我们对data_pd做了删除一行的操作,但是这并没有改变变量data_pd在内存中的值, ...

  5. ES6基础与解构赋值(高颜值弹框小案例!)

    let只作用在当前块级作用域内使用let或者const声明的变量,不能再被重新声明let不存在`变量提升` console.log(dad); var dad = '我是爸爸!';//预定义undef ...

  6. C#常规TcpListener

    1.Xaml <Window x:Class="Server.MainWindow" xmlns="http://schemas.microsoft.com/win ...

  7. 获取redis指定实例中所有的key

    需求:获取redis指定的实例中所有的key的名字. 千万不要使用keys *,可以使用scan命令的递归方式获取. 以下给出自己写的脚本,经过测试效果还可以. db_ip=5.5.5.101 db_ ...

  8. 第十周 11.28 psp0

    课堂测试总结 学生:马小心                                              日期:2017.11.28 作业号 日期 过程 估计数据 实际数据 累计数据 时间 ...

  9. Vue中axios有关请求头的几点小结

    在Vue前端中向后端发起http请求会有着两种写法:一种是在vue文件中直接导入axios模板,另外一种是使用Vue的属性$http. 1.在第一种方式中,在同一个工程中所添加的vue文件直接使用ax ...

  10. Python入门10 —— for循环

    1.字符串依次取值 students = ['egon', 'lxx', 'alex'] i = 0 while i < 3: print(students[i]) i += 1 2.针对循环取 ...