BZOJ 4237 稻草人

Problem

BZOJ

Solution

不知道为什么网上的题解每一篇都要强调一句“此题很简单”，我觉得不是很简单啊qwqqq

按x坐标排序，然后进行分治，每次只统计跨立两个区间的矩形个数，为了避免重复统计，我们仅统计从左下角到右下角的矩形。枚举左边的点，显然我们只需要考虑y坐标比它大的点(当然，x坐标也比它大)，那么我们可以这样按照y坐标从上到下进行归并排序。

由于矩形中不允许有其他的稻草人，那么考虑左边其他点对矩形的影响。假设有点p，点q，其中q在p的右下方，即满足$x_p<x_q,y_p>y_q$，那么我们会发现如果一个跨立的矩形如果包含p点，那么就一定会包含q点，那么加入q点后我们就可以不考虑p点的影响了，可以画个图理解一下。而且又因为我们是按照y轴从大到小加点的，所以这其实就相当于要维护一个x逐渐减小的单调栈。这个单调栈里的点就是对当前点有限制的点。
再来考虑右边的话，就会发现如果q在p的左下方，那么包含p的矩形一定会包含q，同样的我们认为加入q之后就不需要考虑p的影响了，那么右边是一个x逐渐增大的单调栈。这个单调栈里的点是可能与当前点匹配的点。

当我们去考虑当前点A能匹配的点B时，点B的y不能超过A的前一个点，而右边又是单调的，二分这个位置即可。

时间复杂度$O(n\log^2 n)$

Code

#include <algorithm>
#include <cstdio>
#define rg register
using namespace std;
typedef long long ll;
const int maxn=200010;
template <typename Tp> inline void getmin(Tp &x,Tp y){if(y<x) x=y;}
template <typename Tp> inline void getmax(Tp &x,Tp y){if(y>x) x=y;}
template <typename Tp> inline void read(Tp &x)
{
    x=0;int f=0;char ch=getchar();
    while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
    if(ch=='-') f=1,ch=getchar();
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    if(f) x=-x;
}
struct point{
    int x,y;
    bool operator < (const point &t)const{return x<t.x;}
}p[maxn],tmp[maxn];
struct Stack{
    int tp;point a[maxn];
    void clear(){tp=0;}
    point &top(){return a[tp];}
    void push(point x){a[++tp]=x;}
    void pop(){tp--;}
}ls,rs;
int n,tp;
ll ans;
int find(int k)
{
    int res=rs.tp+1,l=1,r=rs.tp+1,m;
    while(l<=r)
    {
        m=(l+r)>>1;
        if(rs.a[m].y<k) res=m,r=m-1;
        else l=m+1;
    }
    return res;
}
void cdq(int l,int r)
{
    int temp;
    if(l>=r) return ;
    int m=(l+r)>>1,i=l,j=m+1;
    cdq(l,m);cdq(m+1,r);
    ls.clear();rs.clear();
    for(;i<=m;i++)
    {
        while(ls.tp&&ls.top().x<p[i].x) ls.pop();
        ls.push(p[i]);
        while(j<=r&&p[j].y>=p[i].y)
        {
            while(rs.tp&&rs.top().x>p[j].x) rs.pop();
            rs.push(p[j++]);
        }
        if(ls.tp==1) ans+=(temp=rs.tp);
        else ans+=(temp=rs.tp-find(ls.a[ls.tp-1].y)+1);
    }
    i=l;j=m+1;tp=0;
    while(i<=m&&j<=r)
    {
        if(p[i].y>p[j].y) tmp[++tp]=p[i++];
        else tmp[++tp]=p[j++];
    }
    while(i<=m) tmp[++tp]=p[i++];
    while(j<=r) tmp[++tp]=p[j++];
    for(rg int i=1;i<=tp;i++) p[l+i-1]=tmp[i];
}
int main()
{
    #ifndef ONLINE_JUDGE
    freopen("in.txt","r",stdin);
    #endif
    read(n);
    for(rg int i=1;i<=n;i++){read(p[i].x);read(p[i].y);}
    sort(p+1,p+n+1);
    cdq(1,n);
    printf("%lld\n",ans);
    return 0;
}