cdq分治学习小结

最新推荐文章于 2024-02-26 17:43:38 发布
最新推荐文章于 2024-02-26 17:43:38 发布
阅读量416 收藏 1
点赞数
分类专栏： 分治 cdq分治
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Lefz_Bamboo/article/details/81944017
版权

一脸懵逼的学了好久，应该算学会了吧。基本了解cdq分治的思想，然而代码....一如既往的丑。

cdq分治只支持离线操作，主要作用是降维，反正就是大大优化了时间和空间。（感觉不扯几句作用非常不好）。

 

二维偏序

题目描述

输入长度为n的序列，并进行m次操作，操作分两种

1.将编号为的数加上y（1 x y）

2.求出某区间每一个数的和（2 x y）

典型的树状数组题，但我们可以把它看做一个二维偏序问题。这个问题一维是时间二维是每个查询或者修改的数的位置。

因为我们是按时间顺序输入的，所以我们默认时间有序。

对于左右两个区间合并，我们只执行左区间的修改以及右区间的查询。

#include<cstdio>
#include<cstring>
using namespace std;
#define N 500000+2000
int n,m,tot,ai;
int a[2*N];
struct node
{
	int pos,opr,val;
}q[3*N],tmp[3*N];
void cdq(int l,int r)
{
	if(l==r)
		return;
	int mid=(l+r)/2,ti=0,i=l,j=mid+1,cnt=0;
	cdq(l,mid),cdq(mid+1,r);
	while(i<=mid&&j<=r)
	{
		if(q[i].pos<=q[j].pos)
		{
			if(q[i].opr==1)
				cnt+=q[i].val;
			tmp[++ti]=q[i];
			i++;
		}
		else
		{
			if(q[j].opr==2)
				a[q[j].val]-=cnt;
			if(q[j].opr==3)
				a[q[j].val]+=cnt;
			tmp[++ti]=q[j];
			j++;
		}
	}
	while(i<=mid)
		tmp[++ti]=q[i],i++;
	while(j<=r)
	{
		if(q[j].opr==2)
			a[q[j].val]-=cnt;
		if(q[j].opr==3)
			a[q[j].val]+=cnt;
		tmp[++ti]=q[j];
		j++;
	}
	for(int k=1;k<=ti;k++)
		q[k+l-1]=tmp[k];
}
int main()
{

	scanf("%d %d",&n,&m);
	for(int i=1;i<=n;i++)
	{
		int x;
		scanf("%d",&x);
		q[++tot].pos=i,q[tot].opr=1,q[tot].val=x;
	}
	for(int i=1;i<=m;i++)
	{
		int x,y,op;
		scanf("%d %d %d",&op,&x,&y);
		if(op==1)
		{
			q[++tot].pos=x,q[tot].opr=1,q[tot].val=y;
		}
		else
		{
			q[++tot].pos=x-1,q[tot].opr=2,q[tot].val=++ai;
			q[++tot].pos=y,q[tot].opr=3,q[tot].val=ai;
		}
	}
	cdq(1,tot);
	for(int i=1;i<=ai;i++)
		printf("%d\n",a[i]);
}

三维偏序

题目描述

有 n 个元素，第 i 个元素有 ai​ 、 bi​ 、 ci​ 三个属性，设 f(i) 表示满足 aj​≤ai​ 且 bj​≤bi​ 且 cj​≤ci​ 的 j 的数量。

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

输入输出格式

输入格式：

第一行两个整数 n 、k ，分别表示元素数量和最大属性值。

之后 n 行，每行三个整数ai​ 、bi​ 、ci​ ，分别表示三个属性值。

输出格式：

输出 n 行，第d+1 行表示f(i)=d 的 i 的数量。

样例输入

10 3
3 3 3
2 3 3
2 3 1
3 1 1
3 1 2
1 3 1
1 1 2
1 2 2
1 3 2
1 2 1		3
1
3
0
1
0
1
0
0
1

