DAY_11（强连通分量）

于 2024-01-27 01:04:03 发布
阅读量618 收藏 15
点赞数 21
文章标签： 图论 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/hycccccch/article/details/135877212
版权

求解强连通分量有两种算法，tarjan和kosaraju

首先来俩概念：

强连通：在有向图G中，若两个点u、v是互相可达的，则称u和v是连通的。如果G中任意两个点都是互相可达的，称G是强连通图。

强连通分量：如果一个有向图G不是强连通图，那么可以把它分成多个子图，其中每个子图的内部是强连通的，而且这些子图已经扩展到最大，不能与子图外的任意点强连通，称这样一个“极大强连通”子图是G的一个强连通分量（Strongly Connected Component）

有个常见的问题：求G中有多少个scc？

我们先来尝试着把每个scc抽象成一个个的群岛，群岛由有着很多或很少的小岛组成，群岛内部连通；群岛与其他群岛之间有单向道路连接，且不会形成环路。然后将每一个群岛虚拟成一个个点，最后得到一个有向无环图GG，GG中点的数量即为scc数量（还能推出：每个群岛挖掉后，并不会影响其他岛内部的连通性）

那么一开始想到的一定是暴力！！！（暴力不可取。。。）

暴力思路：对每个点求连通性，然后进行比较，互相连通的点就组成了scc（dfs或bfs可实现）

。。。时间复杂度：O(n^2+m)

下面介绍Kosaraju和Tarjan算法（时间复杂度：O(n+m)）

Kosaraju

原理：

1、将有向图G的所有边反向构建反图rG，rG不会改变原图G的强连通性。（人话：G的scc数量和rG的数量相同）

2、对G和rG各做一次dfs，可以确定scc数量

有一个我之前特别困惑的点：dfs干啥用的？

——对虚拟点进行拓扑排序，排序后将优先级最高的虚拟点进行反dfs，就可以求出被群岛与外界隔离的那些小岛了

模版在下面\forall \forall \forall \forall

Tarjan

Kosaraju是在图中一个一个将scc给挖出来，而Tarjan能在一次dfs就将所有点按照scc分开

其实理解了Kosaraju，来理解Tarjan也不是很难（个人意见）

有几个变量需要注意：

low[i]：表示能返回到的最远祖先

num[i]：表示第几次循环遍历到i（不太重要的样子）

步骤：首先dfs，赋值每个点的low，然后判断low是否与该点的num值相等，如果相等，那么这个点就是自己的祖先点。额。就没了。。（是不是挺简单的~~）

下面上题目：：

B3609 [图论与代数结构 701] 强连通分量

一道模版题，下面依次提供Kosaraju和Tarjan两种做法：

Kosaraju：：

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+5;
vector<int>g[N],scc[N];
stack<int>st;
int n,m;
int f[N];
int low[N];
int num[N];
int dfn,cnt;
void dfs(int u)
{
	dfn++;
	low[u]=num[u]=dfn;
	st.push(u);
	for(int i=0;i<g[u].size();i++)
	{
		int v=g[u][i];
		if(!num[v])
		{
			dfs(v);
			low[u]=min(low[u],low[v]);
		}
		else if(!f[v])
		low[u]=min(low[u],num[v]);
	}
	if(low[u]==num[u])
	{
		cnt++;
		while(!st.empty())
		{
			int k=st.top();
			f[k]=cnt;
			st.pop();
			scc[cnt].push_back(k);
			if(u==k)
			break;
		}
	}
}
int main()
{
	cin>>n>>m;
	for(int i=1;i<=m;i++)
	{
		int u,v;
		cin>>u>>v;
		g[u].push_back(v);
	}
	memset(low,0,sizeof(low));
	memset(num,0,sizeof(num));
	memset(f,0,sizeof(f));
	cnt=dfn=0;
	for(int i=1;i<=n;i++)
	if(!num[i])
	dfs(i);
	cout<<cnt<<endl;
	for(int i=1;i<=n;i++)
	sort(scc[i].begin(),scc[i].end());
	for(int i=1;i<=n;i++)
	if(num[f[i]])
	{
		for(int j=0;j<scc[f[i]].size();j++)
		cout<<scc[f[i]][j]<<" ";
		cout<<endl;
		num[f[i]]=0;
	}
	return 0;
}

