bzoj2683简单题
#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>
#define maxn 500005
#define maxm 800005
using namespace std; int n,m,cnt,sum[maxn],pos[maxm],ans[maxm];
struct date{
int op,x,y,v,id;
}qs[maxm]; bool comp(date x,date y){
return x.x<y.x;
} int lowbit(int x){
return x&(-x);
} void add(int x,int y){
for (int i=x;i<=n;i+=lowbit(i)){
sum[i]+=y;
}
} int query(int x){
int temp=;
for (int i=x;i>;i-=lowbit(i)){
temp+=sum[i];
}
return temp;
} void cdq_solve(int l,int r){
if (l==r) return;
int mid=(l+r)/,temp=;
cdq_solve(l,mid),cdq_solve(mid+,r);
sort(qs+l,qs+mid+,comp),sort(qs+mid+,qs+r+,comp);
int i=l,j=mid+;
while (j<=r){
while (qs[i].op==&&i<=mid) i++;
while (qs[j].op==&&j<=r) j++;
if (i<=mid&&qs[i].x<=qs[j].x) add(qs[i].y,qs[i].v),i++,temp=i-;
else if (j<=r) ans[qs[j].id]+=query(qs[j].y),j++;
}
for (int t=l;t<=temp;t++) if (qs[t].op==) add(qs[t].y,-qs[t].v);
} int main(){
memset(ans,,sizeof(ans));
memset(sum,,sizeof(sum));
int op,x1,x2,y1,y2;
scanf("%d",&n),m=cnt=;
for (;;){
scanf("%d",&op);
if (op==){
qs[++m].op=op,qs[m].id=m;
scanf("%d%d%d",&qs[m].x,&qs[m].y,&qs[m].v);
}else{
if (op==){
scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
pos[++cnt]=m;
qs[++m].op=op,qs[m].x=x1-,qs[m].y=y1-,qs[m].id=m;
qs[++m].op=op,qs[m].x=x2,qs[m].y=y2,qs[m].id=m;
qs[++m].op=op,qs[m].x=x1-,qs[m].y=y2,qs[m].id=m;
qs[++m].op=op,qs[m].x=x2,qs[m].y=y1-,qs[m].id=m;
}else break;
}
}
cdq_solve(,m);
for (int i=;i<=cnt;i++) printf("%d\n",ans[pos[i]+]+ans[pos[i]+]-ans[pos[i]+]-ans[pos[i]+]);
}
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=2683
题目大意:给定一个n*n的矩形,以及若干个操作,操作有如下两种:
1.给矩形的(x,y)加上一个v;
2.询问某个子矩阵的权值和。
数据范围:n<=5*10^5,操作数<=2*10^5;
初看这题,尼玛这不是二维树状数组(二维线段树也行,不过是作死)傻逼题吗,期望复杂度O(nlog^2 n),空间复杂度n^2,再一看数据范围,GG。
在ls神犇的教导下,我学会了四种通过此题的方法:1.二维线段树动态开点,空间复杂度nlog^2 n),有点卡空间。
2.二维树状数组+hash表,可惜我没掌握这种正确姿势-.-;
3.线段树套平衡树,线段树维护x坐标,每一层套一个splay,以y坐标为关键字,对于一个询问,拆成四个矩阵,再在对应的线段树中在splay中去找,更新答案即可,据说这种写法空间复杂度只有O(n)的复杂度,可以过。
4.进入正题吧,这题的最佳姿势是cdq分治+树状数组维护区间和;
前三种做法略过(在线做法),说说第四种做法(离线做法)吧:
我先来说说cdq_分治的流程吧,cdq分治一般是解决具有一下两种性质的问题(只有修改和查询操作):
1.只能用离线做法。
2.每个修改操作对询问的影响是独立的,即与之前的修改操作无关。
cdq分治与一般的分治不同,一般的分治分出来的子区间是独立的,个个击破即可,而cdq分治分出来的两个子区间是相互联系的。(以下的分治都是指cdq分治)
由于在该问题中,每个询问只与在此之前的修改操作有关。
对于区间l,mid,记为区间1,区间mid+1,r,记为区间2,
过程:
1.递归处理区间(1),
2.递归处理区间(2),
3.统计区间1中的修改操作对区间2中的询问操作的影响。
4.算法结束。
关键在于第3步,对于这个题,我们将两个区间各自按x坐标排序,然后用两个指针分别去扫一遍,遇到区间1中的修改操作,就把纵坐标加入树状数组,然后及时统计区间2中询问的答案即可。具体细节可见代码。
cdq分治+树状数组
bzoj2683简单题的更多相关文章
- [BZOJ2683]简单题/[BZOJ1176][BalkanOI2007]Mokia
[BZOJ2683]简单题 题目大意: 一个\(n\times n(n\le5\times10^5)\)的矩阵,初始时每个格子里的数全为\(0\).\(m(m\le2\times10^5)\)次操作, ...
- Bzoj2683 简单题
Time Limit: 50 Sec Memory Limit: 128 MBSubmit: 1071 Solved: 428 Description 你有一个N*N的棋盘,每个格子内有一个整数, ...
- bzoj2683简单题 cdq分治
2683: 简单题 Time Limit: 50 Sec Memory Limit: 128 MBSubmit: 1803 Solved: 731[Submit][Status][Discuss] ...
- BZOJ2683: 简单题(cdq分治 树状数组)
Time Limit: 50 Sec Memory Limit: 128 MBSubmit: 2142 Solved: 874[Submit][Status][Discuss] Descripti ...
- BZOJ2683 简单题(CDQ分治)
传送门 之前听别人说CDQ分治不难学,今天才知道果真如此.之前一直为自己想不到CDQ的方法二很不爽,今天终于是想出来了一道了,太弱-- cdq分治主要就是把整段区间分成两半,然后用左区间的值去更新右区 ...
- Bzoj2683 简单题 [CDQ分治]
Time Limit: 50 Sec Memory Limit: 128 MBSubmit: 1071 Solved: 428 Description 你有一个N*N的棋盘,每个格子内有一个整数, ...
- 【对询问分块】【主席树】bzoj2683 简单题
对操作序列分块,每S次暴力重建主席树. 当S=sqrt(n*log(n))时,复杂度为O(m*sqrt(n*log(n))). 在线的. #include<cstdio> #include ...
- cdq分治——bzoj2683简单题
https://www.lydsy.com/JudgeOnline/problem.php?id=2683 知识点:1.以操作的顺序进行分治 2.cdq分治维护矩阵 3.计算比mid小的给比mid大 ...
- [BZOJ2683][BZOJ4066]简单题
[BZOJ2683][BZOJ4066]简单题 试题描述 你有一个N*N的棋盘,每个格子内有一个整数,初始时的时候全部为0,现在需要维护两种操作: 命令 参数限制 内容 1 x y A 1<=x ...
随机推荐
- Html5 Egret游戏开发 成语大挑战(八)一般性二级页面处理
在游戏中,我们一般会有各种各样的二级页面,比如游戏暂停界面或者游戏结束界面,这些界面组成了对玩家交互主要手段,在游戏开发中,对于这些界面的coding组织是非常有学问的,如果倒退到十年前,游戏开发的老 ...
- OSWatcher Black Box
Linux监控工具介绍系列--OSWatcher Black Box OSWatcher Balck Box简介 OSWatcher Black Box (oswbb)是Oracle开发.提供的一个小 ...
- 关于NOIP2016与NOI2018
NOIP2016惨淡收场了,距离省一还有相当一大段距离,省队更是差了十条街去了,不过没关系. 既然已经对信息学产生了兴趣,竞赛无疑是最好的锻炼场所. 路是自己选择的,伤痕累累也要走下去. 还有一年,事 ...
- XMLHTTPRequest/Ajax请求 和普通请求的区别
Ajax请求头会多一个x-requested-with参数,值为XMLHttpRequest 详情:http://blog.csdn.net/zhangdaiscott/article/details ...
- linux ISO/IMG make
sudo dd if=/PATH/*.ISO of=/dev/sdb 1.制作启动U盘需要sdb,不能sdb1,否则会提示isolinux.bin文件丢失 2.TF卡,设置sdb1?忘了 /* sy ...
- 链路层的简介和MTU
链路层杂谈(凭个人理解瞎说的,欢迎拍砖) 链路层,说白了就是把网络层的IP数据处理一下,加点东西,放到物理层上去. 加的东西:源.目的地址和CRC校验值,有的还有类型这个字段,用来区分协议. ...
- 图片加载框架Picasso解析
picasso是Square公司开源的一个Android图形缓存库 主要有以下一些特性: 在adapter中回收和取消当前的下载: 使用最少的内存完成复杂的图形转换操作: 自动的内存和硬盘缓存: 图形 ...
- TF400916错误修复办法
在使用TFS作为研发过程管理工具的时候,如果调整了工作项的状态信息,可能会出现下面的错误: 要解决此问题非常简单: 1.找一台安装了VS2015程序的环境.因为我们使用的是TFS2015,所以需要对应 ...
- 自己留存:小经验在asp.net 4.5或者asp.net mvc 5解决A potentially dangerous Request.Form value was detected from the client
以前的解决办法是 <configuration> <system.web> <pages validateRequest="false&q ...
- AMDJS编译工具
amdjs-build可以将amdjs标准的js代码编译成原生的js代码,从而在发布时去除多余的模块加载器代码. 比如像这样的代码: define("a/b", ["b& ...