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 ...
随机推荐
- XXE--XML外部实体注入漏洞
XXE漏洞原理 XXE漏洞全称XML External Entity Injection 即xml外部实体注入漏洞,XXE漏洞发生在应用程序解析XML输入时,没有禁止外部实体的加载,导致可加载恶意外部 ...
- 新iPhone泄密12人被捕,苹果这是下狠手的节奏
一直以来,苹果在保密这件事儿上就秉持着强硬态度.还记得当年乔老爷子在的时候,苹果的保密工作在科技行业算得上是首屈一指.每款iPhone及其他新品在正式发布前,几乎不会被曝出什么消息.而这,或许也是&q ...
- Centos7 死循环登录问题
问题:用户名和密码输入正确,登录后屏幕闪一下又回到初始的登录界面.不知道具体什么原因引起的,先记录下不知道是否正确的解决方案,网上找了些相关的方案有的也实现不了,可能这个问题跟装的虚拟机的版本也有关系 ...
- Ajax请求文件下载操作失败的原因和解决办法
使用Poi做excel表格导出功能,第一个想到的就是用Ajax来发送请求,但是Ajax和后台代码都执行了,就是无法下载文件. 前台代码 function exportExl(){ var form = ...
- opencv摄像头捕获视频
1.ord()函数:它以一个字符(长度为1的字符串)作为参数,返回对应的 ASCII 数值,或者 Unicode 数值,如果所给的 Unicode 字符超出了你的 Python 定义范围,则会引发一个 ...
- Linux(CENTOS7) Nginx负载均衡简单配置
负载均衡的作用 1.转发功能 按照一定的算法[权重.轮询],将客户端请求转发到不同应用服务器上,减轻单个服务器压力,提高系统并发量. 2.故障移除 通过心跳检测的方式,判断应用服务器当前是否可以正常工 ...
- h5-应用级缓存
1.html代码及css代码 <!DOCTYPE html> <!--manifest="应用程序缓存清单文件的路径 建议文件的扩展名是appcacje,这个文件的本质就是 ...
- sudo apt-get update数字签名错误解决方法
lzb@lzb:~/projects/curl-master$ sudo apt-get update 命中: http://mirrors.aliyun.com/ubuntu xenial InRe ...
- os简介
1. 操作系统(Operation System,OS) 操作系统作为接口的示意图  没有安装操作系统的计算机,通常被称为 裸机 如果想在 裸机 上运行自己所编写的程序,就必须用机器语言书写程序 如 ...
- Xcode中SVN不能提交.a及其他文件
Xcode默认忽略的.a 文件.所以无法提交到svn服务器,但是很多第三方的库都有.a文件.所以还是必须提交到服务器. 搜索了一下解决方案: http://wpt205.blog.163.com/bl ...