Openjudge — 7624 山区建小学

时间:2021-06-28 18:41:46

问题描述###

*在某山区修建了一条道路,恰好穿越总共m个村庄的每个村庄一次,没有回路或交叉,任意两个村庄只能通过这条路来往。已知任意两个相邻的村庄之间的距离为di(为正整数),其中,0 < i < m。为了提高山区的文化素质,*又决定从m个村中选择n个村建小学(设 0 < n < = m < 500 )。请根据给定的m、n以及所有相邻村庄的距离,选择在哪些村庄建小学,才使得所有村到最近小学的距离总和最小,计算最小值。

输入###

第1行为m和n,其间用空格间隔

第2行为(m-1) 个整数,依次表示从一端到另一端的相邻村庄的距离,整数之间以空格间隔。

例如

10 3

2 4 6 5 2 4 3 1 3

表示在10个村庄建3所学校。第1个村庄与第2个村庄距离为2,第2个村庄与第3个村庄距离为4,第3个村庄与第4个村庄距离为6,…,第9个村庄到第10个村庄的距离为3。

输出###

各村庄到最近学校的距离之和的最小值。

样例输入###

10 2

3 1 3 1 1 1 1 1 3

样例输出###

18

来源

元培-From Whf

解题思路:###

动态规划来解决,这里需要一个辅助的数组dist,dist[i][j]表示在从i到j这一段区间建一所小学,i到j的村庄都到这个学校来上学的路程和。

状态表达:f[i][j]表示前i个村庄建j所学校,到里那个村庄最近的学校上学的路程和。

状态转移:f[i][j] = std::min(f[i][j],f[k][j-1]+dist[k+1][i]),也就是:

for(int i = 1;i<=n;i++)
for(int j = 1;j<i;j++)
for(int k = 1;k<i;k++)
f[i][j] = std::min(f[i][j],f[k][j-1]+dist[k+1][i]);

状态数量:n^2

转移代价:O(n)

时间复杂度:O(n^3)

空间复杂度:O(n^2)

还有一个关键点,就是有关如何求dist数组的,其实,我们可以发现,在i到j村庄里建小学,选i到j村庄的路程的中间位置村庄,一定是最优的。画个图你就能发现了。

如何快速求dist呢,其实,dist[i][j]等于dist[i][j-1]+num[j]-num[i+j]/2的,建议画个图。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring> int num[1001],dist[1001][1001],f[1001][1001];
int main()
{
int n,m;
std::cin>>n>>m;
for(int i = 2;i<=n;i++){std::cin>>num[i];num[i] += num[i-1];}
for(int i = 1;i<=n;i++)
for(int j = i;j<=n;j++)
dist[i][j] = dist[i][j-1]+num[j]-num[(i+j)/2];
memset(f,0x3f,sizeof(f));
for(int i = 1;i<=n;i++)f[i][1] = dist[1][i];
for(int i = 1;i<=n;i++)
for(int j = 1;j<i;j++)
for(int k = 1;k<i;k++)
f[i][j] = std::min(f[i][j],f[k][j-1]+dist[k+1][i]);
std::cout<<f[n][m];
return 0;
}