CDQ分治学习笔记

最新推荐文章于 2023-07-02 10:25:53 发布
最新推荐文章于 2023-07-02 10:25:53 发布
阅读量376 收藏 1
点赞数 1
文章标签： CDQ分治
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_46258139/article/details/124079307
版权

upd:整体二分也顺便学了，比较简单不写。

概述

CDQ 分治主要用于解决偏序问题。

通常是先按某一维排序，再递归处理分出来的左子问题对右子问题的答案，最后合并。

经典问题-逆序对

同时也是经典的二维偏序问题。

记属性组 ( a , b ) (a,b) (a,b)，其中 a a a 表示位置， b b b 表示值。

那么就是求 a i < a j a_i<a_j ai​<aj​ 且 b i > b j b_i>b_j bi​>bj​ 的个数。

然后就是对 a a a 排序，对 b b b 树状数组从前往后扫一遍统计即可。

时间复杂度 O ( n log ⁡ n ) \mathcal O(n \log n) O(nlogn)。

三维偏序

陌上花开

有 n n n 个元素，第 i i i 个元素有 a i , b i , c i a_i,b_i,c_i ai​,bi​,ci​ 三个属性，设 f ( i ) f(i) f(i) 表示满足 a j ≤ a i a_j \leq a_i aj​≤ai​ 且 b j ≤ b i b_j \leq b_i bj​≤bi​ 且 c j ≤ c i c_j \leq c_i cj​≤ci​ 且 j ≠ i j \ne i j​=i 的 j j j 的数量。

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

1 ≤ n ≤ 1 0 5 1 \leq n \leq 10^5 1≤n≤105， 1 ≤ a i , b i , c i ≤ k ≤ 2 × 1 0 5 1 \leq a_i,b_i,c_i \leq k \leq 2\times 10^5 1≤ai​,bi​,ci​≤k≤2×105。

先按 a a a 维排序，再将左、右子区间按 b b b 维排序。

虽然现在 a a a 的顺序被打乱了，但是前半边还是都小于后半边的，所以要是只计算前半边对后半边的偏序关系，是不会受到 a a a 的影响的。

考虑到 b b b 已经排序了，直接双指针枚举，树状数组统计 c c c 维贡献即可。

时间复杂度 O ( n log ⁡ 2 n ) \mathcal O(n \log^2 n) O(nlog2n)。

#include <bits/stdc++.h>

using namespace std;

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

const int _ = 5e5 + 10;

int n, n_, k, tot, cnt[_], c[_];

struct node
{
	int x, y, z, ans, w;
} a[_], b[_];

inline bool cmp(node u, node v)
{
	if(u.x == v.x)
	{
		if(u.y == v.y) return u.z < v.z;
		return u.y < v.y;
	}
	return u.x < v.x;
}

inline void update(int x, int val)
{
	for(int i = x; i <= k; i += i & -i) c[i] += val;
}

inline int query(int x)
{
	int res = 0;
	for(int i = x; i; i -= i & -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);
	int pl = l, pr = mid + 1, k = l;
	while(pl <= mid && pr <= r)
	{
		if(a[pl].y <= a[pr].y)
		{
			update(a[pl].z, a[pl].w);
			b[k++] = a[pl++];
		}
		else
		{
			a[pr].ans += query(a[pr].z);
			b[k++] = a[pr++];
		}
	}
	while(pl <= mid)
	{
		update(a[pl].z, a[pl].w);
		b[k++] = a[pl++];
	}
	while(pr <= r)
	{
		a[pr].ans += query(a[pr].z);
		b[k++] = a[pr++];
	}
	for(int i = l; i <= mid; ++i) update(a[i].z, -a[i].w);
	for(int i = l; i <= r; ++i) a[i] = b[i];
}

signed main()
{
	n_ = read(), k = read();
	for(int i = 1; i <= n_; ++i)
		b[i].x = read(), b[i].y = read(), b[i].z = read();
	sort(b + 1, b + n_ + 1, cmp);
	int c = 0;
	for(int i = 1; i <= n_; ++i)
	{
		c++;
		if(b[i].x != b[i + 1].x || b[i].y != b[i + 1].y || b[i].z != b[i + 1].z)
			a[++n] = b[i], a[n].w = c, c = 0;
	}
	cdq(1, n);
	for(int i = 1; i <= n; ++i)
		cnt[a[i].ans + a[i].w - 1] += a[i].w;
	for(int i = 0; i < n_; ++i)
		printf("%d\n", cnt[i]);
	return 0;
}

例题

P3157 [CQOI2011]动态逆序对

现在给出 1 ∼ n 1\sim n 1∼n 的一个排列 a a a，按照某种顺序依次删除 m m m 个元素，你的任务是在每次删除一个元素之前统计整个序列的逆序对数。

