欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页  >  移动技术

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

程序员文章站 2022-09-14 14:52:01
简单说一下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).....

简单说一下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;
}

 

本文地址:https://blog.csdn.net/bjfu170203101/article/details/107208862