HDU5126---stars (CDQ套CDQ套 树状数组)
题意:Q次操作,三维空间内 每个星星对应一个坐标,查询以(x1,y1,z1) (x2,y2,z2)为左下顶点 、右上顶点的立方体内的星星的个数。
注意Q的范围为50000,显然离散化之后用三维BIT会MLE。 我们可以用一次CDQ把三维变成二维,变成二维之后就有很多做法了,树套树,不会树套树的话还可以继续CDQ由二维变成一维,,变成一维了就好做了,,最基本的数据结构题目了。。
不得不说、CDQ真的很神奇。
下面做法就是CDQ套CDQ套树状数组。
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
using namespace std;
const int maxn = ;
inline int lowbit (int x)
{
return x & -x;
}
int c[maxn*],MAX;
void add(int x,int d)
{
while (x <= MAX)
{
c[x] += d;
x += lowbit(x);
}
}
int sum(int x)
{
int ans = ;
while (x)
{
ans += c[x];
x -= lowbit (x);
}
return ans;
}
struct Point
{
int x,y,z;
int kind,idx,delt;
Point() {}
Point(int _x,int _y,int _z,int _delt,int _kind,int _idx):
x(_x), y(_y), z(_z), delt(_delt), kind(_kind), idx(_idx) {} } star[maxn << ],star3[maxn << ];
int ans[maxn<<];
bool cmp1(const Point &p1,const Point &p2)
{
return p1.x < p2.x || ((p1.x == p2.x) && (p1.idx < p2.idx) );
}
bool cmp2(const Point &p1,const Point &p2)
{
return p1.y < p2.y || ((p1.y == p2.y) && (p1.idx < p2.idx) );
}
void CDQ2(int l,int r)
{
if (l >= r)
return;
int mid = (l + r) >> ;
CDQ2(l,mid);
CDQ2(mid+,r);
int j = l;
for (int i = mid + ; i <= r; i++)
{
if (star3[i].kind == )
{
for ( ; j <= mid && (star3[j].y <= star3[i].y); j++)
{
if (star3[j].kind == )
add(star3[j].z,star3[j].delt);
}
ans[star3[i].idx] += sum(star3[i].z) * star3[i].delt;
}
}
for (int i = l; i < j; i++)
{
if (star3[i].kind == )
add(star3[i].z,-star3[i].delt);
}
inplace_merge(star3+l,star3+mid+,star3+r+,cmp2);
}
void CDQ1(int l,int r)
{
if (l == r)
return;
int mid = (l + r) >> ;
CDQ1(l, mid);
CDQ1(mid+, r);
int tot = ;
for (int j = l; j <= mid ; j++)
if (star[j].kind == )
star3[tot++] = star[j];
for (int i = mid + ; i <= r; i++)
{
if (star[i].kind == )
star3[tot++] = star[i];
}
sort(star3+,star3+tot,cmp1);
CDQ2(,tot-);
}
int vec[maxn << ],idx;
void hash_(int tot)
{
sort(vec,vec+idx);
idx = unique(vec,vec+idx) - vec;
MAX = idx + ;
for (int i = ; i <= tot; i++)
star[i].z = lower_bound(vec,vec+idx,star[i].z) - vec + ;
}
int main(void)
{
#ifndef ONLINE_JUDGE
freopen("in.txt","r",stdin);
#endif // ONLINE_JUDGE
int T,Q;
scanf ("%d",&T);
while (T--)
{
scanf ("%d",&Q);
int tot = ;
int totq = ;
idx = ;
memset(c,,sizeof (c));
memset(ans,,sizeof(ans));
for (int i = ; i <= Q; i++)
{
int op,x1,y1,z1,x2,y2,z2;
scanf ("%d",&op);
if (op == )
{
scanf ("%d%d%d",&x1,&y1,&z1);
star[tot] = Point(x1,y1,z1,,,totq);
vec[idx++] = z1;
ans[totq] = -;
tot++;
totq++;
}
if (op == )
{
scanf ("%d%d%d%d%d%d",&x1,&y1,&z1,&x2,&y2,&z2);
star[tot] = Point(x1-, y1-, z1-, -, , totq), vec[idx++] = z1-, tot++;
star[tot] = Point(x2, y1-, z1-, , , totq), vec[idx++] = z1-, tot++;
star[tot] = Point(x2 , y2 , z1-, -, , totq), vec[idx++] = z1-, tot++;
star[tot] = Point(x1-, y2, z1-, , , totq), vec[idx++] = z1-, tot++;
star[tot] = Point(x1-, y2, z2 , -, , totq), vec[idx++] = z2 , tot++;
star[tot] = Point(x2 , y2, z2 , , , totq), vec[idx++] = z2 , tot++;
star[tot] = Point(x2 , y1-, z2 , -, , totq), vec[idx++] = z2 , tot++;
star[tot] = Point(x1-, y1-, z2 , , , totq), vec[idx++] = z2 , tot++;
totq++;
}
}
hash_(tot);
CDQ1(,tot-);
for (int i = ; i < totq; i++)
if (~ans[i])
printf("%d\n",ans[i]);
}
return ;
}
HDU5126---stars (CDQ套CDQ套 树状数组)的更多相关文章
- HDU 5618 Jam's problem again (cdq分治+BIT 或 树状数组套Treap)
题意:给n个点,求每一个点的满足 x y z 都小于等于它的其他点的个数. 析:三维的,第一维直接排序就好按下标来,第二维按值来,第三维用数状数组维即可. 代码如下: cdq 分治: #pragma ...
- HDU 4247 Pinball Game 3D(cdq 分治+树状数组+动态规划)
Pinball Game 3D Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) ...
- HDU 5618 Jam's problem again(三维偏序,CDQ分治,树状数组,线段树)
Jam's problem again Time Limit: 5000/2500 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Othe ...
- BZOJ2683: 简单题(cdq分治 树状数组)
Time Limit: 50 Sec Memory Limit: 128 MBSubmit: 2142 Solved: 874[Submit][Status][Discuss] Descripti ...
- 【二维偏序】【树状数组】【权值分块】【分块】poj2352 Stars
经典问题:二维偏序.给定平面中的n个点,求每个点左下方的点的个数. 因为 所有点已经以y为第一关键字,x为第二关键字排好序,所以我们按读入顺序处理,仅仅需要计算x坐标小于<=某个点的点有多少个就 ...
- 【BZOJ1146】网络管理(主席树,树状数组)
[BZOJ1146]网络管理(主席树,树状数组) 题面 BZOJ权限题,洛谷题面 题解 树上带修改主席树 貌似和\(Count\ On\ A\ Tree\)那题很相似呀 只需要套上一个树状数组来维护修 ...
- hdu 5126 stars cdq分治套cdq分治+树状数组
题目链接 给n个操作, 第一种是在x, y, z这个点+1. 第二种询问(x1, y1, z1). (x2, y2, z2)之间的总值. 用一次cdq分治可以将三维变两维, 两次的话就变成一维了, 然 ...
- hdu 5126 stars (四维偏序,离线,CDQ套CDQ套树状数组)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5126 思路:支持离线,那么我们可以用两次CDQ分治使四维降为二维,降成二维后排个序用树状数组维护下就好 ...
- [Bzoj3262]陌上花开(CDQ分治&&树状数组||树套树)
题目链接 题目就是赤裸裸的三维偏序,所以用CDQ+树状数组可以比较轻松的解决,但是还是树套树好想QAQ CDQ+树状数组 #include<bits/stdc++.h> using nam ...
随机推荐
- 项目总结——深入浅出socket网络编程
前言: 为什么会有如题的概念呢,我想对于没有主动听说过socket网络编程的人来说读到题目可能就已经蒙头了,为了很好的让大家进入场景,首先说一下一个需要用到这点东西的业务需求. 首先大家应该明确的是s ...
- android在广播接收器BroadcastReceiver里面再进行发送广播,造成当前广播接收器不断循环执行问题
最近在公司处理项目时,用到锁屏状态弹出activity进行提示,类似QQ消息弹屏提示的功能.当中用到了,假如该弹出activity已经位于锁屏界面外时,将不进行再次弹窗,而是发送广播进行通知数据更新, ...
- 7、第七节课,js逻辑运算
1.条件语句 function MyClick() { ; ; ; if (iNum1==iNum2) { iNum2=iNum3; } alert(iNum2); } 2. iNum3=iNum1& ...
- js添加onclick函数
document.getElementById('Add').setAttribute("onclick",AddNum()); 相当于不停的调用Addnum函数 应改成docum ...
- JavaScript 类、构造函数、原型
类.构造函数.原型 :本质均为函数 利用的原理是:词法作用域,调用对象及作用域链 闭包 属性查找方式 设计和new运算符一起使用的函数叫做构造函数. 构造函数的工作:初始化一个新创建的对象 ...
- [转]C++强制类型转换
dynamic_cast: 通常在基类和派生类之间转换时使用,run-time castconst_cast: 主要针对const和volatile的转换. static_cast: ...
- 带格式导出数据到Excel
导出按钮事件 protected void Button1_Click(object sender, EventArgs e) { Response.Clear(); ...
- 解决百度Ueditor编辑器表格不显示边框问题
一.主要内容 CMS使用百度Ueditor编辑器中的表格功能,在编辑模式下可以正常显示边框,而文章发布之后表格不能显示边框.本博文经过查阅相关资料,最终解决了该问题. 二.使用平台 1. dedecm ...
- linux下实现ls()函数遍历目录
转载请注明原创:http://www.cnblogs.com/StartoverX/p/4600794.html 需求:在linux下遍历目录,输出目录中各文件名. 在linux下遍历目录的相关函数有 ...
- IE10的bug?disabled button如何触发事件
最近碰到一个问题,IE10下的,貌似是个bug,求助! 问题表现为:在内部有dom元素的情况下,disabled状态的button会因为子元素的mouseover产生事件冒泡而触发,也就是说,disa ...