1 ≤ n ≤ 1 0 5 , 1 ≤ m ≤ 5 × 1 0 4 1 \leq n \leq 10^5,1 \leq m \leq 5 \times 10^4 1≤n≤105,1≤m≤5×104。

对于每一个被删的三元组 ( i , t i , a i ) (i,t_i,a_i) (i,ti​,ai​)（分别表示第 i i i 个数的位置，删除时间及权值），消失的逆序对 ( i , j ) (i,j) (i,j) 为：

  • 满足 j < i , a j > a i , t j > t i j<i,a_j>a_i,t_j>t_i j<i,aj​>ai​,tj​>ti​ 的 j j j。
  • 满足 j > i , a j < a i , t j > t i j > i,a_j<a_i,t_j>t_i j>i,aj​<ai​,tj​>ti​ 的 j j j。

明显的三维偏序问题，做一次 CDQ 即可。

时间复杂度 O ( n log ⁡ 2 n ) \mathcal O(n \log^2 n) O(nlog2n)。

#include <bits/stdc++.h>

using namespace std;

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

const int _ = 5e5 + 10;

int n, m, tot, cnt[_], r[_], c[_];

struct node
{
	int x, y, z, ans; // x：位置 y: 时间 z: 权值 
} a[_], b[_];

inline bool cmp(node u, node v)
{
	if(u.x == v.x)
	{
		if(u.y == v.y) return u.z < v.z;
		return u.y < v.y;
	}
	return u.x < v.x;
}

inline void update(int x, int val)
{
	for(int i = x; i <= 1e5 + 10; i += i & -i) c[i] += val;
}

inline int query(int x)
{
	int res = 0;
	for(int i = x; i; i -= i & -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);
	int ql = l, qr = mid + 1, k = l;
	while(ql <= mid && qr <= r)
	{
		if(a[ql].y > a[qr].y)
		{
			update(a[ql].z, 1);
			b[k++] = a[ql++];
		}
		else
		{
			a[qr].ans += query(1e5 + 9) - query(a[qr].z);
			b[k++] = a[qr++];
		}
	}
	while(ql <= mid)
	{
		update(a[ql].z, 1);
		b[k++] = a[ql++];
	}
	while(qr <= r)
	{
		a[qr].ans += query(1e5 + 9) - query(a[qr].z);
		b[k++] = a[qr++];
	}
	for(int i = l; i <= mid; ++i) update(a[i].z, -1);
	ql = l, qr = mid + 1, k = l;
	while(ql <= mid && qr <= r)
	{
		if(a[ql].y < a[qr].y)
		{
			update(a[qr].z, 1);
			b[k++] = a[qr++];
		}
		else
		{
			a[ql].ans += query(a[ql].z);
			b[k++] = a[ql++];
		}
	}
	while(qr <= r)
	{
		update(a[qr].z, 1);
		b[k++] = a[qr++];
	}
	while(ql <= mid)
	{
		a[ql].ans += query(a[ql].z);
		b[k++] = a[ql++];
	}
	for(int i = mid + 1; i <= r; ++i) update(a[i].z, -1);
	for(int i = l; i <= r; ++i) a[i] = b[i];
}

bool cmp2(node u, node v)
{
	return u.y < v.y;
}

signed main()
{
	n = read(), m = read();
	for(int i = 1; i <= n; ++i)
	{
		a[i].x = i;
		a[i].z = read();
		r[a[i].z] = i;
	}
	for(int i = 1, p; i <= m; ++i)
		a[r[read()]].y = i;
	for(int i = 1; i <= n; ++i)
		if(a[i].y == 0) a[i].y = m + 1;
	long long res = 0;
	for(int i = 1; i <= n; ++i)
	{
		res += query(n + 1) - query(a[i].z);
		update(a[i].z, 1);
	}
	for(int i = 1; i <= n; ++i)
		update(a[i].z, -1);
	sort(a + 1, a + n + 1, cmp);
	cdq(1, n);
	sort(a + 1, a + n + 1, cmp2);
	for(int i = 1; i <= m; ++i)
	{
		printf("%lld\n", res);
		res -= a[i].ans;
	}
	return 0;
}

习题：

P4093 [HEOI2016/TJOI2016]序列

#include <bits/stdc++.h>

using namespace std;

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

const int _ = 1e5 + 10;

int n, m, c[_], a[_], mx[_], mn[_], p[_], f[_], ans;

inline void update(int x, int val)
{
	for(int i = x; i <= 1e5; i += i & -i)
		c[i] = max(c[i], val);
}

inline void clear(int x)
{
	for(int i = x; i <= 1e5; i += i & -i)
		c[i] = 0;
}

inline int query(int x)
{
	int res = 0;
	for(int i = x; i; i -= i & -i)
		res = max(res, c[i]);
	return res;
}

