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

\(N,M\le50000\)

\(1\) 操作中 \(|c|\le N\)

\(2\) 操作中 \(|c|\le \text{maxlongint}\)

Solution

首先,操作一是给 \([a,b]\) 的每个位置插入一个数 \(c\),每个位置可以有多个数。(我一直以为是区间加QAQ。)

Solve(l, r, L, R) 表示第 \([L,R]\) 个询问可能的答案在区间 \([l,r]\) 中。每次把 \(c\le mid\) 的插入操作放到左边,把 \(c>mid\) 的放到右边;把答案在 \([l,mid]\) 的询问放到左边,把答案在 \([mid+1,r]\) 的询问放到右边。时间复杂度 \(O(n\log^2n)\)。

Code

#include <cstdio>

const int N = 50005;
typedef long long LL;
struct Node { int l, r, x, f; LL k; } a[N], b[N], c[N];
int ans[N], n, m, cnt; LL tag[N << 3], sum[N << 3]; int read() {
int x = 0, f = 1; char c = getchar();
while (c < '0' || c > '9') { if (c == '-') f = -1; c = getchar(); }
while (c >= '0' && c <= '9') x = (x << 3) + (x << 1) + (c ^ 48), c = getchar();
return x * f;
}
void pushdown(int cur, int l, int r) {
int mid = (l + r) >> 1;
tag[cur << 1] += tag[cur], tag[cur << 1 | 1] += tag[cur];
sum[cur << 1] += tag[cur] * (mid - l + 1), sum[cur << 1 | 1] += tag[cur] * (r - mid);
tag[cur] = 0;
}
void update(int cur, int l, int r, int L, int R, int x) {
if (L <= l && r <= R) { tag[cur] += x, sum[cur] += x * (r - l + 1); return; }
if (tag[cur]) pushdown(cur, l, r);
int mid = (l + r) >> 1;
if (L <= mid) update(cur << 1, l, mid, L, R, x);
if (mid < R) update(cur << 1 | 1, mid + 1, r, L, R, x);
sum[cur] = sum[cur << 1] + sum[cur << 1 | 1];
}
LL query(int cur, int l, int r, int L, int R) {
if (L <= l && r <= R) return sum[cur];
if (tag[cur]) pushdown(cur, l, r);
int mid = (l + r) >> 1; LL res = 0;
if (L <= mid) res = query(cur << 1, l, mid, L, R);
if (mid < R) res += query(cur << 1 | 1, mid + 1, r, L, R);
return res;
}
void solve(int l, int r, int L, int R) {
if (l > r || L > R) return;
if (l == r) {
for (int i = L; i <= R; ++i) if (a[i].f) ans[a[i].x] = n - l + 1;
return;
}
int mid = (l + r) >> 1, p = 0, q = 0;
for (int i = L; i <= R; ++i)
if (a[i].f) {
LL tmp = query(1, 1, n, a[i].l, a[i].r);
if (tmp >= a[i].k) b[++p] = a[i];
else a[i].k -= tmp, c[++q] = a[i];
} else {
if (a[i].k <= mid) update(1, 1, n, a[i].l, a[i].r, 1), b[++p] = a[i];
else c[++q] = a[i];
}
for (int i = 1; i <= p; ++i) if (!b[i].f) update(1, 1, n, b[i].l, b[i].r, -1);
for (int i = 1; i <= p; ++i) a[L + i - 1] = b[i];
for (int i = 1; i <= q; ++i) a[L + p + i - 1] = c[i];
solve(l, mid, L, L + p - 1), solve(mid + 1, r, L + p, R);
}
int main() {
n = read(), m = read();
for (int i = 1; i <= m; ++i) {
int op = read(), x = read(), y = read(); LL z; scanf("%lld", &z);
if (op == 1) a[i] = (Node){x, y, i, 0, n - z + 1};
else a[i] = (Node){x, y, ++cnt, 1, z};
}
solve(1, n, 1, m);
for (int i = 1; i <= cnt; ++i) printf("%d\n", ans[i]);
return 0;
}

