值得一提的是:有向图找最小环可以用dij n^2logn做:枚举起点,把起点入堆并处理完起点连向的变后。
把dis[s]设成INF, 最后dis[s]就是答案。
下面考虑这一题:
考虑floyd:dp[i][j][k],只经过前k个点,i到j的最短距离。
考虑循环处理到k时,dp更新到dp[i][j][k-1]。
即只经过了前k-1个点,此时我们让i ->j -> k -> i 形成一个环,满足环的点数大于3,且无重复点。
处理完k等于1-n,就能遍历完所有环。
然后只需要证明这样遍历环能遍历完所有换即可:
对于一个环 x-->……->y->z ->z其中z的编号最大。
以上述遍历,一定会有k==z时,dp[i][j]表示经过k-1个点i->j的最短路, 我们一定会遍历到 x->y 然后经过k回到x形成的环。
即:一定会遍历到这个环,得证!
所以只需要处理floyd的同时,记录最小环即可。
然后加个pos[x][y]表示 x->y最后一次更新经过的点的id。dfs记录即可,
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define ls (o<<1)
#define rs (o<<1|1)
#define pb push_back
const double PI= acos(-1.0);
const int M = 100+7;
/*
int head[M],cnt;
void init(){cnt=0,memset(head,0,sizeof(head));}
struct EDGE{int to,nxt,val;}ee[M*2];
void add(int x,int y){ee[++cnt].nxt=head[x],ee[cnt].to=y,head[x]=cnt;}
*/
int a[M][M];
int d[M][M],pos[M][M];
vector<int>path;
void get_path(int x,int y)
{
if(pos[x][y]==0)//x直达y是当前k的最短路
{
path.pb(y);
return ;
}
get_path(x,pos[x][y]);
get_path(pos[x][y],y);
}
int main()
{
ios::sync_with_stdio(false);
cin.tie(0);
int n,m;
cin>>n>>m;
memset(d,0x3f,sizeof(d));
memset(a,0x3f,sizeof(a));
for(int i=1;i<=m;i++)
{
int u,v,w;
cin>>u>>v>>w;
if(d[u][v]>w)d[u][v]=d[v][u]=a[u][v]=a[v][u]=w;
}
ll ans=0x3f3f3f3f;//最小环的大小
for(int k=1;k<=n;k++)//只经过小于k的点做中转
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
{
if(i==j||i==k||j==k)continue;
if((ll)d[i][j]+a[i][k]+a[k][j]<ans)//比当前最小权值环还小
{
// cout<<"=== "<<i<<" "<<j<<" "<<k<<endl;
ans=(ll)d[i][j]+a[i][k]+a[k][j];
path.clear();
path.pb(i);
get_path(i,j);
path.pb(k);
}
if(d[i][j]>d[i][k]+d[k][j])//i->a->j不如 i->x->k->y->j
{
d[i][j]=d[i][k]+d[k][j];
pos[i][j]=k;//i->j 中间经过k ,在已经经过k-1个点的基础上
}
}
if(path.size()==0)cout<<"No solution."<<endl;
else
{
for(int i=0;i<path.size()-1;i++)
cout<<path[i]<<" ";
cout<<path[path.size()-1]<<endl;
}
return 0;
}