inline bool cmp1(int x, int y)
{
	return mx[x] < mx[y];
}

inline bool cmp2(int x, int y)
{
	return a[x] < a[y];
}

void cdq(int l, int r)
{
	if(l == r)
	{
		f[l] = max(f[l], 1);
		return;
	}
	int mid = l + r >> 1;
	cdq(l, mid);
	for(int i = l; i <= r; ++i)
		p[i] = i;
	sort(p + l, p + mid + 1, cmp1);
	sort(p + mid + 1, p + r + 1, cmp2);
	int j = l;
	for(int i = mid + 1; i <= r; ++i)
	{
		while(j <= mid && mx[p[j]] <= a[p[i]])
		{
			update(a[p[j]], f[p[j]]);
			++j;
		}
		f[p[i]] = max(f[p[i]], query(mn[p[i]]) + 1);
	}
	for(int i = l; i <= mid; ++i)
		clear(a[i]);
	cdq(mid + 1, r);
}

signed main()
{
	n = read(), m = read();
	for(int i = 1; i <= n; ++i)
	{
		a[i] = read();
		mx[i] = mn[i] = a[i];
	}
	for(int i = 1, x, y; i <= m; ++i)
	{
		x = read(), y = read();
		mx[x] = max(mx[x], y);
		mn[x] = min(mn[x], y);
	}
	cdq(1, n);
	for(int i = 1; i <= n; ++i)
		ans = max(ans, f[i]);
	printf("%d\n", ans);
	return 0;
}

P8253 [NOI Online 2022 提高组] 如何正确地排序

#include <bits/stdc++.h>

using namespace std;

#define int long long

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

const int _ = 2e5 + 10, base = 2e5 + 1;

int n, m, a[4][_], ans, c[_ << 1], tot;

struct abc
{
	int op, x, y, val;
	inline bool operator < (const abc &t) const
	{
		if(x != t.x) return x < t.x;
		return op < t.op;
	}
} q[_ << 1];

inline void update(int x)
{
	for(int i = x; i < 2 * base; i += i & -i) c[i]++;
}

inline int query(int x)
{
	int res = 0;
	for(int i = x; i; i -= i & -i) res += c[i];
	return res;
}

int calc(int x, int y, int z)
{
	int res = 0;
	memset(c, 0, sizeof c), tot = 0;
	for(int i = 1; i <= n; ++i)
	{
		q[++tot] = {0ll, a[x][i] - a[y][i] + (x > y), a[y][i] - a[z][i] + (y > z), 0ll};
		q[++tot] = {1ll, a[y][i] - a[x][i], a[z][i] - a[y][i], a[y][i]};
	}
	sort(q + 1, q + tot + 1);
	for(int i = 1; i <= tot; ++i)
	{
		if(!q[i].op) update(q[i].y + base);
		else res += q[i].val * query(q[i].y + base);
	}
	return res;
}

signed main()
{
	m = read(), n = read();
	for(int i = 0; i < m; ++i)
		for(int j = 1; j <= n; ++j) a[i][j] = read();
	for(int i = m; i <= 3; ++i)
		for(int j = 1; j <= n; ++j) a[i][j] = a[i - m][j];
	for(int i = 0; i <= 3; ++i)
		for(int j = 1; j <= n; ++j) ans += 2 * n * a[i][j];
	for(int i = 0; i <= 3; ++i)
		for(int j = 0; j <= 3; ++j)
			for(int k = 0; k <= 3; ++k)
				if(i != j && j != k && k != i) ans -= calc(i, j, k);
	printf("%lld\n", ans);
	return 0;
}
心怀凉月
关注
  • 1
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
CDQ分治
zixiangzhan的博客
03-21 6275
CDQ分治 CDQ分治，又称基于时间的分治算法，常用于解决多维偏序问题。该算法可以通过增加log(n)的代价将偏序问题降掉一维，从而转化成更易解决的多维偏序问题。事实上，CDQ分治能解决的题目很多都可以用支持动态查询的高级数据结构完成，但是CDQ分治的思维难度和代码实现难度较于高级数据结构减小了很多，并且空间更小。 有些题目要求诸如：只有修改操作的属性值小于（或大于）某一询问操作的答案，该修改操作才能对询问产生影响。这里，我们视操作顺序为第一维，属性值为第二维，修改/询问的值为第三维，发现这就是三维偏序问题
cdq分治 学习笔记
glorious_dream的博客
10-17 664
cdq分治 学习笔记
CDQ分治入门经典.pdf
08-23
CDQ分治入门经典.pdf
从《Cash》谈一类分治算法的应用
10-24
《从《Cash》谈一类分治算法的应用》这篇文章主要探讨了如何运用分治算法来解决一个与货币兑换相关的动态规划问题。在这个问题中，小Y需要计算在已知未来N天金券价值的情况下，从一定初始资金S元出发，经过一系列...
SFDC CDQ-crx插件
04-04
**SFDC CDQ-crx插件详解** SFDC CDQ-crx是一款专为 Salesforce（SFDC）用户设计的浏览器扩展程序，主要服务于日本语环境。它提供了强大的功能，帮助用户更方便地检查和管理SFDC工具中的信息。下面将详细阐述这款...
cdq：RubyMotion的核心数据查询
02-05
核心数据查询（CDQ）是一个库，可帮助您在使用RubyMotion时管理您的核心数据堆栈。 它使用一个数据模型文件，您可以在XCode中生成该文件，也可以使用 。 CDQ由位于俄勒冈州波特兰市和加利福尼亚州旧金山市的网络和...
CDQ分治详解
m0_57344422的博客
05-22 1495
CDQ分治
CBW、CWDE、CDQ、CWD: 扩展指令
whalefall
04-08 3583
 CDQ 是一个让很多人感到困惑的指令。  这个指令把 EAX 的第 31 bit 复制到 EDX 的每一个 bit 上。 它大多出现在除法运算之前。它实际的作用只是把EDX的所有位都设成EAX最高位的值。也就是说，当EAX &lt;80000000, EDX 为00000000；当EAX &gt;= 80000000， EDX 则为FFFFFFFF。例如 :         假设 EAX 是 F...