[BZOJ 3110] [ZJOI 2013] K大数查询的更多相关文章

  1. BZOJ 3110 ZJOI 2013 K大数查询 树套树(权值线段树套区间线段树)

    题目大意:有一些位置.这些位置上能够放若干个数字. 如今有两种操作. 1.在区间l到r上加入一个数字x 2.求出l到r上的第k大的数字是什么 思路:这样的题一看就是树套树,关键是怎么套,怎么写.(话说 ...

  2. [BZOJ 3110] [luogu 3332] [ZJOI 2013]k大数查询(权值线段树套线段树)

    [BZOJ 3110] [luogu 3332] [ZJOI 2013]k大数查询(权值线段树套线段树) 题面 原题面有点歧义,不过从样例可以看出来真正的意思 有n个位置,每个位置可以看做一个集合. ...

  3. 数据结构(树套树):ZJOI 2013 K大数查询

    有几个点卡常数…… 发现若第一维为位置,第二维为大小,那么修改时第一维修改区间,查询时第一维查询区间,必须挂标记.而这种情况下标记很抽象,而且Push_down不是O(1)的,并不可行. 那要怎么做呢 ...

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

    http://www.lydsy.com/JudgeOnline/problem.php?id=3110 题意:-- 思路:其实和之前POJ那道题差不多,只不过是换成区间更新,而且是第k大不是第k小, ...

  5. 解题:ZJOI 2013 K大数查询

    题面 树套树,权值线段树套序列线段树,每次在在权值线段树上的每棵子树上做区间加,查询的时候左右子树二分 本来想两个都动态开点的,这样能体现树套树在线的优越性.但是常数太大惹,所以外层直接固定建树了QA ...

  6. [ZJOI 2013] K大数查询

    [题目链接] https://www.lydsy.com/JudgeOnline/problem.php?id=3110 [算法] 整体二分 + 线段树 时间复杂度 : O(NlogN ^ 2) [代 ...

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

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

  8. 【bzoj 3110】[Zjoi2013]K大数查询

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

  9. 【34.14%】【BZOJ 3110】 [Zjoi2013]K大数查询

    Time Limit: 20 Sec Memory Limit: 512 MB Submit: 5375 Solved: 1835 [Submit][Status][Discuss] Descript ...

随机推荐

  1. SDOI 2018划水记

    Day0 最后一天啦,此时不颓更待何时? 上午10:15坐车从gryz出发,在一路颓废中到了农大 不得不说,农大的宾馆真的好高档啊,壁橱里面居然有保险柜!电视柜厨子里居然有冰箱!!冰箱里居然有饮料!! ...

  2. 微信小程序发红包

    背景: 近期一个朋友公司要做活动,活动放在小程序上.小程序开发倒是不难,不过要使用小程序给微信用户发红包,这个就有点麻烦 确定模式: 小程序目前没有发红包接口,要实现的话,只能是模拟红包,即小程序上做 ...

  3. c/c++ 多线程 unique_lock的使用

    多线程 unique_lock的使用 unique_lock的特点: 1,灵活.可以在创建unique_lock的实例时,不锁,然后手动调用lock_a.lock()函数,或者std::lock(lo ...

  4. 简单用数组模拟顺序栈(c++)

    **栈是一种操作受限制的线性表,太多官方的话我也不说了,我们都知道栈元素是先进后出的,它有两种存储结构,分别是顺序存储结构和链式存储结构. **今天我先记一下顺序存储结构,后面我会加上链式存储结构的. ...

  5. Python简单的多线程demo:装逼写法

    用面向对象来写多线程: import threading class MyThread(threading.Thread): def __init__(self, n): super(MyThread ...

  6. 爬虫系列----scrapy爬取网页初始

    一 基本流程 创建工程,工程名称为(cmd):firstblood: scrapy startproject firstblood 进入工程目录中(cmd):cd :./firstblood 创建爬虫 ...

  7. Eclipse安装Gradle插件

    1.window下安装Gradle见https://www.cnblogs.com/felixzh/p/9203271.html2.eclipse中依次打开help>Install new so ...

  8. NSSM安装服务

    NSSM是一个服务封装程序,它可以将普通exe程序封装成服务,使之像windows服务一样运行.同类型的工具还有微软自己的srvany,不过nssm更加简单易用,并且功能强大.它的特点如下: 支持普通 ...

  9. linux目录1

    1.linux常用快捷键 2.linux命令之vmstat 3.linux命令之df 4.linux命令之netstat 5.linux命令之 tar 6.The authenticity of ho ...

  10. Re:Exgcd解二元不定方程

    模拟又炸了,我死亡 $exgcd$(扩展欧几里德算法)用于求$ax+by=gcd(a,b)$中$x,y$的一组解,它有很多应用,比如解二元不定方程.求逆元等等,这里详细讲解一下$exgcd$的原理. ...