[传智杯 #2 决赛] 传送门（松弛最短路+flord）

[传智杯 #2 决赛] 传送门

题目描述

传智专修学院里有 $n$ 栋教学楼，有 $m$ 条双向通行道路连接这些教学楼，不存在重边和自环。每条道路都有一定的长度，而且所有教学楼之间都可以直接或者间接的通过道路到达。我们可以很容易的求出这些教学楼之间的最短路。

为了使交通更为顺畅，校方决定在两个教学楼里增设一对传送门。传送门可以将这对教学楼的距离直接缩短为 0。利用传送门，某些教学楼之间的最短路的距离就变短了。

由于预算有限，学校里只能安装一对传送门。但是校长希望尽可能方便学生，使任意两点之间的最短路长度的总和最小。当然啦，从 $x$ 教学楼到 $y$ 教学楼的长度和从 $y$ 教学楼到 $x$ 教学楼的长度只需要统计一次就可以了。

输入格式

输入第 1 行两个正整数 $n,m(n\le 100,m\le \frac{1}{2}n(n-1))$，代表教学楼和道路数量。

接下来 $m$ 行，每行三个正整数 $x_i,y_i,w_i(0<w_i\le 10^4)$，表示在教学楼 $x_i$ 和 $y_i$ 之间，有一条长度为 $w_i$ 的道路。

输出格式

输出一行，在最优方案下的任意点对的最短道路之和。

样例 #1

样例输入 #1

4 5
1 2 3
1 3 6
2 3 4
2 4 7
3 4 2

样例输出 #1

14

提示

样例如图。当在 1 和 4 号教学楼架设一对传送门时，1 → 2 的最短路是 3，1 → 3 的最短路是 0+2，1 → 4 的最短路是 0，2 → 3 的最短路是 4，2 → 4 的最短路是 3+0，3 → 4 的最短路是 2，最短路之和是 14，是最佳方案。

思路

这道题如果没有传送门的话，它就是很简单的最短路。

• 但它加上了传送门，即从一点到另一点建了一条边权为 0 的边，但题目更恶心，叫我们自己建一个传送门。
• 还好，题目只要一个传送门，而且数据规模比较小，此时我们直接用flord先预处理出每点直接的最短路径，最后就枚举在每两点间建立传送门，看看最短路的变化如何。（但这边得注意的是，有可能这个方案不是最优的，因此我们得开一个临时数组，不要让它的数据丢失）

代码

//flord
//去你妈xrf，hjp
#include<iostream>
#include<algorithm>
#include<cstring>

using namespace std;

const int N = 5010;

int dist[N][N];
int g[N][N];
int n,m;

void flord(){
    for(int k=1;k<=n;k++){
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                dist[i][j]=min(dist[i][j],dist[i][k]+dist[k][j]);
            }
        }
    }
}

void init(){
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            g[i][j]=dist[i][j];
        }
    }
}

int main(){
    cin>>n>>m;
    
    memset(dist,0x3f,sizeof dist);
    
    for(int i=1;i<=m;i++){
        dist[i][i]=0;
        int a,b,c;
        cin>>a>>b>>c;
        dist[a][b]=dist[b][a]=c;
    }
    
    flord();
    
    int ans=2e9;
    
    //遍历以每点为核心
    
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            init();
            g[i][j]=g[j][i]=0;
            for(int k1=1;k1<=n;k1++){
                for(int k2=1;k2<=n;k2++){
                    g[k1][k2]=min(g[k1][k2],g[k1][i]+g[i][k2]);
                }
            }
            
            for(int k1=1;k1<=n;k1++){
                for(int k2=1;k2<=n;k2++){
                    g[k1][k2]=min(g[k1][k2],g[k1][j]+g[j][k2]);
                }
            }
            
            int res=0;
            
            for(int k1=1;k1<=n;k1++){
                for(int k2=1;k2<k1;k2++){
                    res+=g[k1][k2];
                }
            }
            ans=min(ans,res);
        }
    }
    
    cout<<ans;
    
    return 0;
}