算法竞赛——进阶指南——acwing344. 观光之旅 froyd 处理无向图最小环

值得一提的是：有向图找最小环可以用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;
}