poj 1236 （强连通）

#include<cstdio>
#include<cstring>
#define MAX_N 110
#define MAX(x,y) ((x)>(y)?(x):(y))
using namespace std;
int map[MAX_N][MAX_N],rmap[MAX_N][MAX_N],vis[MAX_N],stack[MAX_N],cmp[MAX_N];
int tot=0,n;
void dfs(int i)
{
	vis[i]=1;
	for(int j=1;j<=n;j++)
	{
		if(!vis[j]&&map[i][j])
		dfs(j);
	}
	stack[tot++]=i;
}
void rdfs(int i,int k)
{
	vis[i]=1;
	cmp[i]=k;
	for(int j=1;j<=n;j++)
	{
		if(!vis[j]&&rmap[i][j])
		rdfs(j,k);
	}
}
int scc()
{
	memset(vis,0,sizeof(vis));
	for(int i=1;i<=n;i++)
	{
		if(!vis[i])
		dfs(i);
	}
	int k=1;
	memset(vis,0,sizeof(vis));
	for(int i=tot-1;i>=0;i--)
	{
		if(!vis[stack[i]])
		{
		vis[stack[i]]=1; 
		rdfs(stack[i],k++);
	    }
	}
	return k;
}
int main()
{
	int j,k,in[MAX_N],out[MAX_N];
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	{
	    while(~scanf("%d",&j)&&j)
		{
			map[i][j]=1;
			rmap[j][i]=1;
		}	
	}
	k=scc();
	memset(in,0,sizeof(in));
	memset(out,0,sizeof(out)) ;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			if(map[i][j]&&cmp[i]!=cmp[j])
			{
				in[cmp[j]]++;
				out[cmp[i]]++;
			}
		}
	}
	int t1=0,t2=0;
	for(int i=1;i<k;i++)
	{
		if(!in[i])
		t1++;
		if(!out[i])
		t2++;
	}
	if(k==2)
	printf("1\n0\n");
	else 
	printf("%d\n%d\n",t1,MAX(t1,t2));
}