POJ1734 AcWing 344. 观光之旅(任意两点间最短路)

题干：

给定一张无向图，求图中一个至少包含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

思路：

题干要求我们求出一个最小环，因为权值不能为负，所以我们可以理解为求从i到j的最短路再加上从j到i的最短路(其中点不能重复)，这样就变成了求任意两点间最短路问题。
首先跑Floyd得出任意两点间的最短路，然后又因为我们需要求这个环的路径。
考虑到Floyd算法是通过中间点来求最短路的(即i->k->j的距离小于i->j的距离)所以我们将每次更新的中间点记录下来，这样每两个点就可能存在一个中间点，我们只需要求记录的时候递归以下就行。
然后再考虑如何求从j->i，方法与Floyd类似(因为本质是求j到i的最短路)，我们还是根据中间点的思想，这里需要注意因为d[i][j]实质是d[k][i][j]表示的是编号不超过k-1的点，从i到j的最短路。所以我们再更新环的时候，只有通过k这个点作为媒介且所以节点的编号都不超过k-1时才能保证不重复尽管只考虑了编号不超过k的节点构成的环，但由于对称性，对结果并不会产生影响。

#include <cstdio>
#include <vector>
#include <cmath>
#include <algorithm>
#include <iostream>
#include <cstring>
using namespace std;
int d[310][310],x[310][310],pos[310][310];
vector<int> q;
void check(int a,int b){
    if(pos[a][b]==0)
        return;
    check(a,pos[a][b]);
    q.push_back(pos[a][b]);
    check(pos[a][b],b);
}
int main()
{
    int n,m,a,b,c;
    scanf("%d%d",&n,&m);
    memset(x,0x3f,sizeof(x));
    memset(pos,0,sizeof(pos));
    for(int i=0;i<=n;i++){
        x[i][i]=0;
    }
    for(int i=0;i<m;i++){
        scanf("%d%d%d",&a,&b,&c);
        x[a][b]=x[b][a]=min(x[a][b],c);
    }
    memcpy(d,x,sizeof(x));
    long long ans=0x3f3f3f3f;
    for(int k=1;k<=n;k++){
        for(int i=1;i<k;i++){ //注意范围需要小于k
            for(int j=i+1;j<k;j++){
                if(ans>(long long)d[i][j]+x[j][k]+x[k][i]){
                    ans=(long long)d[i][j]+x[j][k]+x[k][i];
                    q.clear();
                    q.push_back(i);  
                    check(i,j);  //将i->j中间的点记录
                    q.push_back(j);
                    q.push_back(k);
                }
            }
        }
        for(int i=1;i<=n;i++){  //Floyd算法
            for(int j=1;j<=n;j++){
                if(d[i][j]>d[i][k]+d[k][j]){
                    d[i][j]=d[i][k]+d[k][j];
                    pos[i][j]=k;  //记录中间节点的位置
                }
            }
        }
    }
    if(ans==0x3f3f3f3f)
        printf("No solution.\n");
    else{
        for(int i=0;i<q.size();i++){
            printf("%d ",q[i]);
        }
    }
    return 0;
}