典型的cdq三维分治问题，直接判断每一维的数的大小累加答案即可。
但要注意，相同的元素取位置最往后的那个，以确保答案全部累加。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define N 100000+200
#define lowbit(x) x&(-x)
struct node
{
	int xi,yi,zi,pos,val;
}q[N],tmp[N],tmpp[N];
int n,mxt;
int f[N],ans[N],same[N],tree[2*N];
void add(int x,int ci)
{
	while(x<=mxt)
		tree[x]+=ci,x+=lowbit(x);
}
int query(int x)
{
	int cnt=0;
	while(x)
		cnt+=tree[x],x-=lowbit(x);
	return cnt;
}
int cmp(node a,node b)
{
	if(a.xi==b.xi)
	{
		if(a.yi==b.yi)
			return a.zi<b.zi;
		return a.yi<b.yi;
	}
	return a.xi<b.xi;
}
void cdq1(int l,int r)
{
	for(int i=l;i<=r;i++)
	{
		if(tmp[i].val==0)
			add(tmp[i].zi,1);
		if(tmp[i].val==1)
			f[tmp[i].pos]+=query(tmp[i].zi);
	}
	for(int i=l;i<=r;i++)
		if(tmp[i].val==0)
			add(tmp[i].zi,-1);
}
void cdq(int l,int r)
{
	if(l==r)
		return;
	int mid=(l+r)/2,i=l,j=mid+1,ti=l-1;
	cdq(l,mid),cdq(mid+1,r);
	while(i<=mid&&j<=r)
	{
		if(q[i].yi<=q[j].yi)
			tmp[++ti]=q[i++],tmp[ti].val=0;
		else
			tmp[++ti]=q[j++],tmp[ti].val=1;
	}
	while(i<=mid)
		tmp[++ti]=q[i++],tmp[ti].val=0;
	while(j<=r)
		tmp[++ti]=q[j++],tmp[ti].val=1;
	for(int k=l;k<=r;k++)
		q[k]=tmp[k];
	cdq1(l,r);
}
int main()
{
	
	scanf("%d %d",&n,&mxt);
	for(int i=1;i<=n;i++)
	{
		scanf("%d %d %d",&q[i].xi,&q[i].yi,&q[i].zi);
		q[i].pos=i;
	}
	sort(q+1,q+n+1,cmp);
	for(int i=1;i<=n;)
	{
		int j=i+1;
		while(j<=n&&q[i].xi==q[j].xi&&q[i].yi==q[j].yi&&q[i].zi==q[j].zi)
			j++;
		while(i<j)
            same[q[i++].pos]=q[j-1].pos;
	}
	cdq(1,n);
	for(int i=1;i<=n;i++)
		ans[f[same[q[i].pos]]+1]++;
	for(int i=1;i<=n;i++)
		printf("%d\n",ans[i]);
}

四维偏序 

jzoj4419. hole

Description

GFS打算去郊外建所别墅，享受生活，于是他耗费巨资买下了一块风水宝地，但令他震惊的是，一群DSJ对GFS的富贵生活深恶痛绝，决定打洞以搞破坏。
现在我们简化一下这个问题，在这片土地上会按顺序发生一系列事件。
①一只DSJ在(x,y) 这个点打了一个洞。
②有着高雅品味GFS想建一个等腰直角三角形的别墅，即由(x,y) ,(x+d,y) ,(x,y+d) 三点围成的三角形，但为了地基的牢固，他想知道当前这块三角形土地内的洞的个数。
GFS现在对DSJ已经忍无可忍了，请你帮他回答这些询问。
初始土地上没有洞。GFS毕竟是GFS，你可以认为土地无限大。

Input

第一行一个整数 n，表示事件数。接下来n行，每行3个非负整数x ，y ，d 。
若d=0 表示DSJ打洞的事件。否则表示GFS建房的询问。

Output

对每个询问输出一个整数，表示当时询问的三角形内的洞的个数。

Sample Input/Sample Output

8
1 3 0
1 5 0
3 6 0
4 4 0
2 6 0
1 5 3
1 5 4
1 1 1

 

3
3
0
4
1 5 0
3 7 0
2 5 6
2 3 4
 
1
0

