HDU5917 RAMSEY定理

参考自：http://blog.csdn.net/Miracle_ma/article/details/52737597?locationNum=1
给你n个点，m条边，然后告诉你选择一个点集S
如果里面有一个子集A，A里面的点都不相连，或者都相连，则这个点集不稳定
求不稳定的个数
子集A的大小是大于等于3，所以考虑到6个点的图，里面肯定有3个点，互相有边，或者互相没边
所以用ramsey定理优化，6以上都可以直接求
剩下3，4，5的情况，搜索或者循环
#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cmath>
#include <memory.h>
#include <queue>
#include <map>
#include <set>
using namespace std;
typedef long long ll;
const int maxn = 1e5 + 10;
const int mod = 1e9 + 7;
int n,m,mp[60][60];
ll fact[60];
int ok3(int a,int b,int c){
    if(mp[a][b] == 0&&mp[a][c] == 0&&mp[b][c] == 0) return 1;
    if(mp[a][b] == 1&&mp[a][c] == 1&&mp[b][c] == 1) return 1;
    return 0;
}
int ok4(int a,int b,int c,int d){
    if(ok3(a,b,c)||ok3(a,b,d)||ok3(a,c,d)||ok3(b,c,d)) return 1;
    else return 0;
}
int ok5(int a,int b,int c,int d,int e){
    if(ok4(a,b,c,d)||ok4(a,b,c,e)||ok4(a,b,d,e)||ok4(a,c,d,e)||ok4(b,c,d,e)) return 1;
    else return 0;
}
void init(){
    fact[0]=1;
    for(int i=1;i<=50;i++) fact[i]=fact[i-1]*i%mod;
}
ll qpow(ll a,ll n){
    ll ans=1;
    while(n){
        if(n&1) ans=ans*a%mod;
        a=a*a%mod;
        n>>=1;
    }
    return ans;
}

ll C(int n,int m){
    return fact[n]*qpow(fact[m],mod-2)%mod*qpow(fact[n-m],mod-2)%mod;
}

int main() {
   // freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int T;
    init();
    int cas = 0;
    cin>>T;
    while(T--){
        scanf("%d%d",&n,&m);
        memset(mp,0,sizeof(mp));
        int x,y;
        for(int i = 1;i <= m;i++){
            scanf("%d%d",&x,&y);
            mp[x][y] = 1;
            mp[y][x] = 1;
        }
        ll ans = 0;
        for(int i = 1;i <= n;i++)
            for(int j = i + 1;j <= n;j++)
                for(int k = j + 1;k <= n;k++)
                    if(ok3(i,j,k)) ans++;
        for(int i = 1;i <= n;i++)
            for(int j = i + 1;j <= n;j++)
                for(int k = j + 1;k <= n;k++)
                    for(int l = k + 1;l <= n;l++)
                        if(ok4(i,j,k,l)) ans++;
        for(int i = 1;i <= n;i++)
            for(int j = i + 1;j <= n;j++)
                for(int k = j + 1;k <= n;k++)
                    for(int l = k + 1;l <= n;l++)
                        for(int m = l + 1;m <= n;m++)
                            if(ok5(i,j,k,l,m)) ans++;
        ans %= mod;
         if(n>=6){
            for(int i=6;i<=n;i++){
                ans+=C(n,i);
                if(ans>=mod) ans-=mod;
            }
        }
        printf("Case #%d: %lld
",++cas,ans);
    }

    return 0;
}