离线分治算法——CDQ分治详解
zsyzClb的博客
01-17 1443
以前做过的许多题目，都是在线算法，也就是对与一系列的询问，利用某种数据结构逐个求出答案。 而今天学了一种离线算法——CDQ分治，是将全部的询问放在一起，利用分治一同处理。 CDQ分治，由2008年国际信息学奥林匹克竞赛（IOI）金牌女选手陈丹琦在国家集训队中引入而得名，为算法竞赛界中的一个广泛称呼。 CDQ分治有两种，分别是基于时间的分治和基于值域的整体分治。 下面看就开始吧。 1.基于时间的分...
「学习笔记」CDQ分治
m0_62396648的博客
07-02 400
CDQ 分治的思想最早由 IOI2008 金牌得主陈丹琦在高中时整理并总结，目前这个思想的拓展十分广泛。
CDQ 总结
anchengao1503的博客
07-27 310
tip：　　CDQ 分治主要处理三维偏序问题，解题时主要是找出比较量（三个或两个），并找出适当排序顺序（有时不需））。　　 　　CDQ 分治常常与树状数组搭配，树状数组主要用来统计前缀和（权值前缀和 / 排名）、最值、逆序对。 实战： T1：陌上花开 题干： 　　有 n 朵花,每朵花有三个属性:花形 (s)、颜色 (c)、气味 (m),用三个整数表示。现在要对每朵花评级,一朵花...
汇编指令---CDQ
热门推荐
oBuYiSeng的博客
12-18 1万+
CDQ 是一个让很多人感到困惑的指令。  这个指令把 EAX 的第 31 bit 复制到 EDX 的每一个 bit 上。 它大多出现在除法运算之前。它实际的作用只是把EDX的所有位都设成EAX最高位的值。也就是说，当EAX = 80000000， EDX 则为FFFFFFFF。 例如 :          假设 EAX 是 FFFFFFFB (-5) ，它的第 31 bit (最左边) 是
CDQ分治（二维CDQ 、三维CDQ+树状数组、四维CDQ+CDQ+树状数组）
jkchen's Haven
12-04 1047
CDQ分治 CDQ分治相较于普通分治，多了左区间处理后对于右区间的影响。 利用这一点，CDQ分治可以用来做很多数据结构的题目（树状数组、线段树），加一个log的时间复杂度来优化一维 &amp;amp;amp;amp;amp;ThickSpace;\\\; &amp;amp;amp;amp;amp;ThickSpace;\\\; 二维分治：最长递增子序列的数量 题意： n（5e4）的数组，有多少个LIS（最长递增子序列） 解析： 在CDQ左区间后，我们需要做的...
UVA11990　``Dynamic'' Inversion (树状数组套平衡树)
u013268685的专栏(有一种幸福叫AC，有一种期待叫AK)
08-24 1350
``Dynamic'' Inversion Time Limit: 3000ms Memory Limit: 131072KB This problem will be judged on UVA. Original ID: 11990 64-bit integer IO format: %lld      Java class name: Main Prev Sub
任务安排 cdq分治
最新发布
12-23
任务安排CDQ分治是一种常用的算法设计技巧，主要用于解决一些复杂的问题。CDQ分治算法通常用于解决一些与区间相关的问题，如区间覆盖、区间查询等。 CDQ分治算法的思想是将原始问题分解成若干个子问题，然后分别...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值