Data Constraint

30%的数据n<=3333 。
另30% 的数据 GFS只会在DSJ打完洞后才开始询问,xi,yi<=333333 。
100%的数据 1<=n<=88888，xi,yi<=3333333 。

solution

这道题好像还有其他做法，我用的是最蠢的cdq分治，并且由于代码丑开了O2才卡过去。

首先看大佬题解的一张图。

假设蓝色区域是我们要求的地方，那么这个区域可以用夹在两条深蓝色的线的所有点减去两个平行四边形来表示。

很明显这道题的四维分别是时间t，横坐标x，纵坐标y，以及第四维(x+y)，表示每一个等腰三角形的底，也就是图上的斜线。

时间默认有序，第二维cdq，第三位cdq，第四维用树状数组维护。

#pragma GCC optimize (2)
#pragma G++ optimize (2)
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define N 400000
#define M 3333333+200
#define lowbit(x) x&(-x)
struct node
{
	long long xi,yi,di,frm,val,vall,opr;
}q[2*N],tmp[2*N],tmpp[2*N];
long long n,tot,ast,mtx;
long long tree[5*M],ans[N];
void add(long long x,long long ci)
{
	while(x<=mtx)
		tree[x]+=ci,x+=lowbit(x);
}
long long query(long long x)
{
	long long cnt=0;
	while(x)
		cnt+=tree[x],x-=lowbit(x);
	return cnt;
}
void cdq2(long long l,long long r)
{
	for(long long i=l;i<=r;i++)
	{
		if((tmpp[i].vall==0&&tmpp[i].val==0)&&tmpp[i].opr==1)
			add(tmpp[i].di,1);
		if(tmpp[i].vall==1&&tmpp[i].val==1&&tmpp[i].opr==2)
			ans[tmpp[i].frm]-=query(tmpp[i].di);
		if(tmpp[i].vall==1&&tmpp[i].val==1&&tmpp[i].opr==3)
			ans[tmpp[i].frm]+=query(tmpp[i].di);
	}
	for(long long i=l;i<=r;i++)
		if((tmpp[i].vall==0&&tmpp[i].val==0)&&tmpp[i].opr==1)
			add(tmpp[i].di,-1);
}
void cdq1(long long l,long long r)
{
	if(l==r)
		return;
	long long mid=(l+r)/2,i=l,j=mid+1,ti=l-1,cnt=0;
	cdq1(l,mid),cdq1(mid+1,r);
	while(i<=mid&&j<=r)
	{
		if(tmp[i].yi<=tmp[j].yi)
			tmpp[++ti]=tmp[i++],tmpp[ti].vall=0;
		else
			tmpp[++ti]=tmp[j++],tmpp[ti].vall=1;
	}
	while(i<=mid)
		tmpp[++ti]=tmp[i++],tmpp[ti].vall=0;
	while(j<=r)
		tmpp[++ti]=tmp[j++],tmpp[ti].vall=1;
	for(long long k=l;k<=r;k++)
		tmp[k]=tmpp[k];
	cdq2(l,r);
}
void cdq(long long l,long long r)
{
	if(l==r)
		return;
	long long mid=(l+r)/2,i=l,j=mid+1,ti=l-1,cnt=0;
	cdq(l,mid),cdq(mid+1,r);
	while(i<=mid&&j<=r)
	{
		if(q[i].xi<=q[j].xi)
			tmp[++ti]=q[i++],tmp[ti].val=0;
		else
			tmp[++ti]=q[j++],tmp[ti].val=1;
	}
	while(i<=mid)
		tmp[++ti]=q[i++],tmp[ti].val=0;
	while(j<=r)
		tmp[++ti]=q[j++],tmp[ti].val=1;
	for(long long k=l;k<=r;k++)
		q[k]=tmp[k];
	cdq1(l,r);
}
int main()
{
	freopen("test.in","r",stdin);
	freopen("test.out","w",stdout);
	scanf("%lld",&n);
	for(long long i=1;i<=n;i++)
	{
		long long x,y,d;
		scanf("%lld %lld %lld",&x,&y,&d);
		if(d==0)
		{
			q[++tot].xi=x;
			q[tot].yi=y;
			q[tot].di=x+y;
			q[tot].opr=1;
		}
		else
		{
			
			q[++tot].xi=x-1;
			q[tot].yi=x+y+d;
			q[tot].frm=++ast;
			q[tot].opr=2;
			q[tot].di=x+y+d;
			
			q[++tot].xi=x+y+d;
			q[tot].yi=y-1;
			q[tot].frm=ast;
			q[tot].opr=2;
			q[tot].di=x+y+d;
			
			q[++tot].xi=x+y+d;
			q[tot].yi=x+y+d;
			q[tot].frm=ast;
			q[tot].opr=3;
			q[tot].di=x+y+d;
			
			q[++tot].xi=x-1;
			q[tot].yi=y-1;
			q[tot].frm=ast;
			q[tot].opr=3;
			q[tot].di=x+y-2;
		}
		mtx=max(mtx,x+y+d);
	}
	cdq(1,tot);
	for(long long i=1;i<=ast;i++)
		printf("%lld\n",ans[i]);
}

 

