poj2942 Knights of the Round Table 双连通分支 tarjan

题解：http://blog.csdn.net/lyy289065406/article/details/6756821

讲的很详细我就不多说了。

题目连接：http://poj.org/problem?id=2942

代码（完全白书上的）

 #include <stdio.h>
 #include <string.h>
 #include <iostream>
 #include <algorithm>
 #include <stdlib.h>
 #include <vector>
 #include <queue>
 #include <stack>
 #define loop(s,i,n) for(i = s;i < n;i++)
 #define cl(a,b) memset(a,b,sizeof(a))
 const int maxn = ;
 using namespace std;
 int pre[maxn],iscut[maxn],bccno[maxn],dfsclock,bcc_cnt;
 vector<int>g[maxn],bcc[maxn];
 struct edge
 {
     int u,v,w;
 };
 stack<edge> st;

 int dfs(int u,int fa)
 {
    int lowu= pre[u]=++dfsclock;
    int child=;
     for(int i=;i<g[u].size();i++)
     {
         int v=g[u][i];
         edge e;
         e.u=u,e.v=v;
         if(!pre[v])
         {
             st.push(e);
             child++;
             int lowv=dfs(v,u);
             lowu=min(lowu,lowv);
             if(lowv>=pre[u])
             {
                 iscut[u]=;
                 bcc_cnt++;
                 bcc[bcc_cnt].clear();
                 for(;;)
                 {
                     edge x=st.top();st.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)
         {
             st.push(e);
             lowu=min(lowu,pre[v]);
         }
     }
     if(fa<&&child==)iscut[u]=;
     return lowu;
 }
 void find_bcc(int n)
 {
     int i;
     memset(pre,,sizeof(pre));
     cl(iscut,);
     cl(bccno,);
     dfsclock = bcc_cnt = ;
     loop(,i,n)
     if(!pre[i]) dfs(i,-);
    // puts("yes");
    // printf("%d  """"""\n",bcc_cnt);
 }
 int odd[maxn],color[maxn];
 bool bipartite(int u,int b)
 {
     int i;
     loop(,i,g[u].size())
     {
         int v = g[u][i];
         if(bccno[v] != b)
         continue;
         if(color[v] == color[u])
         return false;
         if(!color[v])
         {
             color[v] = -color[u];
             if(!bipartite(v,b))
             return false;
         }
     }
     return true;
 }
 int map[maxn][maxn];
 int main()
 {
     int i,kase = ,n,m;

     while(scanf("%d %d",&n,&m))
     {
         if(!n)
         break;
         memset(map,,sizeof(map));

         loop(,i,n) g[i].clear();

         loop(,i,m)
         {
             int u,v;
             scanf("%d %d",&u,&v);
             u--,v--;
             map[u][v] = map[v][u] =;
         }
         int j;
         loop(,i,n)
         {
             loop(i+,j,n)
             {
                 if(!map[i][j])
                 g[i].push_back(j),g[j].push_back(i);
             }
         }

         find_bcc(n);

        // cout<<bcc_cnt<<endl;

         memset(odd,,sizeof(odd)) ;
         for(int i=;i<=bcc_cnt;i++)
         {
         memset(color,,sizeof(color));
         for(int j=;j<bcc[i].size();j++)
             bccno[bcc[i][j]]=i;
         int u=bcc[i][];
         color[u]=;
         if(!bipartite(u,i))//判二分图的原因是如果偶数圈会可能是二分图，但是奇数圈不会是~
             for(int j=;j<bcc[i].size();j++)
                 odd[bcc[i][j]]=;
         }
         int ans=n;
         for(int i=;i<n;i++)
             if(odd[i])ans--;

         printf("%d\n",ans);
     }
     return ;
 }