【动态规划】常见基础DP优化技巧
Little_Cake_qwq · · 算法·理论
前言
这是我的第一篇文章,有建议的可以私信我或评论 QwQ。
update:2026.1.31 感谢 @Jacken 指出错误。
正文
废话不多说直接开讲
1. 前缀和优化
前缀和优化 DP 常出现在计数 DP 中,此时的状态通过一段区间和来转移。
例题:AtCoder-dp_m Candies
这个题真的是典中典了。
设状态
时间复杂度是
不难发现状态转移是个区间和,所以设
时间复杂度为
代码
int main()
{
n=read(),m=read();
for(int i=0;i<=m;i++) s[i]=1; //滚掉一维
for(int i=1;i<=n;i++)
{
int a;a=read();
for(int j=0;j<=m;j++)
{
if(j<=a) dp[i][j]=(dp[i][j]+s[j])%mod;
else dp[i][j]=(dp[i][j]+s[j]-s[j-a-1])%mod;
}
s[0]=dp[i][0];
for(int j=1;j<=m;j++) s[j]=s[j-1]+dp[i][j];
}
cout<<dp[n][m];
return 0;
}
2. mx/mn 数组优化
这种优化方法一般出现在求最优解 DP 中。
例题:P1107 雷涛的小猫
也是一个经典题目了。
设状态
在位置
时间复杂度
不难发现柿子式子等价于:
其中的
时间复杂度
代码
int main(){
n=read(),h=read(),d=read();
for(int i=1;i<=n;i++)
{
int t,x;t=read();
while(t--) x=read(),a[x][i]++;
}
for(int i=h;i>=0;i--)
{
for(int j=1;j<=n;j++)
{
dp[i][j]=max(dp[i+1][j],mx[i+t])+a[i][j];
mx[i]=max(mx[i],dp[i][j]);
ans=max(ans,dp[i][j]);
}
}
cout<<ans;
return 0;
}
3. 单调队列(单调栈)优化
重头戏来了!
不会模板的可以去 P1886【模板】单调队列 / 滑动窗口 和 P5788 【模板】单调栈 学习。
这种一般出现在通过区间最大值转移的 DP 优化中。
例题:P1725 琪露诺
设
时间复杂度
不难发现
代码
int main()
{
n=read(),l=read(),r=read();
for(int i=0;i<=n;i++) a[i]=read(),dp[i]=-INF;
dp[0]=0;
for(int i=l;i<=n;i++)
{
while(!q.empty()&&dp[q.back()]<dp[i-l]) q.pop_back();
q.push_back(i-l);
while(!q.empty()&&q.front()<i-r) q.pop_front();
dp[i]=dp[q.front()]+a[i];
}
int ans=-INF;
for(int i=n-r+1;i<=n;i++) ans=max(ans,dp[i]);//注意答案
cout<<ans;
return 0;
}
4. 贪心(数学)优化
这类优化一般在求最优解 DP 中出现,我感觉这类优化挺神奇的。
例题:P13957 [ICPC 2023 Nanjing R] 背包
设
时间复杂度
通过贪心,不难(其实很难)发现选取
时间复杂度
代码
int main()
{
n=read(),m=read(),k=read();
for(int i=1;i<=n;i++) a[i].w=read(),a[i].v=read;
sort(a+1,a+n+1,cmp);
for(int i=1;i<=n;i++)
{
for(int j=m;j>=a[i].v;j--)
dp[j]=max(dp[j],dp[j-a[i].w]+a[i].v);
priority_queue<int> q;
for(int j=i+1;j<=n;j++) q.push(a[j].v);
int qwq=dp[m];
for(int j=1;j<=k;j++)
if(!q.empty()) qwq+=q.top(),q.pop();
ans=max(ans,qwq);
}
cout<<ans;
return 0;
}
最后
你可以拿一些题练手咯:
- P2513 [HAOI2009] 逆序对数列
- CF2049D Shift + Esc
- P3800 Power 收集
- P9119 [春季测试 2023] 圣诞树