「数据结构详解·十二」有向无环图 & 拓扑排序

Abstrummerain

于 2024-04-05 07:51:19 发布

阅读量786

点赞数 18

分类专栏：「数据结构详解」文章标签：数据结构算法 c++ 图论

本文链接：https://blog.csdn.net/leo_chenjy/article/details/137364164

版权

「数据结构详解」专栏收录该内容

12 篇文章 1 订阅

订阅专栏

如果你还没有学过图的基本概念，前往「数据结构详解·九」图的初步。

1. 有向无环图的概念

有向无环图(Directed acyclic graph, DAG)，字面意思，就是一个没有环的有向图。
比如下图就是一个 DAG：

而下图就不是一个 DAG，因为它存在环 $1 - 3 - 5 - 1$ ：

于是，我们就发现，只有父亲指向儿子的边的树是一个 DAG。

DAG 作为特殊的有向图，在各方面都有很重要的作用。
比如，凡是可以进行 dp（动态规划）的数据，其承接关系一定是一个 DAG。

DAG 有一个非常重要的性质：一个 DAG 中，必然存在至少一个顶点的入度为 $0$ ，至少一个顶点的出度为 $0$ 。
下面的拓扑排序便利用了这个性质。
（附：下文中， $\text{din}_i$ 表示节点 $i$ 的入度， $\text{dout}_i$ 表示节点 i 的出度）

2. 拓扑排序的概念及实现

拓扑排序(topological-sort, toposort, topsort)，是对一个 DAG $G$ ，将 $G$ 中所有所有顶点排成一个线性序列，使得图中任意一对顶点 $u, v$ ，若边 $u\rightarrow v∈E(G)$ ，则 $u$ 在线性序列中出现在 $v$ 之前。所以说，它不是对数据的一种排序算法。

我们以 B3644 【模板】拓扑排序 / 家谱树为例。
我们画出样例：

下面我们需要一个队列 $q$ 。
我们首先找到所有入度为 $0$ 的节点，依次加入 $q$ 。比如上述例子中此时 $q=\{2\}$ 。
然后进行循环的操作：

取出队首节点 $x$ ，将所有以它为起点的边删去，也就是对于所有存在的边 $x\rightarrow i$ ，进行 din[i]--；
对于每个进行 din[i]-- 的 $i$ ，看它操作之后的 $\text{din}_i$ 是否为 $0$ ，是则加入队列，然后重复上面的操作。

可以看出来，操作其实就是对每一条有向边 $u\rightarrow v$ ，保证 $u$ 在前 $v$ 在后即可。
最后我们排序出来的一种可能的方案：

参考代码：

#include<bits/stdc++.h>
using namespace std;

vector<int>g[105];
queue<int>q;
int n,din[105];

void toposort()
{
	for(int i=1;i<=n;i++)
	{
		if(din[i]) continue;
		q.push(i);
	}
	while(!q.empty())
	{
		int x=q.front();
		q.pop();
		cout<<x<<' ';
		for(auto i:g[x])
		{
			if(!--din[i]) q.push(i);
		}
	}
}

int main()
{
	cin>>n;
	for(int i=1;i<=n;i++)
	{
		int x;
		cin>>x;
		while(x)
		{
			din[x]++;
			g[i].push_back(x);
			cin>>x;
		}
	}
	toposort();
 	return 0;
}

3. 例题详解

3-1. P4017 最大食物链计数

容易发现，本题需要求的就是一个 DAG 中节点对 $(u, v)$ 满足存在从 $u$ 到 $v$ 的路径且 $\text{din}_u=\text{dout}_v=0$ 。
而根据我们小学学过的加法原理，DAG 中到节点 $p$ 的路径数量一定是其入边上的节点路径数之和。
于是我们就可以拓扑排序的时候 dp。记 $f_i$ 表示到节点 $i$ 的路径数量，状态转移方程为 $f_i=\sum\limits_{x\rightarrow i\in E(G)} f_x$ ，而最后答案就是所有出度为 $0$ 的节点路径数之和，即 $\text{ans}=\sum\limits_{\text{dout}_i=0}f_i$ 。
参考代码：

#include<bits/stdc++.h>
#define mod 80112002
using namespace std;

vector<int>a[5005];
queue<int>q;
int n,m,din[5005],dout[5005],f[5005]; 

void toposort()
{
	while(!q.empty())
	{
		int x=q.front();
		q.pop();
		for(auto i:a[x])
		{
			f[i]=(f[i]+f[x])%mod;
			if(!--din[i]) q.push(i);
		}
	}
}

int main()
{
	cin>>n>>m;
	while(m--)
	{
		int x,y;
		cin>>x>>y;
		a[x].push_back(y);
		dout[x]++;
		din[y]++;
	}
	for(int i=1;i<=n;i++)
	{
		if(!din[i])
		{
			f[i]=1;
			q.push(i);
		}
	}
	toposort();
	int ans=0;
	for(int i=1;i<=n;i++)
	{
		if(!dout[i]) ans=(ans+f[i])%mod;
	}
	cout<<ans;
 	return 0;
}