floyd求最短环路径 及其路径保存 //1494：【例 1】Sightseeing Trip

传送门

再次刷新了我对floyd的认知，因为涉及到了可模板部分，单独放一篇。

//floyd求最短环路径 及其路径保存 
#include<bits/stdc++.h>
using namespace std;
const int N = 310; 
const int INF = 0x3f3f3f3f;
typedef long long ll;
int n,m,ans = INF;
int g[N][N],dis[N][N],pos[N][N];
vector<int> path;
void read(){
	memset(g,0x3f,sizeof g);
	scanf("%d%d",&n,&m);
	for(int i=1;i <= n;i++) g[i][i] = 0;
	for(int i = 1;i <= m;i ++){
		int x,y,w;scanf("%d%d%d",&x,&y,&w);
		//g[x][y] = g[y][x] = min(g[x][y],w);
		g[x][y] = g[y][x] = w;
	}	
	memcpy(dis,g,sizeof g);//这里因为求环，所以需要原路径 
}
void get_path(int i,int j){
	if(pos[i][j] == 0) return;
	int k = pos[i][j]; //注意k保存的只是i,j的中间节点，递归插入 
	get_path(i,k);
	path.push_back(k);
	get_path(k,j);
}
void solve(){
	for(int k = 1; k <= n; k++){
		for(int i = 1;i <k;i++) //因为是先左判断，此时第k条路径还没有加入判断，所以是k-1 
		for(int j = i+1;j < k;j++)//只用判断前k节点，后面根据对称性可知无影响 
			if(ans > (ll)dis[i][j] + g[j][k] + g[k][i]){
				ans = dis[i][j] + g[j][k] + g[k][i];
				path.clear(); //找到新的最小环 
				path.push_back(i);
				get_path(i,j);
				path.push_back(j);
				path.push_back(k);
			}
		for(int i=1;i <= n;i ++)
		for(int j=1;j <= n;j ++)
		if(dis[i][j] > dis[i][k] + dis[k][j]){ //floyd的路径存贮方式 
			dis[i][j] = dis[i][k] + dis[k][j];
			pos[i][j] = k;
		}
 	}
}
int main(){
	read();
	solve();
	if(ans == INF){
		printf("No solution.");
		return 0;
	}
	for(int i = 0;i < path.size();i++){
		printf("%d ",path[i]);
	}
	return 0;
}