floyd求最小环+记录路径（两种）

问题描述： 给你一张无向图，定义环为从i出发到达j然后从j返回i并且所有点都只经过一次（最少为3个点），求所有环当中经过路径最小的环

算法描述： 首先容易想到的是暴力来枚举环，当删除其中一条边ij后再跑一遍从i到j的最短路，然后加上边i j的值就是含有边i j的最小环的值，这样最坏的时间复杂度可以达到O(n^4)，显然复杂度有点大考虑降低时间复杂度，分析下floyd的实现过程，当枚举顶点k之前我们已经求得了顶点为1  -  k-1 的最短路，所以可以在跟新k之前枚举k之前的i和j的组合，我们可以知道dis[i][j]没有经过k点,所以我们就可以知道如果dis[i][j]+mp[i][k]+mp[k][j] != inf(mp[i][j]为没有更新边值) 时就存在一条经过ijk的最小环,所以我们要求的是所有环当中最小的那个。

Sightseeing trip POJ - 1734

第一种常规记录路径：

#include<iostream>
#include<cstdio>

using namespace std;

const int inf=500;
const int N=105;

int n,m,ans,num;
int dis[N][N],mapp[N][N],pre[N][N],path[N];

inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}

void floyd()
{
    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||k==i||k==j){
                    continue;
                }
                if(dis[i][j]+mapp[i][k]+mapp[k][j]<ans){
                    ans=dis[i][j]+mapp[i][k]+mapp[k][j];
                    num=0;
                    path[num++]=k;
                    int temp=j;
                    while(temp!=i){
                        path[num++]=temp;
                        temp=pre[i][temp];
                    }
                    path[num++]=i;
                }
            }
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                if(i==j||k==i||k==j){
                    continue;
                }
                if(dis[i][j]>dis[i][k]+dis[k][j]){
                    dis[i][j]=dis[i][k]+dis[k][j];
                    pre[i][j]=pre[k][j];
                }
            }
        }
    }
}

int main()
{
    n=read();m=read();
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            pre[i][j]=i;
            if(i==j){
                mapp[i][j]=dis[i][j]=0;
            }else{
                mapp[i][j]=dis[i][j]=inf;
            }
        }
    }
    int u,v,w;
    for(int i=1;i<=m;i++){
        u=read();v=read();w=read();
        mapp[u][v]=mapp[v][u]=dis[u][v]=dis[v][u]=min(w,mapp[u][v]);
    }
    floyd();
    if(ans==inf){
        printf("No solution.\n");
    }else{
        for(int i=num-1;i>=0;i--){
            printf("%d",path[i]);
            if(i==0){
                printf("\n");
            }else{
                printf(" ");
            }
        }
    }
    return 0;
}

第二种：通过dfs回溯出路径

#include<iostream>
#include<cstdio>

using namespace std;

const int inf=500;
const int N=105;

int n,m,ans,num;
int dis[N][N],mapp[N][N],pre[N][N],path[N];

inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}

void get(int x,int y)
{
    int mid=pre[x][y];
    if(!mid){
        path[num++]=y;
        return ;
    }
    get(x,mid);
    get(mid,y);
}

void floyd()
{
    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||k==i||k==j){
                    continue;
                }
                if(dis[i][j]+mapp[i][k]+mapp[k][j]<ans){
                    ans=dis[i][j]+mapp[i][k]+mapp[k][j];
                    num=0;
                    path[num++]=i;
                    get(i,j);
                    path[num++]=k;
                }
            }
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                if(i==j||k==i||k==j){
                    continue;
                }
                if(dis[i][j]>dis[i][k]+dis[k][j]){
                    dis[i][j]=dis[i][k]+dis[k][j];
                    pre[i][j]=k;
                }
            }
        }
    }
    if(ans==inf){
        printf("No solution.\n");
    }else{
        for(int i=0;i<num;i++){
            if(i!=0){
                printf(" ");
            }
            printf("%d",path[i]);
        }
        printf("\n");
    }
}

int main()
{
    n=read();m=read();
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            if(i==j){
                mapp[i][j]=dis[i][j]=0;
            }else{
                mapp[i][j]=dis[i][j]=inf;
            }
        }
    }
    int u,v,w;
    for(int i=1;i<=m;i++){
        u=read();v=read();w=read();
        mapp[u][v]=mapp[v][u]=dis[u][v]=dis[v][u]=min(w,mapp[u][v]);
    }
    floyd();
    return 0;
}