1096: [ZJOI2007]仓库建设

1096: [ZJOI2007]仓库建设

思路

　　斜率优化。

代码

#include<cstdio>
#include<iostream>
 
using namespace std;
typedef long long LL;

const int N = 1000100;
LL f[N],s[N],x[N],p[N],c[N];
int q[N],L,R;

inline int read() {
    int x = 0,f = 1;char ch = getchar();
    for (; !isdigit(ch); ch=getchar()) if(ch=='-') f=-1;
    for (; isdigit(ch); ch=getchar()) x = x*10+ch-'0';
    return x * f; 
}
double Slope(int i,int j) {
    return 1.0*(f[i]+s[i]-f[j]-s[j])/(1.0*(p[i]-p[j]));
}
int main() {
    int n = read();
    for (int i=1; i<=n; ++i) 
        x[i] = read(),p[i] = read(),c[i] = read();
    for (int i=1; i<=n; ++i) {
        s[i] = s[i-1] + 1ll * x[i] * p[i];
        p[i] += p[i-1];
    }
    int L = 1,R = 1;
    q[1] = 0;
    for (int i=1; i<=n; ++i) {
        while (L<R && Slope(q[L],q[L+1]) < (double)x[i]) L++;
        int j = q[L];
        f[i] = f[j] + x[i] * (p[i-1] - p[j]) - (s[i-1] - s[j]) + c[i];
        while (L<R && Slope(q[R-1],q[R]) > Slope(q[R],i)) R--;
        q[++R] = i; 
    }
    cout << f[n];
    return 0;
}