标签:初中 color 规划 最小花费 表示 strong 为什么 tco ble
前言
昨天XY讲课!讲到这题!还是IOI的题!不过据说00年的时候DP还不流行。
题面
http://poj.org/problem?id=1160
分析
§ 1 中位数
首先我们考虑,若有x1 < x2 < ... < xn,则当∑abs(x - xi)最小时,x为x1 , x2 , ... , xn这n个数的中位数。
证明如下:我们把x1 , x2 , ... , xn看作数轴上n个点,我们先考虑两端,要使abs(x - x1) + abs(x - xn)最小,那么x必定在x1和xn中间,这是个很显然的初中数学题。
好,那么我们就可以置这两个点不管了;剩下的点我们也这样考虑,如此类推,最后要么剩下一个点,要么就是剩下两个点中间的一段区间,答案就显而易见了。证毕。
(在July 1, 2018的AtCoder Contest 100 Task C Linear Approximation就是用这个思路,但是不知道为什么到比赛结束我还是WA,到时候也搞个这题的题解)
§ 2 给定一段村庄怎么放邮局
我们要求每个村庄到邮局距离总和最近,那么给定一段区间里的村庄,我们要找到一个满足要求的村庄建邮局,由上面关于中位数的结论可知,建邮局的村庄就是这段村庄最中间的那个。
§ 3 动态规划
既然我们知道,在一段给定的村庄下放邮局的最小花费,那么很容易想到一个区间DP的思路。
设DP[i][j]表示在编号为1~i的村庄中,建了j个邮局,所需的最小花费。
这时转移方程已经非常明显了:
DP[i][j] = min{ DP[k][j - 1] + cost[k + 1][i] }
其中j - 1 <= k <= i,cost[l][r]表示在l到r的区间里建一个邮局所需花费。
而cost是可以在O(n2)的时间内预处理出来的,那么时限是绝对绰绰有余了。
参考代码
1 // POJ1160 2 // IOI2000 3 // Post Office 4 #include <cstdio> 5 #include <algorithm> 6 const int MAXV = 305, MAXP = 35, INF = 0x3f3f3f3f; 7 8 int V, P, X[MAXV], DP[MAXV][MAXP], sum[MAXV][MAXV]; 9 10 int main() { 11 scanf("%d%d", &V, &P); 12 int i, j, k, lim; 13 for (i = 1; i <= V; i++) scanf("%d", &X[i]); 14 for (i = 1; i < V; i++) 15 for (j = i + 1; j <= V; j++) sum[i][j] = sum[i][j - 1] + X[j] - X[i + j >> 1]; // 预处理cost,这里用sum写了 16 // DP[0][0] = INF; // 这里注意DP[0][0]=0,但是DP[i][0]=INF,当i > 0时,因为这时候这个状态是没有意义的 17 for (i = 1; i <= V; i++) { 18 DP[i][0] = INF; 19 lim = std::min(i, P); 20 for (j = 1; j <= lim; j++) { 21 DP[i][j] = INF; 22 for (k = j - 1; k < i; k++) DP[i][j] = std::min(DP[i][j], DP[k][j - 1] + sum[k + 1][i]); 23 } 24 } 25 printf("%d\n", DP[V][P]); 26 return 0; 27 }
总结
如果放在现在,这大概就够个联赛题难度了。但是当年呢DP在算法竞赛中还不常用,所以当年要想到这个思路,是非常不容易了。
附
讲题是zzd大佬还讲了一个很奇怪的方程:
DP[i][j]表示前i个村庄,建了j个邮局,其中第i个村庄一定是邮局的最小花费。
当然这个状态转移也不难写,但是要注意答案不是DP[V][P],而是DP[V + 1][P + 1]。(即加一个假村庄,假村庄上建一个假邮局,这样就会把前面n个村庄的花费也正常纳入了,当然要通过一些手段不计算假村庄的花费)
标签:初中 color 规划 最小花费 表示 strong 为什么 tco ble
原文地址:https://www.cnblogs.com/CaptainSlow/p/9262312.html