hdu 3836 tarjain 求强连通分量个数

 //  给你一个有向图,问你最少加几条边能使得该图强连通
 #include <iostream>
 #include <cstdio>
 #include <cstring>
 #include <string>
 #include <utility>
 #include <algorithm>
 #include <vector>
 #include <queue>
 #include <stack>
 using namespace std;
 #define max(x,y) x>=y?x:y
 #define lowbit(x) x&(-x)
 typedef  long long ll;
 const int  N=;
 int n,m;
 int head[N],low[N],num[N];
 int  in[N],out[N],be[N];
 int  ans,cnt,ans1,ans2;
 int  index;
 bool vis[N];
 struct Edge
 {
     int from,to,nex;
 }e[N*];
 stack<int>s;
 void init()
 {
     memset(head,-,sizeof(head));
     memset(num,-,sizeof(num));
     memset(vis,,sizeof(vis));
     memset(low,,sizeof(low));
     memset(in,,sizeof(in));
     memset(out,,sizeof(out));
     cnt=index=;
 }
 void add(int u,int v)
 {
     e[cnt].from=u;
     e[cnt].to=v;
     e[cnt].nex=head[u];
     head[u]=cnt++;
 }
 void tarjan(int u)
 {
     low[u]=num[u]=++index;
     vis[u]=;
     s.push(u);
     for (int i=head[u];i!=-;i=e[i].nex)
     {
         int v=e[i].to;
         if(num[v]==-)
         {
             tarjan(v);
             low[u]=min(low[u],low[v]);
         }
         else if(vis[v])
         {
             low[u]=min(low[u],num[v]);
         }
     }
     if(num[u]==low[u])
     {
         int x;
         ans++;
         do{
             x=s.top();
             s.pop();
             vis[x]=;
             be[x]=ans;
         }while(x!=u);
     }
 }
 int main()
 {
     while(~scanf("%d %d",&n,&m))
     {
         init();
         ans=;
         int a,b;
         for (int i=;i<m;i++)
         {
             scanf("%d %d",&a,&b);
             add(a,b);
 
         }
         for (int i=;i<=n;i++)
         {
             if(num[i]==-)
             {
                 tarjan(i);
             }
         }
         for (int i=;i<cnt;i++)
         {
             int a=e[i].from;
             int b=e[i].to;
             if(be[a]!=be[b])
             {
                 in[be[b]]++;
                 out[be[a]]++;
             }
         }
         ans1=ans2=;
         for (int i=;i<=ans;i++)//ans 从1开始。ans为原图强联通分量的个数、
         {
             if(!in[i])
             {
                 ans1++;
             }
             if(!out[i])
             {
                 ans2++;
             }
         }
         printf("%d\n",ans==?:max(ans1,ans2));
     }
     return ;
 }
 
 /*
 
 3 2
 1 2
 2 3
 
 有3个联通快
 
 */