[CEOI2004]锯木厂选址

题面描述

传送门

思路

状态转移方程

$tot$为所有树都运到山下锯木厂的费用，$dis$为到山下锯木厂的距离，$s$为树的重量前缀和。

$$F_i=\min (tot-di{s}_j\ast s_j-di{s}_i\ast (s_i-s_j))(j<i)$$

不难理解,$[1,j]$的木头运到$j$处锯木厂，$[j+1,i]$运到$i$处锯木厂，剩余的运往山下锯木厂。

决策单调性

设

$$tot-di{s}_j\ast s_j-di{s}_i\ast (s_i-s_j)>tot-di{s}_k\ast s_k-di{s}_i\ast (s_i-s_k)(j<k<i)$$

化简

$$-di{s}_j\ast s_j+di{s}_i\ast s_j>-di{s}_k\ast s_k+di{s}_i\ast s_k$$

在未来状态$t$，证明:

$$-di{s}_j\ast s_j+di{s}_t\ast s_j>-di{s}_k\ast s_k+di{s}_t\ast s_k$$

由$di{s}_t=di{s}_i+val(val<0)$，

$$-di{s}_j\ast s_j+(di{s}_i+val)\ast s_j>-di{s}_k\ast s_k+(di{s}_i+val)\ast s_k$$

仅需证明:

$$val\ast s_j>val\ast s_k$$

$$s_j<s_k$$

证毕。

踢队头

$$-di{s}_j\ast s_j+di{s}_i\ast s_j>-di{s}_k\ast s_k+di{s}_i\ast s_k$$

$$di{s}_j\ast s_j-di{s}_k\ast s_k<di{s}_i\ast s_j-di{s}_i\ast s_k$$

$$calc(j,k)=\frac{di{s}_j\ast s_j-di{s}_k\ast s_k}{s_j-s_k}>di{s}_i$$

$k$优于$j$

即

$$calc(q_{head},q_{head+1})>di{s}_i$$

$q_{head+1}$优于$q_{head}$，踢去$head$

维护一个$calc(q_{head},q_{head+1})<di{s}_i$队列

同时我们可以发现，$di{s}_i$随$i$增大而减小，那么$calc(q_{head},q_{head+1})$随$di{s}_i$的减小，$head$不断增大，而不断减小。

因此，队列维护的应该是一个上凸壳。

踢队尾

由于维护一个上凸壳，

所以

$$calc(q_{tail-1},q_{tail})>calc(q_{tail},i)$$

反之，

当

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

踢去队尾。

AC code

#include<cstdio>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#include<cmath>
#define ll long long 
#define gc getchar()
using namespace std;
const int N=2e4+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;
}
inline void qw(ll x)
{
	if(x<0)x=-x,putchar('-');
	if(x/10)qw(x/10);
	putchar(x%10+48);
}
ll dis[N],s[N],tot,f[N];
int q[N],l,r;
inline double calc(int j,int k)
{
	return (double)(dis[j]*s[j]-dis[k]*s[k])/(s[j]-s[k]);
}
int main()
{
	int n;scanf("%d",&n);ll tot=0;
	for(int i=1;i<=n;i++)qr(s[i]),qr(dis[i]),s[i]+=s[i-1],tot+=s[i]*dis[i];
	for(int i=n;i>=1;i--)dis[i]+=dis[i+1];
	l=1;r=1;q[1]=0;ll ans=1e10;
	for(int i=1;i<=n;i++)
	{
		while(l<r&&calc(q[l],q[l+1])>=dis[i])++l;
		ans=min(ans,tot-dis[q[l]]*s[q[l]]-dis[i]*(s[i]-s[q[l]]));
		while(l<r&&calc(q[r],i)>=calc(q[r-1],q[r]))--r;
		q[++r]=i;
	}
	qw(ans);puts("");
	return 0;
}