FlowsMple
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
CDQ分治入门经典.pdf
08-23
CDQ分治入门经典.pdf
CDQ分治
zixiangzhan的博客
03-21 6378
CDQ分治 CDQ分治，又称基于时间的分治算法，常用于解决多维偏序问题。该算法可以通过增加log(n)的代价将偏序问题降掉一维，从而转化成更易解决的多维偏序问题。事实上，CDQ分治能解决的题目很多都可以用支持动态查询的高级数据结构完成，但是CDQ分治的思维难度和代码实现难度较于高级数据结构减小了很多，并且空间更小。 有些题目要求诸如：只有修改操作的属性值小于（或大于）某一询问操作的答案，该修改操作才能对询问产生影响。这里，我们视操作顺序为第一维，属性值为第二维，修改/询问的值为第三维，发现这就是三维偏序问题
CDQ分治详解
m0_57344422的博客
05-22 1508
CDQ分治
CDQ分治详解，一维、二维、三维偏序
最新发布
EQUINOX1的博客
02-26 1708
如何利用CDQ分治解决偏序问题
从《Cash》谈一类分治算法的应用
10-24
总结起来，分治算法是一种强大的解决问题的工具，尤其适用于那些可以通过递归分解来解决的问题。在实际应用中，结合其他数据结构和算法技术，如平衡树和动态规划，可以进一步提高算法的效率。这篇文章通过《Cash》...
SFDC CDQ-crx插件
04-04
**SFDC CDQ-crx插件详解** SFDC CDQ-crx是一款专为 Salesforce（SFDC）用户设计的浏览器扩展程序，主要服务于日本语环境。它提供了强大的功能，帮助用户更方便地检查和管理SFDC工具中的信息。下面将详细阐述这款...
德力西CDQ3 系列双电源自动切换开关.pdf
09-21
德力西CDQ3 系列双电源自动切换开关pdf,德力西CDQ3系列双电源自动切换开关：实现常用电源（N）与备用电源（R）之间的自动切换（也可设置为手动切换），因一路电源出现断电或异常（欠压、过压、失压）而进行电源之间...
cdq.rar_充电器_Ｃ５１
09-23
标题中的"cdq.rar_充电器_Ｃ51"揭示了这个压缩包文件与充电器设计有关，且采用了C51编程。C51是针对8051系列单片机的专用编程语言，它是标准C语言的一个扩展，增加了对微控制器硬件的直接访问功能。在这个项目中，...
cdq分治
tanjunming2020的博客
09-28 1103
### 介绍 $cdq$分治是一种分治思想的算法，一般用来求二维偏序的问题，可以将$O(n^2)$的时间复杂度优化到$O(nlogn)$ $cdq$分治的主要思想就是将要求的一段区间分为两个区间，先处理两个区间内的贡献，再处理两个区间相互的影响。因为左右区间相互独立，左区间在前，右区间在后，所以只有左区间影响右区间。只要在每次处理区间的时候处理左区间对右区间的影响，就能求出最后的答案。
「学习笔记」CDQ分治
m0_62396648的博客
07-02 409
CDQ 分治的思想最早由 IOI2008 金牌得主陈丹琦在高中时整理并总结，目前这个思想的拓展十分广泛。
【CDQ分治总结】CDQ分治 | 偏序 | N
Keyu Tian's blog
03-30 558
CDQ分治是一种通过降维来解决高维偏序问题的分治算法，可避免使用某些复杂数据结构。 算法考虑的关键在于：将两个答案计算完毕的区间合并时，前段区间将会对后段区间造成怎样的影响。 范例 例1：（二维偏序）求逆序对 求逆序对就是一个很典型的二维偏序问题。设二元组 &lt;a,b&gt;&lt;a, b&gt;<a,b> ...
cdq分治(知识整理+板子总结)
Code92007的博客
09-11 876
思路来源 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分治，解决二维偏序或三维偏序类的...
详解CDQ算法
zhuziyu_edward的博客
09-30 1756
详解CDQ算法 在这篇文章中， 我将带大家走进CDQ算法。 什么是CDQ算法 从二分说起 在解释cdq算法之前， 首先回顾算法中非常经典的分而治之。 Divide And Conquer(分而治之) 分而治之是计算机科学中非常经典的思想， 通常会把一个比较复杂的大问题， 分成两个（或多个）子问题， 分别递归分而治之去解决这些子问题， 最后把子问题的结果合并为原问题的解。 有了以上分而治之的思想后， 我们可以看一个非常经典的算法例子归并排序。 归并排序： 归并排序有以下几个步骤 把数组分成做左区间和右
CDQ分治【分治（真得头疼）
热门推荐
Coco_T的博客
12-13 2万+
前言： 说实话我对于这种没有固定板子，变化多端的算法实在是非常头疼的 但是不学不行，这也是一种很重要的思伟方式 所以趁着这几天的心情比较好（快要放大周），赶紧学一波~鸣谢：tham，stdcallCDQ分治，传说中是一个♀神犇创造的算法 在了解这种算法之前，我们有必要了解一下一种基本的思想：分治知识储备：分治 分治介绍 分而治之，将原问题不断划分成若干个子问题，直到子问题规模小到足以直接解
【总结】CDQ分治
BraketBN
04-19 3560
总的来说，CDQ分治与普通分治不一样的地方在于，CDQ分治的对象是时间。即对于一个时间段[L, R]，我们取mid = (L + R) / 2，（用数据结构题举例）分治的每层只考虑mid之前的修改对mid之后的查询的贡献，然后递归到[L,mid]，(mid,R]。 显然，CDQ分治是一种离线算法，我们需要将所有的修改/查询存下来，一起进行操作。 同时，CDQ分治还需要满足：操作之间相互独立，即
【总结】CDQ分治总结
氧化钠的博客
02-23 652
前言： CDQ分治，严格意义上说并不算一种算法，而是一种思想：将问题分为两部分，先解决左半部分，根据左半部分的信息更新右半部分。我的博客主要是借助三位偏序的模型，来介绍这种算法思想。 一维偏序： 一维偏序问题非常经典：其实就是我们常说的排序。 那么，排序有哪些方法？ 1、归并排序 2、快排 3、堆排序（借助数据结构） …… 我们常用的大约就是以上三种，如果不记得了请自行复习，这几...
离线分治算法——CDQ分治详解
zsyzClb的博客
01-17 1490
以前做过的许多题目，都是在线算法，也就是对与一系列的询问，利用某种数据结构逐个求出答案。 而今天学了一种离线算法——CDQ分治，是将全部的询问放在一起，利用分治一同处理。 CDQ分治，由2008年国际信息学奥林匹克竞赛（IOI）金牌女选手陈丹琦在国家集训队中引入而得名，为算法竞赛界中的一个广泛称呼。 CDQ分治有两种，分别是基于时间的分治和基于值域的整体分治。 下面看就开始吧。 1.基于时间的分...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值