cdq分治(知识整理+板子总结)

思路来源

https://blog.csdn.net/wu_tongtong/article/details/78785836

https://blog.csdn.net/K_R_forever/article/details/81702934

https://www.cnblogs.com/flashhu/p/9381075.html

知识点整理

cdq分治，解决二维偏序或三维偏序类的问题，起到降维的作用

常数小，往往比树套树或kdtree更优，是一种分治

①cdq(l,mid)，cdq(mid+1,r)，解决子问题

②类似归并排序，考虑左子s对右子t的影响，分别对s和t对应处理

考虑到有贡献的对(l,r)，只在分治的一个区间中，

会被划分到左右不同区里，故每对答案只会被计一次

 

例如，统计逆序对(l,r)，

如果左子s大，说明[s,mid]都对此t有贡献，ans+=mid-s+1

 

使用时，先对第一维排升序，

倘若第一维排升序后是1 2 3 4 5 6 7 8，那么[l,mid]和[mid+1,r]在归并的时候

[l,mid]段的第一维，一定比[mid+1,r]小，

虽然合并之后会打乱第一维升序(因为合并是按第二维升序)，

但合并过程中时第一维保持增序，故舍掉了一维

 

第二维类似归并，归并的时候，第二维是有序的，如果二维偏序就计二维的贡献

如果有第三维，就BIT维护第三维，在树状数组中，再统计第三维偏序的相对贡献

 

题目

洛谷 P3374 【模板】树状数组 1

树状数组 n个数 m次单点加修改 和 询问区间求和 

 

题解

树状数组裸题实现单点加和区间求和，

把区间[x,y]求和，前缀和作差，看成是两个独立的[1,y]*1和[1,x-1]*(-1)，

先插前面位置的，后面端点就能统计到

考虑二维偏序(时间，位置)，默认时间增序排序，去掉一维

剩下一维cdq，重载小于号，认为位置靠前的小，位置相同的先改后插

对于有贡献的(l,r)，左子s对右子t的影响(l<=s<=mid,mid+1<=t<=r)

如果s<t且s是单点加，说明s对于[1,t]到[1,r]的前缀和的都有影响，sum+=s.v

如果t<s，说明后面的s都比t大，即此区间内比t小的都统计完，直接计算t的答案

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
const int N=5e5+10;
typedef long long ll;
int n,m,cnt,tot;
int op,x,y;
ll ans[N];
struct node
{
	int op,pos,x;
	bool operator<(const node &a)
	{
		if(pos!=a.pos)return pos<a.pos;
		else return op<a.op;
	}
}e[N*3],b[N*3];

void cdq(node *a,int l,int r)
{
	if(l==r)return;
	int m=(l+r)/2;
	cdq(a,l,m);
	cdq(a,m+1,r);
	int s=l,t=m+1;
	ll sum=0;//sum是本区间的贡献 
	for(int i=l;i<=r;++i)
	{
		if(t>r||s<=m&&a[s]<a[t])
		{
			//s小 可能还有比t大的 s为t统计答案 
			//时间序在前 且位置在前(位置相同 先加后统计) 
			if(a[s].op==1)sum+=a[s].x;
			b[i]=a[s++]; 
		}
		else
		{
			//在t位置前的 在该区间内都已被统计 故计算t的答案
			if(a[t].op==2)ans[a[t].x]-=sum;
			else if(a[t].op==3)ans[a[t].x]+=sum; 
			b[i]=a[t++];
		}
	} 
	for(int i=l;i<=r;++i)
	a[i]=b[i];
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;++i)
	{
		scanf("%d",&y);
		e[++cnt]=node{1,i,y};
	}
	for(int i=1;i<=m;++i)
	{
		scanf("%d%d%d",&op,&x,&y);
		if(op==1)e[++cnt]=node{1,x,y};
		else if(op==2)
		{
			e[++cnt]=node{2,x-1,++tot};
			e[++cnt]=node{3,y,tot};
		}
	}
	cdq(e,1,cnt);
	for(int i=1;i<=tot;++i)
	printf("%lld\n",ans[i]);
	return 0;
}

题目

洛谷 P3810 【模板】三维偏序（陌上花开）

