题目链接:https://vjudge.net/problem/ACdream-1157

题意:

Problem Description

由3钟类型操作:
1)D L R(1 <= L <= R <= 1000000000) 增加一条线段[L,R]
2)C i (1-base) 删除第i条增加的线段,保证每条插入线段最多插入一次,且这次删除操作一定合法
3) Q L R(1 <= L <= R <= 1000000000) 查询目前存在的线段中有多少条线段完全包含[L,R]这个线段,线段X被线段Y完全包含即LY <= LX

<= RX <= RY)
给出N,接下来N行,每行是3种类型之一

Input

多组数据,每组数据N

接下来N行,每行是三种操作之一(1 <= N  <= 10^5)

Output

对于每个Q操作,输出一行,答案

Sample Input

6
D 1 100
D 3 8
D 4 10
Q 3 8
C 1
Q 3 8

Sample Output

2
1

Hint

注意,删除第i条增加的线段,不是说第i行,而是说第i次增加。

比如

D 1 10

Q 1 10

D 2 3

D 3 4

Q 5 6

D 5 6

C 2是删除D 2 3

C 4是删除D 5 6

解法:CDQ分治,将删除操作看作插入一条数量为-1的线段,查询操作看作插入一条数量为0的线段,用cnt[i]表示第i次插入的线段被之前插入的线段包含的次数,按操作顺序进行分治,每次统计[l,mid+1]中有多少元素j满足j.y>=i.y,j.x<=i.x,其中mid+1<=i<=r,这个过程可以通过对两个区间都以x为第一关键字降序排,以y为第二关键字升序排,对于[mid+1,r]中的每个i(i为数量为0的元素,即为查询),将[l,mid]中所有满足j.y>=i.y的j以j.x为下标,cnt[j]为键值插入到树状数组中,那么每次只需统计树状数组中下标小于等于i.x的元素键值之和累加到cnt[i]中即可。

#include <bits/stdc++.h>
using namespace std;
const int maxn = 2e5+10;
int n, h[maxn], tot;
struct node{
int x,y,cnt,id,ans;
bool operator<(const node &rhs) const{
if(y!=rhs.y)return y>rhs.y;
return x<rhs.x;
}
}p[maxn],q[maxn];
bool cmp(node a, node b){
return a.id<b.id;
}
struct BIT{
int b[maxn];
void init(){
memset(b, 0, sizeof(b));
}
inline int lowbit(int x){
return (x&(-x));
}
void add(int x, int v){
while(x<=tot){
b[x]+=v;
x+=lowbit(x);
}
}
int query(int x){
int ret = 0;
while(x){
ret += b[x];
x -= lowbit(x);
}
return ret;
}
}bit;
void CDQ(int l, int r){
if(l == r) return;
int mid = (l+r)/2;
CDQ(l,mid);
CDQ(mid+1,r);
sort(p+l,p+mid+1);
sort(p+mid+1,p+r+1);
int j=l;
for(int i=mid+1; i<=r; i++){
for(;j<=mid&&p[j].y>=p[i].y;j++) bit.add(p[j].x, p[j].cnt);
if(!p[i].cnt) p[i].ans+=bit.query(p[i].x);
}
for(int i=l; i<j; i++) bit.add(p[i].x, -p[i].cnt);
merge(p+l,p+mid+1,p+mid+1,p+r+1,q);
for(int i=0; i<r-l+1; i++) p[l+i]=q[i];
}
int res, l[maxn], r[maxn]; int main()
{
while(~scanf("%d", &n))
{
bit.init();
tot=0,res=1;
for(int i=1; i<=n; i++){
p[i].id = i, p[i].ans = 0;
char op[3];
scanf("%s", op);
if(op[0] == 'D'){
scanf("%d%d",&p[i].x,&p[i].y);
p[i].cnt=1;
l[res]=p[i].x,r[res++]=p[i].y;
h[++tot]=p[i].x,h[++tot]=p[i].y;
}
else if(op[0]=='Q'){
scanf("%d%d",&p[i].x,&p[i].y);
p[i].cnt=0;
h[++tot]=p[i].x,h[++tot]=p[i].y;
}
else{
int temp;
scanf("%d", &temp);
p[i].x=l[temp],p[i].y=r[temp];
p[i].cnt=-1;
}
}
sort(h+1,h+tot+1);
tot = unique(h+1,h+tot+1)-h-1;
for(int i=1; i<=n; i++){
p[i].x=lower_bound(h+1,h+tot+1,p[i].x)-h;
p[i].y=lower_bound(h+1,h+tot+1,p[i].y)-h;
}
CDQ(1, n);
sort(p+1,p+n+1,cmp);
for(int i=1; i<=n; i++){
if(!p[i].cnt) printf("%d\n", p[i].ans);
}
}
return 0;
}

