Floyd算法核心片段如下:
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]);
}
}
}
最初Floyd算法是用来求全源最短路,复杂度为n3,但是如果通过选点来执行n此Dijkstra算法,可以在接近n2longn的时间复杂度下完成全源最短路的求解,但是Floyd算法还具有一些应用的空间,比如求解传递闭包问题。
该题是为判断,给出的一系列二元关系中,判断是否能得到一个有序列,或者这些二元关系是矛盾的。在Floyd算法中,其实状态转移的方程就是一种关系的传递,即如果从i到k,再从k到j,那么必然能够从i到j,所以二者加起来就是一条从i到j的路径,值为其长度。所以同理也能够处理传递闭包的问题,只需要将关系式改为:
d[i][j]|=d[i][k]&d[k][j];
这样,如果d[i][j]的值为1,代表i和j存在这种关系,否则不存在。对于该题本身,需要运用二分法先判断在出现矛盾之前是否存在可推出的序列,如果没有则通过二分寻找是否存在矛盾的关系式,如果仍没有则代表关系式不足不能推出一个序列。在判断时,如果存在i和j使d[i][j]和d[j][i]同时为1,则说明存在矛盾,因为该关系单向的,如果使同时为0,则说明关系不清,因为不能确定二者顺序。还需要注意数据的初始化等等。
题目要求我们在无向图中寻找一个最小环,我们在Floyd算法的基础上,增加了一个新的式子:
d[i][j]+a[j][k]+a[k][i]
其中a数组为点和点之间距离的记录,d数组则是i到j的最短路径,所以该式子代表的是一个环,且起码具有i,j,k三个点,通过枚举变量,最后得到的该式子的最小值就是要求的答案。比较值得注意的是求取最小环的部分,用到了一个分治。代码如下:
#include<cstdio>
#include<vector>
#include<cstring>
#define maxn 305
using namespace std;
int a[maxn][maxn],d[maxn][maxn],pos[maxn][maxn];
int n,m,ans=0x3f3f3f3f;
vector<int>path;
void get_path(int x,int y)
{
if(pos[x][y]==0)return ;
get_path(x,pos[x][y]);
path.push_back(pos[x][y]);
get_path(pos[x][y],y);
}
int main()
{
int n,m;
scanf("%d %d",&n,&m);
memset(a,0x3f,sizeof(a));
for(int i=1;i<=n;i++)a[i][i]=0;
for(int i=1;i<=m;i++)
{
int x,y,z;
scanf("%d %d %d",&x,&y,&z);
a[y][x]=a[x][y]=min(a[x][y],z);
}
memcpy(d,a,sizeof(a));
for(int k=1;k<=n;k++)
{
for(int i=1;i<k;i++)
{
for(int j=i+1;j<k;j++)
{
if((long long)d[i][j]+a[j][k]+a[k][i]<ans)
{
ans=d[i][j]+a[j][k]+a[k][i];
path.clear();
path.push_back(i);
get_path(i,j);
path.push_back(j);
path.push_back(k);
}
}
}
for(int i=1;i<=n;i++)
{
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.");
}
else
{
for(int i=0;i<path.size();i++)
{
printf("%d ",path[i]);
}
}
}