动态规划(算法竞赛、蓝桥杯)--斜率优化DP打印文章

1、B站视频链接:E51【模板】斜率优化DP 打印文章_哔哩哔哩_bilibili

题目链接:Problem - 3507

cpp 复制代码
#include <bits/stdc++.h> 
using namespace std;
typedef long long LL;
const int N=500010;
int n,m,q[N];
LL s[N],f[N];

double sl(int i,int j){
	return (double)(f[i]+s[i]*s[i]-f[j]-s[j]*s[j])
				/(s[i]==s[j]?1e-9:s[i]-s[j]);
}
int main(){
	while(~scanf("%d%d",&n,&m)){
		for(int i=1;i<=n;i++){
			scanf("%lld",&s[i]);
			s[i]+=s[i-1];
		}
	}
	int h=1,t=0;
	for(int i=1;i<=n;i++){
		while(h<t&&sl(i-1,q[t])<=sl(q[t],q[t-1]))t--;
		q[++t]=i-1;
		while(h<t&&sl(q[h+1],q[h])<=2*s[i])h++;
		int j=q[h];
		f[i]=f[j]+(s[i]-s[j])*(s[i]-s[j])+m;
	}
	printf("%lld\n",f[n]);
	return 0;
}
相关推荐
fie88892 小时前
NSCT(非下采样轮廓波变换)的分解和重建程序
算法
晨晖22 小时前
单链表逆转,c语言
c语言·数据结构·算法
im_AMBER4 小时前
Leetcode 78 识别数组中的最大异常值 | 镜像对之间最小绝对距离
笔记·学习·算法·leetcode
鼾声鼾语4 小时前
matlab的ros2发布的消息,局域网内其他设备收不到情况吗?但是matlab可以订阅其他局域网的ros2发布的消息(问题总结)
开发语言·人工智能·深度学习·算法·matlab·isaaclab
LYFlied5 小时前
【每日算法】LeetCode 25. K 个一组翻转链表
算法·leetcode·链表
Swizard5 小时前
别再迷信“准确率”了!一文读懂 AI 图像分割的黄金标尺 —— Dice 系数
python·算法·训练
s09071365 小时前
紧凑型3D成像声纳实现路径
算法·3d·声呐·前视多波束
可爱的小小小狼5 小时前
算法:二叉树遍历
算法
d111111111d6 小时前
在STM32函数指针是什么,怎么使用还有典型应用场景。
笔记·stm32·单片机·嵌入式硬件·学习·算法
AI科技星6 小时前
质量定义方程常数k = 4π m_p的来源、推导与意义
服务器·数据结构·人工智能·科技·算法·机器学习·生活