【强连通分量缩点】【拓扑排序】【dp预处理】CDOJ1640 花自飘零水自流，一种相思，两处闲愁。

#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
#include<vector>
using namespace std;
typedef long long ll;
vector<int>G[100010],rG[100010],vs;
bool used[100010];
int cmp[100010],x[100010],y[100010],a[100010],f[100010],g[100010],ru[100010],ru2[100010];
int n,m,K;
void dfs(int U){
	used[U]=1;
	for(int i=0;i<G[U].size();++i){
		if(!used[G[U][i]]){
			dfs(G[U][i]);
		}
	}
	vs.push_back(U);
}
void rdfs(int U){
	used[U]=1;
	cmp[U]=K;
	for(int i=0;i<rG[U].size();++i){
		if(!used[rG[U][i]]){
			rdfs(rG[U][i]);
		}
	}
}
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;++i){
		scanf("%d%d",&x[i],&y[i]);
		G[x[i]].push_back(y[i]);
		rG[y[i]].push_back(x[i]);
	}
	for(int i=1;i<=n;++i){
		if(!used[i]){
			dfs(i);
		}
	}
	memset(used,0,sizeof(used));
	for(int i=vs.size()-1;i>=0;--i){
		if(!used[vs[i]]){
			++K;
			rdfs(vs[i]);
		}
	}
	for(int i=1;i<=n;++i){
		G[i].clear();
		rG[i].clear();
	}
	for(int i=1;i<=m;++i){
		if(cmp[x[i]]!=cmp[y[i]]){
			G[cmp[x[i]]].push_back(cmp[y[i]]);
			++ru[cmp[y[i]]];
			rG[cmp[y[i]]].push_back(cmp[x[i]]);
			++ru2[cmp[x[i]]];
		}
	}
	for(int i=1;i<=n;++i){
		++a[cmp[i]];
	}
	queue<int>q;
	for(int i=1;i<=K;++i){
		if(!ru[i]){
			q.push(i);
		}
	}
	memset(f,0xaf,sizeof(f));
	f[cmp[1]]=a[cmp[1]];
	while(!q.empty()){
		int U=q.front(); q.pop();
		for(int i=0;i<G[U].size();++i){
			f[G[U][i]]=max(f[G[U][i]],f[U]+a[G[U][i]]);
			--ru[G[U][i]];
			if(!ru[G[U][i]]){
				q.push(G[U][i]);
			}
		}
	}
	
	for(int i=1;i<=K;++i){
		if(!ru2[i]){
			q.push(i);
		}
	}
	memset(g,0xaf,sizeof(g));
	g[cmp[1]]=a[cmp[1]];
	while(!q.empty()){
		int U=q.front(); q.pop();
		for(int i=0;i<rG[U].size();++i){
			g[rG[U][i]]=max(g[rG[U][i]],g[U]+a[rG[U][i]]);
			--ru2[rG[U][i]];
			if(!ru2[rG[U][i]]){
				q.push(rG[U][i]);
			}
		}
	}
	
	int ans=a[cmp[1]];
	for(int i=1;i<=m;++i){
		if(cmp[x[i]]!=cmp[y[i]]){
			ans=(int)max((ll)ans,(ll)f[cmp[y[i]]]+(ll)g[cmp[x[i]]]-(ll)a[cmp[1]]);
		}
	}
	printf("%d\n",ans);
	return 0;
}