ACdream 1157 Segments CDQ分治的更多相关文章

  1. 【ACdream】1157 Segments cdq分治

    Segments   Problem Description 由3钟类型操作:1)D L R(1 <= L <= R <= 1000000000) 增加一条线段[L,R]2)C i ...

  2. ACdream 1157 (cdq分治)

    题目链接 Segments Time Limit: 4000/2000MS (Java/Others)Memory Limit: 20000/10000KB (Java/Others) Problem ...

  3. ACdream 1157 Segments(CDQ分治)

    题目链接:http://acdream.info/problem?pid=1157 Problem Description 由3钟类型操作:1)D L R(1 <= L <= R < ...

  4. ACdream 1157 Segments

    Segments Time Limit: 2000ms Memory Limit: 10000KB This problem will be judged on ACdream. Original I ...

  5. ACdream1157 Segments(CDQ分治 + 线段树)

    题目这么说的: 进行如下3种类型操作:1)D L R(1 <= L <= R <= 1000000000) 增加一条线段[L,R]2)C i (1-base) 删除第i条增加的线段, ...

  6. 【教程】简易CDQ分治教程&学习笔记

    前言 辣鸡蒟蒻__stdcall终于会CDQ分治啦!       CDQ分治是我们处理各类问题的重要武器.它的优势在于可以顶替复杂的高级数据结构,而且常数比较小:缺点在于必须离线操作. CDQ分治的基 ...

  7. BZOJ 2683 简单题 ——CDQ分治

    [题目分析] 感觉CDQ分治和整体二分有着很本质的区别. 为什么还有许多人把他们放在一起,也许是因为代码很像吧. CDQ分治最重要的是加入了时间对答案的影响,x,y,t三个条件. 排序解决了x ,分治 ...

  8. HDU5618 & CDQ分治

    Description: 三维数点 Solution: 第一道cdq分治...感觉还是很显然的虽然题目不能再傻逼了... Code: /*=============================== ...

  9. 初识CDQ分治

    [BZOJ 1176:单点修改,查询子矩阵和]: 1176: [Balkan2007]Mokia Time Limit: 30 Sec  Memory Limit: 162 MBSubmit: 200 ...

随机推荐

  1. CentOS expr和let

    1.expr,用于计算变量等 用法:expr 表达式 用例1: #运算符号和参数之间要有空格分开:[es@bigdata-senior01 ~]$ expr 2 + 3 5 #乘号(*)需要用 \ , ...

  2. BZOJ5322:[JXOI2018]排序问题——题解

    https://www.lydsy.com/JudgeOnline/problem.php?id=5322 https://loj.ac/problem/2543 <-可以看数据,要没有这数据我 ...

  3. Ubuntu 16.04安装NVIDIA驱动后循环登录问题

    问题描述 最近买了两块NVIDIA Titan X Pascal显卡装到了服务器(运行Ubuntu 16.04)上.为了使用这两块GPU显卡,首先需要安装显卡驱动,安装方式为 #安装一个依赖文件,并更 ...

  4. 树莓派搭建LAMP,然后更改根目录

    参考网页: http://shumeipai.nxez.com/2013/10/13/install-and-config-lamp.html http://blog.csdn.net/zzuzadz ...

  5. *和&的使用

    给变量起一个别名: int a = 2; int &b = a; 取a的地址,实参是一个指针: void chage(int *data) { } void main() { int a = ...

  6. [C#] 类型学习笔记一:CLR中的类型,装箱和拆箱

    在学习.NET的时候,因为一些疑问,让我打算把.NET的类型篇做一个总结.总结以三篇博文的形式呈现. 这篇博文,作为三篇博文的第一篇,主要探讨了.NET Framework中的基本类型,以及这些类型一 ...

  7. Linux修改用户密码

    1. root修改自己 # passwd 2. root修改别人 # passwd oracle //修改oracle的密码

  8. 注意for循环中变量的作用域

    for e in collections: pass 在for 循环里, 最后一个对象e一直存在在上下文中.就是在循环外面,接下来对e的引用仍然有效. 这里有个问题容易被忽略,如果在循环之前已经有一个 ...

  9. Item 11 谨慎地覆盖Clone

    1.进行浅拷贝时,只是复制原始数据类型的值,则可以通过覆盖Clone方法来达到.另外,在进行浅拷贝的时候,还要注意,成员对象中不应该要有引用类型,如果有引用类型,那么,进行了浅拷贝之后,两个对象将会共 ...

  10. 【LibreOJ】#541. 「LibreOJ NOIP Round #1」七曜圣贤

    [题意]一开始车上有编号为0~a的红茶,过程中出现的红茶编号仅有[0,b),有三种操作: 1.买进编号未在车上出现过的红茶. 2.丢掉车上指定编号的红茶. 3.将最早丢出去的红茶捡回来. 每次操作后求 ...