CDQ分治

已于 2022-10-10 12:19:13 修改
阅读量1.4k 收藏
点赞数
分类专栏： 数据结构 文章标签： 算法 c++
于 2022-10-09 21:13:30 首次发布
;
本文链接：https://blog.csdn.net/m0_60630928/article/details/127233927
版权

CDQ

1.使用cdq分治的条件： 修改操作对询问的贡献独立，修改操作相互不影响 题目可以使用离线算法，不必强制在线（询问次数可以保存在数组）
2.cdq分治的性质： cdq分治通过对时间复杂度增加一个log来降维 cdq可以用来代替复杂的数据结构 在cdq分治中，对于划分出来的两个区间，前一个子问题需要用来解决后一个子问题。
3. cdq使用步骤： 将整个操作序列分为两个长度相等的部分。 递归处理前一部分的子问题（治1） 计算前一部分子问题的修改操作对后一部分子问题的影响（治2） 递归处理后一部分的子问题（治3）

离线，处理点对关系，常用于降维度，
归并排序（合并左右区间，逐级向上）：

#include<bits/stdc++.h>
using namespace std;
 
int temp[500005],a[500005],n;
long long ans;
inline int read()
{
	int x=0,f=1;
	char ch=getchar();
	while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
	while (ch>='0'&&ch<='9'){x=x*10+ch-48;ch=getchar();}
	return x*f;
}
 
void merge(int l,int r)
{
	if(l==r) return;
	int mid=(l+r)/2;
	merge(l,mid);
	merge(mid+1,r);
	int i=l,j=mid+1,p=l;
	while(i<=mid&&j<=r)
	{
		if(a[i]<=a[j])
			temp[p++]=a[i++];
		else
			ans+=(mid-i+1),//计算逆序数 
			temp[p++]=a[j++];
	}
	while(i<=mid)
		temp[p++]=a[i++];
	while(j<=r)
		temp[p++]=a[j++];
	for(int k=l;k<=r;k++)
		a[k]=temp[k];
}
 
int main()
{
	n=read();
	for(int i=1;i<=n;i++)
		a[i]=read();
	merge(1,n);
	cout<<ans;
	return 0;
}

二维偏序

P2717 寒假作业
在这里插入图片描述

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+10;

int a[N],s[N],n,k;

long long cdq(int l,int r)
{
	long long res=0;
	if(l==r) return a[l]>=0;
	int mid=(l+r)>>1;
	res+=cdq(l,mid);
	res+=cdq(mid+1,r);
	s[mid]=a[mid],s[mid+1]=a[mid+1];
	for(int i=mid-1;i>=l;i--) s[i]=s[i+1]+a[i];//左半部分求后缀和
	for(int i=mid+2;i<=r;i++) s[i]=s[i-1]+a[i];//右半部分求前缀和
	sort(s+l,s+mid+1),sort(s+mid+1,s+r+1);
	int p=l,q=r;
	while(p<=mid)
	{
		while(q>=mid+1&&s[p]+s[q]>=0) q--;
		res+=r-q;//左端点固定，符合条件的右端点数
		p++;
	}
	return res;
}

int main()
{
	cin>>n>>k;
	for(int i=1;i<=n;i++)
	{
		scanf("%d",&a[i]);
		a[i]-=k;
	}
	cout<<cdq(1,n);
	return 0;
}

P5459 [BJOI2016]回转寿司

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+10;

long long a[N];
int n,L,R;

long long cdq(int l,int r)
{
	long long res=0;
	if(l==r) return a[l]>=L&&a[l]<=R;
	int mid=(l+r)>>1;
	res+=cdq(l,mid);
	res+=cdq(mid+1,r);
	int head=l,tail=l-1;
	for(int i=mid+1;i<=r;i++)
	{
		while(tail+1<=mid&&a[i]-a[tail+1]>=L) tail++;//最终a[i]-a[tail]>=L
		while(head<=mid&&a[i]-a[head]>R) head++;//最终a[i]-a[head]<=R
		res+=tail-head+1;//对于每个区间右端点，可行的左端点个数
	}
	sort(a+l,a+r+1);
	return res;
}

int main()
{
	cin>>n>>L>>R;
	for(int i=1;i<=n;i++)
		scanf("%ld",&a[i]),a[i]+=a[i-1];//前缀和数组
	cout<<cdq(1,n);
	return 0;
}

三维偏序

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

