HDU5727 Necklace（枚举 + 二分图最大匹配）

题目大概说有n个yang珠子n个yin珠子，要交替串成一个环形项链，有些yang珠子和某个yin珠子相邻这个yang珠子会不高兴，问最少有几个yang珠子不高兴。

自然会想到直接用状压DP去解，转移很烦，也没写出来。标程是搜索不明觉厉。。听闻了可以枚举一边的顺序，8!，然后用最大匹配解决。

然后想到的是枚举yang的顺序，然后对于每一个yang去其匹配下一个的yin，即X部是yang，而Y部是yin。不过这样开头那个yang可能出现少算的情况。。这个搞了好久都不行。。

其实，枚举yin的顺序，X部是各个yang，而Y部是位置！然后问题就引刃而解。我太菜了。。

另外用最大流超时了，然后用了个匈牙利过了，O((n-1)!*n³)的时间复杂度。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

int n,map[22][22],mat[22];
bool used[22];
bool crosspath(int k){
    for(int i=1; i<=map[k][0]; ++i){
        int j=map[k][i];
        if(!used[j]){
            used[j]=true;
            if(mat[j]==0 || crosspath(mat[j])){
                mat[j]=k;
                return true;
            }
        }
    }
    return false;
}

int hungary(){
    int res=0;
    for(int i=1; i<=n; ++i){
        memset(used,0,sizeof(used));
        if(crosspath(i)) ++res;
    }
    return res;
}

bool rel[11][11];
int main(){
    int m;
    while(~scanf("%d%d",&n,&m)){
        memset(rel,0,sizeof(rel));
        int a,b;
        while(m--){
            scanf("%d%d",&a,&b);
            rel[a][b]=1;
        }

        if(n==1 && rel[1][1]){
            puts("1");
            continue;
        }
        if(n<=1){
            puts("0");
            continue;
        }

        int seq[11];
        for(int i=1; i<=n; ++i){
            seq[i]=i;
        }

        int res=0;
        do{
            for(int i=1; i<=2*n; ++i) map[i][0]=0;
            memset(mat,0,sizeof(mat));
            for(int i=1; i<=n; ++i){
                for(int j=1; j<=n; ++j){
                    if(rel[i][seq[j]] || rel[i][seq[j%n+1]]) continue;
                    map[i][++map[i][0]]=j+n;
                    map[j+n][++map[j+n][0]]=i;
                }
            }
            res=max(res,hungary());
        }while(next_permutation(seq+2,seq+1+n));
        printf("%d
",n-res);
    }
    return 0;
}