蓝桥杯 算法训练 安慰奶牛

资源限制
时间限制：1.0s 内存限制：256.0MB
问题描述
Farmer John变得非常懒，他不想再继续维护供奶牛之间供通行的道路。道路被用来连接N个牧场，牧场被连续地编号为1到N。每一个牧场都是一个奶牛的家。FJ计划除去P条道路中尽可能多的道路，但是还要保持牧场之间 的连通性。你首先要决定那些道路是需要保留的N-1条道路。第j条双向道路连接了牧场Sj和Ej(1 <= Sj <= N; 1 <= Ej <= N; Sj != Ej)，而且走完它需要Lj的时间。没有两个牧场是被一条以上的道路所连接。奶牛们非常伤心，因为她们的交通系统被削减了。你需要到每一个奶牛的住处去安慰她们。每次你到达第i个牧场的时候(即使你已经到过)，你必须花去Ci的时间和奶牛交谈。你每个晚上都会在同一个牧场(这是供你选择的)过夜，直到奶牛们都从悲伤中缓过神来。在早上 起来和晚上回去睡觉的时候，你都需要和在你睡觉的牧场的奶牛交谈一次。这样你才能完成你的 交谈任务。假设Farmer John采纳了你的建议，请计算出使所有奶牛都被安慰的最少时间。

输入格式
第1行包含两个整数N和P。

接下来N行，每行包含一个整数Ci。

接下来P行，每行包含三个整数Sj, Ej和Lj。

输出格式
输出一个整数, 所需要的总时间(包含和在你所在的牧场的奶牛的两次谈话时间)。
样例输入
5 6
10
10
20
6
30
1 2 5
2 3 5
2 4 12
3 4 17
2 5 15
3 5 6
样例输出
178
数据规模与约定
5 <= N <= 10000，N-1 <= P <= 100000，0 <= Lj <= 1000，1 <= Ci <= 1,000。

题目的意思是你先选择出连接各牧场的连通图，使路达到最少。然后你需要选择一个牧场作为你睡觉的地方，这个地方的奶牛需要你走的时候安慰一次，晚上睡觉的时候安慰一次（如果你中间经过他，也是要安慰），他会比其他奶牛都多一次安慰，所以想要达到安慰时间最少，你选择的一定是权值最小的结点。对于每条路径而言，由于是连通图而且你还要回到原处，所以每条路径至少走两遍。所以我们就会得到一个公式：

其中边i是你从p中选择出来的，结点经过几次是根据有几条边和他相连决定的，最后还要加上一个全职最小的结点。
于是，现在的问题就在于如何寻找构成连通图的n-1条边了。我们会想到克鲁斯卡尔算法求最小生成树。但是对于这道题而言，以上的想法还不准确，他并不是简单的克鲁斯卡尔算法，因为他每个端点也是有权值的。于是我们可以转变一下思路，可以将每个路径的权值转化为权值的二倍加上左右端点的权值，这样就将端点考虑进去了。
代码注释如下：

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
int n,p; //结点个数 、边数
int node[10001]; // 结点的权值
int f[10001]; // 判断连通性

struct E{  //定义结构体来实现路径的输入
    int x,y,d; //x、y结点间距离
}e[100001];

bool cmp(E e1, E e2){    //最小生成树算法中，对边的权值进行排序
    return e1.d < e2.d;
}

int find(int x){       //查询两个端点是否已经连接
    return f[x] == x ? x : f[x] = find(f[x]);   //路径压缩，防止超时
}

int Kruscal(){    //最小生成树算法
    int ans = 0;
    sort(e,e+p,cmp);
    for(int i = 0; i < p; i++){    //按照边的权值进行遍历
        int xx = find(e[i].x);
        int yy = find(e[i].y);
        if(xx != yy){
            f[xx] = yy; // 将2点连通
            ans += e[i].d;
        }
    }
    return ans;
}
// 数据初始化
void init(){
    scanf("%d%d",&n,&p);    
    for(int i = 0; i < n; i++){
        scanf("%d",&node[i]);
    }
    for(int i = 0; i < p; i++){
        scanf("%d%d%d",&e[i].x, &e[i].y, &e[i].d);
        e[i].d = e[i].d * 2 + (node[e[i].x-1] + node[e[i].y-1]);    //最重要的一步操作  使用更新的权值
    }
    for(int i = 0; i < n; i++){
        f[i] = i;
    }
}

int main(){
    init();
    int min = 100000;
    for(int i = 0; i < n; i++){
        if(node[i] < min) min = node[i];
    }
    printf("%d\n", min+Kruscal());  //要加入最小权值的一个点

    return 0;
}