Tarjan：：

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+5;
int m,n;
vector<int> g[N];
stack<int> st;
int low[N],num[N],f[N];
int dfn,cnt;
int ans;
int out[N];
int mmp[N];
void dfs(int u)
{
	dfn++;
	low[u]=num[u]=dfn;
	st.push(u);
	for(int i=0;i<g[u].size();i++)
	{
		int v=g[u][i];
		if(!num[v])
		{
			dfs(v);
			low[u]=min(low[u],low[v]);
		}
		else if(!f[v])
		low[u]=min(low[u],num[v]);
	}
	if(low[u]==num[u])
	{
		cnt++;
		while(!st.empty())
		{
			int k=st.top();
			st.pop();
			f[k]=cnt;
			mmp[cnt]++;
			if(u==k)
			break;
		}
	}
}
int main()
{
	cin>>n>>m;
	for(int i=1;i<=m;i++)
	{
		int u,v;
		cin>>u>>v;
		g[u].push_back(v);
	}
	memset(num,0,sizeof(num));
	memset(low,0,sizeof(low));
	memset(f,0,sizeof(f));
	cnt=dfn=0;
	for(int i=1;i<=n;i++)
	if(!num[i])
	dfs(i);
	for(int i=1;i<=n;i++)
	{
		for(int j=0;j<g[i].size();j++)
		if(f[i]!=f[g[i][j]])
		out[f[g[i][j]]]++;
	}
	for(int i=1;i<=cnt;i++)
	if(!out[i])
	ans++;
	cout<<ans<<endl;
}

P2341 [USACO03FALL / HAOI2006] 受欢迎的牛 G

稍微带了点变化，最终是求在只有一个点入度为0的强连通图内点的数目，如果有多个点入度为0，则直接输出0（说明有某点不能连通）

代码：（Kosaraju做法）

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+5;
vector<int> g1[N],g2[N];
stack<int> st;
int n,m;
int in[N];
int f[N];
int cnt;
int mmp[N];
int ans;
int num;
bool vis[N];
void dfs1(int u)
{
	vis[u]=true;
	for(int i=0;i<g1[u].size();i++)
	if(!vis[g1[u][i]])
	dfs1(g1[u][i]);
	st.push(u);
}
void dfs2(int u)
{
	if(vis[u])
	return ;
	f[u]=cnt;
	num++;
	vis[u]=true;
	for(int i=0;i<g2[u].size();i++)
	if(!vis[g2[u][i]])
	dfs2(g2[u][i]);
}
int main()
{
	cin>>n>>m;
	for(int i=1;i<=m;i++)
	{
		int u,v;
		cin>>u>>v;
		g1[u].push_back(v);
		g2[v].push_back(u);
	}
	for(int i=1;i<=n;i++)
	if(!vis[i])
	dfs1(i);
	memset(vis,false,sizeof(vis));
	while(!st.empty())
	{
		int k=st.top();
		st.pop();
		if(!vis[k])
		{
			num=0;
			cnt++;
			dfs2(k);
			mmp[cnt]=num;
		}
	}
	for(int i=1;i<=n;i++)
	for(int j=0;j<g2[i].size();j++)
	if(f[i]!=f[g2[i][j]])
	in[f[g2[i][j]]]++;
	for(int i=1;i<=cnt;i++)
	{
		if(in[i]==0)
		//入度为0的点，f[i]代表编号，需要找f[i]相同的累加 
		{
			if(ans)
			{
				cout<<0<<endl;
				return 0;
			}
			ans=mmp[i];
		}
	}
	cout<<ans<<endl;
	return 0;
}

