BZOJ P1096: [ZJOI2007]仓库建设

第二发斜率优化

t达成了h刷了两三次RE，掉了1%的AC率，好气嗷

下面是代码

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<cstdio>
using namespace std;
long long f[1000003]; 
struct hehe{
	int x,p,c;
}a[1000003];
long long G[1000003],sum[1000003];
int n;
int q[1000003];
bool cmp(hehe a,hehe b){
	return a.x<b.x;
}
double xie(long long j,long long k){
	return(f[k]-f[j]+sum[j]-sum[k]+G[k]*a[k].x*1.0-G[j]*a[j].x)/(G[k]-G[j]);
}
int main(){
	scanf("%lld",&n);
	for(int i=1;i<=n;i++){
		scanf("%lld%lld%lld",&a[i].x,&a[i].p,&a[i].c);
	}
	sort(a+1,a+n+1,cmp);
	for(int i=1;i<=n;i++){
		sum[i]=sum[i-1]+G[i-1]*(a[i].x-a[i-1].x);
		G[i]=G[i-1]+a[i].p;
	}
	f[0]=0;
	int h=1,t=1;
	q[1]=0;
	for(int i=1;i<=n;i++){
		while(h<t&&xie(q[h],q[h+1])<a[i].x){
			h++;
		}
		int j=q[h];
		f[i]=f[j]+sum[i]-sum[j]-G[j]*(a[i].x-a[j].x)+a[i].c;
		while(h<t&&xie(q[t],q[t-1])>xie(i,q[t])){
			t--;
		}
		q[++t]=i;
	}
	printf("%lld",f[n]);
	return 0;
}
/*
in:
3
0 5 10
5 3 100
9 6 10
out:
32
*/