Description

戰狂也在玩《魔方王國》。他只會徵兵而不會建城市，因此他決定對小奇的城市進行轟炸。
小奇有n 座城市，城市之間建立了m 條有向的地下通道。戰狂會發起若干輪轟炸，每輪可以轟炸任意多個城市。
每座城市裡都有戰狂部署的間諜，在城市遭遇轟炸時，它們會通過地下通道撤離至其它城市。非常不幸的是，在地道里無法得知其它城市是否被轟炸，如果存在兩個不同的城市i，j，它們在同一輪被轟炸，並且可以通過地道從城市i 到達城市j，那麼城市i 的間諜可能因為撤離到城市j 而被炸死。為了避免這一情況，戰狂不會在同一輪轟炸城市i 和城市j。
你需要求出戰狂最少需要多少輪可以對每座城市都進行至少一次轟炸。

Solution

語文題不解釋，全世界應該就我一個看漏了每輪可以轟炸任意座城市。
很明顯屬於同一個強聯通分塊的點是一定不會同時選到的，我們先縮一下點，圖就變成的一個DAG。因為在同一條鏈不能同時選，不在同一條鏈的可以選任意個，很明顯就是DAG的最長鏈了。

Code

#include<iostream>
#include<algorithm>
#include<math.h>
#include<stdio.h>
#include<string.h>
using namespace std;
#define fo(i,a,b) for(i=a;i<=b;i++)
 

#define rep(i,x) for(i=la[x];i;i=ne[i])
const int N=1e6+5;
struct arr{int x,y;}B[N];
int da[N],ne[N],la[N],dfn[N],dis[N],blo[N],S[N],du[N],T[N],ST[N];
int n,m,i,x,y,sum,DFN,cnt,ans;
bool P[N];
int read(){
    char c=getchar();
    for(;'0'>c||c>'9';c=getchar());
    int x=0;
    for(;'0'<=c&&c<='9'
 
;c=getchar())x=x*10+(c-'0');
    return x;
}
bool cmp(arr x,arr y){return x.x<y.x||(x.x==y.x&&x.y<=y.y);}
void ins(int x,int y){da[++sum]=y;ne[sum]=la[x];la[x]=sum;}
void tarjan(int x){
    int i,bz;ST[0]=1;ST[1]=x;
    while(ST[0]){
        x=ST[ST[0]];bz=0;
        if(!P[x]) dis[x]=dfn[x]=++DFN,S[++S[0]]=x,P[x]=1;
        rep(i,x){
            if(!dfn[da[i]]){
                ST[++ST[0]]=da[i];
                bz=1;la[x]=i;
                break;
            }
            else if(P[da[i]]) dis[x]=min(dis[x],dis[da[i]]);
        }
        if(bz) continue;
        if(dfn[x]==dis[x]){
            cnt++;
            while(S[S[0]]!=x){
                P[S[S[0]]]=0;
                blo[S[S[0]]]=cnt;
                T[cnt]++;S[0]--;
            }
            P[x]=0;blo[x]=cnt;T[cnt]++;S[0]--;
        }
        ST[0]--;
    }
}
int main(){
    n=read(),m=read();
    fo(i,1,m) x=read(),y=read(),ins(x,y),B[i]=(arr){x,y};
    fo(i,1,n)if(!dfn[i]) tarjan(i);
    fo(i,1,m) B[i].x=blo[B[i].x],B[i].y=blo[B[i].y];
    memset(la,0,sizeof(la));sum=0;
    fo(i,1,m){
        if(B[i].x==B[i-1].x&&B[i].y==B[i-1].y) continue;
        if(B[i].x==B[i].y) continue;
        ins(B[i].x,B[i].y);du[B[i].y]++;
    }
    memset(dis,0,sizeof(dis));
    S[0]=0;fo(i,1,cnt)if(!du[i]) S[++S[0]]=i,dis[i]=T[i];
    fo(i,1,S[0]){
        x=S[i];
        rep(y,x){
            dis[da[y]]=max(dis[da[y]],dis[x]+T[da[y]]);
            du[da[y]]--;if(!du[da[y]]) S[++S[0]]=da[y];
        }
    }
    fo(i,1,cnt) ans=max(ans,dis[i]);
    printf("%d\n",ans);
}