河南省第八届ACM省赛---引水工程

时间:2021-01-07 12:04:21

引水工程

时间限制: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个区域引水管道的费用
输出
对于每组测试数据,输出占一行,即建立整个引水工程的最小费用。
样例输入
1 5 
5 4 4 3 6 
0 2 2 2 2
2 0 3 3 3
2 3 0 4 5
2 3 4 0 1
2 3 5 1 0
样例输出
10
来源

第八届河南省程序设计大赛

prim算法+建立虚点  从下标为0的点开始跑

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<stack>
#include<queue>
#include<vector>
#include<deque>
#include<map>
#include<set>
#include<algorithm>
#include<string>
#include<cstdlib>
#include<cmath>
#include<sstream>
#include<ctime>
using namespace std;
#define eps 1e-6
#define ee exp(1.0)
#define pi acos(-1.0)
const int N = 305;
const int INF = 0x3f3f3f3f;
const int JNF = -0x3f3f3f3f;
int g[N][N];
int ans[N];
int low[N];
int nearvex[N];
int n;
bool check()
{
int i;
for (i = 0; i <= n; i++)
{
if (nearvex[i] != -1)
{
return false;
}
}
return true;
}
int prim(int v0)
{
int i, j;
int sumweight = 0;
for (i = 0; i <= n; i++)
{
low[i] = g[v0][i];
nearvex[i] = v0;
}
nearvex[v0] = -1;
while (!check())
{
int minnum = INF;
int v = -1;
for (j = 0; j <= n; j++)
{
if (nearvex[j] != -1 && low[j]<minnum)
{
v = j;
minnum = low[j];
}
}
if (v != -1)
{
nearvex[v] = -1;
sumweight += low[v];
for (j = 0; j <= n; j++)
{
if (nearvex[j] != -1 && g[v][j]<low[j])
{
low[j] = g[v][j];
nearvex[j] = v;
}
}
}
}
return sumweight;
}
int main()
{
int k;
int i, j;
int pos;
int temp;
int result;
scanf("%d", &k);
while (k--)
{
scanf("%d", &n);
memset(ans, INF, sizeof(ans));
for (i = 1; i <= n; i++)
{
scanf("%d", &g[0][i]);
}
g[0][0] = 0;
for (i = 1; i <= n; i++)
{
g[i][0] = g[0][i];
}
temp = INF;
for (i = 1; i <= n; i++)
{
if (ans[i]<temp)
{
temp = ans[i];
pos = i;
}
}
for (i = 1; i <= n; i++)
{
for (j = 1; j <= n; j++)
{
scanf("%d", &g[i][j]);
}
}
result = prim(0);
printf("%d\n", result);
}
return 0;
}