题目链接:Click here~~
题意:
给一个无向图,求出图中的最小环。
解题思路:
算法分析摘自别人的博客。
floyd求最小环:
抛开Dijkstra算法,进而我们想到用Floyd算法。我们知道,Floyd算法在进行时会不断更新矩阵dist(k)。设dist[k,i,j]表示从结点i到结点j且满足所有中间结点,它们均属于集合{1,2,⋯ ,k}的一条最短路径的权。其中dist[0,i,j ]即为初始状态i到j的直接距离。对于一个给定的赋权有向图, 求出其中权值和最小的一个环。我们可以将任意一个环化成如下形式:u->k->v ->(x1-> x2-> ⋯ xm1)-> u(u与k、k与v都是直接相连的),其中v ->(x1-> 2-> ⋯ m)-> u是指v到u不经过k的一种路径。
在u,k,v确定的情况下,要使环权值最小, 则要求 (x1一>x2->⋯一>xm)->u路径权值最小.即要求其为v到u不经过k的最短路径,则这个经过u,k,v的环的最短路径就是:[v到u不包含k的最短距离]+dist[O,u,k]+dist[O,k,v]。我们用Floyd只能求出任意2点间满足中间结点均属于集合{1,2,⋯ ,k}的最短路径,可是我们如何求出v到u不包含k的最短距离呢?
现在我们给k加一个限制条件:k为当前环中的序号最大的节点(简称最大点)。因为k是最大点,所以当前环中没有任何一个点≥k,即所有点都<k。因为v->(x1->x2->......xm)->u属于当前环,所以x1,x2,⋯ ,xm<k,即x1,x2.⋯。xm≤k一1。这样,v到u的最短距离就可以表示成dist[k一1 ,u,v]。dist[k一1,v,u]表示的是从v到u且满足所有中间结点均属于集合{1,2,⋯ ,k一1}的一条最短路径的权。接下来,我们就可以求出v到u不包含k的最短距离了。这里只是要求不包含k,而上述方法用的是dist[k一1,v,u],求出的路径永远不会包含k+l,k+2,⋯ 。万一所求的最小环中包含k+1,k+2,⋯ 怎么办呢?的确,如果最小环中包含比k大的节点,在当前u,k,v所求出的环显然不是那个最小环。然而我们知道,这个最小环中必定有一个最大点kO,也就是说,虽然当前k没有求出我们所需要的最小环,但是当我们从k做到kO的时候,这个环上的所有点都小于kO了.也就是说在k=kO时一定能求出这个最小环。我们用一个实例来说明:假设最小环为1—3—4—5—6—2—1。的确,在u=l,v=4,k=3时,k<6,dist[3,4,1]的确求出的不是4—5—6—2—1这个环,但是,当u=4,v=6,k=5或u=5,v=2,k=6时,dist[k,v,u]表示的都是这条最短路径.所以我们在Floyd以后,只要枚举u.v,k三个变量即可求出最小环。时间复杂度为O(n3)。我们可以发现,Floyd和最后枚举u,v,k三个变量求最小环的过程都是u,v,k三个变量,所以我们可以将其合并。这样,我们在k变量变化的同时,也就是进行Floyd算法的同时,寻找最大点为k的最小环。
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
const int N = 1e2 + 5;
int w[N][N],d[N][N];
//int floyd(int n) //direct
//{
// memcpy(d,w,sizeof(w));
// for(int k=1;k<=n;k++)
// for(int i=1;i<=n;i++)
// for(int j=1;j<=n;j++)
// d[i][j] = min(d[i][j],d[i][k]+d[k][j]);
// int ans = 2e9;
// for(int i=1;i<=n;i++)
// for(int j=1;j<=n;j++)
// ans = min(ans,d[i][j]+d[j][i]);
// return ans;
//}
int floyd(int n) //un-direct
{
memcpy(d,w,sizeof(w));
int ans = 2e9;
for(int k=1;k<=n;k++)
{
for(int i=1;i<k;i++)
for(int j=i+1;j<k;j++)
if(w[i][k] != w[0][0] && w[k][j] != w[0][0])
ans = min(ans,w[i][k]+w[k][j]+d[j][i]);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
d[i][j] = min(d[i][j],d[i][k]+d[k][j]);
}
return ans==2e9 ? -1 : ans;
}
int main()
{
int n,m;
while(~scanf("%d%d",&n,&m))
{
memset(w,63,sizeof(w));
for(int i=1;i<=n;i++)
w[i][i] = 0;
for(int i=0;i<m;i++)
{
int u,v,ww;
scanf("%d%d%d",&u,&v,&ww);
if(w[u][v] > ww)
w[u][v] = w[v][u] = ww;
}
int ans = floyd(n);
if(ans == -1)
puts("It's impossible.");
else
printf("%d\n",ans);
}
return 0;
}