1096: [ZJOI2007]仓库建设

一看题，DP？方程：f[i]=f[j]+c[i]+（将i到j+1的所有货物搬到i来的费用）

再看数据范围，贪心？

然而证不出来，然后还是DP，观察一下符合斜率优化的条件，就是有一个数据无法在O（1）的时间内求出来（如上）。

遂翻题解，然后。。。。。。。。。我智商拙计是我的错，只想到正着加没想到减法同样可行。

于是令g[n]=sigma(x[i]*p[i])   sum[n]=sigma(p[i])，(sum[i]-sum[j])*x[i]-(g[i]-g[j])即为所求。

然后就很简单了，单调队列维护之。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
const int N=1000000+5;
ll f[N],g[N],sum[N];
ll x[N],p[N],c[N];
int n,q[N];
double slop(int j,int k){
	return (double)(f[j]-f[k]+g[j]-g[k])/(sum[j]-sum[k]);
}
inline int read(){
	int num=0;char ch;
	while(ch<'0'||ch>'9')ch=getchar();
	while(ch>='0'&&ch<='9'){num=num*10+ch-'0';ch=getchar();}
	return num;
}
int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++){
		x[i]=read();p[i]=read();c[i]=read();
	}
	for(int i=1;i<=n;i++){
		sum[i]=sum[i-1]+p[i];
		g[i]=g[i-1]+x[i]*p[i];
	}
	int head=0,tail=0;
	for(int i=1;i<=n;i++){
		while(head<tail&&slop(q[head+1],q[head])<x[i])head++;
		int tmp=q[head];
		f[i]=f[tmp]+(sum[i]-sum[tmp])*x[i]-(g[i]-g[tmp])+c[i];
		while(head<tail&&slop(i,q[tail])<slop(q[tail],q[tail-1]))tail--;
		q[++tail]=i;
	}
	printf("%lld",f[n]);
	return 0;
}