洛谷 P2501
第一问非常经典,构造 bi=ai−i,那么答案就是 n− b 的最长不下降子序列的长度。
考虑第二问,设 fi 表示以 i 结尾的最长不下降子序列的长度,prei 是能转移到 i 的集合,即 ∀j∈prei,fi=fj+1 且 bj≤bi。那么设 gi 表示在改变的数最小的前提下,前 i 个数合法的最小代价,则 gi=minj∈preigj+w(j+1,i−1),其中 w(l,r) 表示把区间 [l,r] 变的合法的最小代价。
结论:对于区间 [l,r],使其单调不降,则存在分界点 k 满足 b′i=bl(i≤k),b′i=br(i>k)。此时代价最小。
证明可以看这位大佬的,我就懒得写了。
具体细节看代码。
Code:
#include <bits/stdc++.h>
using namespace std;
#define pb push_back
typedef long long ll;
const int N = 35005, inf = 0x3f3f3f3f;
int n;
int a[N];
int dp[N], len, g[N];
vector <int> vec[N];
ll f[N], pre[N], suf[N];
int main() {
scanf("%d", &n);
for (int i = 1; i <= n; ++i) scanf("%d", &a[i]), a[i] -= i;
a[n + 1] = inf;
for (int i = 1; i <= n + 1; ++i) {
int l = 0, r = len;
while (l < r) {
int mid = l + r + 1 >> 1;
if (dp[mid] <= a[i]) l = mid;
else r = mid - 1;
}
if (l == len) ++len;
dp[g[i] = l + 1] = a[i], vec[g[i]].pb(i);
}
vec[0].pb(0);
printf("%d\n", n - len + 1);
memset(f, 0x3f, sizeof f), f[0] = 0, a[0] = -inf;
for (int i = 1; i <= n + 1; ++i) {
for (auto j : vec[g[i] - 1]) {
if (j > i || a[j] > a[i]) continue;
pre[j] = 0;
for (int k = j + 1; k < i; ++k) pre[k] = pre[k - 1] + abs(a[k] - a[j]);
suf[i - 1] = 0;
for (int k = i - 2; k >= j; --k) suf[k] = suf[k + 1] + abs(a[k + 1] - a[i]);
for (int k = j; k < i; ++k) f[i] = min(f[i], f[j] + pre[k] + suf[k]);
}
}
printf("%lld", f[n + 1]);
return 0;
}
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列01:轻松3步本地部署deepseek,普通电脑可用
· 按钮权限的设计及实现
· 25岁的心里话