【BZOJ1597】【USACO2008 Mar】土地购买 题解

题面：传送门

---

    我们会发现，长和宽均小于另一块土地的土地是可以直接删去的。如果将剩下来的土地按长从小到大的顺序排序的话，那么宽一定是满足从大到小的。

    这时我们可以考虑$dp$（建立在长从小到大排序的基础上），$dp_i$表示第$1$块土地到第$i$块土地都购买好了的最小总代价，则得出$dp$方程：

dpi=min{dpj+widthj+1⋅lengthi,0≤j<i} d p i = m i n { d p j + w i d t h j + 1 · l e n g t h i , 0 ≤ j < i }

$$dp_i=min\{dp_j+width_{j+1}·length_i,0\le j\lt i\}$$

    上述时间复杂度为$O(n^2)$，对于$n\le 50000$的数据，会超时，我们考虑优化。

    由于长单调增，宽单调减，所以可以对原方程进行斜率优化，$slope_{i,j}=\frac{dp_i-dp_j}{width_{j+1}-width_{i+1}}$，用单调队列维护，最终复杂度为$O(n)$。

    全部代码：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
int n;
pair<int,int>a[50010];
bool nok[50010];
int cnt;
pair<int,int>land[50010];
long long dp[50010];
double slope(int i,int j)
{
    return (double)(dp[i]-dp[j])/(double)(land[j+1].second-land[i+1].second);
}
int q[50010];
int l,r;
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        a[i]=make_pair(x,y);
    }
    sort(a+1,a+n+1);
    int mx=0;
    for(int i=n;i>=1;i--)
    {
        if(a[i].second<=mx)nok[i]=true;
        else mx=a[i].second;
    }
    for(int i=1;i<=n;i++)
    {
        if(!nok[i])land[++cnt]=a[i];
    }
    for(int i=1;i<=cnt;i++)
    {
        while(l<r && slope(q[l],q[l+1])<=land[i].first)l++;
        dp[i]=dp[q[l]]+1LL*land[q[l]+1].second*land[i].first;
        while(l<r && slope(q[r-1],q[r])>=slope(q[r],i))r--;
        q[++r]=i;
    }
    printf("%lld",dp[cnt]);
    return 0;
}