传送门

题意:

先有两种操作,插入和查询,插入操作则插入一个点\((x,y,z)\),查询操作给出两个点\((x_1,y_1,z_1),(x_2,y_2,z_2)\),回答满足\(x_1\leq x\leq x_2,y_1\leq y\leq y_2,z_1\leq z\leq z_2\)的\((x,y,z)\)的个数为多少。

思路:

就是带修改的四维偏序问题。

首先我们可以将每个点看作三维空间中的一个点,那么每次询问就相当于询问立方体中的点。

我们将询问拆成\(8\)个询问,那么每次解决的就是一个前缀和问题,也就是满足\(x\leq x_0,y\leq y_0,z\leq z_0\)的点的个数。

那么就对一维排序,之后\(cdq\)套\(cdq\)再套个树状数组就好啦。

注意我们现在只有左边区间的修改操作才会对右边区间的询问操作有影响~

注意一下空间,要开大点。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 50005; int T, q; struct node{
int x, y, z, op, id, part;
}a[N * 10], b[N * 10], d[N * 10]; int hs[N << 2]; bool isq[N];
int ans[N]; int c[N << 2];
int lowbit(int x) {return x & (-x);} void update(int x, int v) {
for(; x < N << 2; x += lowbit(x)) c[x] += v;
} int query(int x) {
int ans = 0;
for(; x; x -= lowbit(x)) ans += c[x];
return ans;
} void cdq2(int l, int r) {
if(l == r) return ;
int mid = (l + r) >> 1;
cdq2(l, mid); cdq2(mid + 1, r);
int t1 = l, t2 = mid + 1;
for(int i = l; i <= r; i++) {
if(t2 > r || (t1 <= mid && b[t1].y <= b[t2].y)) {
if(b[t1].part == 0 && b[t1].op == 0) {
update(b[t1].z, 1);
}
d[i] = b[t1++];
} else {
if(b[t2].part && b[t2].op != 0) {
ans[b[t2].id] += b[t2].op * query(b[t2].z);
}
d[i] = b[t2++];
}
}
for(int i = l; i <= mid; i++) {
if(b[i].part == 0 && b[i].op == 0) update(b[i].z, -1);
}
for(int i = l; i <= r; i++) b[i] = d[i];
} void cdq(int l, int r) {
if(l == r) return ;
int mid = (l + r) >> 1;
cdq(l, mid); cdq(mid + 1, r);
int t1 = l, t2 = mid + 1;
for(int i = l; i <= r; i++) {
if(t2 > r || (t1 <= mid && a[t1].x <= a[t2].x)) {
b[i] = a[t1++];
b[i].part = 0;
} else {
b[i] = a[t2++];
b[i].part = 1;
}
}
for(int i = l; i <= r; i++) a[i] = b[i];
cdq2(l, r);
} int main() {
// freopen("input.in", "r", stdin);
ios::sync_with_stdio(false); cin.tie(0);
cin >> T;
while(T--) {
cin >> q;
for(int i = 1; i <= q; i++) ans[i] = 0, isq[i] = false;
int cnt = 0; hs[0] = 0;
for(int i = 1; i <= q; i++) {
int op; cin >> op;
if(op == 1) {
int x, y, z; cin >> x >> y >> z;
a[++cnt] = {x, y, z, 0, i, -1};
hs[++hs[0]] = z;
} else {
isq[i] = true;
int x1, y1, z1, x2, y2, z2;
cin >> x1 >> y1 >> z1 >> x2 >> y2 >> z2;
--x1, --y1, --z1;
a[++cnt] = {x1, y1, z1, -1, i, -1};
a[++cnt] = {x1, y2, z1, 1, i, -1};
a[++cnt] = {x2, y2, z1, -1, i, -1};
a[++cnt] = {x2, y1, z1, 1, i, -1};
a[++cnt] = {x1, y1, z2, 1, i, -1};
a[++cnt] = {x1, y2, z2, -1, i, -1};
a[++cnt] = {x2, y2, z2, 1, i, -1};
a[++cnt] = {x2, y1, z2, -1, i, -1};
hs[++hs[0]] = z1, hs[++hs[0]] = z2;
}
}
sort(hs + 1, hs + hs[0] + 1);
hs[0] = unique(hs + 1, hs + hs[0] + 1) - hs - 1;
for(int i = 1; i <= cnt; i++) a[i].z = lower_bound(hs + 1, hs + hs[0] + 1, a[i].z) - hs;
cdq(1, cnt);
for(int i = 1; i <= q; i++) {
if(isq[i]) cout << ans[i] << '\n';
}
}
return 0;
}

