[HDU-1269] SCC tarjan求强连通分量模板题

   日期:2020-07-10     浏览:99    评论:0    
核心提示:简单说一下tarjan 求SCC的思路:dfn[i]:点i的时间戳low[i]:点i通过i的子树点可达(且有用的)点中时间戳最小的点,x -> y 是有用的:1:x经过横叉边到达y,且从y出发存在一条路径z(z是x的祖先节点);这样x->z->y 形成一个环路2:x经过后向边到达其祖先y,每次经过点x把它加入到栈中。x回溯时判断dfn[x]是否等于low[x].若等于,则把栈中元素弹出到x,与x构成SCC(因为此时其他点一定不会与x构成SCC).._ive a simple directed graph

简单说一下tarjan 求SCC的思路:

dfn[i]:点i的时间戳

low[i]:点i通过i的子树点可达(且有用的)点中时间戳最小的点,

x -> y 是有用的: 

1:x经过横叉边到达y,且从y出发存在一条路径z(z是x的祖先节点);这样x->z->y 形成一个环路

2:x经过后向边到达其祖先y,

每次经过点x把它加入到栈中。

x回溯时判断dfn[x]是否等于low[x].

若等于,则把栈中元素弹出到x,与x构成SCC(因为此时其他点一定不会与x构成SCC)

 

#include <cstdio>
#include<iostream>
#include<algorithm>
#include<vector>
#include<cstring>
using namespace std;
typedef long long ll;
#define ls (o<<1)
#define rs (o<<1|1)
#define pb push_back
const double PI= acos(-1.0);
const int M = 1e5+7;
const int N = 1e5+7;
int head[N],cnt=1;

struct EDGE{int to,nxt,w;}ee[M*2];
void add(int x,int y,int w){ee[++cnt].nxt=head[x],ee[cnt].w=w,ee[cnt].to=y,head[x]=cnt;}

int n,m,num,top,ct;
int sk[N],ins[N],c[N];
int dfn[N],low[N];
vector<int>scc[N];
void tarjan(int x)
{
	dfn[x]=low[x]=++num;
	sk[++top]=x,ins[x]=1;
	//cout<<x<<" -  "<<endl;
	for(int i=head[x];i;i=ee[i].nxt)
	{
		int y=ee[i].to;
		if(!dfn[y])
		{
			tarjan(y);
			low[x]=min(low[x],low[y]);
		}
		else if(ins[y])low[x]=min(low[x],dfn[y]);
		
	}
	if(dfn[x]==low[x])
	{
		ct++;int y;
		do{
			y=sk[top--],ins[y]=0;
			c[y]=ct,scc[ct].pb(y);
		}while(x!=y);
	}
}
void init(){
	cnt=1,memset(head,0,sizeof(head));
	memset(dfn,0,sizeof(dfn));
	num=ct=0;
}
int main()
{
	ios::sync_with_stdio(false);
  	cin.tie(0);
  	while(1)
	{
		cin>>n>>m;
		if(n==0)break;
		init();
		for(int i=1;i<=m;i++)
		{
			int x,y;
			cin>>x>>y;
			add(x,y,1);
		}
		for(int i=1;i<=n;i++)
			if(!dfn[i])tarjan(i);
		if(ct==1)cout<<"Yes"<<endl;
		else cout<<"No"<<endl;
	}
	return 0;
}

 

 
打赏
 本文转载自:网络 
所有权利归属于原作者,如文章来源标示错误或侵犯了您的权利请联系微信13520258486
更多>最近资讯中心
更多>最新资讯中心
0相关评论

推荐图文
推荐资讯中心
点击排行
最新信息
新手指南
采购商服务
供应商服务
交易安全
关注我们
手机网站:
新浪微博:
微信关注:

13520258486

周一至周五 9:00-18:00
(其他时间联系在线客服)

24小时在线客服