bzoj 3437 斜率优化DP

/**************************************************************
Problem: 3437
User: BLADEVIL
Language: C++
Result: Accepted
Time:3404 ms
Memory:39872 kb
****************************************************************/

//By BLADEVIL
#include <cstdio>
#define maxn 1000010
#define LL long long

using namespace std;

LL n;
LL a[maxn],sum[maxn],sum_[maxn];
LL que[maxn];
LL w[maxn];

LL get(LL i) {
return (w[i]+sum_[i]);
}

int main() {
scanf("%lld",&n);
for (LL i=1;i<=n;i++) scanf("%lld",&a[i]);
for (LL i=1;i<=n;i++) scanf("%lld",&sum[i]);
for (LL i=1;i<=n;i++) sum_[i]=i*sum[i];
for (LL i=1;i<=n;i++) sum[i]+=sum[i-1],sum_[i]+=sum_[i-1];
LL h(1),t(1);
for (LL i=1;i<=n;i++) {
while ((t-h>0)&&((sum[que[h]]-sum[que[h+1]])*i<=(w[que[h]]+sum_[que[h]]-w[que[h+1]]-sum_[que[h+1]]))) h++;
w[i]=w[que[h]]+(sum[i]-sum[que[h]])*i-(sum_[i]-sum_[que[h]])+a[i];
//w[i]=w[que[h]]+sum_[que[h]]-i*sum[que[h]]+i*sum[i]+a[i]-sum_[i];
while ((t>h)&&((get(que[t-1])-get(i))*(sum[que[t-1]]-sum[que[t]])<=(get(que[t-1])-get(que[t]))*(sum[que[t-1]]-sum[i]))) t--;
/*
while ((t-h>0)&&(
(w[que[t-1]]+sum_[que[t-1]]-w[i]-sum_[i])*(sum[que[t-1]]-sum[que[t]])<=
(w[que[t-1]]+sum_[que[t-1]]-w[que[t]]-sum_[que[t]])*(sum[que[t-1]]-sum[i]))
) t--;
*/
que[++t]=i;
}
printf("%lld\n",w
);
return 0;
}