CDQ分治学习笔记

CDQ分治

解析

CDQ分治主要用来解决子问题之间相互影响的问题

它一共分为四步：

1. 将当前处理区间分为左右两个等大的子区间；
2. 递归处理左子区间；
3. 处理左区间对于右区间的影响，并对于右区间或者答案进行更改与修正；
4. 递归处理右子区间；

上面就是CDQ分治的四个步骤，这四个步骤之中第一、二、四步对于不同的题目来说基本上是相同的，因为毕竟分区间，递归没有什么好更改的。对于不同的题目来说不同点就是第三步，这一步也是CDQ分治的难点，对于这一步的讲解也要借助于例题。

注意CDQ分治不支持强制在线

模板题 P3810

题解

三维偏序

代码

#include<bits/stdc++.h>//陌上花开 
#define M 200009 
using namespace std;
int read(){
	int f=1,re=0;
	char ch;
	for(ch=getchar();!isdigit(ch)&&ch!='-';ch=getchar());
	if(ch=='-'){f=-1,ch=getchar();}
	for(;isdigit(ch);ch=getchar()) re=(re<<3)+(re<<1)+ch-'0';
	return re*f;
}
struct zb{int x,y,z,num,id,sum;}a[M];
int n,k,tr[M],g[M],cnt;  
int lowbit(int i){return i&(-i);}
void update(int x,int y){
	while(x<=k){
		tr[x]+=y;
		x+=lowbit(x);
	}
}
int solve(int x){
	int ans=0;
	while(x>0){
		ans+=tr[x];
		x-=lowbit(x);
	}return ans;
}
bool check(int x,int y){
	if(a[x].x^a[y].x||a[x].y^a[y].y||a[x].z^a[y].z) return false;
	return true;
} 
bool cmp1(const zb &a,const zb &b){
	if(a.x!=b.x) return a.x<b.x;
	if(a.y!=b.y) return a.y<b.y;
	return a.z<b.z;
}
bool cmp2(const zb &a,const zb &b){
	if(a.y!=b.y) return a.y<b.y;
	if(a.z!=b.z) return a.z<b.z;
	return a.x<b.x;
}
void cdq(int l,int r){
	if(l==r) return;
	int mid=(l+r)>>1;
	cdq(l,mid),cdq(mid+1,r);
	sort(a+l,a+r+1,cmp2);
	for(int i=l;i<=r;i++)
		if(a[i].id<=mid) update(a[i].z,a[i].num);
		else a[i].sum+=solve(a[i].z);
	for(int i=l;i<=r;i++)
		if(a[i].id<=mid) update(a[i].z,-a[i].num);
}
int main(){
	n=read(),k=read();
	for(int i=1;i<=n;i++) a[i].x=read(),a[i].y=read(),a[i].z=read();
	sort(a+1,a+n+1,cmp1);
	for(int i=1;i<=n;a[cnt].num++,i++)
		if(!check(i,i-1)) a[++cnt]=a[i],a[cnt].id=cnt,a[i].num=0;
	cdq(1,cnt);
	for(int i=1;i<=cnt;i++) g[a[i].sum+a[i].num-1]+=a[i].num;
	for(int i=0;i<n;i++) printf("%d\n",g[i]);
	return 0;
}