POJ 2553 The Bottom of a Graph

题意：

定义一个点所能到达的所有点都能到达这个点，那么该点就是一个sink点，现给你n个点m条边的有向图，问你图中有多少个sink点，并指出对应的点。

思路：

首先明白一点，强连通分量里的所有点都是相互可达的，里面的点和分量外面的点都是只有一方可达的，所有我们只要用Tarjan缩点，然后对于出度为0的强连通分量，里面的所有点就是我们所求的点。

#include<cstdio>
#include<vector>
#include<stack>
#include<cstring>
#include<queue>
#include<algorithm>
using namespace std;

const int MAXN=5005;
int n,m,time,cnt,dfn[MAXN],low[MAXN],id[MAXN],out[MAXN];
bool vis[MAXN],use[MAXN][MAXN];
vector<int> gr[MAXN],se[MAXN],ans;
stack<int> s;

void init(){
	for(int i=1;i<=n;i++){
		gr[i].clear();
		se[i].clear();
	}
	while(!s.empty()){
		s.pop();
	}
	ans.clear();
	memset(dfn,0,sizeof(dfn));
	memset(low,0,sizeof(low));
	memset(out,0,sizeof(out));
	memset(vis,false,sizeof(vis));
	memset(use,false,sizeof(use));
	time=0;
	cnt=0;
}

void Tarjan(int u){
	int v;
	vis[u]=true;
	dfn[u]=low[u]=++time;
	s.push(u);
	for(int i=0;i<gr[u].size();i++)
	{
		v=gr[u][i];
		if(!dfn[v]){
			Tarjan(v);
			low[u]=min(low[u],low[v]);
		}
		else if(vis[v]){
			low[u]=min(low[u],dfn[v]);
		}
	}
	if(low[u]==dfn[u]){
		cnt++;
		while(!s.empty())
		{
			v=s.top();
			s.pop();
			vis[v]=false;
			se[cnt].push_back(v);
			id[v]=cnt;
			if(u==v) break;
		}
	}
}

int main(){
	while(~scanf("%d",&n),n){
		scanf("%d",&m);
		init();
		int u,v;
		for(int i=1;i<=m;i++){
			scanf("%d%d",&u,&v);
			gr[u].push_back(v);
		}
		for(int i=1;i<=n;i++){
			if(!dfn[i]){
				Tarjan(i);
			}
		}
		for(int i=1;i<=n;i++){
			for(int j=0;j<gr[i].size();j++){
				int v=gr[i][j];
				if(use[id[i]][id[v]]==false&&id[i]!=id[v]){
					use[id[i]][id[v]]=true;
					out[id[i]]++;
				}
			}
		}
		for(int i=1;i<=cnt;i++){
			if(!out[i]){
				for(int j=0;j<se[i].size();j++){
					int v=se[i][j];
					ans.push_back(v);
				}
			}
		}
		sort(ans.begin(),ans.end());
		for(int i=0;i<ans.size();i++){
			if(i!=0) printf(" ");
			printf("%d",ans[i]);
		}
		printf("\n");
	}
	return 0;
}