Time Limit: 20 Sec Memory Limit: 512 MB

Submit: 5375 Solved: 1835

[Submit][Status][Discuss]

Description

有N个位置,M个操作。操作有两种,每次操作如果是1 a b c的形式表示在第a个位置到第b个位置,每个位置加入一个数c

如果是2 a b c形式,表示询问从第a个位置到第b个位置,第C大的数是多少。

Input

第一行N,M

接下来M行,每行形如1 a b c或2 a b c

Output

输出每个询问的结果

Sample Input

2 5

1 1 2 1

1 1 2 2

2 1 1 2

2 1 1 1

2 1 2 3

Sample Output

1

2

1

HINT

【样例说明】

第一个操作 后位置 1 的数只有 1 , 位置 2 的数也只有 1 。 第二个操作 后位置 1

的数有 1 、 2 ,位置 2 的数也有 1 、 2 。 第三次询问 位置 1 到位置 1 第 2 大的数 是

1 。 第四次询问 位置 1 到位置 1 第 1 大的数是 2 。 第五次询问 位置 1 到位置 2 第 3

大的数是 1 。‍

N,M<=50000,N,M<=50000

a<=b<=N

1操作中abs(c)<=N

2操作中c<=Maxlongint

【题解】



树套树。一个权值的线段树套一个区间的线段树。

比如[4,6]这个权值范围->第一层

我们就在第二层里面获得在所有的区间[l,r]内有多少个在[4,6]范围内的数。

即权值范围为[be,en],然后求出在任意的[l,r]内有多少个权值范围在[be,en]内的数。

最后二分获取答案即可。

二分的时候枚举答案;

获取中点m=(l+r)>>1;

看看在m+1..r这个权值范围内,且区间在【L,R】的数字有多少个;

然后做相应的事情

看代码吧

#include <cstdio>

const int MAXN = 50000+100;
const int MAXM = MAXN * 16 * 16; int n, m,L,R,c;
int root[MAXN<<2] = { 0 },cnt,lc[MAXM],rc[MAXM];
long long sum[MAXM],lz[MAXM]; void modify(int &rt, int begin, int end)
{
if (!rt)
rt = ++cnt;
if (L <= begin && end <= R)
{
sum[rt] += end - begin + 1;
lz[rt] ++;
return;
}
int m = (begin + end) >> 1;
if (L <= m)
modify(lc[rt], begin, m);
if (m < R)
modify(rc[rt], m + 1, end);
sum[rt] = sum[lc[rt]] + sum[rc[rt]] + lz[rt] * (end - begin + 1);
} void up_data()
{
int l = 1, r = n, now = 1;//now代表的是l..r这个权值范围所代表的第二层线段树的根节点
while (l != r)
{
int m = (l + r) >> 1;
modify(root[now], 1, n);//虽然修改的是部分区间,但是整个大区间都会变化
//即now这个根节点所代表的权值区间,它的二层范围区间内的数字的个数会发生变化
if (c <= m)
r = m, now <<= 1;
else
l = m + 1, now = now << 1 | 1;
}
modify(root[now], 1, n);
} int min(int a, int b)
{
return a > b ? b : a;
} int max(int a, int b)
{
return a > b ? a : b;
} long long cal(int rt, int begin, int end)
{
if (!rt)
return 0;
if (L <= begin && end <= R)
return sum[rt];
long long temp1 = 0, temp2 = 0;
int m = (begin + end) >> 1;
if (L <= m)
temp1 = cal(lc[rt], begin, m);
if (m < R)
temp2 = cal(rc[rt], m + 1, end);
return temp1 + temp2 + lz[rt] * (min(end, R) - max(L, begin)+1);
} int query()
{
int l = 1, r = n, now = 1;
while (l != r)
{
int m = (l + r) >> 1;
long long temp = cal(root[now << 1|1], 1, n);//计算比m大的数有多少个.(满足在【L,R】区间)
if (c <= temp)
l = m + 1, now = now << 1 | 1;
else//c>temp
r = m, now = now << 1, c -= temp;
}
return l;
} int main()
{
//freopen("F:\\rush.txt", "r", stdin);
scanf("%d%d", &n, &m);
int opt;
for (int i = 1; i <= m; i++)
{
scanf("%d%d%d%d", &opt, &L, &R, &c);
if (opt == 1)
up_data();
else
printf("%d\n", query());
}
return 0;
}

