[无向图最小环][floyd变形]观光之旅 AcWing344

给定一张无向图，求图中一个至少包含 3 个点的环，环上的节点不重复，并且环上的边的长度之和最小。

该问题称为无向图的最小环问题。

你需要输出最小环的方案，若最小环不唯一，输出任意一个均可。

输入格式

第一行包含两个整数 N 和 M，表示无向图有 N 个点，M 条边。

接下来 M 行，每行包含三个整数 u，v，l，表示点 u 和点 v 之间有一条边，边长为 l。

输出格式

输出占一行，包含最小环的所有节点（按顺序输出），如果不存在则输出 No solution.。

数据范围

1≤N≤100,
1≤M≤10000,
1≤l<500

输入样例：

5 7
1 4 1
1 3 300
3 1 10
1 2 16
2 3 100
2 5 15
5 3 20

输出样例：

1 3 5 2

题意： 给出一个无向图，求至少有3个点的边权和最小的环。

分析： 考虑floyd过程中dis[i][j]的含义，当最外层循环为k且还没更新k时，dis[i][j]表示可能途径前k-1个时从点i到点j的最短路，此时如果dis[i][j]+mp[j][k]+mp[k][i]不为无穷大，那就说明此时找到了一个环，可以发现这样一定能遍历完所有的环，首先任取一个环，当k为其中编号最大的点且k两侧的点分别为i和j时，dis[i][j]就是以i、j为起点和终点且不经过大于等于k号点的最短路，将所有的环按照环上最大点划分集合，之后就可以发现遍历完i和j以后得到的就是环上最大点为k时的答案。于是floyd过程中记录一下路径就可以了，之前记录路径是用的pos[i][j] = pos[i][k]，但是这样只有跑完floyd才能得到正确的最短路，这道题目里需要边更新边处理，所以之前的做法会出错，其实可以用pos[i][j] = k，然后统计路径时再加一个递归就可以了。

具体代码如下：

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <string>
#include <vector>
#define int long long
#define inf 0x3f3f3f3f
using namespace std;

int n, m, dis[105][105], mp[105][105], pos[105][105];
vector<int> path;
 
void get_path(int a, int b)
{
	int k = pos[a][b];
	if(k == 0)
		return;
	get_path(a, k);
	path.push_back(k);
	get_path(k, b);
}
 
signed main()
{
	cin >> n >> m;
	for(int i = 1; i <= n; i++)
		for(int j = 1; j <= n; j++)
			(i==j)?mp[i][j]=0:mp[i][j]=inf;
	for(int i = 1; i <= m; i++)
	{
		int u, v, w;
		scanf("%lld%lld%lld", &u, &v, &w);
		mp[u][v] = mp[v][u] = min(mp[u][v], w);
	} 
	memcpy(dis, mp, sizeof dis);
	int ans = inf;
	for(int k = 1; k <= n; k++)
		for(int i = 1; i <= n; i++)
		{
			for(int j = 1; j <= n; j++)
			{
				if(i != j && j != k && k != i)
				{
					if(ans > dis[i][j]+mp[i][k]+mp[k][j])
					{
						path.clear();
						ans = dis[i][j]+mp[i][k]+mp[k][j];
						path.push_back(i);
						get_path(i, j);
						path.push_back(j);
						path.push_back(k);
					}
				}
			}
			for(int j = 1; j <= n; j++)
				if(dis[i][j] > dis[i][k]+dis[k][j])
				{
					dis[i][j] = dis[i][k]+dis[k][j];
					pos[i][j] = k;
				}
		}
	if(!path.size())
		puts("No solution.");
	else
		for(int i = 0; i < path.size(); i++)
			cout << path[i] << " ";
    return 0;
}