HDU tarjan算法模版 强连通分量+最小路径覆盖

注意缩点之后，一个强连通分量内的点，不要在二分图拆点时连边。。

#include<cstdio>
#include<vector>
#include<algorithm>
#include<cstring>
using namespace std;
struct arc
{
    int to,val;
    int next;
    arc(int a=-1,int b=0):to(a),val(b) {}
};
int head[5050];
arc edge[100050];
int e_cnt=0;

int low[5050];//每一点通过儿子,回边可到达的最低深度优先数
int dfn[5050];//每一点的深度优先数
int stack[5050],top=0; //记录深搜的栈
int id[5050]; //每一点所属的强连通分量
int cnt; //记录深度优先数dfn所需递加的变量
int scnt; //强连通分量个数
int n,m; //n为顶点数 m为路径数

bool G[5050][5050]; //二分图
bool vis[5050];
int link[5050];

inline void addedge(int a,int b)
{
    edge[e_cnt].to=b;
    edge[e_cnt].next=head[a];
    head[a]=e_cnt;
    e_cnt++;
}

void dfs(int x)
{
    low[x]=dfn[x]=++cnt;
    stack[++top]=x;
    int v;
    for(int i=head[x];i!=-1;i=edge[i].next)
    {
        v=edge[i].to;
        if(!dfn[v])
        {
            dfs(v);
            low[x]=min(low[v],low[x]);
        }
        else if(!id[v]) low[x]=min(dfn[v],low[x]);
    }
    if(low[x]==dfn[x])
    {
        //printf("%d\n",x);
        scnt++;
        int tmp=0;
        do
        {
            v=stack[top--];
            id[v]=scnt;
        }while(v!=x);
    }
    return;
}

void tarjan()
{
    cnt=scnt=top=0;
    memset(dfn,0,sizeof(dfn));
    for(int i=1;i<=n;i++)
    {
        if(!dfn[i]) dfs(i);
    }
    return;
}

bool find(int u)
{
    for(int i=1;i<=scnt;i++)
    {
        if(!vis[i]&&G[u][i])
        {
            vis[i]=true;
            if(link[i]==-1||find(link[i]))
            {
                link[i]=u;
                return true;
            }
        }
    }
    return false;
}

int slove()
{
    int ans=0;
    memset(link,-1,sizeof(link));
    for(int i=1;i<=scnt;i++)
    {
        memset(vis,0,sizeof(vis));
        if(find(i)) ans++;
    }
    return ans;
}

int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        memset(id,0,sizeof(id));
        memset(head,-1,sizeof(head));
        e_cnt=0;
        scanf("%d%d",&n,&m);
        int a,b;
        for(int i=1;i<=m;i++)
        {
            scanf("%d%d",&a,&b);
            addedge(a,b);
        }
        tarjan();

//        printf("%d\n",scnt);
//        for(int i=1;i<=n;i++)
//        {
//            printf("i=%d id=%d\n",i,id[i]);
//        }

        memset(G,0,sizeof(G));
        for(int i=1;i<=n;i++)
        {
            for(int j=head[i];j!=-1;j=edge[j].next)
            {
                if(id[i]==id[edge[j].to]) continue;
                G[id[i]][id[edge[j].to]]=1;
            }
        }

//        for(int i=1;i<=scnt;i++)
//        {
//            for(int j=1;j<=scnt;j++)
//            {
//                printf("%d ",G[i][j]);
//            }
//            putchar('\n');
//        }
//        putchar('\n');
        printf("%d\n",scnt-slove());
    }
    return 0;
}