HDU-3836 Equivalent Sets（Tarjan+强连通缩点）

题意

给定一个 $n$ 个节点的有向图，求最少加几条边，可以使得整个图强连通。
$1 \leq n \leq 20000$
$1 \leq m \leq 50000$

思路

不难发现出入度和强连通间有一个可否的关系，当一个有向图如果把边都变成双向时是强连通图且图中的每个点入度出度均不为零时，这张有向图图强连通。特别的，只有一个节点的图也是强连通图。
对于这道题，首先要考虑的是强连通子图的问题，不难发现，原本的强连通子图需要缩成一点，比如下图：

不难发现，由于对应的无向图不是强连通，所以这张图不是强连通。所以，$\{1,2,3\}$ 和 $\{4,5,6\}$ 都有必要“缩点”。缩点后只要把对应的 $DAG$ 强连通即可。注意到增加一条边的贡献是将起点出度加一，将终点入度加一。所以，要让 $DAG$ 上所有点的入度出度均不为 $0$ ，所需的边数最少是入度为 $0$ 的点的个数与出度为 $0$ 的点的个数之 $max$ 。仍然需要注意的是只有一个点的缩点图就没有必要加边了。

代码

#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define FOR(i,x,y) for(int i=(x);i<=(y);i++)
#define DOR(i,x,y) for(int i=(x);i>=(y);i--)
#define N 20003
#define M 50003
typedef long long LL;
using namespace std;
template<const int maxn,const int maxm>struct Linked_list
{
    int head[maxn],to[maxm],nxt[maxm],tot;
    void clear(){memset(head,-1,sizeof(head));tot=0;}
    void add(int u,int v){to[++tot]=v,nxt[tot]=head[u];head[u]=tot;}
    #define EOR(i,G,u) for(int i=G.head[u];~i;i=G.nxt[i])
};
Linked_list<N,M>G;
int dfn[N],low[N],bel[N],stk[N],ind[N],oud[N];
int ord,n,m,B,tp,ans;

void tarjan(int u)
{
    dfn[u]=low[u]=++ord;
    stk[++tp]=u;
    EOR(i,G,u)
    {
        int v=G.to[i];
        if(!dfn[v])
        {
            tarjan(v);
            low[u]=min(low[u],low[v]);
        }else if(!bel[v])low[u]=min(low[u],dfn[v]);
    }
    if(dfn[u]==low[u])
    {
        B++;int v;
        do
        {
            v=stk[tp--];
            bel[v]=B;
        }while(u!=v);
    }
}

void clear()
{
    G.clear();
    memset(dfn,0,sizeof(dfn));
    memset(bel,0,sizeof(bel));
    memset(ind,0,sizeof(ind));
    memset(oud,0,sizeof(oud));
    ord=ans=B=tp=0;
}

int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        clear();
        FOR(i,1,m)
        {
            int u,v;
            scanf("%d%d",&u,&v);
            G.add(u,v);
        }
        FOR(i,1,n)if(!dfn[i])tarjan(i);
        int I=0,O=0;
        FOR(u,1,n)
            EOR(i,G,u)
            {
                int v=G.to[i];
                if(bel[u]==bel[v])continue;
                oud[bel[u]]++,ind[bel[v]]++;
            }
        if(B==1)printf("0\n");
        else
        {
            FOR(i,1,B)
            {
                if(!ind[i])I++;
                if(!oud[i])O++;
            }
            printf("%d\n",max(I,O));
        }
    }
    return 0;
}