2018SEERC Points and Rectangles (CDQ分治)
题:http://codeforces.com/gym/101964/problem/K
分析:https://blog.csdn.net/qq_43202683/article/details/98115901
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
inline int read(){
int sum=,x=;
char ch=getchar();
while(ch<''||ch>''){
if(ch=='-')
x=;
ch=getchar();
}
while(ch>=''&&ch<='')
sum=(sum<<)+(sum<<)+(ch^),ch=getchar();
return x?sum:-sum;
}
inline void write(ll x){
if(x<)
putchar('-'),x=-x;
if(x>)
write(x/);
putchar(x%+'');
}
#define pb push_back
#define lowbit(x) (x&(-x))
const int M=5e5+;
int ans[M];
int tree[M];
void add(int i,int x){
while(i<M)
tree[i]+=x,i+=(i&(-i));
return ;
}
int query(int i){
int res=;
while(i){
res+=tree[i];
i-=(i&(-i));
}
return res;
}
struct node{
int x,y,op,id;
bool operator <(const node &b)const{
return id<b.id;
}
}a[M],temp[M];
vector<pair<int, int> > book;
void cdq1(int l, int r)//算点对矩阵的贡献(小于等于)
{
if(r-l<=) return;
int midd = (r+l)>>;
cdq1(l, midd);
cdq1(midd, r);
int i= l,j= midd,n = ;
while(i< midd &&j< r){
if(a[i].x <= a[j].x){
if(a[i].op == ){
add(a[i].y, );
book.push_back(make_pair(a[i].y, ));
}
temp[n++] = a[i++];
}
else{
if(a[j].op==){
ans[a[j].id]+=query(a[j].y);
}
else if(a[j].op==)
ans[a[j].id]-=query(a[j].y);
temp[n++]=a[j++];
}
}
while(i<midd)
temp[n++] = a[i++];
while(j<r){
if(a[j].op == )
ans[a[j].id] += query(a[j].y);
else if(a[j].op == )
ans[a[j].id] -= query(a[j].y);
temp[n++]=a[j++];
}
for(i =l;i<r;++i)
a[i]=temp[i - l];
for(i = ; i <book.size(); ++i)
add(book[i].first, -book[i].second);
book.clear();
return;
}
void cdq2(int l,int r)
{
if(r-l<= )
return;
int midd=(l+r)>>;
cdq2(l,midd);
cdq2(midd,r);
int i= l,j=midd,n=;
while(i<midd&&j<r){
if(a[i].x<a[j].x){
if(a[i].op == ){
add(a[i].y,);
book.pb(make_pair(a[i].y, ));
}
else if(a[i].op==){
add(a[i].y,-);
book.pb(make_pair(a[i].y,-));
}
temp[n++]=a[i++];
}
else{
if(a[j].op == ){
ans[a[j].id] +=query(a[j].y-);
}
temp[n++]=a[j++];
}
}
while(i<midd)
temp[n++]=a[i++];//temp[o++] = e[lp++]?? = temp[lp++],debug???
while(j<r){
if(a[j].op == ) {
ans[a[j].id] +=query(a[j].y - );
}
temp[n++]=a[j++];
}
for(i=l;i<r;i++)
a[i]=temp[i-l];
for(i=;i<book.size();i++)
add(book[i].first,-book[i].second);
book.clear();
return;
}
int lisan[M<<];
int main(){
int n=read();
int tot=,num=;
for(int i=;i<n;i++){
int op=read();
if(op==){
a[tot].x=read(),a[tot].y=read();
a[tot].op=,a[tot].id=i;
lisan[++num]=a[tot].x,lisan[++num]=a[tot].y;
tot++;
}
else{
int x1=read(),y1=read(),x2=read(),y2=read();
a[tot].op=,a[tot].x=x1-,a[tot].y=y1-,a[tot].id=i,tot++;
a[tot].op=,a[tot].x=x2,a[tot].y=y2,a[tot].id=i,tot++;
a[tot].op=,a[tot].x=x1-,a[tot].y=y2,a[tot].id=i,tot++;
a[tot].op=,a[tot].x=x2,a[tot].y=y1-,a[tot].id=i,tot++;
lisan[++num]=x1-,lisan[++num]=x2,lisan[++num]=y1-,lisan[++num]=y2;
} }
sort(lisan+,lisan+num+);
int m=unique(lisan+,lisan++num)-lisan-;
for(int i=;i<tot;i++){
a[i].x=lower_bound(lisan+,lisan++m,a[i].x)-lisan;
a[i].y=lower_bound(lisan+,lisan++m,a[i].y)-lisan;
// cout<<a[i].x<<"~~~~"<<a[i].y<<endl;
}
cdq1(,tot); sort(a,a+tot);
cdq2(,tot);
ll ANS=0ll;
for(int i=;i<n;i++){
ANS+=ans[i];
write(ANS);
puts("");
}
return ;
}
2018SEERC Points and Rectangles (CDQ分治)的更多相关文章
- BZOJ 2716 Violet 3 天使玩偶 CDQ分治
题目大意:初始给定平面上的一个点集.提供两种操作: 1.将一个点增加点集 2.查询距离一个点最小的曼哈顿距离 K-D树是啥...不会写... 我仅仅会CDQ分治 对于一个询问,查询的点与这个点的位置关 ...
- 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 ...
- HDU - 6183:Color it (线段树&动态开点||CDQ分治)
Do you like painting? Little D doesn't like painting, especially messy color paintings. Now Little B ...
- CodeForces - 762E:Radio stations (CDQ分治||排序二分)
In the lattice points of the coordinate line there are n radio stations, the i-th of which is descri ...
- HDU 6183 Color it cdq分治 + 线段树 + 状态压缩
Color it Time Limit: 20000/10000 MS (Java/Others) Memory Limit: 132768/132768 K (Java/Others) Pro ...
- 【教程】简易CDQ分治教程&学习笔记
前言 辣鸡蒟蒻__stdcall终于会CDQ分治啦! CDQ分治是我们处理各类问题的重要武器.它的优势在于可以顶替复杂的高级数据结构,而且常数比较小:缺点在于必须离线操作. CDQ分治的基 ...
- BZOJ 2683 简单题 ——CDQ分治
[题目分析] 感觉CDQ分治和整体二分有着很本质的区别. 为什么还有许多人把他们放在一起,也许是因为代码很像吧. CDQ分治最重要的是加入了时间对答案的影响,x,y,t三个条件. 排序解决了x ,分治 ...
- HDU5618 & CDQ分治
Description: 三维数点 Solution: 第一道cdq分治...感觉还是很显然的虽然题目不能再傻逼了... Code: /*=============================== ...
- 初识CDQ分治
[BZOJ 1176:单点修改,查询子矩阵和]: 1176: [Balkan2007]Mokia Time Limit: 30 Sec Memory Limit: 162 MBSubmit: 200 ...
随机推荐
- 循环(while,break,continue),转义字符
01. 程序的三大流程 在程序开发中,一共有三种流程方式: 顺序 -- 从上向下,顺序执行代码 分支 -- 根据条件判断,决定执行代码的 分支 循环 -- 让 特定代码 重复 执行 02. while ...
- python-IDE使用
集成开发工具 集成开发环境(IDE,Integrated Development Environment ) VIM 经典的linux下的文本编辑器,只有小白和的牛会用这个 Emacs linux文本 ...
- 吴裕雄--天生自然MySQL学习笔记:MySQL 删除数据库
使用普通用户登陆 MySQL 服务器,可能需要特定的权限来创建或者删除 MySQL 数据库,所以使用 root 用户登录,root 用户拥有最高权限. 在删除数据库过程中,务必要十分谨慎,因为在执行删 ...
- 对PHP-GC(垃圾回收)的一点理解
一直对php的垃圾回收机制不明不白故自己开贴记录下. 首先,说到垃圾回收,得先知道什么情况下才能产生垃圾. 如果一个变量refcount在增加,说明在被使用,不是垃圾. 如果一个变量的refcount ...
- 寒假day10
今天开始写论文,同时爬取并清洗了毕设的人才动态相关部分数据
- Physicoochemical|CG content|
NCBI存在的问题: 数据用户的增长 软件开发受限 数据分析缺乏 有些传统束缚,仅用底层语言书写 Pangenome Open gene是随菌株数量增大而增大的gene,Closed gene是随菌株 ...
- awk使用笔记
awk特殊字符打印方法: 1.awk打印双引号: awk '{print "\""}' 2.awk打印单引号: awk '{print "'\''&quo ...
- 收藏基本Java项目开发的书
一.Java项目开发全程实录 第1章 进销存管理系统(Swing+SQL Server2000实现) 第2章企业内部通信系统(Swing+JavaDB实现) 第3章 企业人事管理系统( Swing+H ...
- CEO、COO、CFO、CTO、CIO是什么意思/?
CEO Chief Executive Officer 首席执行官 COO Chief Operating Officer 首席运营官 CFO Chief Financial Officer 首席财务 ...
- 干货 | 京东云应用负载均衡(ALB)多功能实操
应用负载均衡(Application Load Balancer,简称ALB)是京东云自主研发的一款七层负载均衡产品,主要面向HTTP和HTTPS流量的WEB应用程序,提供灵活的功能配置.应用负载均衡 ...