【APIO2010】特别行动队(斜率优化dp)

题目:
我是超链接
题解:
斜率优化dp,维护上凸壳
这里写图片描述
那么橙色的点及其以上部分是要保留的,将下面的直线和点包括在内,意为上凸(但长得像下凸,我也不知道到底叫什么)
多亏了学姐的blog呐
代码:

#include <cstdio>
#define LL long long 
using namespace std;
LL s[1000050],f[1000005];
int a,b,c,n;
int q[1000005];
LL K(int j){return -2*a*s[j];}
LL B(int j){return a*s[j]*s[j]-b*s[j]+f[j];}
LL Y(int i,int j){return K(j)*s[i]+B(j);}
bool fg(int x1,int x2,int x3)//是否完全覆盖一条直线 
{
    LL w1=(LL)(K(x1)-K(x2))*(LL)(B(x3)-B(x1));
    LL w2=(LL)(K(x1)-K(x3))*(LL)(B(x2)-B(x1));
    return w1>=w2;
}
int main()
{
    int i,x,r;
    scanf("%d",&n);
    scanf("%d%d%d",&a,&b,&c);
    for (i=1;i<=n;i++)
      scanf("%lld",&s[i]),s[i]+=s[i-1];
    int l=r=0;
    for (i=1;i<=n;i++)
    {
        while (l<r && Y(i,q[l])<=Y(i,q[l+1])) l++;
        f[i]=Y(i,q[l])+a*s[i]*s[i]+b*s[i]+c;
        while (l<r && fg(i,q[r],q[r-1])) r--;
        q[++r]=i;
    }
    printf("%lld",f[n]);   
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值