【POJ1734】Sightseeing trip

题目链接：https://www.acwing.com/problem/content/346/

题目大意：确定无向带权图上至少包含 3 个节点的最小环

solution

一道无向图上的最小环问题 , 考虑 (Floyd) , 设 (i) 到 (j) 间的道路长为 (f[i][j]) , 最短路径长 (g[i][j]) , 若每次 (Floyd) 最外层 遍历到 (k) 时 , (g[i][j]) 决定了 (i) 到 (j) 间经过节点编号不超过 (k) 的最短路径 , 不妨以 (k) 为环上的节点 , 枚举 (i) 和 (j) , 即是枚举 (k) 在环上的相邻节点 , 那么这个环的最小长度为(f[i][k] + f[k][j] + g[i][j]) , 这样就可以找出环上节点编号不超过 (k) 的最小环 , 由对称性 , 最终可以遍历到每一个最小环

拓展 : 对于有向图上的最小环 , 可以依次枚举起点 (s) , 再在 (dijstra) 中 (s) 向外拓展后 , 把 (shortestpath[s]) 设为极大值 , 最终的 (shortestpath[s]) 即为从 (s) 开始的最小环的长度 , 笔者在这不多加赘述

时间复杂度 : (O(n^3))

code

#include<bits/stdc++.h>
using namespace std;
template <typename T> inline void read(T &FF) {
	int RR = 1; FF = 0; char CH = getchar();
	for(; !isdigit(CH); CH = getchar()) if(CH == '-') RR = -RR;
	for(; isdigit(CH); CH = getchar()) FF = FF * 10 + CH - 48;
	FF *= RR;
}
inline void file(string str) {
	freopen((str + ".in").c_str(), "r", stdin);
	freopen((str + ".out").c_str(), "w", stdout);
}
const int N = 105;
int g[N][N], ans = INT_MAX, f[N][N];
int pos[N][N], n, m;
queue<int> path;
void get_path(int xi, int yi) {
	if(pos[xi][yi] == 0) return;
	get_path(xi, pos[xi][yi]);
	path.push(pos[xi][yi]);
	get_path(pos[xi][yi], yi);
}
int main() {
	//file("");
	int u, v, w;
	read(n), read(m);
	memset(g, 0x3f, sizeof(g));
	memset(f, 0x3f, sizeof(f));
	for(int i = 1; i <= n; i++) f[i][i] = g[i][i] = 0; 
	for(int i = 1; i <= m; i++) {
		read(u), read(v), read(w);
		g[u][v] = g[v][u] = min(g[u][v], w);
		f[u][v] = f[v][u] = g[u][v];
	}
	ans = f[0][0];
	for(int k = 1; k <= n; k++) {
		for(int i = 1; i < k; i++)
			for(int j = i + 1; j < k; j++)
				if((long long)g[i][j] + f[i][k] + f[k][j] < ans) {
					ans = g[i][j] + f[i][k] + f[k][j];
					while(!path.empty()) path.pop();
					path.push(i);
					get_path(i, j);
					path.push(j); path.push(k);
				}
		for(int i = 1; i <= n; i++)
			for(int j = 1; j <= n; j++)
				if(g[i][k] + g[k][j] < g[i][j]) {
					g[i][j] = g[i][k] + g[k][j];
					pos[i][j] = k;
				}
	}
	if(path.empty()) {
		puts("No solution.");
		return 0;
	}
	while(!path.empty()) cout << path.front() << " ", path.pop();
	cout << endl;
	return 0;
}