Tarjan（做法）：

 

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+5;
int m,n;
vector<int> g[N];
stack<int> st;
int low[N],num[N],f[N];
int dfn,cnt;
int ans;
int in[N];
int mmp[N];
void dfs(int u)
{
	dfn++;
	low[u]=num[u]=dfn;
	st.push(u);
	for(int i=0;i<g[u].size();i++)
	{
		int v=g[u][i];
		if(!num[v])
		{
			dfs(v);
			low[u]=min(low[u],low[v]);
		}
		else if(!f[v])
		low[u]=min(low[u],num[v]);
	}
	if(low[u]==num[u])
	{
		cnt++;
		while(!st.empty())
		{
			int k=st.top();
			st.pop();
			f[k]=cnt;
			mmp[cnt]++;
			if(u==k)
			break;
		}
	}
}
int main()
{
	cin>>n>>m;
	for(int i=1;i<=m;i++)
	{
		int u,v;
		cin>>u>>v;
		g[u].push_back(v);
	}
	memset(num,0,sizeof(num));
	memset(low,0,sizeof(low));
	memset(f,0,sizeof(f));
	cnt=dfn=0;
	for(int i=1;i<=n;i++)
	if(!num[i])
	dfs(i);
	for(int i=1;i<=n;i++)
	{
		for(int j=0;j<g[i].size();j++)
		if(f[i]!=f[g[i][j]])
		in[f[i]]++;
	}
	for(int i=1;i<=cnt;i++)
	{
		if(in[i]==0)
		{
			if(ans)
			{
				cout<<0<<endl;
				return 0;
			}
			ans=mmp[i];
		}
	}
	cout<<ans<<endl;
}

P2002 消息扩散

本题最后求出度为0的scc有多少个

代码：（Kosaraju）

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+5;
vector<int> g1[N],g2[N];
stack<int> st;
int n,m;
int in[N],out[N];
int f[N];
int cnt;
int mmp[N];
int ans;
int num;
bool vis[N];
void dfs1(int u)
{
	vis[u]=true;
	for(int i=0;i<g1[u].size();i++)
	if(!vis[g1[u][i]])
	dfs1(g1[u][i]);
	st.push(u);
}
void dfs2(int u)
{
	if(vis[u])
	return ;
	f[u]=cnt;
	num++;
	vis[u]=true;
	for(int i=0;i<g2[u].size();i++)
	if(!vis[g2[u][i]])
	dfs2(g2[u][i]);
}
int main()
{
	cin>>n>>m;
	for(int i=1;i<=m;i++)
	{
		int u,v;
		cin>>u>>v;
		g1[u].push_back(v);
		g2[v].push_back(u);
	}
	for(int i=1;i<=n;i++)
	if(!vis[i])
	dfs1(i);
	memset(vis,false,sizeof(vis));
	while(!st.empty())
	{
		int k=st.top();
		st.pop();
		if(!vis[k])
		{
			num=0;
			cnt++;
			dfs2(k);
			mmp[cnt]=num;
		}
	}
	for(int i=1;i<=n;i++)
	for(int j=0;j<g2[i].size();j++)
	if(f[i]!=f[g2[i][j]])
	out[f[i]]++;
	for(int i=1;i<=cnt;i++)
	if(!out[i])
	ans++;
	cout<<ans<<endl;
	return 0;
}

P2863 [USACO06JAN] The Cow Prom S

群岛内小岛数>1的scc个数

代码：

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+5;
int m,n;
vector<int> g[N];
stack<int> st;
int low[N],num[N],f[N];
int dfn,cnt;
int ans;
int out[N];
int mmp[N];
void dfs(int u)
{
	dfn++;
	low[u]=num[u]=dfn;
	st.push(u);
	for(int i=0;i<g[u].size();i++)
	{
		int v=g[u][i];
		if(!num[v])
		{
			dfs(v);
			low[u]=min(low[u],low[v]);
		}
		else if(!f[v])
		low[u]=min(low[u],num[v]);
	}
	if(low[u]==num[u])
	{
		cnt++;
		while(!st.empty())
		{
			int k=st.top();
			st.pop();
			f[k]=cnt;
			mmp[cnt]++;
			if(u==k)
			break;
		}
	}
}
int main()
{
	cin>>n>>m;
	for(int i=1;i<=m;i++)
	{
		int u,v;
		cin>>u>>v;
		g[u].push_back(v);
	}
	memset(num,0,sizeof(num));
	memset(low,0,sizeof(low));
	memset(f,0,sizeof(f));
	cnt=dfn=0;
	for(int i=1;i<=n;i++)
	if(!num[i])
	dfs(i);
	for(int i=1;i<=cnt;i++)
	if(mmp[i]>1)
	ans++;
	cout<<ans<<endl;
}

