可能严格来讲不算CDQ分治把。。我也不太清楚>_<。
按照y坐标来分治,然后首先处理上下两块(下面的y坐标小)的答案。
然后对于某一段,首先按照x坐标进行排序,考虑上面的一块对下面一块产生的答案。考虑上面一块中的某一个点i,以它为右上角的矩形的左下角的x坐标是有范围的限制的(左下角在下面一块中)。显然左下角的x坐标不能<第一个在上面的离i最近的y坐标比i小的点的x坐标(有点拗口233),然后左下角的点也是有限制的,也就是最为左下角的点如果有两个点i,j使得x[i]<x[j]且y[i]<y[j],那么i是不满足的。
于是直接将上面一块的点维护一个y坐标单调递增的栈,下面的维护一个y坐标单调递减的栈,然后二分查找一下即可。
AC代码如下:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 200005
#define ll long long
using namespace std;
int n,m,p[N],q[N]; struct node{ int x,y; }a[N],b[N]; ll ans;
int read(){
int x=0; char ch=getchar();
while (ch<'0' || ch>'9') ch=getchar();
while (ch>='0' && ch<='9'){ x=x*10+ch-'0'; ch=getchar(); }
return x;
}
bool cmpy(node u,node v){ return u.y<v.y; }
int find(int x,int l,int r){
while (l+1<r){
int mid=(l+r)>>1;
if (a[q[mid]].x<x) l=mid; else r=mid;
}
return l;
}
void solve(int l,int r){
if (l==r) return; int mid=(l+r)>>1;
solve(l,mid); solve(mid+1,r);
int tp1=0,tp2=0,i,j=l,k;
for (i=mid+1; i<=r; i++){
while (tp1 && a[i].y<a[p[tp1]].y) tp1--;
p[++tp1]=i;
for (; a[j].x<a[i].x && j<=mid; j++){
while (tp2 && a[j].y>a[q[tp2]].y) tp2--;
q[++tp2]=j;
}
ans+=tp2-find(a[p[tp1-1]].x,0,tp2+1);
}
j=l; k=mid+1;
for (i=l; i<=r; i++)
b[i]=(j<=mid && a[j].x<a[k].x || k>r)?a[j++]:a[k++];
for (i=l; i<=r; i++) a[i]=b[i];
}
int main(){
n=read(); int i;
for (i=1; i<=n; i++){
a[i].x=read(); a[i].y=read();
}
sort(a+1,a+n+1,cmpy);
solve(1,n); printf("%lld\n",ans);
return 0;
}
by lych
2016.3.18