hdu 3507 斜率优化

#include<iostream>
#include<string>
#include <bits/stdc++.h>
using namespace std;

int dp[500005];
int q[500005];
int sum[500005];
int head, tail, n, m;

int getDP(int i, int j)
{
        return dp[j] + m + (sum[i] - sum[j]) * (sum[i] - sum[j]);
}

int getUP(int j, int k) 
{
        return dp[j] + sum[j] * sum[j] - (dp[k] + sum[k] * sum[k]);
}

int getDOWN(int j, int k) 
{
        return 2 * (sum[j] - sum[k]);
}

int main()
{
        int i;
        while (scanf("%d%d", &n, &m) == 2)
        {
                sum[0] = dp[0] = 0;
                for (i = 1; i <= n; i++)
                {
                        scanf("%d", &sum[i]);
                        sum[i] += sum[i - 1];
                }
                head = tail = 1;
                q[tail] = 0;
                for (i = 1; i <= n; i++)
                {
                        while (head < tail && getUP(q[head + 1], q[head]) <= sum[i]*getDOWN(q[head + 1], q[head]))
                        {
                                head++;
                        }
                        dp[i] = getDP(i, q[head]);
                        while (head < tail && getUP(i, q[tail])*getDOWN(q[tail], q[tail - 1]) <= getUP(q[tail], q[tail - 1])*getDOWN(i, q[tail]))
                        {
                                tail--;
                        }
                        q[++tail] = i;
                }
                printf("%d\n", dp[n]);
        }
        return 0;
}

 

转载于:https://www.cnblogs.com/Aragaki/p/7687468.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值