#include<cstdio>
#include<algorithm>
#define maxn 200005
using namespace std;

struct node
{
	int a,b,c,cnt,ans;
}s1[maxn],s2[maxn];
int n,m,k,mx,top,su[maxn];
int c[maxn];//树状数组

bool cmp1(node x,node y)
{
	if(x.a==y.a)
	{
		if(x.b==y.b)return x.c<y.c;
		else return x.b<y.b;
	}
	else return x.a<y.a;
}//第一维排序

bool cmp2(node x,node y)
{
	if(x.b==y.b)
	return x.c<y.c;
	else return x.b<y.b;
}//第二维排序

int lowbit(int x)
{
	return x&(-x);
}


void add(int pos,int val)
{
	for(int i=pos;i<=mx;i+=lowbit(i))
		c[i]+=val;
}

int query(int pos)
{
	int tot=0;
	for(int i=pos;i;i-=lowbit(i))
		tot+=c[i];
	return tot;
}

void cdq(int l,int r)
{
	if(l==r) return;
	int mid=(l+r)>>1;
	cdq(l,mid);
	cdq(mid+1,r);
	sort(s2+l,s2+mid+1,cmp2);
	sort(s2+mid+1,s2+r+1,cmp2);
	int i=l,j;
	for(int j=mid+1;j<=r;j++)
	{
		while(s2[j].b>=s2[i].b&&i<=mid)
		{
			add(s2[i].c,s2[i].cnt);
			i++;
		}
		s2[j].ans+=query(s2[j].c);
	}
	for(int p=l;p<i;p++)
		add(s2[p].c,-s2[p].cnt);
}

int main()
{
	scanf("%d%d",&n,&k);
	mx=k;//树状数组的区间
	for(int i=1;i<=n;++i)
		scanf("%d%d%d",&s1[i].a,&s1[i].b,&s1[i].c);
	sort(s1+1,s1+1+n,cmp1);//第一维为关键字排序
	for(int i=1;i<=n;++i)
	{
		top++;
		if(s1[i].a!=s1[i+1].a||s1[i].b!=s1[i+1].b||s1[i].c!=s1[i+1].c)
		{
			s2[++m]=s1[i];
			s2[m].cnt=top;
			top=0;
		}
	}//第一维已有序,合并相同节点
	cdq(1,m);//cdq分治
	for(int i=1;i<=m;++i)
		su[s2[i].ans+s2[i].cnt-1]+=s2[i].cnt;
	for(int i=0;i<n;++i)
	printf("%d\n",su[i]);
	return 0;
}

P3157 [CQOI2011]动态逆序对

对答案有贡献的点 （ i , j ） （i,j） （i,j）对满足的条件：
v a l i < v a l j , p o s i > p o s j , t i m e i < = t i m e j val_i<val_j, pos_i>pos_j, time_i<=time_j vali​<valj​,posi​>posj​,timei​<=timej​
那么这个问题就变成了经典的三维偏序问题，可以通过cdq分治来解决。
先对time维排序（不需要额外的处理），在time有序的情况下使得归并区间内的pos有序，将val加到树状数组中。

#include<bits/stdc++.h>
using namespace std;
const int N=2e5+10;
struct Node
{
	int val,pos,time,cnt;
}e[N];
bool cmp1(Node x,Node y)
{
	return x.pos<y.pos;
}
int pos[N],n,m,tot,c[N];
//树状数组在val维上建 
long long ans[N];
int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9') ch=='-'&&(f=-1),ch=getchar();
    while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}

int lowbit(int x)
{
	return x&(-x);
}

void add(int pos,int x)
{
	for(int i=pos;i<=n;i+=lowbit(i))
		c[i]+=x;
}

