HDU1599 find the mincost route(floyd求最小环)

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1599

find the mincost route

Time Limit: 1000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 2435    Accepted Submission(s): 973

Problem Description

杭州有N个景区，景区之间有一些双向的路来连接，现在8600想找一条旅游路线，这个路线从A点出发并且最后回到A点，假设经过的路线为V1,V2,....VK,V1,那么必须满足K>2,就是说至除了出发点以外至少要经过2个其他不同的景区，而且不能重复经过同一个景区。现在8600需要你帮他找一条这样的路线，并且花费越少越好。

 

Input

第一行是2个整数N和M（N <= 100, M <= 1000)，代表景区的个数和道路的条数。
接下来的M行里，每行包括3个整数a,b,c.代表a和b之间有一条通路，并且需要花费c元(c <= 100)。

Output

对于每个测试实例，如果能找到这样一条路线的话，输出花费的最小值。如果找不到的话，输出"It's impossible.".

 

Sample Input

  
  

   
   3 3
1 2 1
2 3 1
1 3 1
3 3
1 2 1
1 2 3
2 3 1
  
  

 

Sample Output

  
  

   
   3
It's impossible.
  
  

floyd求最小环模板。

floyd核心思想：枚举每一个点去松弛其他两个点之间的最短路

思路：用Floyed解决，枚举每一个点K，和连接它的两个点i,j（i,j均<k)，使i,j,k构成环。则ans=min{ans,dis[i][j]+g[k][i]+g[k][j]}，根据Floyed原理，若i,j<k则dis[i][j]与k无关。

结合网上的代码，我理解了好长时间，应该注意的有以下几点：

1，初始化dp和mp数组时全部赋为INF，避免自己和自己形成环；

2，INF不能使用0x7fffffff，利用floyd求最短路的时候可能会超出Int的范围而爆掉

3，模板新增代码是核心代码，应加强理解ans = mmin(ans,dp[i][j]+mp[i][k]+mp[k][j]); 

i,j<k，这样才能使dp[i][j]与k无关，即dp[i][j]不会经过K点；

  我的理解是：两个点之间不经过k点的最短路径+经过k点中转的路径和取最小值为最小环，三点成环

4，此题有重边的数据，应取最小的费用，另外也是双向边

5，此题的dp数组也可以用动态规划的思想理解

floyd与最小环参考：http://blog.sina.com.cn/s/blog_476a25110100mag6.html

#include <iostream>
#include <cstring>
#include <cstdio>
#include <queue>
#define mmin(a,b) a<b? a:b
using namespace std;
const int INF=0xffffff;

int mp[105][105];
int n,m;
int dp[105][105];

//floyd求最小环模板
int floyd()
{
    int i,j,k;
    int ans = INF;

    for(k=1;k<=n;k++)
    {
        //floyd代码新增求最小环部分
        for(i=1;i<k;i++)
            for(j=i+1;j<k;j++)
                //mp[i][k]+mp[k][j]表示经过k中转从而形成环，此环包括i到j的最短路
                //再在这些最小环中找到费用最小的一个赋给ans
                ans = mmin(ans,dp[i][j]+mp[i][k]+mp[k][j]);
                
        //floyd原始代码
        for(i=1;i<=n;i++)
            for(j=1;j<=n;j++)
                dp[i][j]=mmin(dp[i][j],dp[i][k]+dp[k][j]);
    }
    return ans;
}


int main()
{
    while(scanf("%d %d",&n,&m)!=EOF)
    {
        int x,y,z;
        int i,j;

        for(i=1;i<=n;i++)
            for(j=1;j<=n;j++)
                mp[i][j]=dp[i][j]=INF;

        for(i=0;i<m;i++)
        {
            scanf("%d %d %d",&x,&y,&z);
            if(z<mp[x][y])
            {
                mp[x][y]=z;mp[y][x]=z;
                dp[x][y]=z;dp[y][x]=z;
            }

        }
        int result = floyd();
        if(result==INF) printf("It's impossible.\n");
        else printf("%d\n",result);

    }
    return 0;
}