[点双连通分量-奇环判定]poj 2942:Knights of the Round Table
大致题意:
给出一个无向图,求出这个图的补图G。求出在G中存在多少个点使得这些点不属于任何一个奇环中。
大致思路:
这道题从一开始就是个杯具!先是把题目想成边的双连通分量做,狂哇。后来用割点的模版照着网上的代码改造出一个点双连通分量的代码,依然哇!!原因是同一个点可能属于不同的点双连通分量!! 后来把搜索染色的过程加到Tarjan里面,还是哇!!检查出来原因是搜索u的时候,起点并没有标记其在那一块双连通分量!! 前前后后花了我八天才AC……大概得终身难忘这道题了>_< 发图供大牛BS
给出一组可以判边双连通分量死刑的数据
#include<iostream>
#include<cstdio>
#include <algorithm>
#include<cstring>
using namespace std;
const int inf=1<<30;
const int nMax=30015;
const int mMax=5000000;
class edge{
public:
int u,v,nex;
};edge e[mMax];
int k,k1,head[nMax];//head[i]是以点i为起点的链表头部
int dfn[nMax],low[nMax],sta[nMax],top,atype,belon[nMax],dep; //Tarjan需要的数据记录 //atype 强连通分量的个数
bool insta[nMax];
int n,m;
int ans,vis[nMax],num[nMax];
void addedge(int a,int b){
e[k].u=a;
e[k].v=b;
e[k].nex=head[a];
head[a]=k;k++;
}
int change(int a){
if(a==1)return 2;
if(a==2)return 1;
return 0;
}
int clor[nMax],cnt;
bool dfs(int loc,int val){
int i,v;
clor[loc]=val;
for(i=head[loc];i;i=e[i].nex){
v=e[i].v;
if(atype!=belon[v])continue;
if(clor[v]==0){
if(!dfs(v,change(val))){
return 0;
}
}
else{
if(clor[v]!=change(clor[loc])){
return 0;
}
}
}
return 1;
}
void Tarjan(int u,int rt){ //我的Tarjan模版
int i,t;
dfn[u]=low[u]=++dep;
sta[++top]=u;
insta[u]=1;
for(i=head[u];i;i=e[i].nex){
int v=e[i].v;
if(v==rt)continue;
if(!dfn[v]){
Tarjan(v,u);
low[u]=min(low[u],low[v]);
if(dfn[u]<=low[v]){ ///要注意一个点可能同时属于不同的点双连通分量
atype++; //点双通分量个数
do{
t=sta[top--];
insta[t]=0;
belon[t]=atype; //第j个点属于第atype个连通块
num[++cnt]=t;
}while(v!=t);
num[++cnt]=u;
memset(clor,0,sizeof(clor));
if(cnt>=3&&!dfs(u,1)){
while(cnt>0)
vis[num[cnt--]]=1;
}
else cnt=0;
}
}
else{
if(insta[v])low[u]=min(low[u],dfn[v]);
}
}
}
void init(){
k=k1=1;
dep=1;
top=atype=0;
memset(insta,0,sizeof(insta)); //是否在栈中
memset(head,0,sizeof(head)); //静态链表头指针
memset(low,0,sizeof(low)); //Tarjan的low数组
memset(dfn,0,sizeof(dfn)); //Tarjan的dfn数组
memset(belon,0,sizeof(belon)); //记录每个点属于哪一个双连通分量
}
bool map[2000][2000];
int main(){
int i,j,a,b;
while(scanf("%d%d",&n,&m)!=EOF&&(n||m)){
init();
ans=0;
cnt=0;
memset(map,0,sizeof(map));
while(m--){
scanf("%d%d",&a,&b);
map[a][b]=map[b][a]=1;
}
for(i=1;i<=n;i++){
for(j=1;j<=n;j++){
if(!map[i][j]&&i!=j){
addedge(i,j);
}
}
}
memset(vis,0,sizeof(vis));
for(i=1;i<=n;i++){
if(!dfn[i])Tarjan(i,-1);
}
for(i=1;i<=n;i++){
if(!vis[i])ans++;
}
printf("%d\n",ans);
}
return 0;
}