BZOJ1096: [ZJOI2007]仓库建设（洛谷P2120）

最新推荐文章于 2019-06-27 21:56:44 发布

forezxl

最新推荐文章于 2019-06-27 21:56:44 发布

阅读量197

点赞数 1

分类专栏：洛谷 BZOJ DP---斜率优化DP 蒟蒻zxl的Blog专栏文章标签： BZOJ 洛谷斜率优化

本文链接：https://blog.csdn.net/a1799342217/article/details/83477436

版权

蒟蒻zxl的Blog专栏同时被 3 个专栏收录

460 篇文章 6 订阅

订阅专栏

洛谷

292 篇文章 1 订阅

订阅专栏

BZOJ

281 篇文章 1 订阅

订阅专栏

斜率优化

BZOJ题目传送门
 洛谷题目传送门

设 $f [i]$ 为前 $i$ 个工厂且在 $i$ 建仓库的最小代价。那么有 $f[i]=min\{f[j]+c[i]+\sum_{k=j}^ip[k]*(X[i]-X[k])\}$ 。因为第 $n$ 个工厂一定要建仓库，那么答案就是 $f [n]$ 。

记 $s1[i]=\sum p[i],s2[i]=\sum X[i]*p[i]$ 。那么 $f[i]=min\{f[j]+c[i]+(s1[i]-s1[j])X[i]-s2[i]+s2[j]\}$ 。这是 $n^2$ 的，考虑斜率优化。

设 $x < y$ 且 $x$ 优于 $y$ ：
$\begin{aligned} f[x]+c[i]+(s1[i]-s1[x])X[i]-(s2[i]-s2[x])&<f[y]+c[i]+(s1[i]-s1[y])X[i]-(s2[i]-s2[y])\\ f[x]-f[y]+s2[x]-s2[y]&<(s1[x]-s1[y])X[i]\\ \frac{f[x]-f[y]+s2[x]-s2[y]}{s1[x]-s1[y]}&<X[i] \end{aligned}$
代码：

#include<cctype>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1000005
#define F inline
using namespace std;
typedef long long LL;
int n,X[N],c[N],p[N],q[N];
LL f[N],s1[N],s2[N];
F char readc(){
	static char buf[100000],*l=buf,*r=buf;
	if (l==r) r=(l=buf)+fread(buf,1,100000,stdin);
	return l==r?EOF:*l++;
}
F int _read(){
	int x=0,f=1; char ch=readc();
	while (!isdigit(ch)){ if (ch=='-') f=-1; ch=readc(); }
	while (isdigit(ch)) x=(x<<3)+(x<<1)+(ch^48),ch=readc();
	return x*f;
}
#define calc(x,y) (1.0*(f[x]-f[y]+s2[x]-s2[y])/(s1[x]-s1[y]))
int main(){
	n=_read(); int l=0,r=0;
	for (int i=1;i<=n;i++){
		X[i]=_read(),p[i]=_read(),c[i]=_read();
		s1[i]=s1[i-1]+p[i],s2[i]=s2[i-1]+1ll*p[i]*X[i];
	}
	for (int i=1;i<=n;q[++r]=i,i++){
		while (l<r&&calc(q[l],q[l+1])<X[i]) l++;
		f[i]=f[q[l]]+c[i]+(s1[i]-s1[q[l]])*X[i]-s2[i]+s2[q[l]];
		while (l<r&&calc(q[r-1],q[r])>calc(q[r],i)) r--;
	}
	return printf("%lld\n",f[n]),0;
}