POJ(2942)(Knights of the Round T
2018-08-11 本文已影响0人
kimoyami
链接:https://vjudge.net/problem/POJ-2942
思路:本来算是一个多个算法的综合模板题,但是我不熟悉就拿来熟悉模板了,大概就是先用tarjan求出双连通分量,然后利用二分图对每个分量染色,并将能成功染色(必定为奇圈)的顶点标记,最后没标记的点就不能参加会议
代码:
#include<cstdio>
#include<vector>
#include<queue>
#include<algorithm>
#include<cstring>
using namespace std;
const int maxn = 1001;
vector<int> G[1001],bcc[1001];
int odd[1001],color[1001];
int dfn[1001];
int low[1001];
int iscut[1001];
int bccno[1001];
int ntime;
int n,m,bcc_cnt;
struct edge{
int u,v;
edge(){}
edge(int uu,int vv):u(uu),v(vv){}
};
deque<edge> edges;
//tarjan模板
void tarjan(int u,int f){
int i,j,k;
int child = 0;
//low在这里代表非父边所能连到的最早的结点的时间(与求强连通分量时候的low数组意思不一样,一定注意!),dfn表示该结点的开始搜索时间,初始时间dfn等于low
low[u] = dfn[u] = ++ntime;
for(i=0;i<G[u].size();i++){
int v = G[u][i];
if(!dfn[v]){
child++;//统计孩子个数,为1且没有父节点的话则为割顶
edges.push_back(edge(u,v));
tarjan(v,u);
low[u] = min(low[u],low[v]);//更新最早的时间
if(dfn[u]<=low[v]){
edge tmp;
iscut[u] = 1;
//bcc_cnt用来记录双连通分量个数,bcc数组用来记录其中有哪些顶点
bcc_cnt++;
bcc[bcc_cnt].clear();
do{
tmp = edges.back();
edges.pop_back();
if(bccno[tmp.u]!=bcc_cnt){
bcc[bcc_cnt].push_back(tmp.u);
bccno[tmp.u] = bcc_cnt;
}
if(bccno[tmp.v]!=bcc_cnt){
bcc[bcc_cnt].push_back(tmp.v);
bccno[tmp.v] = bcc_cnt;
}
}while(!(tmp.u==u&&tmp.v == v));
}
}
//如果v能连回之前的结点且不是他自己的父节点
else if(dfn[v]<dfn[u]&&v!=f){
edges.push_back(edge(u,v));
low[u] = min(low[u],dfn[v]);
}
}
if(f<0&&child==1)iscut[u] = 0;//标记为割顶
}
void find_bcc(int n){
memset(dfn,0,sizeof(dfn));
memset(iscut,0,sizeof(iscut));
memset(bccno,0,sizeof(bccno));
ntime = bcc_cnt = 0;
for(int i=0;i<n;i++){
if(!dfn[i])tarjan(i,-1);
}
}
//二分图利用dfs染色
bool bipartite(int u,int b){
for(int i=0;i<G[u].size();i++){
int v = G[u][i];
if(bccno[v]!=b)continue;
if(color[v]==color[u])return false;
if(!color[v]){
color[v] = 3-color[u];
if(!bipartite(v,b))return false;
}
}
return true;
}
int A[1001][1001];
int main(){
int kase = 0;
while(scanf("%d%d",&n,&m)&&n){
for(int i=0;i<n;i++)G[i].clear();
memset(A,0,sizeof(A));
for(int i=0;i<m;i++){
int u,v;
scanf("%d%d",&u,&v);
u--;
v--;
A[u][v] = A[v][u] = 1;
}
for(int u=0;u<n;u++){
for(int v=u+1;v<n;v++){
if(!A[u][v]){
G[u].push_back(v);
G[v].push_back(u);
}
}
}
find_bcc(n);
memset(odd,0,sizeof(odd));
for(int i=1;i<=bcc_cnt;i++){
memset(color,0,sizeof(color));
for(int j=0;j<bcc[i].size();j++)
bccno[bcc[i][j]] = i;//处理割顶
int u = bcc[i][0];
color[u] = 1;
if(!bipartite(u,i))
for(int j=0;j<bcc[i].size();j++)
odd[bcc[i][j]] = 1;//可以参加会议
}
int ans = n;
for(int i=0;i<n;i++)if(odd[i])ans--;
printf("%d\n",ans);
}
return 0;
}