HDU 1599 find the mincost route(Floyd求最小环)

题目链接：

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

题意：

给定一个无向图，可能存在重边和自环，问所成环的最小权值，环中最少需3个点；

分析：

$floyd$求最小环，在一遍一遍更新最短路的时候，来找环，那么怎样会形成环；假设有$i<j<k$不相同三点，$i$到$j$之间有若干个点，$j$是$k$之前的一个点，那么这三点所在环，一定有，$i$和$k$直接连接，那么$e[i][k]$与$e[j][k]$一定不等于0(e数组代表直接相连两点间距离，dis代表两点间最短路)，那么最小环的值一定是求最短路过程中更新到的最小的$dis[i][j]+e[j][k]+e[k][i]$；
比较好的floyd求最小环理解博客：https://www.xuebuyuan.com/756949.html
代码中需要注意的是dis数组的初始化的inf值，不要太大，不然会爆int，或者把数组都开成long long

#include<algorithm>
#include<iostream>
#include<cstring>
#include<string>
#include<cstdio>
#include<vector>
#include<queue>
#include<stack>
#include<cmath>
#include<set>
#include<map>

using namespace std;

#define  inf 0x7f7f7f
#define  maxn 100411
#define  N 2
#define mod 100003

typedef long long ll;
typedef struct {
    int u, v, next, w;
} Edge;
Edge e[N];
int cnt, head[N];

inline void add(int u, int v) {
    e[cnt].u = u;
    e[cnt].v = v;
    //e[cnt].w=w;
    // e[cnt].f=f;
    e[cnt].next = head[u];
    head[u] = cnt++;
    //e[cnt].u=v;
//    e[cnt].v=u;
//    e[cnt].w=0;
//    e[cnt].f=-f;
//    e[cnt].next=head[v];
//    head[v]=cnt++;
}

inline void write(int x) {
    if (x < 0)
        putchar('-'), x = -x;
    if (x > 9)
        write(x / 10);
    putchar(x % 10 + '0');
}

inline int read() {
    int x = 0, f = 1;
    char c = getchar();
    while (c < '0' || c > '9') {
        if (c == '-')
            f = -1;
        c = getchar();
    }
    while (c >= '0' && c <= '9') {
        x = x * 10 + c - '0';
        c = getchar();
    }
    return x * f;
}

int n, m;
int dis[110][110],ed[110][110];

int main() {
    while (scanf("%d%d", &n, &m) != EOF) {
        int a, b, c;
        int ans = inf;
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                if(i!=j)
                    ed[i][j]=dis[j][i]=inf;
                else ed[i][j]=dis[i][j]=0;
        for (int i = 0; i < m; i++) {
            scanf("%d%d%d", &a, &b, &c);
            dis[b][a] = dis[a][b] = ed[b][a] = ed[a][b] = min(c, ed[a][b]);
        }
        for (int k = 1; k <= n; k++) {
            for (int i = 1; i < k; i++)
                for (int j = i + 1; j < k; j++) {
                    //cout<<ans<<endl;
                    ans = min(ans, dis[i][j] + ed[i][k] + ed[j][k]);
                }
            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]);
        }
        if (ans == inf)cout << "It's impossible." << endl;
        else cout << ans << endl;
    }
    return 0;
}

我们坚持一件事情，并不是因为这样做了会有效果，而是坚信，这样做是对的。
——哈维尔