第二发斜率优化
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
*/