poj 2942 点的双连通分量
思路:
对于该图,直接用建图貌似没法解,所以也很容易想到建补图,这样存在边的两个点就能再圆桌上做一起。也就将问题转化为对双连通分量中是否存在奇圈了。
我们将每次查询的边保存在stack中,当遇到关键点的时候,stack里面保存的就是一个连通分量。在该连通分量中进行深搜,每次标记一个与父节点相反的颜色。当某次子节点与父节点颜色相同,那么就存在奇圈,且该连通分量中所有的点都在奇圈中。将这些点标记,最后进行遍历就行了。
引用discuss里的话:
一个块若无法做二分图染色,势必存在一个长度为奇数的环
任找一个奇环C,则对于任意一个非环上的点A,一定有两条不相交的路,连向这个奇环,交奇环于两个不同的点P、Q(否则这就不是一个双连通分量)
那么在环C上有两条P-->Q的路径,一条经过奇数条边,一条经过偶数条边
其中一条同PA、AQ相连后,一定是个奇环,所以A一定也在一个奇环上
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<queue> #define Maxn 1010 #define Maxm Maxn*Maxn using namespace std; int index[Maxn],vi[Maxn],dfn[Maxn],col[Maxn],low[Maxn],map[Maxn][Maxn],e,n,lab=0,stack[Maxm],top,odd[Maxn]; void init() { memset(index,-1,sizeof(index)); memset(vi,0,sizeof(vi)); memset(map,1,sizeof(map)); memset(col,0,sizeof(col)); memset(low,0,sizeof(low)); memset(dfn,0,sizeof(dfn)); memset(odd,0,sizeof(odd)); e=lab=top=0; } struct Edge{ int from,to,next,v; }edge[Maxm]; void addedge(int from, int to) { edge[e].v=0; edge[e].from=from; edge[e].to=to; edge[e].next=index[from]; index[from]=e++; edge[e].v=0; edge[e].to=from; edge[e].from=to; edge[e].next=index[to]; index[to]=e++; } int find(int u) { int i,j,temp; for(i=index[u];i!=-1;i=edge[i].next) { temp=edge[i].to; if(vi[temp]) { if(col[temp]==-1) { col[temp]=!col[u]; if(find(temp))//寻找奇圈 return 1; } else if(col[temp]==col[u]) return 1; } } return 0; } int color(int u) { memset(col,-1,sizeof(col)); memset(vi,0,sizeof(vi)); col[u]=0; int i; do{//将该连通分量进行标记 i=stack[--top]; vi[edge[i].from]=1; vi[edge[i].to]=1; } while(edge[i].from!=u); if(find(u))//如果找到就进行标记 { for(i=1;i<=n;i++) { if(vi[i]) odd[i]=1; } } return 0; } int dfs(int u) { dfn[u]=low[u]=++lab; int i,j,temp; for(i=index[u];i!=-1;i=edge[i].next) { temp=edge[i].to; if(edge[i].v) continue;//一开始没加这个判断,一直WA edge[i].v=edge[i^1].v=1; stack[top++]=i; if(!dfn[temp]) { dfs(temp); if(low[temp]>=dfn[u]) color(u); low[u]=min(low[u],low[temp]); } low[u]=min(low[u],dfn[temp]); } return 0; } int main() { int m,i,j,a,b; while(scanf("%d%d",&n,&m),n||m) { init(); for(i=1;i<=m;i++) { scanf("%d%d",&a,&b); map[a][b]=map[b][a]=0; } for(i=1;i<=n;i++) { for(j=i+1;j<=n;j++) { if(map[i][j]) addedge(i,j); } } for(i=1;i<=n;i++) if(!dfn[i]) dfs(i); int ans=0; for(i=1;i<=n;i++) if(!odd[i]) ans++; printf("%d\n",ans); } return 0; }