[USACO08MAR]土地征用Land Acquisition

题面描述

传送门

思路

瞎扯

乍一看，貌似找不到什么有序的东西。

先排序吧(第一关键字为$x$坐标，第二关键字为$y$坐标，都以从小到大排序)。

类似下图

发现相邻的两块土地(排序后)$i,i+1$,如果$y_i<y_{i+1}$，那么买第$i+1$块土地，顺便可以带上第$i$块土地(岂不美哉！，买一送一啊！)，

按照这样的思想，我们可以买一送$n$啊。

之后去掉送的土地，保留有用的(至于如何有用，详见上文)土地。

如下图

之后我们就可以愉快地推状态转移方程啦！

设离散化之后的数组为$b$(保证$x$坐标不下降)

$$F_i=\min (F_j+y_{b_{j+1}}\ast x_{b_i})$$

$y_{b_{j+1}}$保证是$b_{[j+1,i]}$后最大的纵坐标哦，因为后面如果有比它大的纵坐标，$b_{j+1}$根本不存在于数组中。

证明决策单调性

设有
$$F_k+y_{b_{k+1}}\ast x_{b_i}\le F_j+y_{b_{j+1}}\ast x_{b_i}(j<k<i)$$

证明:对于未来状态$t$,都有

$$F_k+y_{b_{k+1}}\ast x_{b_t}\le F_j+y_{b_{j+1}}\ast x_{b_t}(j<k<i<t)$$

由于$x_{b_t},y_{b_{k+1}},y_{b_{j+1}}$不变，证毕。

这怕不是不用证明

踢队头

$$F_k+y_{b_{k+1}}\ast x_{b_i}\le F_j+y_{b_{j+1}}\ast x_{b_i}(j<k<i)$$

$$F_j-F_k\ge y_{b_{k+1}}\ast x_{b_i}-y_{b_{j+1}}\ast x_{b_i}$$

由$y_{b_{k+1}}\le y_{b_{j+1}}$可知，

$$\frac{F_j-F_k}{y_{b_{k+1}}-y_{b_{j+1}}}\le x_{b_i}$$

因此，当上述不等式成立时，$k$优于$j$。

如果$$calc(q_{head},q_{tail})=\frac{F_{head}-F_{head+1}}{y_{b_{head+1}}-y_{b_{head}}}\le x_{b_i}$$

说明$head+1$优于$head$。

踢队尾

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

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

现有一个标准值$x_{b_t}(tail<t)$

$calc(q_{tail-1},q_{tail)}\ge x_{b_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>
#define gc getchar()
#define ll long long 
using namespace std;
const int N=5e4+10;
inline void qr(ll &x)
{
	x=0;int f=1;char c=gc;
	while(c<'0'||c>'9')c=gc;
	while(c>='0'&&c<='9'){x=x*10+(c^48);c=gc;}
	x*=f;
}
void qw(ll x)
{
	if(x/10)qw(x/10);
	putchar(x%10+48);
}
struct node
{
	ll x,y;
}a[N];
bool cmp(node a,node b){return a.x==b.x?a.y<b.y:a.x<b.x;}
ll f[N];int b[N],q[N],l,r;
inline double calc(int j,int k){return (f[j]-f[k])/(double)(a[b[k+1]].y-a[b[j+1]].y);}
int main()
{
	int n;scanf("%d",&n);
	for(int i=1;i<=n;i++)qr(a[i].x),qr(a[i].y);
	sort(a+1,a+n+1,cmp);
	int tot=0;
	for(int i=1;i<=n;i++)
	{
		while(tot&&a[b[tot]].y<=a[i].y)--tot;
		b[++tot]=i;
	}
	l=1;r=1;q[1]=0;
	for(int i=1;i<=tot;i++)
	{
		while(l<r&&calc(q[l],q[l+1])<=1.0*a[b[i]].x)++l;
		f[i]=f[q[l]]+a[b[i]].x*a[b[q[l]+1]].y;
		while(l<r&&calc(q[r],i)<=calc(q[r-1],q[r]))--r;
		q[++r]=i;
	}
	qw(f[tot]);puts("");
	return 0;
}