1353. 滑雪场设计
分析
所有输入的山峰高度都在[0, 100]
最优解里, 修改完后的山峰高度是否会存在<0或者>100的山峰高度
直觉上讲, 是不可能的
证明:
假设在最优解里面, 某个山峰的高度<0, 因为最优解保证所有山峰高度差在17以内
我们可以将这些山峰移动到以0为起始的位置, 在我们假定的最优解里是将后面的值变到<0的位置
将山峰按照棕色线变化, 那么代价会变得更小, 最终会将所有山峰变到0这个位置, 它也是一个合法的方案, 并且也是更好的方案, 与<0是最优解,矛盾
因此, 最优解里不会存在某个山峰最优解<0
如果最优解横跨0, 所以在第2种情况下, 也不会有最优解<0, 因为可以按照如上方法, 将<0的点移动到0, 总的代价会减小
因此, 在最优解中, 所有山峰高度>= 0; 同理可以证明, 最优解所有山峰高度<= 100
有了这个结论后, 我们就可以直接枚举[0, 100]中, 长度为17的区间
如果当前区间是最优解, 怎么修改能够得到最小值呢, 如果山峰在区间内的话, 不需要做任何修改; 如果当前点在区间左边的话, 将当前到修改到左边界; 如果山峰在右侧, 那么将山峰修改到右边界; 这样的话, 对于每个山峰代价都能够取到最小值(原因是每个h[i]
是独立的)
扩展
题目要求:只能改变整数单位的高度, 如果题目可以改变小数单位的高度, 那么刚刚的做法, 还能用吗?
举个反例:
如果两个山峰高度是0, 18, 按照刚刚的做法, 最优的区间[0, 17], 或者[1, 18], 这样的话花费就是1
但是如果可以改变小数位置的高度, 那么可以将区间改成[0.5, 17.5]
0.25 + 0.25 = 0.5, 比刚才的代价更小, 所以说如果是小数的话, 不能用刚才的做法来做
小数的话, 可以用三分来做
code
#include <iostream>
using namespace std;
const int N = 1010;
int h[N], n;
int main(){
scanf("%d", &n);
for (int i = 1; i <= n; i ++ ) scanf("%d", &h[i]);
int res = 1e9;
for (int i = 0; i + 17 <= 100; i ++ ) {
int cost = 0, l = i, r = i + 17;
for (int j = 1; j <= n; j ++ )
if (h[j] < l) cost += (h[j] - l) * (h[j] - l);
else if (h[j] > r) cost += (h[j] - r) * (h[j] - r);
res = min(res, cost);
}
cout << res << endl;
return 0;
}
进阶要求
如果山峰高度可以改动成小数高度
假设当前山峰高度是
h
i
h_i
hi,
当
h
i
h_i
hi在区间右边的时候, 即区间左端点
x
i
+
17
≤
h
i
x_i + 17 \leq h_i
xi+17≤hi, 最小代价是
(
h
i
−
x
i
−
17
)
2
(h_i - x_i - 17)^2
(hi−xi−17)2
当
x
i
≤
h
i
≤
x
i
+
17
x_i \leq h_i \leq x_i + 17
xi≤hi≤xi+17, 代价是0
当
x
i
>
h
i
x_i > h_i
xi>hi, 代价是
(
h
i
−
x
i
)
2
(h_i - x_i)^2
(hi−xi)2
h
i
h_i
hi是固定的, 以
x
i
x_i
xi为自变量, 画图, 可以发现是个凸函数, 所以可以用三分求极值
482. 合唱队形
分析
dp问题
code
#include <iostream>
using namespace std;
const int N = 110;
int a[N], f[N], g[N];
int n;
int main(){
scanf("%d", &n);
for (int i = 1; i <= n; i ++ ) scanf("%d", &a[i]);
int res = 0;
for (int i = 1; i <= n; i ++ ){
f[i] = 1;
for (int j = 1; j < i; j ++ )
if (a[j] < a[i])
f[i] = max(f[i], f[j] + 1);
res = max(res, f[i]);
}
for (int i = n; i >= 1; i -- ){
g[i] = 1;
for (int j = n; j > i; j -- )
if (a[j] < a[i])
g[i] = max(g[i], g[j] + 1);
res = max(res, g[i]);
}
for (int i = 1; i <= n; i ++ ) res = max(res, f[i] + g[i] - 1);
cout << n - res << endl;
return 0;
}