hdu1599——find the mincost route

最新推荐文章于 2024-09-14 20:30:05 发布

lazy_virus

最新推荐文章于 2024-09-14 20:30:05 发布

阅读量885

点赞数

文章标签： graph 算法 disk c

本文链接：https://blog.csdn.net/k1246195917/article/details/5841116

版权

Floyd 的改进写法可以解决最小环问题，时间复杂度依然是 O(n^3)，储存结构也是邻接矩阵

int mincircle = infinity;

Dist = Graph;

for ( int k=0;k<nVertex;++k){

     //新增部分:
     for ( int i=0;i<k;++i)
         for ( int j=0;j<i;++j)
            mincircle = min(mincircle,Dist[i][j]+Graph[j][k]+Graph[k][i]);

//通常的 floyd 部分:

for ( int i=0;i<nVertex;++i)
for ( int j=0;j<nVertex;++j){

int temp = Dist[i][k] + Disk[k][j];

if (temp < Dist[i][j])

Dist[i][j] = Dist[j][i] = temp;

}

上面是对无向图的情况。

Floyd 算法保证了最外层循环到 k 时所有顶点间已求得以 0…k-1 为中间点的最短路径。一个环至少有3个顶点，设某环编号最大的顶点为 L ，在环中直接与之相连的两个顶点编号分别为 M 和 N (M,N < L)，则最大编号为 L 的最小环长度即为 Graph(M,L) + Graph(N,L) + Dist(M,N) ，其中 Dist(M,N) 表示以 0…L-1 号顶点为中间点时的最短路径，刚好符合 Floyd 算法最外层循环到 k=L 时的情况，则此时对 M 和 N 循环所有编号小于 L 的顶点组合即可找到最大编号为 L 的最小环。再经过最外层 k 的循环，即可找到整个图的最小环。

引用：http://dantvt.spaces.live.com/blog/cns!D87988A6CAC0A480!916.entry

#include<stdio.h>
#include<string.h>
#define max 105
#define maxcost 1000000
int g[max][max],a[max][max];
int n,m;
void floyd()
{
    int i,j,k,ans=maxcost;
    for(i=0;i<=n;i++)
        for(j=0;j<=n;j++)
        {
            a[i][j]=g[i][j];
        }
        for(k=1;k<=n;k++)
        {
            for(i=1;i<=k;i++)
                for(j=1;j<i;j++)
                {
                    ans=ans>a[i][j]+g[j][k]+g[k][i]?a[i][j]+g[j][k]+g[k][i]:ans;
                }
            for(i=1;i<=n;i++)
                for(j=1;j<=n;j++)
                {
                    if(a[i][j]>a[i][k]+a[k][j])
                        a[i][j]=a[i][k]+a[k][j];
                }
        }
        if(ans==maxcost)
            printf("It's impossible./n");
        else
            printf("%d/n",ans);
}
main()
{
    int i,j,x,y,c;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        for(i=0;i<=n;i++)
        {
            for(j=0;j<=n;j++)
            {
                g[i][j]=maxcost;
            }
        }
        for(i=0;i<m;i++)
        {
            scanf("%d%d%d",&x,&y,&c);
            if(g[x][y]>c)
            {
                g[x][y]=c;
                g[y][x]=c;
            }
        }
        floyd();
    }
}