P2746 [USACO5.3] 校园网Network of Schools

A任务：入度为0的scc个数

B任务：入度为0和出度为0的scc个数的较大值

代码：

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+5;
int m,n;
vector<int> g[N];
stack<int> st;
int low[N],num[N],f[N];
int dfn,cnt;
int ans1,ans2;
int out[N],in[N];
void dfs(int u)
{
	dfn++;
	low[u]=num[u]=dfn;
	st.push(u);
	for(int i=0;i<g[u].size();i++)
	{
		int v=g[u][i];
		if(!num[v])
		{
			dfs(v);
			low[u]=min(low[u],low[v]);
		}
		else if(!f[v])
		low[u]=min(low[u],num[v]);
	}
	if(low[u]==num[u])
	{
		cnt++;
		while(!st.empty())
		{
			int k=st.top();
			st.pop();
			f[k]=cnt;
			if(u==k)
			break;
		}
	}
}
int main()
{
	cin>>n;
	for(int i=1;i<=n;i++)
	{
		int v;
		while(1)
		{
			cin>>v;
			if(v==0)
			break;
			g[i].push_back(v);
		}
	}
	memset(num,0,sizeof(num));
	memset(low,0,sizeof(low));
	memset(f,0,sizeof(f));
	cnt=dfn=0;
	for(int i=1;i<=n;i++)
	if(!num[i])
	dfs(i);
	for(int i=1;i<=n;i++)
	{
		for(int j=0;j<g[i].size();j++)
		if(f[i]!=f[g[i][j]])
		{
			out[f[i]]++;
			in[f[g[i][j]]]++;
		}
	}
	ans1=0;
	ans2=0;
	for(int i=1;i<=cnt;i++)
	{
		if(!out[i])
		ans2++;
		if(!in[i])
		ans1++;
	}
	cout<<ans1<<endl;
	if(cnt==1)
	cout<<0<<endl;
	else
	cout<<max(ans1,ans2)<<endl;
}

有兴趣建议看看P1073。。。

完结撒花~~~~~~

