题解 洛谷 P3332
权值线段树套线段树板子题
首先观察题目,判断为二维偏序问题
操作1为区间修改,所以一定是外部线段树维护权值,内部线段树维护所在区间,否则时间复杂度爆炸qwq
为方便查找,哈希时我采用哈希每个数的相反数的方法将求第k大转换为求第k小
询问可以直接想到的做法就是二分答案,查询1~ans在区间内的个数,时间复杂度 O(nlog^3n)
尝试去掉一个log,思考发现可以直接在权值线段树上二分,每次查询左子树表示的数在区间内的个数p,若p大于等于当前查询的第k小则直接进入左子树,否则进入右子树并将k减p。时间复杂度O(nlog^2n)
内部线段树采用动态开点,空间复杂度O(nlog^2n) (n、m同阶)
但是我第一次写的树套树貌似常数过大,T了7个点,于是又加上了标记永久化,省去了标记下放的时间和不必要的空间浪费。区间修改时直接在经过的节点上修改权值并在原先打标记的节点上打上永久化标记,查询是将经过的点的标记都加起来乘区间长度再加上询问区间的权值即为答案。
具体实现见代码
#include <iostream>
#include <cstdio>
#include <algorithm>
#define gc getchar
#define re register
using namespace std;
template <typename T> void rd(T& s)
{
s = 0;
bool p = 0;
char ch;
while (ch = gc(), p |= ch == '-', ch < '0' || ch > '9');
while (s = s * 10 + ch - '0', ch = gc(), ch >= '0' && ch <= '9');
s *= (p ? -1 : 1);
}
template <typename T, typename... Args> void rd(T& s, Args&... args)
{
rd(s);
rd(args...);
}
int cnt, tot;
#define LL long long
#define new_node(ls, rs, val) (st[++cnt] = node(ls, rs, val), cnt)
const int MAXN = 50050;
const int MAX = 50000000;
struct node
{
int ls, rs;
LL val;
int tag;
node(int ls, int rs, LL val) : ls(ls), rs(rs), val(val), tag(0) {}
node() {}
}st[MAX];
struct que
{
int opt, a, b;
LL c;
que(int opt, int a, int b, LL c) : opt(opt), a(a), b(b), c(c) {}
que() {}
}qu[MAXN];
int n;
LL hs[MAXN];
int tree[MAXN << 2];
void modify(int& o, re int l, re int r, re int ll, re int rr)
{
if (!o)
o = new_node(0, 0, 0);
st[o].val += (r - l + 1);
if (l == ll && r == rr)
{
st[o].tag += 1;
return;
}
re int mid = (ll + rr) >> 1;
if (r <= mid)
modify(st[o].ls, l, r, ll, mid);
else
if (l > mid)
modify(st[o].rs, l, r, mid + 1, rr);
else
modify(st[o].ls, l, mid, ll, mid),
modify(st[o].rs, mid + 1, r, mid + 1, rr);
}
LL query(re int o, re int l, re int r, re int ll, re int rr, re int tag)
{
if (!o)
return tag * 1ll * (r - l + 1);
if (l == ll && r == rr)
return st[o].val + tag * 1ll * (r - l + 1);
re int mid = (ll + rr) >> 1;
if (r <= mid)
return query(st[o].ls, l, r, ll, mid, tag + st[o].tag);
else
if (l > mid)
return query(st[o].rs, l, r, mid + 1, rr, tag + st[o].tag);
else
return query(st[o].ls, l, mid, ll, mid, tag + st[o].tag) + query(st[o].rs, mid + 1, r, mid + 1, rr, tag + st[o].tag);
}
void add(re int o, re int k, re int l, re int r, re int ll, re int rr)
{
modify(tree[o], l, r, 1, n);
if (ll == rr)
return;
re int mid = (ll + rr) >> 1;
if (k <= mid)
add(o << 1, k, l, r, ll, mid);
else
add((o << 1) | 1, k, l, r, mid + 1, rr);
}
int midsearch(re int o, re int ll, re int rr, re int l, re int r, re int k)
{
if (ll == rr)
return ll;
re LL p = query(tree[o << 1], l, r, 1, n, 0);
if (p >= k)
return midsearch(o << 1, ll, (ll + rr) >> 1, l, r, k);
else
return midsearch((o << 1) | 1, ((ll + rr) >> 1) + 1, rr, l, r, k - p);
}
int main()
{
re int m = 0, opt, a, b;
re LL c;
rd(n, m);
for (re int i = 1; i <= m; ++i)
{
rd(opt, a, b);
rd(c);
qu[i] = que(opt, a, b, c);
if (opt == 1)
hs[++tot] = -c;
}
sort(hs + 1, hs + 1 + tot);
tot = unique(hs + 1, hs + 1 + tot) - hs - 1;
for (re int i = 1; i <= m; ++i)
{
if (qu[i].opt == 1)
{
add(1, lower_bound(hs + 1, hs + 1 + tot, -qu[i].c) - hs, qu[i].a, qu[i].b, 1, tot);
}
else
{
printf("%lld\n", -hs[midsearch(1, 1, tot, qu[i].a, qu[i].b, qu[i].c)]);
}
}
return 0;
}
题解 洛谷 P3332的更多相关文章
- 洛谷 P3332 BZOJ 3110 [ZJOI2013]K大数查询
题目链接 洛谷 bzoj 题解 整体二分 Code #include<bits/stdc++.h> #define LL long long #define RG register usi ...
- 题解 洛谷P5018【对称二叉树】(noip2018T4)
\(noip2018\) \(T4\)题解 其实呢,我是觉得这题比\(T3\)水到不知道哪里去了 毕竟我比较菜,不大会\(dp\) 好了开始讲正事 这题其实考察的其实就是选手对D(大)F(法)S(师) ...
- 题解 洛谷 P3396 【哈希冲突】(根号分治)
根号分治 前言 本题是一道讲解根号分治思想的论文题(然鹅我并没有找到论文),正 如论文中所说,根号算法--不仅是分块,根号分治利用的思想和分块像 似却又不同,某一篇洛谷日报中说过,分块算法实质上是一种 ...
- 题解-洛谷P5410 【模板】扩展 KMP(Z 函数)
题面 洛谷P5410 [模板]扩展 KMP(Z 函数) 给定两个字符串 \(a,b\),要求出两个数组:\(b\) 的 \(z\) 函数数组 \(z\).\(b\) 与 \(a\) 的每一个后缀的 L ...
- 题解-洛谷P4229 某位歌姬的故事
题面 洛谷P4229 某位歌姬的故事 \(T\) 组测试数据.有 \(n\) 个音节,每个音节 \(h_i\in[1,A]\),还有 \(m\) 个限制 \((l_i,r_i,g_i)\) 表示 \( ...
- 题解-洛谷P4724 【模板】三维凸包
洛谷P4724 [模板]三维凸包 给出空间中 \(n\) 个点 \(p_i\),求凸包表面积. 数据范围:\(1\le n\le 2000\). 这篇题解因为是世界上最逊的人写的,所以也会有求凸包体积 ...
- 题解-洛谷P4859 已经没有什么好害怕的了
洛谷P4859 已经没有什么好害怕的了 给定 \(n\) 和 \(k\),\(n\) 个糖果能量 \(a_i\) 和 \(n\) 个药片能量 \(b_i\),每个 \(a_i\) 和 \(b_i\) ...
- 题解-洛谷P5217 贫穷
洛谷P5217 贫穷 给定长度为 \(n\) 的初始文本 \(s\),有 \(m\) 个如下操作: \(\texttt{I x c}\),在第 \(x\) 个字母后面插入一个 \(c\). \(\te ...
- 题解 洛谷 P2010 【回文日期】
By:Soroak 洛谷博客 知识点:模拟+暴力枚举 思路:题目中有提到闰年然后很多人就认为,闰年是需要判断的其实,含有2月29号的回文串,前四位是一个闰年那么我们就可以直接进行暴力枚举 一些小细节: ...
随机推荐
- echarts 实战 : 想让图例颜色和元素颜色对应,怎么办?
首先,在 series 里设置颜色. (我是用js生成 echarts 需要的option对象,所以可能很难看懂) normalData.sData.forEach((item, index) =&g ...
- Java基础之NIO
NIO简介: Java NIO(New IO)是从Java 1.4版本开始引入的一个新的IO API,可以替代标准的Java IO API.NIO与原来的IO有同样的作用和目的,但是使用的方式完全不同 ...
- 如何在CSDN博客开头处加上版权声明?
1.首先在CSDN账号头像处打开"管理博客"选项. 2.然后在管理博客界面左侧找到设置下面的"博客设置"选项. 3.将博客设置里的"版权声明" ...
- Pycharm远程解释器SFTP开发和调试
转载:https://blog.csdn.net/ll641058431/article/details/53049453 使用PyCharm进行远程开发和调试 你是否经常要在Windows 7或MA ...
- Spring Data JPA根据属性名查询
https://blog.csdn.net/chengqiuming/article/details/82528961
- proj0的具体实现 #CS61B-sp18
https://github.com/Centurybbx/sp18-century/tree/master/proj0 proj0的具体实现在上面的Github中. 在proj0中我明显感受到国外大 ...
- Django学习路15_创建一个订单信息,并查询2020年\9月的信息都有哪些
在 app5.models.py 中添加一个 Order 表 class Order(models.Model): o_num = models.CharField(max_length= 16 ,u ...
- PHP getDocNamespaces() 函数
实例 返回 XML 文档的根节点中声明的命名空间: <?php$xml=<<<XML高佣联盟 www.cgewang.com<?xml version="1.0 ...
- ZROI 提高十连测 DAY3
由于我不太会写 觉得从比赛开始就冷静分析.然后看完三道题心态有点爆炸没有紧扣题目的性质. 这个心态是不可取的尽量不要有畏难心理 不要草草的写暴力. LINK:[最长01子序列](http://zhen ...
- luogu 4331 [BalticOI 2004]Sequence 数字序列
LINK:数字序列 这是一道论文题 我去看了一眼论文鸽的论文. 发现讲的还算能懂.可并堆的操作也讲的比较清晰. 对于这道题首先有一个小trick 我们给a数组全部减去其对应的下标这样我们求出来的b数组 ...