第八届河南省程序设计大赛 引水工程 Prim算法

引水工程

时间限制：2000 ms  |  内存限制：65535 KB

难度：3

描述

南水北调工程是优化水资源配置、促进区域协调发展的基础性工程，是新中国成立以来投资额最大、涉及面最广的战略性工程，事关中华民族长远发展。“南水北调工程”，旨在缓解中国华北和西北地区水资源短缺的国家战略性工程。就是把中国长江流域丰盈的水资源抽调一部分送到华北和西北地区。我国南涝北旱，南水北调工程通过跨流域的水资源合理配置，促进南北方经济、社会与人口、资源、环境的协调发展。

整个工程分东线、中线、西线三条调水线。东线工程位于东部，因地势低需抽水北送至华北地区。中线工程从汉水与其最大支流丹江交汇处的丹江口水库引水，自流供水给黄淮海平原大部分地区，20多座大中城市；西线工程在青藏高原上，由长江上游向黄河上游补水。

现在有N个区域需要建设水资源工程，它们可以自建水库解决缺水问题，也可以从已有水源的地区建立管道引水过来。当然，这些建设都需要大量投资。

你能不能给出一个优化水资源配置方案，在保证每个区域都能用上水的前提下，使得整个引水工程费用最低。

输入

第一行： K 表示有多少组测试数据。
接下来对每组测试数据：
第1行: N 表示有N个区域（ 1<=N<=300 ）
第2 行： W1 W2 …. WN Wi表示第i个区域自建水库需要的费用
再有N行： Pi1 Pi2 …. Pin Pij表示建立第i个区域与第j个区域引水管道的费用

输出

对于每组测试数据，输出占一行，即建立整个引水工程的最小费用。

样例输入

155 4 4 3 60 2 2 2 22 0 3 3 32 3 0 4 52 3 4 0 12 3 5 1 0

样例输出

10

 

很简单的prim算法，但是注意一点，需要从建立水库成本最少的点出发，另外如果从一个点到达另一个点成本已经超过自建成本则将路线成本修改为自建成本；

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int main()
{
    int a[303][303]={0},dis[303]={0},book[303]={0},a1[303]={0},sum=0,bs=0,k,n;
    cin>>k;
    while(k--)
    {
        cin>>n;
        int y;
        sum=0;
        bs=0;
        memset(a,0,sizeof(a));
        memset(dis,0,sizeof(dis));
        memset(book,0,sizeof(book));
        memset(a1,0,sizeof(a1));
        int min1=99999999;
        for(int i=1;i<=n;i++)
        {
            cin>>a1[i];
            if(a1[i]<min1){min1=a1[i];y=i;}
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                cin>>a[i][j];
                if(a[i][j]>a1[j])a[i][j]=a1[j];
            }
        }
        for(int i=1;i<=n;i++){
            dis[i]=a[y][i];
        }
        book[y]=1;
        bs++;
        while(bs<n)
        {
            int min2=99999999,j;
            for(int i=1;i<=n;i++){
                if(book[i]==0&&dis[i]<min2){
                    min2=dis[i];j=i;
                }
            }
            book[j]=1;sum+=min2;bs++;
            for(int i=1;i<=n;i++){
                if(book[i]==0&&dis[i]>a[j][i])
                    if(a[j][i]<a1[i])dis[i]=a[j][i];
                    else dis[i]=a1[i];
            }
        }
        int sc=sum+min1;
        printf("%d\n",sc);
    }
    return 0;
}