C++ 强连通分量 - 缩点（Tarjan算法）

最新推荐文章于 2026-01-20 13:32:17 发布

原创最新推荐文章于 2026-01-20 13:32:17 发布 · 1.3k 阅读

2 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#算法 #图论 #c++ #1024程序员节

强连通分量专栏收录该内容

2 篇文章

订阅专栏

本文介绍了Tarjan算法用于求解有向图的强连通分量，强调了算法原理和在洛谷P3387题中的应用。通过深度优先搜索，利用dfn和low值判断强连通分量，并提供了代码实现。

分析：需要找出一条点权最大的路径就行了，不限制点的个数。考虑对于一个环上的点被选择了，那么一整条环都被选择。很关键的是题目还允许重复经过某条边或者某个点，因此整个环实际上可以看成一个点（选了其中一个点就应该选其他的点）。

显然，是一道缩点 + 拓扑排序

既然缩点，就不得不提一下 Tarjan 算法：

Tarjan 算法是基于对图深度优先搜索的算法，每个强连通分量为搜索树中的一棵子树。搜索时，把当前搜索树中未处理的节点加入一个堆栈，回溯时可以判断栈顶到栈中的节点是否为一个强连通分量。
定义 dfn[u] 为节点 u 搜索的次序编号(时间戳)，low[u] 为 u 或 u 的子树能够追溯到的最早的栈中节点的次序号。

当 dfs[u]=low[u] 时，以 u 为根的搜索子树上所有节点是一个强连通分量。

【代码实现】

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<map>
#include<vector>
#include<string>
#include<cstring>
#include<queue>
#include<stack>
using namespace std;
inline int read()//快读
{
    int num=0,w=1;
    char ch=getchar();
    while(ch<'0'||ch>'9') {if(ch=='-') w=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') {num=(num<<1)+(num<<3)+ch-'0';ch=getchar();}
    return num*w;
}
int n,m,ans;
int a[1000001];
int head[1000001],tot;
int dis[1000001];
int h[1000001],lol;
struct node//链式前向星
{
	int from;
	int to;
	int nex;
}g[1000001],mapp[1000001];//分别存储缩点前后的图
int u,v;
int dfn[1000001],cnt;
int low[1000001];
int stac[1000001],top;//记录环中的点
int fa[1000001],sum;
bool vis[1000001];
int ru[1000001];
inline void add_1(int xx,int yy)
{
	g[++tot]=(node){xx,yy,head[xx]};
	head[xx]=tot;
	return;
}
inline void add_2(int xx,int yy)
{
	mapp[++lol]=(node){xx,yy,h[xx]};
	h[xx]=lol;
	return;
}
inline void tanjor(int k)//Tarjan 算符
{
	dfn[k]=low[k]=++cnt;
	stac[++top]=k;
	vis[k]=1;
	for(int i=head[k];i;i=g[i].nex)
	{
		int to=g[i].to;
		if(!dfn[to]) tanjor(to),low[k]=min(low[k],low[to]);
		else if(vis[to]) low[k]=min(low[k],dfn[to]);
	}
	if(low[k]==dfn[k])
	{
		int qwq;
		while((qwq=stac[top--]))
		{
			vis[qwq]=0;
			fa[qwq]=k;
			if(qwq==k) break;
			a[k]+=a[qwq];
		}
	}
	return;
}
queue<int> q;
inline void tuopusort()//拓扑排序
{
	for(int i=1;i<=n;i++) if(fa[i]==i&&ru[i]==0) q.push(i),dis[i]=a[i];
	while(!q.empty())
	{
		int k=q.front();
		q.pop();
		for(int i=h[k];i;i=mapp[i].nex)
		{
			int to=mapp[i].to;
			ru[to]--;
			dis[to]=max(dis[to],dis[k]+a[to]);
			if(ru[to]==0) q.push(to);
		}
	}
	for(int i=1;i<=n;i++) ans=max(ans,dis[i]);
	return;
}
int main()
{
	//freopen(".in","r",stdin);
	//freopen(".out","w",stdout);
	n=read(),m=read();
	for(int i=1;i<=n;i++) a[i]=read();
	for(int i=1;i<=m;i++) u=read(),v=read(),add_1(u,v);
	for(int i=1;i<=n;i++) if(!dfn[i]) tanjor(i);
	for(int i=1;i<=m;i++) {u=fa[g[i].from],v=fa[g[i].to];if(u!=v) add_2(u,v),ru[v]++;}
	tuopusort();
	printf("%d\n",ans);
	//fclose(stdin);
	//fclose(stdout);
	return 0;//完结撒花
}

又又又是一道板子题！！！