【BZOJ3156】—防御准备（斜率优化）

最新推荐文章于 2019-11-05 22:53:27 发布

stargazer.

最新推荐文章于 2019-11-05 22:53:27 发布

阅读量151

点赞数

CC 4.0 BY-SA版权

分类专栏：斜率优化dp

本文链接：https://blog.youkuaiyun.com/qq_42555009/article/details/97100586

斜率优化dp 专栏收录该内容

6 篇文章

订阅专栏

本文通过斜率优化技术解决了一道动态规划问题。利用自定义结构体表示点，并实现点运算，通过比较斜率的方式进行优化。代码中详细展示了如何进行斜率优化，包括查找最优前缀状态和维护凸包等步骤。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

传送门

水题
随便推个式子斜率优化就完了

#include<bits/stdc++.h>
using namespace std;
const int N=1000005;
#define int long long
#define gc getchar
inline int read(){
	char ch=gc();
	int res=0,f=1;
	while(!isdigit(ch))f^=ch=='-',ch=gc();
	while(isdigit(ch))res=(res+(res<<2)<<1)+(ch^48),ch=gc();
	return f?res:-res;
}
#define ll long long
#define cs const
#define pb push_back
#define db double
struct pt{
	db x,y;
	pt(db _x=0,db _y=0):x(_x),y(_y){}
	friend inline pt operator +(cs pt &a,cs pt &b){
		return pt(a.x+b.x,a.y+b.y);
	}
	friend inline pt operator -(cs pt &a,cs pt &b){
		return pt(a.x-b.x,a.y-b.y);
	}
	friend inline db operator *(cs pt &a,cs pt &b){
		return a.x*b.y-a.y*b.x;
	}
	inline double slope(){
		return y/x;
	}
}p[N];
int stk[N],top;
ll f[N];
int a[N],n;
inline void trans(int i,int j){
	f[i]=f[j]+1ll*(i-j)*(i-j-1)/2+a[i];
}
inline bool check(int k,double s){
	return (p[stk[k]]-p[stk[k-1]]).slope()<2*s;
}
signed main(){
	freopen("lx.cpp","r",stdin);
	n=read();
	for(int i=1;i<=n;i++)a[i]=read();
	for(int i=1;i<=n;i++){
		int l=1,r=top,res=0;
		while(l<=r){
			int mid=(l+r)>>1;
			if(check(mid,i))l=mid+1,res=mid;
			else r=mid-1;
		}
		trans(i,stk[res]);
		p[i].x=i,p[i].y=2ll*f[i]+1ll*i*i+i;
		while(top&&(p[stk[top]]-p[stk[top-1]])*(p[i]-p[stk[top-1]])<=0)top--;
		stk[++top]=i;
	}
	cout<<f[n];
}