【34.14%】【BZOJ 3110】 [Zjoi2013]K大数查询的更多相关文章

  1. BZOJ 3110: [Zjoi2013]K大数查询 [树套树]

    3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 6050  Solved: 2007[Submit][Sta ...

  2. 树套树专题——bzoj 3110: [Zjoi2013] K大数查询 &amp; 3236 [Ahoi2013] 作业 题解

    [原题1] 3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec  Memory Limit: 512 MB Submit: 978  Solved: 476 Descri ...

  3. bzoj 3110: [Zjoi2013]K大数查询 树状数组套线段树

    3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 1384  Solved: 629[Submit][Stat ...

  4. BZOJ 3110: [Zjoi2013]K大数查询( 树状数组套主席树 )

    BIT+(可持久化)权值线段树, 用到了BIT的差分技巧. 时间复杂度O(Nlog^2(N)) ---------------------------------------------------- ...

  5. BZOJ 3110([Zjoi2013]K大数查询-区间第k大[段修改,在线]-树状数组套函数式线段树)

    3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec   Memory Limit: 512 MB Submit: 418   Solved: 235 [ Submit][ ...

  6. BZOJ 3110 [Zjoi2013]K大数查询(整体二分)

    3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 11654  Solved: 3505[Submit][St ...

  7. [BZOJ 3110] [Zjoi2013] K大数查询 【树套树】

    题目链接: BZOJ - 3110 题目分析 这道题是一道树套树的典型题目,我们使用线段树套线段树,一层是区间线段树,一层是权值线段树.一般的思路是外层用区间线段树,内层用权值线段树,但是这样貌似会很 ...

  8. BZOJ 3110 [Zjoi2013]K大数查询 (CDQ分治+树状数组)

    题目描述 有N个位置,M个操作.操作有两种,每次操作如果是1 a b c的形式表示在第a个位置到第b个位置,每个位置加入一个数c如果是2 a b c形式,表示询问从第a个位置到第b个位置,第C大的数是 ...

  9. bzoj 3110 [Zjoi2013]K大数查询——线段树套线段树(标记永久化)

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3110 第一道线段树套线段树! 第一道标记永久化! 为什么为什么写了两个半小时啊…… 本想线段 ...

  10. BZOJ.3110.[ZJOI2013]K大数查询(整体二分 树状数组/线段树)

    题目链接 BZOJ 洛谷 整体二分求的是第K小(利用树状数组).求第K大可以转为求第\(n-K+1\)小,但是这样好像得求一个\(n\). 注意到所有数的绝对值\(\leq N\),将所有数的大小关系 ...

随机推荐

  1. Java学习笔记二.2

    5.运算符:变量的意义就在于可以进行运算,运算需要用运算符来表示,分为以下几种 class test { public static void main(String[] args) { //赋值运算 ...

  2. 找不到或无法载入主类 org.jivesoftware.openfire.starter.ServerStarter

    刚接触openfire的配置就出现了这个错误.解决方法非常easy,忘记了将openfire的源文件加入到user entries中了

  3. 设计模式之十二:状态模式(State)

    状态模式: 当一个对象的内部状态发生变化时同意改变它的行为. Allow an object to alter its behavior when its internal state changes ...

  4. BFS模版程序

    本文转自q=bfs&u=cnyali&t=blog">http://so.csdn.net/so/search/s.do?q=bfs&u=cnyali& ...

  5. 【Codeforces Round #447 (Div. 2) C】Marco and GCD Sequence

    [链接] 我是链接,点我呀:) [题意] 在这里输入题意 [题解] 把gcd(a[1..n])放在输入的n个数之间. [代码] /* 1.Shoud it use long long ? 2.Have ...

  6. 目标识别(object detection)中的 IoU(Intersection over Union)

    首先直观上来看 IoU 的计算公式: 由上述图示可知,IoU 的计算综合考虑了交集和并集,如何使得 IoU 最大,需要满足,更大的重叠区域,更小的不重叠的区域. 两个矩形窗格分别表示: 左上点.右下点 ...

  7. Python 实用第三方库

    1. youtube 视频下载 使用:you-get pip install you-get you-get 'https://www.youtube.com/watch?v=jNQXAC9IVRw'

  8. H5移动端IOS/Android兼容性总结,持续更新中…

    H5移动端IOS/Android兼容性总结,持续更新中… 1. IOS不识别日期 new Date("2018-07-01 08:00:00")在Android下正常显示可以直接进 ...

  9. ds1302模块的一个arduino程序

    /* * 读写DS1302 时钟芯片 * @author Yangtf * 很棒的文档 http://www.21ic.com/jichuzhishi/datasheet/DS1302/data/18 ...

  10. 硬件——nrf51822第一篇,GPIO的使用

    未完,待续...... 本实现是基于一个开发箱,包括:综合应用开发系统主板XT-EDU-AK   1套: 手持终端系统 XT-EDU-HK 1套: GPIO操作 工程: 这是一个关于流水灯的程序: 我 ...