【洛谷】P2803 学校选址II

题目地址:

题目描述: 在一条大路一旁有许多栋楼,每栋楼里有许多小学生(哈哈哈一波小学生来袭!)。但是这条路上没有小学!!!!所以唯恐世界不乱的牛A打算在路上(汽车什么的都不敢来这个小学生云集的地方咯,所以不用担心安全问题)任选几点(可以和楼重合,当然也可以不重合)建立小学,且使所有小学生上学走的路程之和最短。牛A发现修建一所小学根本无法满足他唯恐世界不乱的(变态)心理,所以他准备建立 K K K所小学。

输入格式: 第一行 2 2 2个整数,表示楼数 n n n,学校数 k k k( 1 ≤ n , k ≤ 100 1le n,kle 100 1≤n,k≤100) 第二行 n n n个整数,表示每栋楼的学生数( 0 < 0< 0<每栋楼学生数 ≤ 100 ≤100 ≤100) 第三行 n − 1 n-1 n−1个数,分别表示楼 i i i到楼 i + 1 i+1 i+1之间距离( 1 ≤ 1≤ 1≤距离 ≤ 100 ≤100 ≤100, 1 ≤ i ≤ n − 1 1≤i≤n-1 1≤i≤n−1)

输出格式: 即学生走的距离和的最小值

以第 1 1 1个楼的位置作为坐标原点,设第 i i i个楼的位置为 d [ i ] d[i] d[i]。容易证明,如果要在 [ d [ i − 1 ] , d [ i ] ] [d[i-1],d[i]] [d[i−1],d[i]]中选一个地方放学校,则该学校最优位置一定是两个端点其一处。当前 i − 1 i-1 i−1个大楼的总学生数大于第 i i i到 n n n个大楼的总学生数的时候,显然要选 d [ i ] d[i] d[i]处建学校;如果小于,则在 d [ i − 1 ] d[i-1] d[i−1]处建学校;如果等于,则两个点都可以。

我们可以预处理一个二维数组 g [ i ] [ j ] g[i][j] g[i][j],表示只考虑第 i ∼ j isim j i∼j的楼的学生的时候,最优的放置学校的位置。 g g g是很容易求出的。再定义 f [ i ] [ k ] f[i][k] f[i][k]为在前 i i i个楼的范围内建 k k k个学校的情况下,最小的总距离和是多少。那么,如果 k ≥ i kge i k≥i,则 f [ i ] [ k ] = 0 f[i][k]=0 f[i][k]=0;如果 k = 1 k=1 k=1,即只建一个学校,则根据定义, f [ i ] [ k ] = g [ 1 ] [ i ] f[i][k]=g[1][i] f[i][k]=g[1][i];剩余情况,可以按最右边的那个学校服务哪些楼的学生来分类,如果是建在 d [ l ] ∼ d [ i ] d[l]sim d[i] d[l]∼d[i]之间(并且也服务这些学生),那么最小总距离为 f [ l − 1 ] [ k − 1 ] + g [ l ] [ i ] f[l-1][k-1]+g[l][i] f[l−1][k−1]+g[l][i],所以: f [ i ] [ k ] = min ⁡ l = 1 , 2 , . . . , i f [ l − 1 ] [ k − 1 ] + g [ l ] [ i ] f[i][k]=min_{l=1,2,...,i} f[l-1][k-1]+g[l][i] f[i][k]=l=1,2,...,iminf[l−1][k−1]+g[l][i]最后返回 f [ n ] [ k ] f[n][k] f[n][k]即可。代码如下:

#include <iostream>
#include <cstring>
using namespace std;

const int N = 110;
int n, K;
int dis[N];
int g[N][N], f[N][N], w[N];

int main() {
          
   
  scanf("%d%d", &n, &K);
  for (int i = 1; i <= n; i++) {
          
   
    int x;
    scanf("%d", &x);
    w[i] = w[i - 1] + x;
  }
  for (int i = 2; i <= n; i++) {
          
   
    int x;
    scanf("%d", &x);
    dis[i] = dis[i - 1] + x;
  }

  for (int i = 1; i <= n; i++)
    for (int j = 1; j <= i; j++) {
          
   
      // 找到只考虑第j到第i个楼的情况下,学校放哪里总距离和最小
      int k;
      for (k = j; k <= i; k++)
        if (w[k] - w[j - 1] >= w[i] - w[k]) break;
      k = min(k, i);

      for (int l = j; l <= i; l++)
        g[j][i] += abs(dis[k] - dis[l]) * (w[l] - w[l - 1]);
    }

  memset(f, 0x3f, sizeof f);
  for (int i = 1; i <= n; i++)
    for (int j = 1; j <= K; j++) {
          
   
      if (j >= i) f[i][j] = 0;
      else if (j == 1) f[i][j] = g[1][i];
      else {
          
   
        for (int l = 1; l <= i; l++)
          f[i][j] = min(f[i][j], f[l - 1][j - 1] + g[l][i]);
      }
    }

  printf("%d
", f[n][K]);
}

时间复杂度 O ( n 2 ( n + k ) ) O(n^2(n+k)) O(n2(n+k)),空间 O ( n ( n + k ) ) O(n(n+k)) O(n(n+k))。

经验分享 程序员 微信小程序 职场和发展