HDU1269（有向图强连通分量Tarjan算法)

今天学习有向图的强连通分量，学习了Tarjan算法（其实还有Gabow和Kosaraju算法…)，同时，Tarjan对缩点是很有用的。

Tarjan算法是基于dfs的，由于任何一个强连通分量，必定是对原图的深度优先搜索树的子树。那么其实只要确定每个强连通分量的子树的根，然后根据这些根从树的最低层开始，一个一个的拿出强连通分量即可。
两个贯穿算法的数组：
（1)dfn[i]：i结点被首次深度优先搜索到的时间戳；
（2)low[i]：i及其子孙结点中最早搜到的结点的时间;

HDU:1269
题意：是否所有结点强连通；
思路：tarjan跑一遍，看是否强连通分量数为1；
代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<stack>
using namespace std;
#define maxn 10050

int dfn[maxn],low[maxn];
//dfn[i]为i结点第一次搜索到的时间戳(不变);
//low[i]为i及其子孙结点中最早搜到的时间; 
int scc[maxn];//表示结点所对应的强连通分量标号; 
bool vis[maxn];//判断是否在栈中;
stack<int>st;
struct Edge{
	int to;
	int next;
}edge[100050];//链式前向星建图;
int head[maxn];
int cnt,dfn_num,top,scc_num;

//加边操作; 
void add(int u,int v)
{
	edge[++cnt].to=v;
	edge[cnt].next=head[u];
	head[u]=cnt;
}

void tarjan(int u)
{
	//初始化 
	dfn[u]=low[u]=++dfn_num; 
	st.push(u);
	vis[u]=true;//表示在栈中; 
	for(int i=head[u];i!=-1;i=edge[i].next)
	{
		int v=edge[i].to;
		if(!dfn[v])//若该节点还没有被搜索过;
		{
			tarjan(v);
			low[u]=min(low[u],low[v]);//回溯时更新low值; 
		} 
		else if(vis[v])//如果在栈中,更新low值; 
		{
			low[u]=min(low[u],dfn[v]);
		}
	}
	if(low[u]==dfn[u])//为该强连通分量的根结点;
	{
		scc_num++;
		int now;
		do{
			now=st.top();//top的位置为栈中最后一个结点的位置+1;
			st.pop();
			scc[now]=scc_num;
			vis[now]=false;//表示已弹出栈; 
		}while(now!=u);
	 } 	
} 

inline void ini()
{
	memset(vis,false,sizeof(vis));
	memset(head,-1,sizeof(head));
	memset(dfn,0,sizeof(dfn));
	memset(low,0,sizeof(low));
	cnt=0;
	dfn_num=0;
	scc_num=0;
	while(!st.empty())
	{
		st.pop();
	}
}

int main()
{
	int i,j;
	int n,m;
	while(scanf("%d%d",&n,&m)!=EOF)
	{
		if(n==0&&m==0)
		{
			break;
		}
		ini();
		int u,v;
		for(i=1;i<=m;i++)
		{
			scanf("%d%d",&u,&v);
			add(u,v);
		}
		for(i=1;i<=n;i++)//确保在非联通时每个点都能搜到; 
		{
			if(dfn[i]==0)
			{
				tarjan(i);
			}
		}
		//printf("%d\n",scc_num);
		if(scc_num==1)
		{
			printf("Yes\n");
		}
		else
		{
			printf("No\n");
		}
	}
	return 0;
}