强连通+最小路径覆盖 hdu 3861

题意：有 n 座城市和 m 条有向边，现在要把这 n 座城市分成一些洲，要求： 如果城市 u，v 之间可以互达则两座城市要分在一个洲里面，

        每个洲里面的任意两个城市 u，v之间至少需要存在一条边，问做少能够分成多少个洲。

先将可以互相到达的点缩点，然后要求用最小的边将点连接，所以是最小路径覆盖问题。

#include<stdio.h>
#include<string.h>
#include<stack>
#include<queue>
#include<algorithm>
using namespace std;
#define ll long long
#define me(a,b) memset(a,b,sizeof(a))
struct node
{
    int u,val,next;
} e[100005];
int dfn[5005],low[5005],head[5005],scnt[5005],head2[5005],book[5005],f[5005];
int tot,cnt,top,top2;
void add(int u,int v)
{
    e[top].u=v;
    e[top].next=head[u];
    head[u]=top++;
}
node e2[100005];
void add2(int u,int v)
{
    e2[top2].u=v;
    e2[top2].next=head2[u];
    head2[u]=top2++;
}
stack<int>q;
int n,m;
int trijan(int x)
{
    dfn[x]=low[x]=++tot;
    q.push(x);
    for(int i=head[x];i!=-1;i=e[i].next)
    {
        int u=e[i].u;
        if(!dfn[u])
        {
            trijan(u);
            low[x]=min(low[u],low[x]);
        }
        else if(!scnt[u])low[x]=min(low[x],dfn[u]);
    }
    if(low[x]==dfn[x])
    {
        int v=-1;
        cnt++;
        while(v!=x)
        {
            v=q.top();
            q.pop();
            scnt[v]=cnt;
        }
    }
}
int dfs(int x)
{
    for(int i=head2[x];i!=-1;i=e2[i].next)
    {//printf("ok");
        int u=e2[i].u;
        if(!book[u])
        {
            book[u]=1;
            if(f[u]==-1||dfs(f[u]))
            {
                f[u]=x;
                return 1;
            }
        }
    }
    return 0;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        int n,m;
        me(head,-1);
        me(dfn,0);
        me(low,0);
        me(scnt,0);
        cnt=tot=top=top2=0;
        scanf("%d%d",&n,&m);
        for(int i=0;i<m;i++)
        {
            int a,b;
            scanf("%d%d",&a,&b);
            add(a,b);
        }
        for(int i=1;i<=n;i++)
        {
            if(!dfn[i])
            {
                trijan(i);
            }
        }
        me(head2,-1);
        for(int i=1;i<=n;i++)
        {
            for(int j=head[i];j!=-1;j=e[j].next)
            {
                int u=e[j].u;
                if(scnt[i]!=scnt[u])
                {
                    add2(scnt[i],scnt[u]);
                }
            }
        }
        me(f,-1);
        int ans=cnt;
        for(int i=1;i<=cnt;i++)
        {
            me(book,0);
            ans-=dfs(i);
        }
        printf("%d\n",ans);
    }
}