BZOJ 1911 [Apio2010]特别行动队

最新推荐文章于 2018-09-28 15:02:23 发布

原创最新推荐文章于 2018-09-28 15:02:23 发布 · 458 阅读

0 ·

CC 4.0 BY-SA版权

动态规划专栏收录该内容

107 篇文章

订阅专栏

本文介绍了一种使用斜率优化技巧来加速动态规划算法的方法。通过维护一个凸壳队列来保证状态转移方程的单调性，从而实现O(n)的时间复杂度。文章提供了完整的C++代码实现，并通过一个具体的例子展示了如何应用该方法。

【题目分析】
斜率优化

【代码】

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
long long dp[1000001],x[1000001],n;
long long que[1000001];
long long a,b,c,hd=1,tl=0;
inline double ang(long long j,long long k)
{
    return (((double)dp[k]-dp[j])+a*((double)x[k]*x[k]-x[j]*x[j]))/(x[k]-x[j]);
}
int main()
{
    scanf("%lld%lld%lld%lld",&n,&a,&b,&c);
    for (long long i=1;i<=n;++i)
    {
        scanf("%lld",&x[i]);
        x[i]+=x[i-1];
    }
    memset(dp,129,sizeof dp);
    dp[0]=0;
//    for (long long i=1;i<=n;++i)
//        for (long long j=0;j<i;++j)
//            dp[i]=max(dp[i],dp[j]+a*(x[i]-x[j])*(x[i]-x[j])+b*(x[i]-x[j])+c);
    que[++tl]=0;
    for (long long i=1;i<=n;++i)
    {
        while (hd<tl&&ang(que[hd],que[hd+1])>=(double)2*a*x[i]+b) hd++;
        dp[i]=max(dp[i],dp[que[hd]]+a*(x[i]-x[que[hd]])*(x[i]-x[que[hd]])+b*(x[i]-x[que[hd]])+c);
        while (hd<tl&&ang(que[tl-1],que[tl])<(ang(que[tl],i))) tl--;
        que[++tl]=i;
    }
    printf("%lld\n",dp[n]);
//    while (1);
}