有 n个元素，第 i 个元素有 ai​、bi​、ci​ 三个属性，

设 f(i) 表示满足 aj<=ai且bj<=bi且cj<=ci的j的数量，

对于 d∈[0,n)，求 f(i)=d 的数量

思路来源

https://blog.csdn.net/element_hero/article/details/79944523

题解

有k-dtree题解，回头再学，咕

注意到，输入可能重复，cdq处理不了完全相同的三元组的贡献，

所以合并一下，单独考虑相同的值对其内部的贡献，

如有w个完全相同的(a,b,c)，那内部有w-1个对自己产生了贡献，故形如ans+w-1的值多了w个

先按a的增序排序，在cdq中比较s.b与t.b的大小，将其归并

如果s.b<=t.b，说明t.b可能贡献还没记满，即可能还有s.b<=t.b，那先把s.b插入BIT；

如果t.b<=s.b，对于当前t，在当前区间[l,r]中不可能有s再满足s.b<=t.b，直接统计t的答案记在t.ans里

注意，BIT在每段[l,r]求完答案之后，应对[l,mid]中向BIT中插的值进行撤销操作，

这样均摊下来的复杂度才有保证O(nlogn)，不然每次都清全部会被生生卡成O(n^2logn)

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
const int N=1e5+10;
typedef long long ll;
int n,k,x,y,z,cnt;
int res[N]; 

struct node
{
	int a,b,c,w,ans;
	bool operator<(const node &x)
	{
		return a<x.a||(a==x.a&&b<x.b)||(a==x.a&&b==x.b&&c<x.c);
	}
	bool operator==(const node &x)
	{
		return a==x.a&&b==x.b&&c==x.c;
	}
}e[N],b[N];

struct BIT
{
	const static int N=2e5+10;
	int n,tr[N];
	void init(int _n)
	{
		n=_n;
		memset(tr,0,sizeof tr);
	}
	void add(int x,int v)
	{
		for(int i=x;i<=n;i+=i&-i)
		tr[i]+=v;
	}
	int sum(int x)
	{
		int sum=0; 
		for(int i=x;i>0;i-=i&-i)
		sum+=tr[i];
		return sum;
	}
}tr;

void cdq(node *a,int l,int r)
{
	if(l==r)return;
	int m=(l+r)/2;
	//分治处理子区间 
	cdq(a,l,m);
	cdq(a,m+1,r);
	//处理左区间对右区间的影响 
	int s=l,t=m+1;
	for(int i=l;i<=r;++i)
	{
		if(t>r||s<=m&&a[s].b<=a[t].b)//左区间的修改 
		{
			//s小于等于 可能还有比t大的 s为t的答案提供贡献 
			//时间序在前 且位置在前(位置相同 先加后统计) 
			tr.add(a[s].c,a[s].w);
			b[i]=a[s++];
		}
		else//右区间的查询 没有s.b<=t.b的了 统计t的答案 
		{
			//在t位置前的 在该区间内都已被统计 故计算t的答案
			a[t].ans+=tr.sum(a[t].c);
			b[i]=a[t++];
		}
	} 
	for(int i=l;i<=m;++i)
	tr.add(a[i].c,-a[i].w);
	for(int i=l;i<=r;++i)
	a[i]=b[i];
}

int main()
{
	scanf("%d%d",&n,&k);
	tr.init(k);
	for(int i=1;i<=n;++i)
	{
		scanf("%d%d%d",&x,&y,&z);
		e[i]=node{x,y,z,1,0};
	}
	sort(e+1,e+n+1);
	for(int i=1;i<=n;)
	{
		int j=i;
		while(j+1<=n&&e[j+1]==e[j])j++;//去重 完全相同的算一个 [i,j]
		e[++cnt]=e[i];
		e[cnt].w=j-i+1;
		i=j+1; 
	}
	cdq(e,1,cnt);
	for(int i=1;i<=cnt;++i)
	res[e[i].ans+e[i].w-1]+=e[i].w;//w个相同 每个都有w-1个完全相同的额外贡献 
	for(int i=0;i<n;++i)
	printf("%d\n",res[i]);
	return 0;
}