poj 2942 双连通分量+染色判二部图（驱除骑士）

题意：已知骑士有些人之间有矛盾，如果一个骑士能够出现在满足如下条件的任何一个圈中，则可以保留他；否则这个骑士必须要驱除。圈的条件是：1、人数大于1；2、总人数为奇数；3、有仇恨的骑士不能挨着坐。问需要驱除几个骑士（有几个骑士不能和任何人形成合法圈）？

思路：双连通分支+交叉染色。建图方法是将没有矛盾的骑士连边，求双连通分量（关于点的）。对于每个双连通分量，看是否存在奇圈，若存在那么这个双连通分量中的任意两骑士都可以同时出现在一个奇圈里。原因如下：因为每个不在那个奇圈里的点都一定能找到两条连到奇圈上不同节点的路径，然后奇圈被分割为一个总数为奇数的节点串和一个总数为偶数的节点串，这时做出正确的选择即可构成奇圈。

判断奇圈的方法显然交叉染色法即可。

细节：割点需要仔细考虑一下，因为它可以出现在"下面"的双连通分支中（即从它开始遍历）或者"上面"（在某个祖先结点处遍历）的中，不要计算重复。我的做法是在从x遍历的时候记录一下它是否出现在某个奇圈中（cut数组的作用），如果出现，那么它上面的那次（该点需要出栈的那次）就不判断了；否则当作非割点来判断。

再有那个do-while循环的判断条件写错了（一开始写成x!=stack[st-1]）导致WA了一阵(比如样例：5 4     1 2 1 3 2 4 3 4)。

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <algorithm>
#include <queue>
using namespace std;
#define INF 0x3fffffff
#define clr(s,t) memset(s,t,sizeof(s))
#define N 1005
#define M 1000005
struct edge{
    int y,next;
}e[M<<1];
int g[N][N],first[N],dfn[N],low[N],stack[N],flag[N],c[N],res,cut[N];
int top,con,id,st,n,m;
void add(int x,int y){
    e[top].y = y;
    e[top].next = first[x];
    first[x] = top++;
}
int dfs(int x,int col){
    c[x] = col;
    for(int i = first[x];i!=-1;i=e[i].next){
        if(flag[e[i].y] != flag[x])//不在一个双连通分支
            continue;
        if(c[e[i].y]==-1){
            if(!dfs(e[i].y,1-col))
                return 0;
        }else if(c[e[i].y]==col)//出现奇圈
            return 0;
    }
    return 1;
}
void tarjan(int x,int fa){
    int i,j,y;
    dfn[x] = low[x] = ++id;
    stack[st++] = x;
    for(i = first[x];i!=-1;i=e[i].next){
        y = e[i].y;
        if(y == fa)
            continue;
        if(dfn[y] == -1){
            tarjan(y, x);
            low[x] = min(low[x],low[y]);
            if(dfn[x] <= low[y]){
                con++;
                j = 0;
                do{
                    flag[stack[--st]] = con;
                    if(cut[stack[st]] == -1)//否则已经出现在某个圈中了，这里就不判断了
                        j++;
                }while(y != stack[st]);
                /*一开始写的x!=stack[st-1]WA。。原因在于x可能先连了一个能回到更早结点的儿子，才出现在一个双连通分支中，这样x和y之间塞入了其他的结点*/
                flag[x] = con;
                clr(c, -1);
                if(dfs(x,0))
                    res += j;
                else//表示出现在一个奇圈中，标定一下
                    cut[x] = 1;
            }
        }else
            low[x] = min(low[x],dfn[y]);
    }
}
int main(){
    while(scanf("%d %d",&n,&m) && (n+m)){
        int i,j,a,b;
        top = con = id = res = st = 0;
        clr(first, -1);
        memset(g,0,sizeof(g));
        clr(dfn, -1);
        clr(cut, -1);
        clr(flag, -1);
        for(i = 1;i<=m;i++){
            scanf("%d %d",&a,&b);
            g[a][b] = g[b][a] = 1;
        }
        for(i = 1;i<=n;i++)
            for(j = i+1;j<=n;j++)
                if(!g[i][j])
                    add(i,j),add(j,i);
        for(i = 1;i<=n;i++)
            if(dfn[i] == -1){
                tarjan(i,-1);
                if(cut[i] == -1)
                    res++;
            }
        printf("%d\n",res);
    }
    return 0;
}