ACM

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;
}
上一篇下一篇

猜你喜欢

热点阅读