[SDOI2016]征途 —— 斜率优化DP
时隔多年没有碰斜率优化了。。。
想当年被斜率优化虐的死去活来,现在看看。。。也就那样吧。
Pine开始了从S地到T地的征途。
从S地到T地的路可以划分成n段,相邻两段路的分界点设有休息站。
Pine计划用m天到达T地。除第m天外,每一天晚上Pine都必须在休息站过夜。所以,一段路必须在同一天中走完。
Pine希望每一天走的路长度尽可能相近,所以他希望每一天走的路的长度的方差尽可能小。
帮助Pine求出最小方差是多少。
设方差是v,可以证明,v×m2是一个整数。为了避免精度误差,输出结果时输出v×m2。
Input
第一行两个数 n、m。
第二行 n 个数,表示 n 段路的长度
Output
一个数,最小方差乘以 m^2 后的值Sample Input
5 21 2 5 8 6
Sample Output
36
HINT
1≤n≤3000,保证从 S 到 T 的总路程不超过 30000
首先根据方差的定义我们可以将方差的公式转换为
(数组为前缀和)
(具体转换过程可以百度百科方差)
对于式子中的是已经确定的常数,目前不用考虑。于是我们要求的便是。
于是题目就被我们转换为了将段路分为个块,使它们的平方和最小,不难得出DP的状态转移方程是:
只需要枚举k即可,其中状态的第一维定义为第段路,第二维定义为第个块。
然而不难发现时间复杂度过高,于是可以使用斜率优化。详见代码
#include <bits/stdc++.h>
using namespace std;
#define N 3100
#define LL long long
int n,A[N],S[N],m,q[3*N];
LL f[N][N];
int fucky (int a,int j) {
return f[a][j-1]+S[a]*S[a];
}
int fuck (int a,int b,int j) {
return (fucky(a,j)-fucky(b,j))/(S[a]-S[b]);
}
int main() {
cin>>n>>m;
for(int i=1;i<=n;i++) cin>>A[i],S[i]=S[i-1]+A[i];
for(int i=1;i<=n;i++) f[i][1]=S[i]*S[i];
for(int j=2;j<=m;j++) {
int l=1,r=1;
for(int i=1;i<=n;i++) {
while(l<r && fuck(q[l],q[l+1],j)<2*S[i]) l++;
int BestPos=q[l];
f[i][j]=f[BestPos][j-1]+(S[i]-S[BestPos])*(S[i]-S[BestPos]);
while(l<r && fucky(q[r],i,j)<fucky(q[r-1],q[r],j)) r--;
q[++r]=i;
}
}
cout<<f[n][m]*m-S[n]*S[n];
}
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 一个奇形怪状的面试题:Bean中的CHM要不要加volatile?
· [.NET]调用本地 Deepseek 模型
· 一个费力不讨好的项目,让我损失了近一半的绩效!
· .NET Core 托管堆内存泄露/CPU异常的常见思路
· PostgreSQL 和 SQL Server 在统计信息维护中的关键差异
· CSnakes vs Python.NET:高效嵌入与灵活互通的跨语言方案对比
· DeepSeek “源神”启动!「GitHub 热点速览」
· 我与微信审核的“相爱相杀”看个人小程序副业
· Plotly.NET 一个为 .NET 打造的强大开源交互式图表库
· 上周热点回顾(2.17-2.23)