最小环
问题¶
给出一个图,问其中的有
图的最小环也称围长。
暴力解法¶
设
那么最小环是
总时间复杂度
Dijkstra¶
相关链接:最短路/Dijkstra
枚举所有边,每一次求删除一条边之后对这条边的起点跑一次 Dijkstra,道理同上。
时间复杂度
Floyd¶
相关链接:最短路/Floyd
记原图中
我们注意到 Floyd 算法有一个性质:在最外层循环到点
由最小环的定义可知其至少有三个顶点,设其中编号最大的顶点为
故在循环时对于每个
时间复杂度:
下面给出 C++ 的参考实现:
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | // C++ Version
int val[maxn + 1][maxn + 1]; // 原图的邻接矩阵
inline int floyd(const int &n) {
static int dis[maxn + 1][maxn + 1]; // 最短路矩阵
for (int i = 1; i <= n; ++i)
for (int j = 1; j <= n; ++j) dis[i][j] = val[i][j]; // 初始化最短路矩阵
int ans = inf;
for (int k = 1; k <= n; ++k) {
for (int i = 1; i < k; ++i)
for (int j = 1; j < i; ++j)
ans = std::min(ans, dis[i][j] + val[i][k] + val[k][j]); // 更新答案
for (int i = 1; i <= n; ++i)
for (int j = 1; j <= n; ++j)
dis[i][j] = std::min(
dis[i][j], dis[i][k] + dis[k][j]); // 正常的 floyd 更新最短路矩阵
}
return ans;
}
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | # Python Version
val = [[0 for i in range(maxn + 1)] for j in range(maxn + 1)] # 原图的邻接矩阵
def floyd(n):
dis = [[0 for i in range(maxn + 1)] for j in range(maxn + 1)] # 最短路矩阵
for i in range(1, n + 1):
for j in range(1, n + 1):
dis[i][j] = val[i][j] # 初始化最短路矩阵
ans = inf
for k in range(1, n + 1):
for i in range(1, k):
for j in range(1, i):
ans = min(ans, dis[i][j] + val[i][k] + val[k][j]) # 更新答案
for i in range(1, n + 1):
for j in range(1, n + 1):
dis[i][j] = min(dis[i][j], dis[i][k] + dis[k][j]) # 正常的 floyd 更新最短路矩阵
return ans
|
例题¶
GDOI2018 Day2 巡逻
给出一张
- 删除一个图中的点以及与它有关的边
- 恢复一个被删除点以及与它有关的边
- 询问点
x
对于
对于每一个点
那么枚举所有与
或者直接对每次询问跑一遍 Floyd 求最小环,
对于
还是利用 Floyd 求最小环的算法。
若没有删除,删去询问点将简单环裂开成为一条简单路。
然而第二步的求解改用 Floyd 来得出。
那么答案就是要求出不经过询问点
怎么在线?
强行离线,利用离线的方法来避免删除操作。
将询问按照时间顺序排列,对这些询问建立一个线段树。
每个点的出现时间覆盖所有除去询问该点的时刻外的所有询问,假设一个点被询问
完成之后遍历一遍整棵线段树,在经过一个点时存储一个 Floyd 数组的备份,然后加入被插入在这个区间上的所有点,在离开时利用备份数组退回去即可。
这个做法的时间复杂度为
还有一个时间复杂度更优秀的在线做法。
对于一个对点
那么一定能找出一条非树边,满足这条非树边的两个端点在根的不同子树中,使得这条非树边
证明:
显然最小环包含至少两个端点在根的不同子树中一条非树边。
假设这条边为
那么就可以枚举所有非树边,更新答案。
每次询问的复杂度为跑一次单源最短路的复杂度,为
总时间复杂度为
build本页面最近更新:,更新历史
edit发现错误?想一起完善? 在 GitHub 上编辑此页!
people本页面贡献者:OI-wiki
copyright本页面的全部内容在 CC BY-SA 4.0 和 SATA 协议之条款下提供,附加条款亦可能应用