LOJ 10189 仓库建设 ——斜率优化dp

题目：https://loj.ac/problem/10189

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define db double
#define ll long long
using namespace std;
const int N=1e6+;
int n,d[N],p[N],C[N];
ll s[N],c[N],dp[N];
int q[N],he,tl;
int rdn()
{
  int ret=;bool fx=;char ch=getchar();
  while(ch>''||ch<''){if(ch=='-')fx=;ch=getchar();}
  while(ch>=''&&ch<='') ret=(ret<<)+(ret<<)+ch-'',ch=getchar();
  return fx?ret:-ret;
}
ll Y(int a){return dp[a]+c[a]*d[a]-s[a];}
db slp(int u,int v){return (db)(Y(v)-Y(u))/(d[v]-d[u]);}
int main()
{
  n=rdn();
  for(int i=;i<=n;i++)
    d[i]=rdn(),p[i]=rdn(),C[i]=rdn();
  for(int i=;i<=n;i++)
    d[i]=d[n]-d[i];
  for(int i=n;i>=;i--)
    c[i]+=c[i+]+p[i],s[i]=s[i+]+(ll)p[i]*d[i];
  he=;q[tl=]=n;dp[n]=C[n];
  for(int i=n-;i>=;i--)
    {
      while(he<tl&&slp(q[he],q[he+])<=c[i+])he++;
      int j=q[he];
      dp[i]=dp[j]+c[j]*d[j]-s[j]-c[i+]*d[j]+s[i+]+C[i];
      while(he<tl&&slp(q[tl-],q[tl])>=slp(q[tl-],i))tl--;
      q[++tl]=i;
    }
  printf("%lld\n",dp[]);
  return ;
}