POJ 2942—— Knights of the Round Table双连通分量+交叉染色

很经典的题目，可是我做了四天才做出来，题意是：一群骑士，某些骑士之间互相憎恨，如果在一起容易发生争斗事件，因此他们只有满足一定条件才能参加圆桌会议：1.圆桌边上任意相邻的两个骑士不能互相憎恨；2。同一个圆桌边上的骑士数量必须是奇数；

分析： 让不互相憎恨的骑士才能坐在一起，可以如下处理：将各骑士看成点，不互相憎恨的骑士之间存在边，建立无向图。先按照要求1，将所有能坐在一起的骑士分为一组，全部骑士分为若干组，每一组在图中是一个双连通分量。然后根据双连通分量的性质，判断双连通分量中是否存在奇圈，如果存在，则这一组骑士都能参与会议，反则这一组都不能参与会议；

具体的方法是：

1.搜索双连通分量。深度优先搜索过程中，用一个栈保存所有经过的节点，判断割点，碰到割点就标记当前栈顶点的节点并退栈，直到 当前节点停止并标记当前割点。标记过的节点处于同一个双连通分量

2.交叉染色搜索奇圈。在一个节点大于2的双连通分量中，必定存在一个圈经过该连通分量的所有节点；如果这个圈是奇圈，则该连通分量内所有的点都满足条件；若这个圈是圈，如果包含奇圈，则必定还有一个奇圈经过所有剩下的点。因此一个双连通分量中只要存在一个奇圈，那么该双连通分量内所有的点都处于一个奇圈中。根据这个性质，只需要在在一个双连通分量内找奇圈即可判断该联通分量是否满足条件。交叉染色法就是在dfs的过程中反复交换着用两种不同的颜色对为染色过的点染色，若某次dfs中当前节点的子结点和当前节点同色，则找到奇圈。

3.需要注意的地方是：因为同一个点可能在多个双连通分量中，因此标记某个点是否满足条件必须专门用一数组标记

#include <stdio.h>
#include <memory.h>

#define N 1001
#define init(a) memset(a,0,sizeof(a))
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
struct block
{
 int belongto;
 block *next;
 block(int u,int v):belongto(v),next(P[u]){}
}*P[N+1];
int G[N+1][N+1],used[N+1],part[N+1],deep[N+1],anc[N+1],open[N+1],color[N+1];
int n,m,block_n,top;

void readData()
{
 int i,j;
 init(G);
 while(m--)
 {
  scanf("%d%d",&i,&j);
  G[i][j]=G[j][i]=1;
 }
 for (i=1;i<=n;i++)
  G[i][i]=1;
 for (i=1;i<=n;i++)
 {
  for (j=1;j<=n;j++)
  G[0][j]=G[i][j];
  for (j=1;j<=n;j++)
  if(!G[0][j])
  G[i][++G[i][0]]=j;
 }
}

void con_dfs(int s,int father,int d)
{
 int i,j,k;
 anc[s]=deep[s]=d;
 used[s]=1;
 open[top++]=s;
 for(i=1;i<=G[s][0];i++)
 {
  j=G[s][i];
  if(j!=father&&used[j]==1)
  anc[s]=min(anc[s],deep[j]);
  if (!used[j])
  {
  con_dfs(j,s,d+1);
  anc[s]=min(anc[s],anc[j]);
  if(anc[j]>=deep[s])
  {
  block_n++;
  P[s]=new block(s,block_n);
  for(k=open[top];k!=j;P[k]=new block(k,block_n))
  k=open[--top];
  }
  }
 }
 used[s]=2;
}

void connex()
{
 init(used),init(P);
 block_n=top=0;
 for (int i=1;i<=n;i++)
  if(!used[i])
  con_dfs(i,-1,1);
}

int dfs_circle(int s,int col)
{
 color[s]=col;
 for (int i=1;i<=G[s][0];i++)
 {
  int j=G[s][i];
  if (part[j])
  {
  if(color[j]==0&&dfs_circle(j,-col))
  return 1;
  if(color[j]==col)
  return 1;
  }
 }
 return 0;
}

int cal_circle()
{
 int i,j,count=0;
 init(used);
 for (i=1;i<=block_n;i++)
 {
  init(part);
  init(color);
  for (j=1;j<=n;j++)
  for(block *L=P[j];L;L=L->next)
  if (L->belongto==i)
  {
  part[j]=1;
  break;
  }
  for (j=1;j<=n;j++)
  if (part[j])
  {
  if(dfs_circle(j,1))
  for(j=1;j<=n;j++)
  used[j]+=part[j];
  break;
  }
 }
 for (i=1;i<=n;i++)
  if(!used[i])
  count++;
 return count;
}

int main()
{
 while(scanf("%d%d",&n,&m)!=EOF&&n+m)
 {
  readData();
  connex();
  printf("%d/n",cal_circle());
 }
 return 0;
}