[bzoj4237]稻草人

题目大意

给定n个点（横纵坐标是不大于$10^9$的非负整数且两两x,y都不相同）
求有多少个矩形满足：
1. 边平行于横、纵轴
2. 左下角、右上角都是给定的点
3. 内部不包含其它点

n≤200000

分析

如果枚举一个点作为右上角，那么枚举在它左下角的点，求多少个点可以满足条件。显然满足这些条件的点按横坐标排序后，纵坐标是严格下降的。

直接维护这个单调队列显然不好做，可以考虑分治：
首先对所有点按横坐标排序。
假设当前分治区间是[l,r]，先对区间[l,mid],[mid+1,r]求完答案后，求左下角在左区间，右上角在右区间的答案。其中最后一个流程相当于确定一条直线x=k，然后统计跨过这条线的答案。

首先可以对左右区间的点分别按纵坐标排序。
然后按纵坐标从大到小顺序枚举，对左区间维护横坐标单调递减的队列（不满足的点一定会被替代掉），对右区间维护横坐标单调递增的队列，然后当插入一个左下角后，右边队列中会有些不符合条件的。但是由于它们还满足了横坐标越来越小，那么满足二分性质。二分求答案即可。

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;

const int maxn=200005;

typedef long long LL;

int n,s1[maxn],t1,s2[maxn],t2;

LL ans;

struct p
{
    int x,y;
}a[maxn];

bool cmp1(p a,p b)
{
    return a.x<b.x;
}

bool cmp2(p a,p b)
{
    return a.y>b.y;
}

char c;

int read()
{
    for (c=getchar();c<'0' || c>'9';c=getchar());
    int x=c-48;
    for (c=getchar();c>='0' && c<='9';c=getchar()) x=x*10+c-48;
    return x;
}

int find(int x)
{
    int l,r,mid;
    for (l=1,r=t2,mid=l+r>>1;l<r;mid=l+r>>1)
        if (a[s2[mid]].y<x) r=mid;else l=mid+1;
    if (!t2 || a[s2[l]].y<x) l--;
    return l;
}

void work1(int x)
{
    for (;t1>0 && a[s1[t1]].x<a[x].x;t1--);
    ans+=t2;
    if (t1>0) ans-=find(a[s1[t1]].y);
    s1[++t1]=x;
}

void work2(int x)
{
    for (;t2>0 && a[s2[t2]].x>a[x].x;t2--);
    s2[++t2]=x;
}

void solve(int l,int r)
{
    if (l==r) return;
    int mid=l+r>>1;
    solve(l,mid); solve(mid+1,r);
    sort(a+l,a+mid+1,cmp2);
    sort(a+mid+1,a+r+1,cmp2);
    t1=t2=0;
    int i=l,j=mid+1;
    while (i<=mid && j<=r)
        if (a[i].y>a[j].y) work1(i++);else work2(j++);
    while (i<=mid) work1(i++);
}

int main()
{
    freopen("scarecrows.in","r",stdin); freopen("scarecrows.out","w",stdout);
    n=read();
    for (int i=1;i<=n;i++)
    {
        a[i].x=read(); a[i].y=read();
    }
    sort(a+1,a+n+1,cmp1);
    solve(1,n);
    printf("%lld\n",ans);
    fclose(stdin); fclose(stdout);
    return 0;
}