HDU5126 stars(cdq分治)的更多相关文章

  1. HDU5126 stars(CDQ分治)

    传送门 大意: 向三维空间中加点,询问一个三维区间中点的个数. 解题思路: 带修改CDQ,将修改和询问一起插入CDQ分治询问. (询问可以由8个前缀和加减操作实现) 其中第一层CDQ维护x有序. 第二 ...

  2. hdu 5126 stars cdq分治套cdq分治+树状数组

    题目链接 给n个操作, 第一种是在x, y, z这个点+1. 第二种询问(x1, y1, z1). (x2, y2, z2)之间的总值. 用一次cdq分治可以将三维变两维, 两次的话就变成一维了, 然 ...

  3. HDU - 5126 stars (CDQ分治)

    题目链接 题目大意:一共有Q(1<=Q<=50000)组操作,操作分为两种: 1.在x,y,z处添加一颗星星 2.询问以(x1,y1,z1)与(x2,y2,z2)为左上和右下顶点的矩形之间 ...

  4. HDU5126 stars【CDQ分治】*

    HDU5126 stars Problem Description John loves to see the sky. A day has Q times. Each time John will ...

  5. [学习笔记]CDQ分治和整体二分

    序言 \(CDQ\) 分治和整体二分都是基于分治的思想,把复杂的问题拆分成许多可以简单求的解子问题.但是这两种算法必须离线处理,不能解决一些强制在线的题目.不过如果题目允许离线的话,这两种算法能把在线 ...

  6. CDQ分治笔记

    以前一直不会CDQ……然后经常听到dalao们说“这题直接CDQ啊”“CDQ不就秒了吗”的时候我只能瑟瑟发抖QAQ CDQ分治 其实CDQ分治就是二分分治,每次将$[l,r]$的问题划分为$[l,mi ...

  7. CDQ分治--用时间降维的美丽算法

    CDQ分治–用时间降维的美丽算法 CDQ分治,网上的阐述很多,太专业性的文字我就不赘述,这里指谈谈自己的感受 还是%一下CDQ大神的论文 CDQ分治的主要想法就是降维(比如三维问题降维到二维问题),并 ...

  8. CDQ分治的嵌套

    CDQ的嵌套 上一篇博客介绍了一下CDQ的入门思想.这里再介绍一下它的进阶,CDQ套CDQ.其实如果对入门思想掌握的透彻,嵌套也是很容易掌握的,思想是一样的. 什么是嵌套 简单地说,有的问题,如果用一 ...

  9. CDQ分治&整体二分学习个人小结

    目录 小结 CDQ分治 二维LIS 第一道裸题 bzoj1176 Mokia bzoj3262 陌上花开 bzoj 1790 矩形藏宝地 hdu5126四维偏序 P3157 [CQOI2011]动态逆 ...

  10. 技巧专题3(cdq分治、整体二分等)

    cdq分治与整体二分 cdq来源于2008年国家集训队作业陈丹琦(雅礼巨佬),用一个log的代价完成从静态到动态(很多时候是减少时间那一维的). 对于一个时间段[L, R],我们取mid = (L + ...

随机推荐

  1. nowcoder911J 异或的路径

    题目链接 题意 给出一棵树,每条边有边权.求\(\sum\limits_{i=1}^n{f(i,j)}\),\(f(i,j)\)表示从i到j路径的异或和. 思路 \(g_i\)表示从根到\(i\)的异 ...

  2. [LeetCode] 64. Minimum Path Sum 最小路径和

    Given a m x n grid filled with non-negative numbers, find a path from top left to bottom right which ...

  3. 初始socket模块和巧解粘包问题

    1.什么是socket? 两个进程如果需要进行通讯最基本的一个前提能够唯一的标示一个进程,在本地进程通讯中我们可以使用PID来唯一标示一个进程,但PID只在本地唯一,网络中的两个进程PID冲突几率很大 ...

  4. innodb部分内部操作

    tab_node_t* tab_create_graph_create( /*====================*/ dict_table_t* table, /*!< in: table ...

  5. SpringBoot 分环境变量配置

    方式一:通过不同环境的配置文件 1.1 springboot目录下配置不同环境的配置文件 1.2 在application.properties中配置环境 1.3运行springboot的run方法就 ...

  6. tmp/ccdLyHub.o:(.eh_frame+0x12): undefined reference to `__gxx_personality_v0' collect2: ld returned 1 exit status

    其实就是一个问题,gcc只能编译.c文件,你如果取名为.cpp,那么gcc编译就会就会出现这个错误. 这种情况下: 1.用g++编译(.c 或.c++都可以编译) 2.仍用gcc编译,但是文件后缀改为 ...

  7. pandas的使用(7)--分组

    pandas的使用(7)--分组

  8. 我的周记9——"所以快乐才是真谛"

    七月又名:巧月.瓜月.霜月.相月.凉月.初商.肇秋.初秋.首秋.早秋.新秋.上秋.孟秋.夷则.申月等. 你好七月 ,想起大学的时候我在校园的点歌台点了一首歌,歌名是July . 七月一号,新的一月,新 ...

  9. bean的shutdown

    使用@Bean注解,在不配置destroyMethod时,其默认值为: String destroyMethod() default AbstractBeanDefinition.INFER_METH ...

  10. 一个简单 System.Threading.Tasks.Dataflow.BufferBlock 示例

    直接贴代码了: using System; using System.Threading.Tasks; using System.Threading.Tasks.Dataflow; namespace ...