uvaoj 1364 - Knights of the Round Table

直接抄袭训练指南上的代码。。

从这道题学到了一个结论，一个双连通分量，要么含有奇圈，要么是二分图

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<stack>
using namespace std;

const int maxn=1005;

struct Edge{
    int u,v;
    Edge(int u,int v):u(u),v(v){}
};

int pre[maxn],iscut[maxn],bccno[maxn],dfs_clock,bcc_cnt;
vector<int> G[maxn],bcc[maxn];

stack<Edge> S;

int dfs(int u,int fa)
{
    int lowu=pre[u]=++dfs_clock;
    int child=0;
    for(int i=0;i<G[u].size();++i)
    {
        int v=G[u][i];
        Edge e=Edge(u,v);
        if(!pre[v])
        {
            S.push(e);
            child++;
            int lowv=dfs(v,u);
            lowu=min(lowu,lowv);
            if(lowv>=pre[u])
            {
                iscut[u]=1;
                bcc_cnt++;bcc[bcc_cnt].clear();
                for(;;)
                {
                    Edge x=S.top();S.pop();
                    if(bccno[x.u]!=bcc_cnt){bcc[bcc_cnt].push_back(x.u);bccno[x.u]=bcc_cnt;}
                    if(bccno[x.v]!=bcc_cnt){bcc[bcc_cnt].push_back(x.v);bccno[x.v]=bcc_cnt;}
                    if(x.u==u&&x.v==v)break;
                }
            }
        }
        else if(pre[v]<pre[u]&&v!=fa)
        {
            S.push(e);
            lowu=min(lowu,pre[v]);
        }
    }
    if(fa<0&&child==1)iscut[u]=0;
    return lowu;
}

void find_bcc(int n)
{
    memset(pre,0,sizeof(pre));
    memset(iscut,0,sizeof(iscut));
    memset(bccno,0,sizeof(bccno));
    dfs_clock=bcc_cnt=0;
    for(int i=0;i<n;++i)
    {
        if(!pre[i])dfs(i,-1);
    }
}

int A[maxn][maxn];
int odd[maxn],color[maxn];

bool bipartite(int u,int b)
{
    for(int i=0;i<G[u].size();++i)
    {
        int v=G[u][i];if(bccno[v]!=b)continue;
        if(color[v]==color[u])return false;
        if(!color[v])
        {
            color[v]=3-color[u];
            if(!bipartite(v,b))return false;
        }
    }
    return true;
}

int main()
{
//    freopen("data.txt","r",stdin);
    int kase=0,n,m;
    while(scanf("%d%d",&n,&m)==2&&n)
    {
        for(int i=0;i<n;++i)
        {
            G[i].clear();
        }
        memset(A,0,sizeof(A));
        for(int i=0;i<m;++i)
        {
            int u,v;
            scanf("%d%d",&u,&v);
            u--;v--;
            A[u][v]=A[v][u]=1;
        }
        for(int u=0;u<n;++u)
        {
            for(int v=u+1;v<n;++v)
            {
                if(!A[u][v]){G[u].push_back(v);G[v].push_back(u);}
            }
        }

        find_bcc(n);

        memset(odd,0,sizeof(odd));
        for(int i=1;i<=bcc_cnt;++i)
        {
            memset(color,0,sizeof(color));
            for(int j=0;j<bcc[i].size();++j)bccno[bcc[i][j]]=i;//割顶的编号会跟连通分量不一样，这里给改成一样的
            int u=bcc[i][0];//开始的时候一直以为每个bcc保存的第一个都是割顶，就一直没理解这个算法。后来发现u只是随便选的一个点
            color[u]=1;
            if(!bipartite(u,i))
                for(int j=0;j<bcc[i].size();++j)odd[bcc[i][j]]=1;
        }
        int ans=n;
        for(int i=0;i<n;++i)
        {
            if(odd[i])ans--;
        }
        printf("%d\n",ans);
    }
    return 0;
}