HDU3836 Tarjan缩点

/*
      通过强连通分量Tarjan算法后，进行缩点，然后算强连通点的入度跟出度的分量数目即可。
      答案即为入度为0或出度为0的点数目的最大值。  这点可以将缩点后图形堪称有向无环图，仿照树的定义，根为入度为0的点，叶子为出度为0的点。很明显的，要使其成为强连通图，令f为树根数，g为叶子数，则答案为max(f,g)，特别的，当缩点只有一个点时，答案为0.
*/

/*
 * @author ipqhjjybj
 * @date  20130620
 */
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <ctime>

#include <iostream>
#include <cmath>
#include <algorithm>
#include <numeric>
#include <utility>

#include <cstring>
#include <vector>
#include <queue>
#include <map>
#include <stack>
#include <string>
#include <memory.h>
using namespace std;

#define inf 0x3f3f3f3f
#define MAXN 30000
#define clr(x,k) memset((x),(k),sizeof(x))
#define cpy(x,k) memcpy((x),(k),sizeof(x))
#define Base 10000

#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
int n,m;
struct node{
    int to,next,from;
}Edge[50050];
int head[50050];
int tot;
int low[50050],temp[50050],dfn[50050];
stack<int> st;
int scc_num,lay;
void add(int a,int b){
    Edge[tot].to = b;
    Edge[tot].from = a;
    Edge[tot].next = head[a];
    head[a]=tot++;
}
//构图
void init(){
    int a,b;
    tot = 0;
    clr(head,-1);
    for(int i = 0;i < m;i++){
        scanf("%d %d",&a,&b);
        add(a,b);
    }
}
//Targan算法实现
int Targan(int k,int lay,int& scc_num){
    int j;
    temp[k]=1; dfn[k]=low[k]=lay;
    st.push(k);
    for(int i=head[k];i != -1;i=Edge[i].next){
        if(temp[Edge[i].to]==0){
            Targan(Edge[i].to,++lay,scc_num);
        }
        if(temp[Edge[i].to]==1)
            low[k] = min(low[k],low[Edge[i].to]);
    }
    if(dfn[k]==low[k]){
        ++scc_num;
        do{
            j = st.top();
            low[j]=scc_num;
            temp[j]=2;
            st.pop();
        }while(j!=k);
    }
    return 0;
}
int in[MAXN],out[MAXN];
//缩点
void degreeSet(){
    clr(in,0),clr(out,0);
    for(int i = 0;i < tot;i++){
        if(low[Edge[i].to]!=low[Edge[i].from]){
            //缩点操作
            out[low[Edge[i].to]]++;
            in[low[Edge[i].from]]++;
        }
    }
}
int main(){
    while(scanf("%d %d",&n,&m)!=EOF){
        init();
        clr(low,0);
        clr(temp,0);
        scc_num=0,lay=1;
        st.empty();
        for(int i = 1;i<=n;i++){
            if(temp[i]==0)
                Targan(i,lay,scc_num);
        }

        degreeSet();
        int inNum,outNum;
        inNum = outNum=0;
        for(int i=1;i <= scc_num;i++){
            if(!in[i])
                inNum++;
            if(!out[i])
                outNum++;
        }
        if(scc_num==1){
            printf("0\n");
        }else printf("%d\n",max(inNum,outNum));
    }
    return 0;
}