强连通 HDU 3861

t个样例

n个点m条边

分成一些区

2个点互相能到达必须分在一个区

一个区中任何2个点可以u->v 或者v->u

任何点都要有自己的区

求最小的区的数目

强联通缩点

成新图

二分匹配 求最大匹配

最小路径覆盖=点数-最大匹配数

 #include<stdio.h>
 #include<algorithm>
 #include<string.h>
 #include<stack>

 using namespace std;

 #define MAXN 5010
 #define MAXN1 100010
 int head[MAXN],dfn[MAXN],low[MAXN],cou[MAXN],fa[MAXN];
 int cnt,k,num;
 bool vis[MAXN],mark[MAXN];

 struct edg
 {
     int next,to,fr;
 }x[MAXN1];
 void add(int u,int v)
 {
     x[cnt].next=head[u];
     x[cnt].fr=u;
     x[cnt].to=v;
     head[u]=cnt++;
 }
 stack<int>s;

 void dfs(int u)
 {
     low[u]=dfn[u]=k++;
     vis[u]=;
     s.push(u);
     int i;
     for(i=head[u];i!=-;i=x[i].next)
     {
         int v=x[i].to;
         if(!dfn[v])
         {
             dfs(v);
             low[u]=min(low[u],low[v]);
         }
         else if(vis[v])
             low[u]=min(low[u],dfn[v]);
     }
     if(dfn[u]==low[u])
     {
         num++;
         while(!s.empty())
         {
             int now=s.top();
             s.pop();
             vis[now]=;
             fa[now]=num;
             if(now==u)break;
         }
     }
 }
 int pa[MAXN];

 bool dfs1(int u)
 {
     int i;
     for(i=head[u];i!=-;i=x[i].next)
     {
         if(mark[x[i].to])
             continue;
         mark[x[i].to]=;
         if(pa[x[i].to]==-||dfs1(pa[x[i].to]))
         {
             pa[x[i].to]=u;
             return ;
         }
     }
     return ; //很重要
 }
 int main()
 {
     int t;
     scanf("%d",&t);

     while(t--)
     {
         int n,m,i;
         scanf("%d%d",&n,&m);
         cnt=;
         memset(head,-,sizeof(head));
         memset(dfn,,sizeof(dfn));
         memset(low,,sizeof(low));

         for(i=;i<=m;i++)
         {
             int a,b;
             scanf("%d%d",&a,&b);
             add(a,b);
         }
         k=;
         num=;
         for(i=;i<=n;i++) //强联通
             if(!dfn[i])
                 dfs(i);
         memset(head,-,sizeof(head));
         memset(pa,-,sizeof(pa));
         int en=cnt;
         cnt=;
         for(i=;i<en;i++)
         {
             int u,v;
             u=fa[x[i].fr];
             v=fa[x[i].to];
             if(u!=v)
             {
                 add(u,v);
             }
         }
         int ans=;

         for(i=;i<=num;i++) //二分匹配
         {
             memset(mark,,sizeof(mark));
             if(dfs1(i))
                 ans++;
         }
         printf("%d\n",num-ans);
     }

     return ;
 }