[ZJOI2007]仓库建设

题面描述

传送门

思路

决策单调性

其实很容易想到方程，

$F_i=\min (F_j+c_i$$+$$[j+1,i-1]$的物品运到$i$的花费)

关键在于如何快速求到$[j+1,i-1]$的物品运到$i$的花费

运用前缀和思想，我们将$[j+1,i-1]$的物品放在原点，运到$i$的花费为(此时$p$为物品数的前缀和，下同)

$$di{s}_i\ast (p_{i-1}-p_j)$$

但是物品$[j+1,i-1]$并不在原点啊，而是在$[j+1,i-1]$点上啊，

那么我们需要再来一个前缀和$s$，来求从原点出发，物品数$j+1$到点$j+1$的花费，物品数$j+2$到点$j+2$的花费，$\cdots$，物品数$i-1$到点$i-1$的花费。

那么$[j+1,i-1]$的物品运到$i$的花费为

$$di{s}_i\ast (p_{i-1}-p_j)-s_{i-1}+s_j$$

好了，真正的方程

$$F_i=\min (F_j+c_i+di{s}_i\ast (p_{i-1}-p_j)-s_{i-1}+s_j)$$

假设有$$F_j+c_i+di{s}_i\ast (p_{i-1}-p_j)-s_{i-1}+s_j\ge F_k+c_i+di{s}_i\ast (p_{i-1}-p_j)-s_{i-1}+s_k(j<k<i)$$

化简为

$$F_j-di{s}_i\ast p_j+s_j\ge F_k-di{s}_i\ast p_k+s_k$$

①

那么对于未来状态$t$，是否也有$k$优于$j$呢?

证明:$$F_j+c_t+di{s}_t\ast (p_{t-1}-p_j)-s_{t-1}+s_j\ge F_k+c_t+di{s}_t\ast (p_{t-1}-p_j)-s_{t-1}+s_k(j<k<i<t)$$

现将原式化简

$$F_j-di{s}_t\ast p_j+s_j,F_k-di{s}_t\ast p_k+s_k$$

一定有$di{s}_t>di{s}_i$，则$di{s}_t=di{s}_i+val(val>0)$，

因此

$$F_j-(di{s}_i+val)\ast p_j+s_j,F_k-(di{s}_i+val)\ast p_k+s_k$$

$$F_j-di{s}_i\ast p_j-val\ast p_j+s_j,F_k-di{s}_i\ast p_j-val\ast p_k+s_k$$

②

将①代入②，证明$-val\ast p_j\ge -val\ast p_k$即可

$$val\ast p_j\le val\ast p_k$$

$$p_j\le p_k$$

由$p_k>p_j$得证。

因此决策单调性得出

($F_j+c_i+di{s}_i\ast (p_{i-1}-p_j)-s_{i-1}+s_j\ge F_k+c_i+di{s}_i\ast (p_{i-1}-p_j)-s_{i-1}+s_k(j<k<i)$时，$k$优于$j$，$j$没有用了)

踢队头

我们用上述①式来求解:

$$F_j-di{s}_i\ast p_j+s_j\ge F_k-di{s}_i\ast p_k+s_k$$

$$F_j-F_k+s_j-s_k\ge di{s}_i\ast (p_j-p_k)$$

由$p_j<p_k$

$$\frac{F_j-F_k+s_j-s_k}{p_j-p_k}\le di{s}_i$$

时，$k$优于$j$.

用$calc(j,k)$表示$\frac{F_j-F_k+s_j-s_k}{p_j-p_k}$,就有

$$calc(j,k)\le di{s}_i$$

如果$$calc(q_{head},q_{head+1})\le di{s}_i$$

则证明$q_{head+1}$优于$q_{head}$

踢队尾

我们首先一定要确定$i$这个状态，要优于$q_{tail}$

要维护的队列由于$q_{tail-1}$一定优于$q_{tail}$

现有一个标准值$di{s}_t(tail<t)$

$calc(q_{tail-1},q_{tail)}\ge di{s}_t$

如果$q_{tail}$要优于$i$，

则$calc(q_{tail},i)\ge calc(q_{tail-1},q_{tail})$

才能确保$q_{tail}$优于$i$

反之，当$calc(q_{tail},i)\le calc(q_{tail-1},q_{tail})$时，需要删去队尾。

貌似就没了哦。

AC code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<cmath>
#include<iostream>
#define gc getchar()
#define ll long long 
using namespace std;
const int N=1e6+10;
inline void qr(ll &x)
{
	x=0;int f=1;char c=gc;
	while(c<'0'||c>'9'){if(c=='-')f=-1;c=gc;}
	while(c>='0'&&c<='9'){x=x*10+(c^48);c=gc;}
	x*=f;
}
void qw(ll x)
{
	if(x<0)x=-x,putchar('-');
	if(x/10)qw(x/10);
	putchar(x%10+48);
}
ll s[N],p[N],c[N],f[N],dis[N];
inline double calc(int j,int k)
{
	return (f[j]-f[k]+s[j]-s[k])/(double)(p[j]-p[k]);
}
int q[N],l,r,n;
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)qr(dis[i]),qr(p[i]),qr(c[i]);
	for(int i=1;i<=n;i++)
	{
		s[i]=s[i-1]+dis[i]*p[i];
		p[i]+=p[i-1];
	}
	l=1;r=1;q[1]=0;
	for(int i=1;i<=n;i++)
	{
		while(l<r&&calc(q[l],q[l+1])<=1.0*dis[i])++l;
		f[i]=f[q[l]]+(p[i-1]-p[q[l]])*dis[i]-(s[i-1]-s[q[l]])+c[i];
		while(l<r&&calc(q[r],i)<=calc(q[r-1],q[r]))--r;
		q[++r]=i;
	}
	qw(f[n]);puts("");
	return 0;
}