题意
现在有一个长度为n的序列a,你最多可以进行k次单点修改,最小化max(|ai−ai−1|)
n,k≤2×103,−109≤ai≤109
样例
In:
63
123789
out:
1
Solution
看到最小化一个最大值这条信息,马上想到的应该就是二分。
所以我们二分答案这个最大值,问题就被我们转化成了:判断是否可以通过不超过 k 次操作使得序列任意相邻两个数差的绝对值小于mid
而这个东西我们是可以DP的。但如何设计状态?
我们发现如果仅仅这样设计:”设fi表示1到i中要满足条件最少要修改多少次”,这个东西是很难进行转移的,因为序列是会进行修改的,而这个DP状态无法考虑修改的情况。
我举个例子可能会更好理解一些:
正常地,按照这个例子,若ai进行了修改,状态转移是:fi=min(fi,fi−1+1),条件是:|ai−ai−1|>mid
但我们发现:fi−1这个值不一定是在ai−1没被修改的情况下得到的。而若其恰好会被修改,换句话说,ai−1并不是原来的值,那么|ai−ai−1|>mid这个条件就出现错误了
正确设计状态
所以我们针对ai−1并不是原来的值这个问题来思考如何改进状态。
因为在原先的状态中,我们并不知道ai−1是否被修改,所以我们不妨强制其不被修改,然后把DP值中” 最少要修改多少次 “改为” 最多有多少个数不被修改 “
那么状态就变为了:设fi表示1到i中要满足条件最多有多少个数不被修改,且强制ai不被修改。
而这里状态转移略微有些变化:由于我们不知道上一个不修改的地方是哪里,所以我们枚举j(1≤j<i)
那么[j+1,i−1]这段区间进行修改需要满足的条件是|ai−aj|≤mid×(i−j) (应该比较好理解)
直接n2枚举进行状态转移,与LIS问题的转移方程有点类似
而最后判断是否可行时,我们需要枚举序列中最后一个没被修改的位置x,看看fx是否≥n−k。而如果整个序列找不到这样的x说明不可行
那么这道题就做完了。
时间复杂度:O(n2)
Code
#include <bits/stdc++.h>
const int N = 2010;
#define LL long long
int n, K, dp[N];
LL a[N];
LL get_abs(LL a, LL b)
{
if(a < b) std::swap(a, b);
return a - b;
}
bool check(LL mid)
{
for(int i=1;i<=n;i++) dp[i] = 1;
for(int i=1;i<=n;i++)
for(int j=1;j<i;j++)
{
if(get_abs(a[i], a[j]) <= 1LL * (i - j) * mid) dp[i] = std::max(dp[i], dp[j] + 1);
}
for(int i=1;i<=n;i++) if(dp[i] >= (long long)n - K) return true;
return false;
}
int main()
{
scanf("%d%d", &n, &K);
for(int i=1;i<=n;i++) scanf("%lld", &a[i]);
LL l = 0, r = 2000000000;
while(l < r)
{
LL mid = (l + r) >> 1;
if(check(mid)) r = mid;
else l = mid + 1;
}
printf("%lld", l);
return 0;
}