int query(int pos)
{
	int res=0;
	for(int i=pos;i;i-=lowbit(i))
		res+=c[i];
	return res;
}
void cdq(int l,int r)
{
	if(l==r) return;
	int mid=(l+r)>>1;
	cdq(l,mid);
	cdq(mid+1,r);
	sort(e+l,e+mid+1,cmp1);
	sort(e+mid+1,e+r+1,cmp1);//按照pos排序
	int i=l,j;//i为左部分指针，j为右部分指针 
	for(j=mid+1;j<=r;j++)
	{
		while(i<=mid&&e[i].pos<=e[j].pos)//time_i<=time_j&&pos_i<=pos_j
		{
			add(e[i].val,e[i].cnt);
			i++;
		} 
		ans[e[j].time]+=1ll*e[j].cnt*(query(n)-query(e[j].val));
		//(query(n)-query(e[j].val)):pos_i<=pos_j&&val_i>val_j形成的逆序对 
	}
	for(int p=l;p<i;p++)
		add(e[p].val,-e[p].cnt);//还原 
	i=mid;
	for(int j=r;j>mid;j--)
	{
		while(i>=l&&e[i].pos>=e[j].pos)//time_i<=time_j&&pos_i>=pos_j
		{
			add(e[i].val,e[i].cnt);
			i--;
		}
		ans[e[j].time]+=1ll*e[j].cnt*query(e[j].val-1); 
		// query(e[j].val-1):pos_i>=pos_j&&val_i<val_j形成的逆序对 
	}
	for(int p=mid;p>i;p--)
		add(e[p].val,-e[p].cnt);//还原 
}

int main()
{
	cin>>n>>m;
	for(int i=1;i<=n;i++)
	{
		int x=read();
		pos[x]=i;
		e[++tot]={x,i,0,1};
	}
	for(int i=1;i<=m;i++)
	{
		int x=read();
		e[++tot]={x,pos[x],i,-1};
	}
	//此时 time维 有序 
	cdq(1,tot);
	//此时ans[0]是原始数列逆序对数，之后的 ans[1]~ans[m]都<=0 
	for(int i=1;i<=m;i++) ans[i]+=ans[i-1];
	for(int i=0;i<m;i++) printf("%lld\n",ans[i]);
	return 0;
}

P4169 [Violet]天使玩偶/SJY摆棋子

只考虑左下角的所有点，对答案有贡献的点 （ i , j ） （i,j） （i,j）对满足的条件：
x i < x j , y i < y j , t i < t j x_i<x_j, y_i<y_j, t_i<t_j xi​<xj​,yi​<yj​,ti​<tj​
那么这个问题就变成了经典的三维偏序问题，可以通过cdq分治来解决。
先对time维排序（不需要额外的处理），在time有序的情况下使得归并区间内的x有序，将(x+y)加到y轴的树状数组中。

#include<bits/stdc++.h>
using namespace std;
const int N=3e5+10;
const int MAX=2e6+10;
int n,m,tot,num;
struct Node
{
	int x,y,t,type;
}e[N];

int c[MAX],maxn;
int ques[N],ans[N];

int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9') ch=='-'&&(f=-1),ch=getchar();
    while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}

bool cmp(Node a,Node b)
{
	return a.x<b.x;
}

int lowbit(int x)
{
	return x&(-x);
}

void add(int pos,int val)
{
	for(int i=pos;i<=maxn;i+=lowbit(i))
		c[i]=max(c[i],val);
}

void cls(int pos)
{
	for(int i=pos;i<=maxn;i+=lowbit(i))
		c[i]=0;
}

int query(int pos)
{
	int res=0;
	for(int i=pos;i;i-=lowbit(i))
		res=max(c[i],res);
	return res;
}

void cdq(int l,int r)
{
	if(l==r) return;
	int mid=(l+r)>>1;
	cdq(l,mid);
	cdq(mid+1,r);
	sort(e+l,e+mid+1,cmp);
	sort(e+mid+1,e+r+1,cmp);
	int i=l,j;
	for(j=mid+1;j<=r;j++)
	{
		while(i<=mid&&e[i].x<=e[j].x)
		{
			if(e[i].type==1)
			{
				//cout<<"i:"<<i<<' ';
				//cout<<e[i].x<<' '<<e[i].x+e[i].y<<endl;
				add(e[i].y,e[i].x+e[i].y);
			}
			i++;
		}
		if(e[j].type==2)
		{
			for(int k=1;k<=maxn;k++)
				cout<<" c["<<k<<"]:"<<c[k];
			cout<<endl;
			int temp=query(e[j].y);
			cout<<"query("<<e[j].y<<"):"<<temp<<endl;
			ans[e[j].t]=min(e[j].x+e[j].y-temp,ans[e[j].t]);
		}
	}
	for(int p=l;p<i;p++)
		cls(e[p].y);
}
int main()
{
	cin>>n>>m;
	memset(ans,0x3f,sizeof(ans));
	for(int i=1;i<=n;i++)
	{
		int x=read()+1,y=read()+1;
		e[++tot]={x,y,0,1};
		maxn=max(maxn,x+y);
	}
	for(int i=1;i<=m;i++)
	{
		int op=read(),x=read()+1,y=read()+1;
		e[++tot]={x,y,i,op};
		if(op==2) ques[++num]=i;
		maxn=max(maxn,x+y);
	}
	maxn++;
	cdq(1,tot);
	for(int i=1;i<=num;i++)
		printf("%d\n",ans[ques[i]]);
	return 0;
}

