hdu1269--求强连通分量个数--tarjan--图的静态邻接表

#include<stdio.h>
#include<string.h>
#include<stack>
using namespace std;
stack<int>s;
int ret;
struct node
{
	int v,next;
} e[100010];
int ins[10010],fang[10010],head[10010],low[10010];
int n,m,yong;
void tarjan(int k)
{
	int j,u;
	fang[k]=low[k]=yong++;
	ins[k]=1;
	s.push(k);//强连通为什么要用栈？栈中其往上的都是可能和他是同一连通分量的，不可以访问非栈上节点是他们不可能在同一连通分量上（访问过而不在栈上，还可能是同一连通分量么），若访问 访问过而不在栈上 的节点(不是同一连通分量，他们早已出栈，说明他们的根还是比较高的，对他没有更新)
	for(j=head[k];j;j=e[j].next)
	{
		u=e[j].v;
		if(fang[u]==0)
		{
			tarjan(u);
			if(low[u]<low[k])
				low[k]=low[u];
		}else if(ins[u]&&fang[u]<low[k])//u可能在上一路上已被访问过，或者u就是k的祖先，他们一定是同一连通分量;与双连通不同的是：他要求的是自己不越过自己，所以，即使他有回边，早晚更新都是一样；双连通担心的是他提前更新，导致他的孩子也提前更新，从而fang[本节点]>low[孩子]=low[本节点]=low[本节点的父节点]
			low[k]=fang[u];
	}
	if(low[k]==fang[k])//其能探到的最低的不越过他自己，说明其以下是一个连通分量
	{
		ret++;
		do{
		j=s.top();
		s.pop();
		ins[j]=0;
		}while(j!=k);
	}
}
int main()
{
	int i,a,b;
	while(1)
	{
		ret=0;
		yong=1;
		scanf("%d%d",&n,&m);
		if(n==0&&m==0)
			break;
		memset(head,0,sizeof(head));
		for(i=1;i<=m;i++)
		{
			scanf("%d%d",&a,&b);
			e[yong].v=b;
			e[yong].next=head[a];
			head[a]=yong;
			yong++;
		}
		yong=1;
		memset(ins,0,sizeof(ins));
		memset(fang,0,sizeof(fang));
		for(i=1;i<=n;i++)
			if(fang[i]==0)
				tarjan(i);
		if(ret==1)
			printf("Yes\n");
		else printf("No\n");
	}
	return 0;
}