【题解】洛谷P2120 [ZJOI2007]仓库建设 斜率优化

给定长为$n$(1e6)的三个数列X（单调递增），P，C，现在将要其分成若干连续段。
每段$[l,r]$的代价是$C_r+{\Sigma }_{i=l}^r{P}_i\ast (X_r-X_i)$。
求代价之和的最小值。

---

又是分段emm

为了快速求得段代价，我们维护两个前缀和数组：
$s_i={\Sigma }_{j=1}^i{P}_k$
$t_i={\Sigma }_{j=1}^i{P}_k{X}_k$

设$d_i$表示把前$i$个数分好的最小代价，$d_0=0$
d i = m i n { d j + C i + X i ( s i − s j ) − ( t i − t j ) } ， 0 < = j < i d_i=min\{d_j+C_i+X_i(s_i-s_j) -(t_i-t_j)\}，0<=j<i di​=min{dj​+Ci​+Xi​(si​−sj​)−(ti​−tj​)}，0<=j<i

整理一下： d i = ( C i + X i s i − t i ) + m i n { d j + t j − X i s j } d_i=(C_i+X_is_i-t_i)+min\{d_j+t_j-X_is_j\} di​=(Ci​+Xi​si​−ti​)+min{dj​+tj​−Xi​sj​}

设$a<b<i$并且从$b$转移优于从$a$转移，那么：
$d_b+t_b-X_i{s}_b<d_a+t_a-X_i{s}_a$

整理得：$$\frac{(d_b+t_b)-(d_a+t_a)}{s_b-s_a}<X_i$$

由斜率优化原理：

1. $y_i=d_i+t_i，x_i=s_i$；
2. 维护下凸包；
3. 目标斜率是$X_i$，单调递增；
4. 双端队列。

---

比较基础的一道题目，可以作为模板题使用

WA点：

1. 计算dp时i和j搞混
2. 计算斜率时把i写成q[r]

看来我真的困了。。。

---

/* LittleFall : Hello! */
#include <bits/stdc++.h>
using namespace std; using ll = long long; inline int read();
const int M = 1000016, MOD = 1000000007;

ll X[M], P[M], C[M], s[M], t[M], dp[M];
inline ll subx(int a, int b){return s[b]-s[a];}
inline ll suby(int a, int b){return dp[b]+t[b]-dp[a]-t[a];}
inline ll cal(int i, int j){
	return (C[i] + X[i]*s[i] - t[i]) + (dp[j] + t[j] - X[i]*s[j]);
}
int q[M], l, r;
int main(void)
{
	#ifdef _LITTLEFALL_
	freopen("in.txt","r",stdin);
    #endif

	int n = read();
	for(int i=1; i<=n; ++i)
	{
		X[i] = read(), P[i] = read(), C[i] = read();
		s[i] = s[i-1] + P[i];
		t[i] = t[i-1] + P[i]*X[i];
	}

	q[r++] = 0;
	for(int i=1; i<=n; ++i)
	{
		// dp[i] = 1e18;
		// for(int j=0; j<i; ++j)
		// 	dp[i] = min(dp[i], cal(i,j));

		while(r-l>=2 && suby(q[l], q[l+1])<=subx(q[l], q[l+1])*X[i]) ++l;
		dp[i] = cal(i, q[l]);
		while(r-l>=2 &&
			suby(q[r-2], q[r-1])*subx(q[r-1], i)>=subx(q[r-2], q[r-1])*suby(q[r-1], i)
		) --r;
		q[r++] = i;
	}

	cout << dp[n] << "\n";

    return 0;
}


inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9') {if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}