Algo_Tarjan,Scc强联通分量

supimo

已于 2022-03-24 21:57:44 修改

阅读量383

点赞数

分类专栏：计算机知识文章标签：深度优先图论算法

于 2022-03-13 12:01:02 首次发布

本文链接：https://blog.csdn.net/weixin_42712593/article/details/123455940

版权

计算机知识专栏收录该内容

226 篇文章

订阅专栏

catalog

Scc

Scc: Strong connected component 强联通分量: 针对有向图, 将这(有向图), 按照点集进行划分, 可以得到若干个 scc

在这里插入图片描述

该有向图有4个scc

scc的定义为: 一个scc里的任意两个点ab, 均可达 (a可以达到b, b也可以到达a)

若将原有向图, 按照 scc, 进行缩点, 则原图会变成一个 DAG( directed acyclic graph)

在这里插入图片描述

因为, 不同的scc之间, 要么是不可达的, 要么是单向可达; 这正好与scc的定义相悖; 所以, 缩点后是一个dag

我们的目的是: 设置一个scc_id[ ], 对于同一个scc里的点, 他的scc_id都是一样的; 表示, 该点所在的scc 的编号

tarjan算法, 是O(n)的时间

int Scc_count = 0; ' 总共的scc个数; 且每个scc的编号是[0, Scc_count) '
int Scc_id[ PointCount]; ' 每个点所在的scc的 编号 '
int Scc_size[ PointCount]; ' 每个scc里面, 有多少个点 ' 

int Vis_ind = 0;
int Vis_id[ PointCount]; '  tarjan是dfs, 每个点只会进行1次dfs; '
' 每访问一个点, 他的Vis_id就设置为 Vis_ind; 然后Vis_ind ++; 最终, 所有点的Vis_id 均为 [0, n)之间 '

int Stk[ PointCount]; ' 因为是dfs, dfs本身就符合 栈stack结构 '  
int Stk_size = 0;
bool Is_inStk[ PointCount];

void Tarjan_scc( int _node){
	Vis_id[ _node] = Vis_ind;
	Vis_ind ++;
	
	Stk[ Stk_size] = _node;
	Stk_size ++;
	Is_inStk[ _node] = true;
	
	Scc_id[ _node] = Vis_id[ _node];
	
	for( int i = Head[ _node]; i != -1; i = Nex[ i]){
		int j = Ver[ i];
		if( -1 == Vis_id[ j]){
			Tarjan_scc( j);
			if( Is_inStk[ j]){
				Scc_id[ _node] = min( Scc_id[ _node], Scc_id[ j]);
			}
		}
		else{
			if( Is_inStk[ j]){
				Scc_id[ _node] = min( Scc_id[ _node], Scc_id[ j]);
			}
		}
	}
	
	if( Scc_id[ _node] == Vis_id[ _node]){
		while( Stk_size > 0){
			int i = Stk[ Stk_size - 1];
			-- Stk_size;
			Is_inStk[ i] = false;
			
			Scc_id[ i] = Scc_count;
			Scc_size[ Scc_count] ++;
			
			if( i == _node){
				break;
			}
		}
		Scc_count ++;
	}
}

int main(){

	memset( Is_inStk, false, sizeof( Is_inStk));
	memset( Vis_id, -1, sizeof( Vis_id));
	memset( Scc_size, 0, sizeof( Scc_size));
	for( int i = 1; i <= n; ++i){
		if( Vis_id[ i] == -1){
			Tarjan_scc( i);
		}
	}
}

注意点:

外界每次调用的: Tarjan_scc( )之后, 整个Stk栈会是空的!!!
即, 每次Tarjan_scc()完后, 所访问到点, 他们的scc就确立下来了;
而且, 比如, 之前的scc 已经有: 0, 1, 2, 3, 本次Tarjan_scc() 新产生的scc有: 4, 5, 6;
那么, 本次Tarjan_scc()所遍历到的点, 他们所在的scc 一定都是4, 5, 6; 不会是其他的scc
而且, 本次Tarjan_scc()完后, 对于没有图中还没有访问到的点, 他们所的scc, 一定是: > 6的;
dfs顺序是未知的

虽然A所在的scc, 是DAG的起点, 但, 他并不是第一个的tarjan; 第一次的tarjan, 是(D, E, F)
scc的编号

共有3个scc, 我们可以将这3个scc的编号, 设置为: (3) -> (0) -> (1)
即, 每个scc的编号, 为该scc里所有点的最小的vis_id;
这是很容易做到的, 因为根据tarjan算法, 一个scc的确立是取决于 (该scc里最小的vis_id这个点)
但是, 我们通常不会这样做!!! 因为不方便, 虽然知道有3个scc, 但是, 他们的编号, 我们很不方便得到
我们会把每个scc的编号, 设置为: [0, 1, 2]; 即有n个scc, 则每个scc编号为: [0, n)
Is_inStk, 这个标识, 不仅表示: 元素是否在stk里
还表示一个重要含义: 若Is_inStk[ x] == true, 则表示: x已经dfs访问过了, 但还没有求出scc
但是, 如果Is_inStk == false, (可能: x已经求出scc了) (也可能: x还没有dfs访问过)

dfs回溯更新

if( -1 == Vis_id[ j]){
		Tarjan_scc( j);
		if( Is_inStk[ j]){
			Scc_id[ _node] = min( Scc_id[ _node], Scc_id[ j]);
		}
	}

正在进行tarjan算法时, scc_id表示什么呢? 所谓正在进行tarjan, 即if( Scc_id[ _node] == Vis_id[ _node]) 还没有执行
我们知道, 最终的scc_id, 代表该scc的编号; (最先求出来的scc, 编号是0; 依次从0开始递增)
但是, 这是在 (已经求出scc, scc已经确立了)的前提下; 那么, 在还没求出scc, scc_id表示什么呢?
这真的是, tarjan算法的核心;
1. 当访问一个点时, 会有: Scc_id[ _node] = Vis_id[ _node];, 即该点的 scc_id 和其vis_id 是相同的;
2. 然后, 执行dfs; 这个过程, scc_id 可能会变小, 因为只有min操作
3. 当一个scc确立时, 即在执行: if( Scc_id[ _node] == Vis_id[ _node])时, 此时会有:
  比如一个scc, 他所有的点的编号, 是: 3, 5, 6, 7, 假设, vis_id[ 5]最小; 即5相对3/6/7, 是最先dfs访问到的
  根据tarjan, 我们知道, 一个scc 由(vis_id最小的点) 确立, 即当dfs回溯到5时, 会确立这个 scc(3, 5, 6, 7)
  此时当回溯到5, 要确立scc时,会有: scc_id[ 5] = vis_id[ 5]; scc_id[ other] = {>= vis_id[5]}
  从上面的1, 我们知道, scc_id 最初= vis_id; 因为vis_id[ 5]最小, 即vis_id[ other] > vis_id[ 5]; 再根据1, 有: scc_id[ other] > vis_id[ 5]
  但是, 由于2, scc_id 可能会变小; 会变小多少呢? 最小为vis_id[ 5]
  
  当回溯到5时, 此时这个scc就可以确立下来的; 该scc里点的 scc_id, 有0 也是{>0}的; 是挺混乱的
  所以, 确立一个scc, 不是根据一个点的scc_id来确立的!!!; 而是根据stk栈里元素, 具体讲: 当回溯到5时 (知道, 5是scc里最小的vis_id), 则此时 stk栈里, 所有在5前面的元素, 都是该scc里的元素 (这是判定scc的依据)