BZOJ1176: [Balkan2007]Mokia CDQ分治
最近很不对啊=w= 写程序全是bug啊
ans数组开小了竟然一直不知道,小数据没问题大数据拍不过,交上去RE
蛋疼半天
这个主要把每次询问拆成3个询问。
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<iostream>
#define dout printf
using namespace std; const int Maxw=+,Maxm=+,Maxq=+;
typedef long long ll;
int S,W,tot,C[Maxw],ans[Maxm+Maxq*]; void Add(int x,const int&d){
for(;<x&&x<=W;x+=x&-x)C[x]+=d;
}
int Query(int x){
int ret=;
for(;<x&&x<=W;x-=x&-x)ret+=C[x];
return ret;
} const int Hmod=1e5+;
struct node{
bool tp;//0 for modify and 1 for query
int x,y,ans,id;
node(int x,int y,int tp,int ans,int id):tp(tp),x(x),y(y),ans(ans),id(id){}
node(){ans=;}
}da[Maxq*+Maxm],q[Maxq*+Maxm]; struct prob{
int s1,s2,s3,s4,t;
void init(int x1,int y1,int x2,int y2){
++tot,da[s1=tot]=node(x2,y2,,,tot);
if(x1==||y1==)s2=;
else ++tot,da[s2=tot]=node(x1-,y1-,,,tot);
if(x1==)s3=;
else ++tot,da[s3=tot]=node(x1-,y2,,,tot);
if(y1==)s4=;
else ++tot,da[s4=tot]=node(x2,y1-,,,tot);
t=(x2-x1+)*(y2-y1+);
}
int calc(){
return ans[s1]+ans[s2]-(ans[s3]+ans[s4])+t*S;
}
}pr[Maxq];int totpr=; inline bool cmp(const node&a,const node&b){
if(a.x!=b.x)return a.x<b.x;
return a.y<b.y;
}
void init(){
scanf("%d%d",&S,&W);
for(int x1,y1,x2,y2,d,opt;~scanf("%d",&opt)&&opt!=;){
if(opt==)scanf("%d%d%d",&x1,&y1,&d),++tot,da[tot]=node(x1,y1,,d,tot);
else scanf("%d%d%d%d",&x1,&y1,&x2,&y2),pr[++totpr].init(x1,y1,x2,y2);
}
} void CDQ(int l,int r){
if(l==r)return;
int mid=(l+r)>>;
CDQ(l,mid);
CDQ(mid+,r);
int i=l;
for(int j=mid+;j<=r;j++){
for(;i<=mid&&da[i].x<=da[j].x;i++)if(!da[i].tp)Add(da[i].y,da[i].ans);
if(da[j].tp)da[j].ans+=Query(da[j].y);
}
for(i--;i>=l;i--)if(!da[i].tp)Add(da[i].y,-da[i].ans);
merge(da+l,da+mid+,da+mid+,da+r+,q,cmp);
memcpy(da+l,q,sizeof(da[])*(r-l+));
}
int main(){
freopen("in.txt","r",stdin);
freopen("out.txt","w",stdout); init();
CDQ(,tot);
for(int i=;i<=tot;i++)ans[da[i].id]=da[i].ans;
for(int i=;i<=totpr;i++)printf("%d\n",pr[i].calc()); return ;
}
BZOJ1176: [Balkan2007]Mokia CDQ分治的更多相关文章
- BZOJ 1176: [Balkan2007]Mokia( CDQ分治 + 树状数组 )
考虑cdq分治, 对于[l, r)递归[l, m), [m, r); 然后计算[l, m)的操作对[m, r)中询问的影响就可以了. 具体就是差分答案+排序+离散化然后树状数组维护.操作数为M的话时间 ...
- BZOJ 1176[Balkan2007]Mokia(CDQ分治)
1176: [Balkan2007]Mokia Time Limit: 30 Sec Memory Limit: 162 MBSubmit: 3381 Solved: 1520[Submit][S ...
- [BZOJ1176][Balkan2007]Mokia cdq+树状数组
1176: [Balkan2007]Mokia Time Limit: 30 Sec Memory Limit: 162 MBSubmit: 3134 Solved: 1395[Submit][S ...
- BZOJ 1176: [Balkan2007]Mokia [CDQ分治]
题意: 有一个n * n的棋盘,每个格子内有一个数,初始的时候全部为0.现在要求维护两种操作: 1)Add:将格子(x, y)内的数加上A. 2)Query:询问矩阵(x0, y0, x1, y1)内 ...
- BZOJ 1176 [Balkan2007]Mokia ——CDQ分治
[题目分析] 同BZOJ2683,只需要提前处理s对结果的影响即可. CDQ的思路还是很清晰的. 排序解决一维, 分治时间, 树状数组解决一维. 复杂度是两个log [代码] #include < ...
- bzoj1176: [Balkan2007]Mokia cdq
链接 bzoj 思路 cdq入门题,拆成4个矩阵,然后cdq. 代码 /************************************************************** P ...
- bzoj1176: [Balkan2007]Mokia【cdq分治】
把询问搞成4个,cdq分治. #include <bits/stdc++.h> #define rep(i, a, b) for (int i = a;i <= b; i++) #d ...
- BZOJ1176 [Balkan2007]Mokia 【CDQ分治】
题目 维护一个W*W的矩阵,初始值均为S.每次操作可以增加某格子的权值,或询问某子矩阵的总权值.修改操作数M<=160000,询问数Q<=10000,W<=2000000. 输入格式 ...
- [bzoj1176]Mokia[CDQ分治]
啃了一天论文,发现CDQ分治的原理其实很简单,大概就是这样的一类分治:将左右区间按一定规律排序后分开处理,递归到底时直接计算答案,对于一个区间,按照第二关键字split成两个区间,先处理左区间,之后因 ...
随机推荐
- boost::thread boost库线程
一.boost::thread的创建 1.线程创建方法一: boost::shared_ptr<boost::thread> writeThread_; boost::function0& ...
- javascript 学习笔记之面向对象编程(一):类的实现
~~想是一回事,做是一回事,写出来又是一回事~~一直以来,从事C++更多的是VC++多一些,从面向过程到面向对象的转变,让我对OO的编程思想有些偏爱,将一个客观存在的规律抽象出来总是让人比较兴奋,通过 ...
- thinkphp3.2 namespace及use用法
PHP 5.3中的namespace其实是个不错的东西,可以简化编程,下面介绍三类在代码中 访问namespace中类的方法 1 引用namespace和类 假设namespace的程序为name ...
- 移动web页面使用微软雅黑字体的问题
很多前端工程师在开发手机页面的时候,发现视觉设计师们喜欢用微软雅黑作为中文字体进行设计,于是写页面的时候也定义 font-family 为微软雅黑,后来发到线上后,细心的产品经理发现页面的字体不是微软 ...
- Python快速排序
快排,取一个key值,一般取第一个即可,将小于key的放到左边,大于key的放到右边,递归实现 import random def quicksort(data, low = 0, high = No ...
- Spire.XLS
又一款Excel处理神器Spire.XLS,你值得拥有(二) 前言:上篇 C#组件系列——又一款Excel处理神器Spire.XLS,你值得拥有 介绍了下组件的两个功能,说不上特色,但确实能解决我 ...
- objective-C: nonatomic retain copy assgin 等属性详解
http://my.oschina.net/u/728866/blog/90798 property,可以提供的功能有:提供成员变量的访问方法的声明.控制成员变量的访问权限.控制多线程时成员变量的访问 ...
- Hdu 3177 Crixalis's Equipment
Crixalis's Equipment Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Oth ...
- 微软HoloLens技术解谜
HoloLens 是什么? HoloLens 是微软发布的可穿戴式增强现实计算设备,它拥有这么几个关键要素: 它是增强现实产品,即 Augmented Reality(AR),AR 技术将计算机生成的 ...
- MySQL中group_concat函数,用符号连接查询分组里字段值
http://blog.csdn.net/my_yang/article/details/7186844