【bzoj4237】【分治】稻草人

稻草人

Time Limit: 40 Sec Memory Limit: 256 MB

Description

JOI村有一片荒地，上面竖着N个稻草人，村民们每年多次在稻草人们的周围举行祭典。
有一次，JOI村的村长听到了稻草人们的启示，计划在荒地中开垦一片田地。和启示中的一样，田地需要满足以下条件：
田地的形状是边平行于坐标轴的长方形；
左下角和右上角各有一个稻草人；
田地的内部(不包括边界)没有稻草人。
给出每个稻草人的坐标，请你求出有多少遵从启示的田地的个数

Input

第一行一个正整数N，代表稻草人的个数
接下来N行，第i行(1<=i<=N)包含2个由空格分隔的整数Xi和Yi，表示第i个稻草人的坐标

Output

输出一行一个正整数，代表遵从启示的田地的个数

Sample Input

4
0 0
2 2
3 4
4 3

Sample Output

3

HINT

所有满足要求的田地由下图所示：
1<=N<=2*10^5
0<=Xi<=10^9(1<=i<=N)
0<=Yi<=10^9(1<=i<=N)
Xi(1<=i<=N)互不相同。
Yi(1<=i<=N)互不相同。

Source

JOI 2013~2014 春季training合宿 竞技3 By PoPoQQQ

网上说的CDQ分治，然而我并不知道什么是CDQ分治，也不会CDQ分治，感觉这道题实际上和归并求逆序对差不多，大概说一下

---

先按y排序，为了把图分开，好进行分治，然后再归并使x有序，然后利用栈的思想，每次维护mid+1~r y单增，l~mid单减，然后再枚举上面每两个点之间有多少个点，这个就是答案

图什么的画得太丑了，不管了
代码见下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<set>
#include<map>
#include<queue>
#include<algorithm>
#include<vector>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<stack>
#define INF 2100000000
#define ll long long
#define clr(x)  memset(x,0,sizeof(x))
#define maxclr(x)  memset(x,127,sizeof(x))
#ifdef WIN32
#define AUTO "%I64d"
#else
#define AUTO "%lld"
#endif

using namespace std;

inline int read()
{
    char c;
    int ret=0;
    while(!(c>='0'&&c<='9'))
        c=getchar();
    while(c>='0'&&c<='9')
    {
        ret=(c-'0')+(ret<<1)+(ret<<3);
        c=getchar();
    }
    return ret;
}

#define M 200005

struct node
{
    int x,y;
}a[M],b[M];

bool com(node A,node B)
{
    return A.y<B.y;
}

int n,t1[M],t2[M];
ll ans;

int fid(int x,int l,int r)
{
    int mid;
    while(r>l+1)
    {
        mid=(l+r)>>1;
        if(a[t2[mid]].x>=x)r=mid;
        else l=mid;
    }
    return l;
}

void work(int l,int r)
{
    if(l==r)return ;
    int mid=(l+r)>>1;
    work(l,mid);
    work(mid+1,r);
    int j=l,p1=0,p2=0,k;
    for(int i=mid+1;i<=r;i++)
    {
        while(p1&&a[i].y<a[t1[p1]].y)p1--;
        t1[++p1]=i;
        for(;a[j].x<a[i].x&&j<=mid;j++)
        {
            while(p2&&a[j].y>a[t2[p2]].y)p2--;
            t2[++p2]=j;
        }
        ans+=p2-fid(a[t1[p1-1]].x,0,p2+1);
    }
    j=l;k=mid+1;
    for(int i=l;i<=r;i++)
        b[i]=((j<=mid&&a[j].x<a[k].x)||k>r)?a[j++]:a[k++];
    for(int i=l;i<=r;i++)
        a[i]=b[i];
}

int main()
{
    freopen("in.txt","r",stdin);
    freopen("out.txt","w",stdout);
    n=read();
    for(int i=1;i<=n;i++)
        a[i].x=read(),a[i].y=read();
    sort(a+1,a+1+n,com);
    work(1,n);
    printf(AUTO,ans);
    return 0;
}

大概就是这个样子，如果有什么问题，或错误，请在评论区提出，谢谢。