将点坐标对称，每次更新答案，可以求出最小值。

#include<bits/stdc++.h>
using namespace std;
const int N=6e5+10;
const int MAX=3e6+10;
int n,m,tot,num,cnt;
struct Node
{
	int x,y,t,type;
}e[N],a[N];

int c[MAX],maxn=-1e9;
//在y维上建树状数组 
int ques[N],ans[N];

int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9') ch=='-'&&(f=-1),ch=getchar();
    while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}

void init(int d)
{
    int mxx=0,myy=0;cnt=0;
    for(int i=1;i<=n+m;++i)
    {
        if(d==1) a[i].x=maxn-a[i].x;//横坐标全部对称过去 
        else if(d==2) a[i].y=maxn-a[i].y;//纵坐标全部对称过去 
        if(a[i].type==2) mxx=max(mxx,a[i].x),myy=max(myy,a[i].y);//更新边界 
    }
    for(int i=1;i<=n+m;++i)
    	if(a[i].x<=mxx&&a[i].y<=myy)//可以被询问到的点 
    		e[++cnt]=a[i];
}

bool cmp(Node a,Node b)
{
	return a.x<b.x;
}

int lowbit(int x)
{
	return x&(-x);
}

void add(int pos,int val)
{
	for(int i=pos;i<=maxn;i+=lowbit(i))
		c[i]=max(c[i],val);
}

void cls(int pos)
{
	for(int i=pos;i<=maxn;i+=lowbit(i))
		c[i]=0;
}

int query(int pos)
{
	int res=0;
	for(int i=pos;i;i-=lowbit(i))
		res=max(c[i],res);
	return res;
}

void cdq(int l,int r)
{
	if(l==r) return;
	int mid=(l+r)>>1;
	cdq(l,mid);
	cdq(mid+1,r);
	sort(e+l,e+mid+1,cmp);//左半部分按x排序 
	sort(e+mid+1,e+r+1,cmp);//右半部分按x排序 
	int i=l,j;
	for(j=mid+1;j<=r;j++)
	{
		while(i<=mid&&e[i].x<=e[j].x)
		{
			if(e[i].type==1)
				add(e[i].y,e[i].x+e[i].y);//更新区域最大值 
			i++;
		}
		if(e[j].type==2)
		{
			//for(int k=1;k<=maxn;k++)
			//	cout<<" c["<<k<<"]:"<<c[k];
			//cout<<endl;
			int temp=query(e[j].y);//查询0~y之间的区域最大值 
			//cout<<"query("<<e[j].y<<"):"<<temp<<endl;
			if(temp) ans[e[j].t]=min(e[j].x+e[j].y-temp,ans[e[j].t]);
		}
	}
	for(int p=l;p<i;p++)
		if(e[p].type==1) cls(e[p].y);//还原，清零 
}

