hdu 2767 Proving Equivalences(强连通入门题）

/*************************************************
     Proving Equivalences(hdu 2767)
     强连通入门题
     给个有向图，求至少加多少条边使得图是所有点都是强连通的
     由a->b->c->a易知n个点至少要n条边，每个出度和入度都要大
     于1。先求所有所有强连通分量，把每个强连通分量看成一个点
     在找每个点的出度和入度，最后还差的出度和入度的最大值就是
     答案。

*************************************************/


#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<stack>
#include<vector>
using namespace std;

const int mx=20005;
stack<int>s;
vector<int>g[mx];
int dfn[mx],low[mx];
int sccno[mx],vs[mx];
int cut_dfs,cut_scc;
int in[mx],out[mx];

void Init(int n)
{
    cut_dfs=cut_scc=0;
    memset(sccno,0,sizeof(sccno));
    memset(dfn,0,sizeof(dfn));
    memset(low,0,sizeof(low));
    memset(vs,0,sizeof(vs));
    for (int i=1;i<=n;i++) g[i].clear();
    while (!s.empty()) s.pop();
}

void dfs(int u)
{

    s.push(u);
    low[u]=dfn[u]=++cut_dfs;
    vs[u]=1;
    for (int i=0;i<g[u].size();i++)
    {
        int v=g[u][i];
        if (!vs[v])
        {
            dfs(v);
            low[u]=min(low[u],low[v]);
        }
        else if (!sccno[v]) low[u]=min(low[u],dfn[v]);///一定要加if (!sccno[v])
    }

    if (low[u]==dfn[u])
    {
        cut_scc++;
        while (1)
        {
            int v=s.top();
            s.pop();
            sccno[v]=cut_scc;
            if (v==u) break;
        }
    }
}

int main()
{
    int t,n,m;
    scanf("%d",&t);
    while (t--)
    {
        scanf("%d%d",&n,&m);
        Init(n);
        int a,b;
        while (m--)
        {
            scanf("%d%d",&a,&b);
            g[a].push_back(b);
        }
        for (int i=1;i<=n;i++)
        {
            if (!dfn[i]) dfs(i);
        }

        if (cut_scc==1)
        {
            printf("0
");
            continue;
        }
        for (int i=1;i<=cut_scc;i++) in[i]=out[i]=1;
        for (int i=1;i<=n;i++)
        {
            for (int j=0;j<g[i].size();j++)
            {
                int v=g[i][j];
                if (sccno[i]!=sccno[v]) out[sccno[i]]=in[sccno[v]]=0;
            }
        }
        a=b=0;
        for (int i=1;i<=cut_scc;i++)
        {
            if (out[i]) a++;
            if (in[i]) b++;
        }
        printf("%d
",max(a,b));
    }
}