像这题,有的时候正向的dp式子不好进行斜率优化,
就是 f【i】=min{f【j】+cost【j+1~i】}样的式子,一般好多都可以斜率优化,但是有一些cost【j+1~i】不能在常数时间内差分之类求出来,
反正这题是正向不好差分,但是反向可以差分,就是锯木厂选址那样的思路,一个s【i】表示sigema j b【j】*dis[j] (距离和值相乘这一类的)
1:反向来做easy
2:求cost【j~i】的就是s【i】表示值与距离相乘的前缀和(也就是到第一个点的数值),然后再存距离和值的前缀和,O(1)求出cost【j~i】区间的值
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
using namespace std;
typedef long long ll;
const int N=1000005;
ll f[N],s[N],sum[N],a[N],b[N];
int n,q[N],head,tail;
ll getup(int j,int k)
{
return (f[j]-s[j]+sum[j]*j)-
(f[k]-s[k]+sum[k]*k);
}
ll getdown(int j,int k)
{
return j-k;
}
ll getdp(int i,int j)
{
return f[j]+s[i]-s[j]-(sum[i]-sum[j])*j+a[i+1];
}
int main()
{
scanf("%d",&n);
for (int i=n;i>=1;i--) scanf("%lld",&a[i]);
for (int i=n;i>=1;i--) scanf("%lld",&b[i]);
n++;
for (int i=1;i<=n;i++) sum[i]=sum[i-1]+b[i];
for (int i=1;i<=n;i++) s[i]=s[i-1]+b[i]*(i-1);
head=tail=1;
q[head]=0;
f[0]=a[1];
for (int i=1;i<=n-1;i++)
{
while (head<tail &&
getup(q[head+1],q[head])<=sum[i]*getdown(q[head+1],q[head])
) head++;
f[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("%lld",f[n-1]);
return 0;
}