HDU 3836 Equivalent Sets(强连通分量)

题目链接：

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

题意：

给定一个有向图，问你最少需要添加多少条边，使得整个图成为一个连通图；

分析：

       首先求出所有的强连通图，然后缩点，求每个强连通分量的出度入度，需要加的边数即为出度为零强连通分量的个数、入度为零强连通分量个数，这两个数中的较大者；(下面所说的点都特指缩点后的强连通分量)假如出度为0的点个数大，那么显然，至少需要从这些每个出度为0的点延伸出一条边，来保证从这些出度为0的点能够返回至图的其他位置；如果入度为0的点的数量大，那么同理。

代码：

#include<iostream>
#include<string>
#include<queue>
#include<set>
#include<vector>
#include<stack>
#include<map>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;

const int inf=0x7f7f7f7f;
const int maxn=5e4+50;
const int N=5e4+50;
typedef long long ll;
typedef struct{
    ll u,v,next,w;
}Edge;
Edge e[N];

typedef struct B{
    int l,r;
    ll sum,lazy1,lazy2;
    void update(ll value1,ll value2){
        sum=sum*value1;
        lazy2=lazy2*value1;
        sum=(sum+(r-l+1)*value2);
        lazy1=(lazy1*value1+value2);
    }
}Tree;
Tree tree[4];
int cnt,head[N];

inline void add(int u,int v){
    e[cnt].u=u;
    e[cnt].v=v;
    //e[cnt].w=w;
    // e[cnt].f=f;
    e[cnt].next=head[u];
    head[u]=cnt++;
    // e[cnt].u=v;
    // e[cnt].v=u;
    // e[cnt].w=0;
    // e[cnt].f=-f;
    // e[cnt].next=head[v];
    // head[v]=cnt++;
}

inline int read()
{
    int x = 0;
    int f = 1;
    char c = getchar();
    while (c<'0' || c>'9')
    {    
        if (c == '-')
            f = -1;
        c = getchar();
    }
    while (c >= '0'&&c <= '9')
    {
        x = x * 10 + c - '0';
        c = getchar();
    }
    return x*f;
}

int num,res,n ,m,t,low[N],dfn[N],_time,vis[N],_scc[N],in[N],out[N];
stack<int >s;
void tarjan(int u){
    low[u]=dfn[u]=++_time;
    s.push(u);
    vis[u]=1;
    for(int i=head[u];i!=-1;i=e[i].next){
        int v=e[i].v;
        if(!dfn[v]){
            tarjan(v);
            low[u]=min(low[v],low[u]);
        }
        else if(vis[v])low[u]=min(low[u],dfn[v]);
    }
    if(dfn[u]==low[u]){
        res++;
        do{
            t=s.top();
            s.pop();
            vis[t]=0;
            _scc[t]=res;
        }while(t!=u);
    }

}

int main() {
    while(scanf("%d%d",&n,&m)!=EOF){
        cnt=num=res=_time=0;
        int a,b;
        memset(dfn,0,sizeof(dfn));
        memset(head,-1,sizeof(head));
        memset(in,0,sizeof(in));
        memset(out,0,sizeof(out));
        memset(vis,0,sizeof(vis));
        for(int i=0;i<m;i++){
            cin>>a>>b;
            add(a,b);
        }
        for(int i=1;i<=n;i++)
            if(!dfn[i])
                tarjan(i);
        if(res==1){
            cout<<0<<endl;
            continue;
        }
        for(int i=0;i<cnt;i++){
            int u=e[i].u,v=e[i].v;
            if(_scc[u]!=_scc[v]){
                in[_scc[v]]++;
                out[_scc[u]]++;
            }      
        }
        int cnt1=0,cnt2=0;
        for(int i=1;i<=res;i++){
            if(in[i]==0)cnt1++;
            if(out[i]==0)cnt2++;
        }
        cout<<max(cnt1,cnt2)<<endl;
    }
}   

(仅供个人理解)