LA3523 / poj2942 圆桌骑士 (二着色+双连通分量)

可以将所有人的连坐关系表示为一张无向图。每个人为一个节点,任意二人不憎恨则可以连坐,用一条边相连;憎恨则不相连。
因为参会者坐成一个圆圈,所以有:
“不参加任何会议的人数” == “不在任何一个简单奇圈的节点数”
插入知识:

  1. 简单圈:圈中每个节点只出现一次
  2. 奇数圈:圈中节点总数为奇数
  3. 二分图:一个图为二分图,当且仅当它不含长度为奇数的圈

由上可知,简单圈上的点一定属于同一个双连通分量,因此需要事先找到所有的双连通分量,而二分图中是没有奇圈的,因此我们只需要关注那些不是二分图的双连通分量即可。

综上可得思路:
1.根据不同人间的关系建立图
2.求出所有BCC
3.二着色法检查每个BCC,若其不为二分图,则把它的所有节点标记为”奇点“
4.输出”奇点“总数

第一次忘了在处理第二组数据的时候重置数组,WA了一发

//LA3523 / poj2942
#include <cstdio>
#include <stdio.h>
#include <algorithm>
#include <vector>
#include <cstring>
#include <stack>
using namespace std;
const int maxn=1000+10;


const int BLACK=-1;
const int WHITE=1;
const int NONE=0;
int color[maxn];
int noedge[maxn][maxn];
bool isodd[maxn];

struct Edge{
	int u,v;
};
int pre[maxn];
bool iscut[maxn];
int bccno[maxn];
int dfs_clock;
int bcc_cnt;
vector<int> g[maxn];
vector<int> 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]=true;
				bcc_cnt++;
				bcc[bcc_cnt].clear();
				while(true){
					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]=false;
	return lowu;
}
void find_bcc(int n){
	memset(pre,0,sizeof(pre));
	memset(iscut,false,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);
	}
}

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

int main(){
	int n,m,a,b;
	while(~scanf("%d%d",&n,&m)&&(n||m)){
		memset(noedge,0,sizeof(noedge));
		for(int i=0;i<n;i++){
			g[i].clear();
			bcc[i].clear();
		}
		while(!s.empty()) s.pop();
		for(int i=0;i<m;i++){
			scanf("%d%d",&a,&b);
			a--;
			b--;
			noedge[a][b]=1;
			noedge[b][a]=1;
		}
		for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
		if(i!=j&&(!noedge[i][j])){ //注意不要把自己和自己连一条边
			g[i].push_back(j);
		}
		find_bcc(n);
		memset(isodd,false,sizeof(isodd));
		for(int i=1;i<=bcc_cnt;i++){
			memset(color,NONE,sizeof(color));
			for(int j=0;j<bcc[i].size();j++) bccno[bcc[i][j]]=i; 
			//因为割顶可以是不同双连通分量的公共点,所以要将割顶编号与双联通分量统一
			int u=bcc[i][0]; //双连通分量的第一个几点用于启动染色
			color[u]=BLACK; //预染色
			if(!bipartite(u,i)){
				for(int j=0;j<bcc[i].size();j++) isodd[bcc[i][j]]=true;
			}
		}
		int ans=n;
		for(int i=0;i<n;i++) if(isodd[i]) ans--;
		printf("%d\n",ans);
	}
	return 0;
}

参考资料:
1.刘汝佳的柏树
2.例题5.5 圆桌骑士 LA3523

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值