int main()
{
	n=read(),m=read();
	for(int i=1;i<=n;i++)
	{
		int x=read()+1,y=read()+1;//题给数据有0 
		a[++tot]={x,y,0,1};//初始t==0 
		maxn=max(maxn,x);
		maxn=max(maxn,y);
	}
	for(int i=1;i<=m;i++)
	{
		int op=read(),x=read()+1,y=read()+1;
		a[++tot]={x,y,i,op};
		if(op==2) ques[++num]=i,ans[i]=1e9;//储存下询问 
		maxn=max(maxn,x);
		maxn=max(maxn,y);
	}
	maxn++;
	//此时t维已经有序 
	init(0); cdq(1,cnt);
    init(1); cdq(1,cnt);
    init(2); cdq(1,cnt);
    init(1); cdq(1,cnt);
	for(int i=1;i<=num;i++)
		printf("%d\n",ans[ques[i]]);
	return 0;
}
春弦_
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
CDQ分治详解
m0_57344422的博客
05-22 1522
CDQ分治
CDQ分治总结（CDQ，树状数组，归并排序）
a83229442的博客
07-28 305
闲话 CDQ是什么？ 是一个巨佬，和莫队、HJT（不是我这个蒟蒻）一样，都发明出了在OI中越来越流行的算法/数据结构。 CDQ分治思想 分治就是分治，“分而治之”的思想。 那为什么会有CDQ分治这样的称呼呢？ 这一类分治有一个重要的思想——用一个子问题来计算对另一个子问题的贡献。 有了这种思想，就可以方便地解决更复杂的问题。 这样一句话怎样理解好呢？还是做做题目吧。 例题1 三维偏序问题...
[cdq分治习题练习]
MekakuCityActors的博客
08-20 539
bzoj1935 题意：给出n棵树的位置(xi,yi)，有m次询问，每次询问一个以(Xi,Yi)为左下角，(Li,Ri)为右上角的矩形内树的个数 题解：显然的三维偏序问题，一维是询问的时间T，一维是横坐标，一维是纵坐标。把初始的n棵树的位置当作插入，把询问当成4个二维前缀和相加减，由于时间是按照输入的顺序，所以第一维不需要排序，直接cdq分治处理第二维，树状数组维护第三维即可 /*********...
CDQ分治与整体二分
hipamp的博客
04-01 357
前言 CDQ分治和整体二分都是离线算法，并且都是基于分治的思想。大部分非强制在线的树套树题都可以使用整体二分或CDQ分治解决，当然也有一些是不行的，比如动态求区间种类数或者动态求区间 mexmexmex。 CDQ 分治 分治 首先要了解什么是分治，分治的全称为分而治之，英文名为 divide and conquerdivide\ and\ conquerdivide and conquer，和分治相关的算法有很多，譬如线段树，点分治，归并排序等。下面以归并排序为例。
BZOJ 1176[Balkan2007]Mokia （cdq分治，矩阵加矩阵求和）
Dch19990825的博客
09-25 353
BZOJ 1176[Balkan2007]Mokia （cdq分治解三维偏序） Description 维护一个W*W的矩阵，初始值均为S.每次操作可以增加某格子的权值,或询问某子矩阵的总权值.修改操作数M<=160000,询问数Q<=10000,W<=2000000. Input 第一行两个整数,S,W;其中S为矩阵初始值;W为矩阵大小 接下来每行为一下三种输入之一(不包含引号...
【教程】简易CDQ分治教程&学习笔记
weixin_34040079的博客
12-25 162
前言 　　辣鸡蒟蒻__stdcall终于会CDQ分治啦！ CDQ分治是我们处理各类问题的重要武器。它的优势在于可以顶替复杂的高级数据结构，而且常数比较小；缺点在于必须离线操作。 　　CDQ分治的基本思想和实现都很简单，但是因为没有人给本蒟蒻详讲，所以我对着几篇论文头疼了一个下午，最终在menci和sxysxy大佬的帮助下学会了CDQ分治。本文介绍一些非常simple的CDQ分治问...
cdq分治
tanjunming2020的博客
09-28 1116
### 介绍 $cdq$分治是一种分治思想的算法，一般用来求二维偏序的问题，可以将$O(n^2)$的时间复杂度优化到$O(nlogn)$ $cdq$分治的主要思想就是将要求的一段区间分为两个区间，先处理两个区间内的贡献，再处理两个区间相互的影响。因为左右区间相互独立，左区间在前，右区间在后，所以只有左区间影响右区间。只要在每次处理区间的时候处理左区间对右区间的影响，就能求出最后的答案。
任务安排 cdq分治
最新发布
12-23
任务安排CDQ分治是一种常用的算法设计技巧，主要用于解决一些复杂的问题。CDQ分治算法通常用于解决一些与区间相关的问题，如区间覆盖、区间查询等。 CDQ分治算法的思想是将原始问题分解成若干个子问题，然后分别解决这些子问题，最后将子问题的结果合并得到原始问题的解。在这个过程中，CDQ分治算法通常会利用一些优化技巧，如树状数组、线段树等数据结构，以提高算法的效率。 CDQ分治算法的具体步骤通常包括三个阶段：分解阶段、求解阶段和合并阶段。在分解阶段，原始问题会被分解成若干个子问题；在求解阶段，每个子问题会被独立求解；最后，在合并阶段，子问题的结果会被合并得到原始问题的解。 总的来说，CDQ分治算法是一种高效的算法设计技巧，能够有效地解决一些复杂的问题。通过合理地分解和求解问题，并利用一些优化技巧，CDQ分治算法能够在较短的时间内得到问题的解，因此在实际的问题求解中得到了广泛的应用。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

春弦_

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值