HDU 1599 find the mincost route 最小环

题目：

http://acm.hdu.edu.cn/showproblem.php?pid=1599

题意：

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算法最外层循环执行到k时已求出以1~k-1为中间点的最短路径，无向图中一个环最少有三个点，设环中标号最大的顶点为w，跟w直接相连的两个点为v,u，则此最小环长度为dis[v][u] +g[v][w] + g[w][u]，dis[v][u]表示以1~k-1为中间节点时的最短路径，符合最外层循环k=w时的情况，之所以只求到以k-1为中间节点的最短路，是因为v到u的最短路中不能包含k，那么接下来枚举小于标号小于w的v、u两点的组合，即可求出最大标号为w的最小环的长度，当最外层循环结束时，就可以求出整个图的最小环。

#include <bits/stdc++.h>
using namespace std;

const int N = 110, INF = 1<<28;
int n, m;
int g[N][N], dis[N][N];

int floyd()
{
    int res = INF;
    for(int k = 1; k <= n; k++)
    {
        for(int i = 1; i < k; i++)
            for(int j = i+1; j < k; j++)
                res = min(res, dis[i][j] + g[i][k] + g[k][j]);
        for(int i = 1; i <= n; i++)
            for(int j = 1; j <= n; j++)
                dis[i][j] = min(dis[i][j], dis[i][k] + dis[k][j]);
    }
    return res;
}
int main()
{
    while(~ scanf("%d%d", &n, &m))
    {
        for(int i = 1; i <= n; i++)
            for(int j = 1; j <= n; j++)
                if(i == j) g[i][j] = dis[i][j] = 0;
                else g[i][j] = dis[i][j] = INF;
        for(int i = 1; i <= m; i++)
        {
            int v, u, c;
            scanf("%d%d%d", &v, &u, &c);
            g[v][u] = g[u][v] = dis[v][u] = dis[u][v] = min(g[v][u], c);
        }
        int res = floyd();
        if(res == INF) printf("It's impossible.\n");
        else printf("%d\n", res);
    }
    return 0;
}

另外这题还可以枚举两点，删除这两点的边，然后求两点间的最短路，环的长度即为两点间的最短路加上删除边的值，取最小的环长度即可，不过效率不高

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
const int N = 110, INF = 0x3f3f3f3f;

struct edge
{
    int to, cost, next;
}g[N*N*2];

int dis[N], mpa[N][N];
bool vis[N];
int cnt, head[N];
inline void add_edge(int v, int u, int cost)
{
    g[cnt].to = u, g[cnt].cost = cost, g[cnt].next = head[v], head[v] = cnt++;
}
void spfa(int s, int t)
{
    queue<int> que;
    memset(dis, 0x3f, sizeof dis);
    memset(vis, 0, sizeof vis);
    que.push(s), dis[s] = 0, vis[s] = true;
    while(! que.empty())
    {
        int v = que.front(); que.pop();
        vis[v] = false;
        for(int i = head[v]; i != -1; i = g[i].next)
        {
            int u = g[i].to;
            if((v == s && u == t) || (v == t && u == s)) continue;
            if(dis[u] > dis[v] + g[i].cost)
            {
                dis[u] = dis[v] + g[i].cost;
                if(! vis[u]) que.push(u), vis[u] = true;
            }
        }
    }
}
int main()
{
    int n, m;
    while(~scanf("%d%d", &n, &m))
    {
        cnt = 0;
        memset(head, -1, sizeof head);
        memset(mpa, 0x3f, sizeof mpa);
        int v[N*10], u[N*10], c[N*10];
        for(int i = 1; i <= m; i++)
        {
            scanf("%d%d%d", &v[i], &u[i], &c[i]);
            mpa[v[i]][u[i]] = mpa[u[i]][v[i]] = min(mpa[v[i]][u[i]], c[i]);
        }
        for(int i = 1; i <= n; i++)
            for(int j = 1; j <= n; j++)
                if(mpa[i][j] != INF) add_edge(i, j, mpa[i][j]);
        int res = INF;
        for(int i = 1; i <= n; i++)
            for(int j = i+1; j <= n; j++)
                if(mpa[i][j] != INF)
                {
                    spfa(i, j);
                    res = min(res, dis[j] + mpa[i][j]);
                }
        if(res == INF) printf("It's impossible.\n");
        else printf("%d\n", res);
    }
    return 0;
}