HDU 3836 Equivalent Sets(强连通分量)

HDU 3836 Equivalent Sets(强连通分量)

http://acm.hdu.edu.cn/showproblem.php?pid=3836

题意:给你一个有向图,问你最少加几条边能使得该图强连通

分析:本题与HDU 2767 基本一样:

http://blog.csdn.net/u013480600/article/details/31805017

首先求出图的所有强连通分量,然后缩点成DAG图.求出新图中所有点的出度=0与入度=0的节点数分别为a与b.

max(a,b)即为所求需要添加的边.

AC代码:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<stack>
using namespace std;
const int maxn= 20000+10;
int n,m;
vector<int> G[maxn];
stack<int> S;
int dfs_clock,scc_cnt;
int pre[maxn],low[maxn],sccno[maxn];
int in[maxn],out[maxn];
void dfs(int u)
{
    pre[u]=low[u]=++dfs_clock;
    S.push(u);
    for(int i=0;i<G[u].size();i++)
    {
        int v=G[u][i];
        if(!pre[v])
        {
            dfs(v);
            low[u]=min(low[u],low[v]);
        }
        else if(!sccno[v])
            low[u]=min(low[u],pre[v]);
    }
    if(low[u]==pre[u])
    {
        scc_cnt++;
        while(true)
        {
            int x=S.top(); S.pop();
            sccno[x]=scc_cnt;
            if(x==u) break;
        }
    }
}
void find_scc(int n)
{
    dfs_clock=scc_cnt=0;
    memset(pre,0,sizeof(pre));
    memset(sccno,0,sizeof(sccno));
    for(int i=1;i<=n;i++)
        if(!pre[i]) dfs(i);
}
int main()
{
    while(scanf("%d%d",&n,&m)==2)
    {
        for(int i=1;i<=n;i++) G[i].clear();
        while(m--)
        {
            int u,v;
            scanf("%d%d",&u,&v);
            G[u].push_back(v);
        }
        find_scc(n);

        for(int i=1;i<=scc_cnt;i++) in[i]=out[i]=0;
        for(int u=1;u<=n;u++)
        for(int i=0;i<G[u].size();i++)
        {
            int v=G[u][i];
            int x=sccno[u], y=sccno[v];
            if(x!=y)
            {
                in[y]++;
                out[x]++;
            }
        }
        int a=0,b=0;
        for(int i=1;i<=scc_cnt;i++)
        {
            if(!in[i]) a++;
            if(!out[i]) b++;
        }
        printf("%d\n",scc_cnt==1?0:max(a,b));
    }
    return 0;
}