bzoj1010

Posted mikuknight

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj1010相关的知识,希望对你有一定的参考价值。

1.HDU3507

数学分析

2.bzoj1010

原理解释

题解

根据动规方程状态i从状态j转化而来,
y只与j有关,k与i,j有关,b只与i或常数有关,可得直线方程y=kx+b,本题让b尽量小,故维护下凸包,使k单调递增

则已确认了当前直线的斜率,取得的j即为当前直线向左碰到的第一个点

无论如何,一定不会碰到斜率更小的,更新head;插入当前点,同时维护下凸包性质,更新tail

#include<cstdio>
#include<algorithm>
using namespace std;
#define maxn 50002
long long s[maxn],q[maxn],f[maxn];
double Y(long long j){return f[j]+s[j]*s[j];}
double X(long long j){return s[j];}
double slop(long long j,long long k){return(Y(k)-Y(j))/(X(k)-X(j));}
int main(){
    long long n,L,i,x,head,tail;
    scanf("%lld%lld",&n,&L);L++;s[0]=0;
    for(i=1;i<=n;i++){scanf("%lld",&s[i]);s[i]+=s[i-1];}
    for(i=1;i<=n;i++)s[i]+=i;
    head=tail=1;q[1]=0;
    for(i=1;i<=n;i++){
        while(head<tail && slop(q[head],q[head+1])<=2*(s[i]-L))head++;
        int j=q[head];
        f[i]=f[j]+(s[i]-s[j]-L)*(s[i]-s[j]-L);
        while(head<tail && slop(q[tail-1],q[tail])>slop(q[tail],i))tail--;
        q[++tail]=i;
    }
    printf("%lld
",f[n]);
    return 0;
}

 


以上是关于bzoj1010的主要内容,如果未能解决你的问题,请参考以下文章

bzoj1010

bzoj1010HNOI2008玩具装箱 toy

1010bzoj

BZOJ - 1010斜率优化DP

BZOJ1010单调性DP优化

BZOJ 1010 玩具装箱