hycccccch
关注
day 111
BoneInscri的博客
01-11 149
它在以太网媒介系统中的应用属于数据通信中的两种位同步方法里的自同步法（另一种是外同步法），即接收方利用包含有同步信号的特殊编码从信号自身提取同步信号来锁定自己的时钟脉冲频率，达到同步目的。曼彻斯特编码将时钟和数据包含在数据流中，在传输代码信息的同时，也将时钟同步信号一起传输到对方，每位编码中有一跳变，不存在直流分量，因此具有自同步能力和良好的抗干扰性能。例如：传输二进制信息0，若将0看作一位，我们以0为中心，在两边用虚线界定这一位的范围，然后在这一位的中间画出一个电平由高到低的跳变。
day 97
BoneInscri的博客
01-10 249
由于在此指令流水线上执行第1条指令是流水线的建立过程，是顺序执行的，因此从开始执行第1条指令到获得其执行结果所需时间为3个基本操作时间之和，即(1+1+1)△t=3△t。UDP在一个低的水平上完成以上的功能，使用端口号完成进程到进程的通信，但在收到用户数据报时没有流量控制机制，也没有确认，而且只提供有限的差错控制。因此，UDP是一个无连接、不可靠的传输层协议，如果用户应用程序使用UDP协议进行数据传输，必须传输层的上层即应用层提供可靠性方面的全部工作。D.含n个顶点的连通图是一个环，则它有n棵生成树。
求解有向图的强连通分量：Kosaraju算法
个人工作中学到的专业技术与职业发展知识笔记
05-09 1564
1、强连通分量介绍： 在有向图G中，如果两个顶点vi,vj间有一条从vi到vj的有向路径，同时还有一条从vj到vi的有向路径，则称两个顶点强连通(strongly connected)。如果有向图G的每两个顶点都强连通，称G是一个强连通图。有向图的极大强连通子图，称为强连通分量(strongly connected components)。 求解强连通分量的方法如下： 2、Kosaraju算法： 算法基于的事实： 逆图G^T（同图中的每边的方向相反）和原图G有着完全相同的连通分支.也就是说，如果顶点s和t在
强连通分量之Kosaraju算法
weixin_42222609的博客
02-09 657
今天学习了强连通分量的Kosaraju算法，网上写的人也不多，但是跟着视频教程讲解，还有去网上搜了博客，感觉他们的讲解都存在一定的问题，我在学习的时候碰到的一些困惑，他们并没有讲的清楚明白，当然，他们说的大致的思路是正确的。 算法思想： 先说说什么是强连通分量 强连通分量是针对有向图定义的，为了区别无向图的联通分量的概念。在一个强连通分量中，任意两点都相互可达。 这个图中就存在三个强连通分量，分...
强连通分量 Kosaraju算法
白菜
04-23 688
#include "stdio.h" #include "string.h" #define M 20 int map[M][M]; int n, m; int mark[M]; int rank[M]; int pos; int scc; int all[M]; void init(){ int i; int u, v; memset(map, 0, sizeof(map));
【图论】Kosaraju算法详解
360°顺滑
04-11 5483
讲 Kosaraju 算法之前要先知道什么是强连通分量（SCC） 强连通分量 对于一个有向图顶点的子集 S ，如果在 S 内取两个顶点 u 和 v ，都能找到一条从 u 到 v 的路径，那么就称 S 是强连通的。 如果在强连通的顶点集合 S 中加入其他任意顶点集合后，它都不再是强连通的，那么就称 S 是原图的一个强连通分量（SCC:strongly connected componen...
强连通分量Kosaraju 算法（笔记）
2301_79286317的博客
01-26 213
kosaraju
Day11 轮廓线 dp 插头 dp
William_Wang_的博客
08-08 312
插头 dp
【ZSTU4213 2015年12月浙理工校赛 D】【双连通分量tarjan算法】One-Way Roads 无向连通图确定边的方向使得全图任意两点间可达
snowy_smile的博客
12-23 1493
4213: One-Way Roads Time Limit: 1 Sec  Memory Limit: 128 MB  Special Judge Submit: 133  Solved: 45 Description In the ACM kingdom, there are N cities connected by M two-way roads. These citi
python刷题day6.1.rar
04-23
问题可能是找出所有相互认识的人形成的群组，即查找图中的连通分量。可以使用深度优先搜索或者并查集（Disjoint Set）数据结构来解决此类问题。 4. **LRU Cache（最近最久未使用缓存）**： LRU Cache是一种常用的...
Kosaraju算法解决强连通问题
CandyLien
08-29 2163
Kosaraju算法是求解有向图强连通分量（strong connected component）的三个著名算法之一，能在线性时间求解出一个图的强分量。用Sedgewick爷爷的话说，就是“现代算法设计的胜利”。 什么是强连通分量？在这之前先定义一个强连通性（strong
强连通分量——[Kosaraju算法]
20164225的博客
01-19 333
思路：[Kosaraju算法] Kosaraju是基于对有向图及其逆图两次DFS的方法，其时间复杂度也是O(N+M)。kosaraju算法的步骤如下：(1)在有向图G上，从某个顶点出发进行DFS，并按其所有邻接点的搜索都完成(即退出DFS函数)的顺序将顶点排列起来。(2)在有向图G上，从最后完成搜索的顶点出发，在图G的逆图G’上进行DFS，若此次遍历不能访问到有向图中所有顶点，则从余下的顶点中最后
强连通分量（kosaraju）
booyoungxu的专栏
07-30 528
强连通分量特点 1.有向图 2.强连通分量中任意2个点可以互达 3.在一个非强连通图中极大的强连通子图就是该图的强连通分量 4.一个非强连通图可以有多个强连通分量 kosaraju算法 事实：置换图（将图中的边的方向换向）和原始图有相同的强连通分量 描述： 1.首先使G为一个有向图，并且S是一个空栈 2.while S非空：任意选择一个不在栈S中的顶点。从这个顶点v出
强连通分量——Kosaraju算法
Good_Luck_Kevin2018的博客
05-07 423
    要不是一开始看不懂Tarjan算法，上网搜了搜其他算法，发现有好多同学把Kosaraju算法当成Tarjan算法了。看来还是K算法更让人印象深刻啊（除了名字）    另一种Tarjan算法的入口：https://blog.csdn.net/qq_34374664/article/details/77488976    借用郭伟老师的原话：在有向图G中，如果任意两个不同的顶点相互可达，则称该...
