bzoj3437(斜率优化

像这题，有的时候正向的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;
}