【算法】强连通分量分解——Kosaraju算法
tjj1998的博客
05-09 815
算法描述： 在有向图中，我们可能会想知道是否存在一个环，这个环上的任意一个顶点都可以访问同一个环上的其他顶点，这个环上的顶点所组成的集合就叫做一个强连通分量。而将有向图分解成若干个不想交的强连通分量，则叫做强连通分量分解。 算法原理： 1.将有向图存储为G[V]和其反向图rG[V]。 void add_edge(int from,int to){ G[from].push_...
Kosaraju算法求强连通分量
m0_54123999的博客
05-21 357
Kosaraju算法 该算法旨在得到深度优先后续排列后的递归探索中，每次调用DFS的所有顶点都属于同一强联通分量。所以可以这么理解：当递归进入一个强联通分量时，把他锁死在这个强联通分量中（即不能从该强联通分量中的任意顶点达到外部顶点） ...
强连通分量Kosaraju算法
lllllan
08-31 238
优秀博客 [C++图论] 强连通 强连通图的算法 概念   强连通： 在有向图G中，如果两个点u、v是互相可达的，即从u出发可以到达v，从v出发也能到达u，则称u和v是强连通的。如果G中的任意两个点都是互相可达的，称G是强连通图。   强连通分量： 如果一个有向图G不是强连通图，那么可以把他分成多个子图，其中每个子图的内部是强连通的，而且这些子图已经扩展到最大，不能与子图外的任意点强连通，称这样的一个 “极大强连通” 子图是G的一个强连通分量（SCC）。 算法原理 一个有向图G...
Kosaraju算法(强连通分量分解)
肘子的博客
08-12 9639
对于一个有向图顶点的子集S，如果在S内取两个顶点u和v，都能找到一条从u到v的路径，那么就称S是强连通的。如果在强连通的顶点集合S中加入其他任意顶点集合后，它都不再是强连通的，那么就称S是原图的一个强连通分量(SCC:strongly connected component)。任意有向图都可以分解成若不相干的强连通分量，这就是强连通分量分解。把分解后的强连通分量缩成一个顶点，就得到了一个DAG(有...
DAY_OF_WEEK,
最新发布
03-28
DAY_OF_WEEK是一个常见的日期操作，用于获取给定日期的星期几。在不同的编程语言和日期库中，DAY_OF_WEEK的实现方式可能会有所不同，但基本原理是一致的。 在大多数编程语言中，日期通常以特定的格式表示，例如年、月、日。DAY_OF_WEEK函数可以接受一个日期作为输入，并返回该日期对应的星期几。通常，星期几用数字表示，例如1表示星期一，2表示星期二，以此类推。 以下是一个示例，展示了如何使用DAY_OF_WEEK函数获取给定日期的星期几： ```python import datetime def get_day_of_week(date): # 将日期字符串转换为datetime对象 date_obj = datetime.datetime.strptime(date, "%Y-%m-%d") # 使用weekday()方法获取星期几（0表示星期一，1表示星期二，以此类推） day_of_week = date_obj.weekday() return day_of_week # 示例调用 date = "2022-01-01" day_of_week = get_day_of_week(date) print(day_of_week) # 输出：5（表示星期六） ``` 请注意，具体的实现方式可能因编程语言和日期库而异。上述示例是使用Python中的datetime库来实现的。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值