嘟嘟嘟

有谁能想到这题会用到主席树呢?(不愧是WJMZBMR出的题)

首先考虑如果区间是固定的话,中位数该怎么求。
没错,二分。如果大于当前二分值\(mid\)的数比小于\(mid\)的数多,说明\(mid\)还可以再变大,向右二分;否则向左二分。
如果我们把小于\(mid\)的数都标记成\(-1\),大于的标记成\(1\),那么只用判断这个区间的和是否\(\geqslant 0\)就行了。

但现在区间不固定。首先\([b + 1, c - 1]\)是一定要选的。对于\([a, b]\)和\([c, d]\),因为要让中位数尽量大,所以应该选\([a, b]\)的最大后缀和以及\([c, d]\)的最大前缀和。

主要思路就是这些,但单次查询的复杂度是\(O(n \log{n})\)的,过不了。得想办法优化。

如果对每一个二分的值建一棵区间线段树(这里的二分在序列中的值中进行就行,而不是\(1\)到\(1e9\),所以只用建\(n\)棵),把小于他的都标记成\(1\),大于标记成\(-1\),那么每一次查询就能达到\(O(\log ^ 2{n})\)了。
但是很显然这样空间开不下,而且预处理复杂度过高。所以现在得想办法减少预处理的时间。
如果把序列中的数排一个序,会发现对于相邻的两个不一样的数(因为数字可能有重),建的线段树只有一处不一样,而这一处不一样只会导致线段树中的一条链改变。所以我们只要单独把这条链提建出来就行了。
然后就会发现这其实就是一棵主席树呀。

于是这题就写完了。

#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<cctype>
#include<vector>
#include<stack>
#include<queue>
using namespace std;
#define enter puts("")
#define space putchar(' ')
#define Mem(a, x) memset(a, x, sizeof(a))
#define rg register
typedef long long ll;
typedef double db;
const int INF = 0x3f3f3f3f;
const db eps = 1e-8;
const int maxn = 2e4 + 5;
const int maxt = 2e6 + 5;
inline ll read()
{ll ans = 0;char ch = getchar(), last = ' ';while(!isdigit(ch)) {last = ch; ch = getchar();}while(isdigit(ch)) {ans = (ans << 1) + (ans << 3) + ch - '0'; ch = getchar();}if(last == '-') ans = -ans;return ans;
}
inline void write(ll x)
{if(x < 0) x = -x, putchar('-');if(x >= 10) write(x / 10);putchar(x % 10 + '0');
}int n, _n, m, a[maxn], b[maxn], q[4];
vector<int> v[maxn];struct Tree
{int ls, rs;int sum, lmax, rmax;
}t[maxt];
int root[maxn], cnt = 0;
void pushup(int now)
{t[now].sum = t[t[now].ls].sum + t[t[now].rs].sum;t[now].lmax = max(t[t[now].ls].lmax, t[t[now].ls].sum + t[t[now].rs].lmax);t[now].rmax = max(t[t[now].rs].rmax, t[t[now].rs].sum + t[t[now].ls].rmax);
}
void build(int& now, int l, int r)
{if(!now) now = ++cnt;if(l == r) {t[now].sum = t[now].lmax = t[now].rmax = 1; return;}int mid = (l + r) >> 1;build(t[now].ls, l, mid);build(t[now].rs, mid + 1, r);pushup(now);
}
void insert(int old, int& now, int l, int r, int id)
{t[now = ++cnt] = t[old];if(l == r) {t[now].sum = t[now].lmax = t[now].rmax = -1; return;}int mid = (l + r) >> 1;if(id <= mid) insert(t[old].ls, t[now].ls, l, mid, id);else insert(t[old].rs, t[now].rs, mid + 1, r, id);pushup(now);
}
int querySum(int now, int l, int r, int L, int R)
{if(R < L) return 0;if(l == L && r == R) return t[now].sum;int mid = (l + r) >> 1;if(R <= mid) return querySum(t[now].ls, l, mid, L, R);else if(L > mid) return querySum(t[now].rs, mid + 1, r, L, R);else return querySum(t[now].ls, l, mid, L, mid) + querySum(t[now].rs, mid + 1, r, mid + 1, R);
}
int queryL(int now, int l, int r, int L, int R)
{if(l == L && r == R) return t[now].lmax;int mid = (l + r) >> 1;if(R <= mid) return queryL(t[now].ls, l, mid, L, R);else if(L > mid) return queryL(t[now].rs, mid + 1, r, L, R);else{int ret1 = queryL(t[now].ls, l, mid, L, mid);int ret2 = querySum(t[now].ls, l, mid, L, mid) + queryL(t[now].rs, mid + 1, r, mid + 1, R);return max(ret1, ret2);}
}
int queryR(int now, int l, int r, int L, int R)
{if(l == L && r == R) return t[now].rmax;int mid = (l + r) >> 1;if(R <= mid) return queryR(t[now].ls, l, mid, L, R);else if(L > mid) return queryR(t[now].rs, mid + 1, r, L, R);else{int ret1 = queryR(t[now].rs, mid + 1, r, mid + 1, R);int ret2 = querySum(t[now].rs, mid + 1, r, mid + 1, R) + queryR(t[now].ls, l, mid, L, mid);return max(ret1, ret2);}
}bool judge(int x)
{int ans1 = querySum(root[x], 1, n, q[1] + 1, q[2] - 1);int ans2 = queryR(root[x], 1, n, q[0], q[1]);int ans3 = queryL(root[x], 1, n, q[2], q[3]);return ans1 + ans2 + ans3 >= 0;
}
int solve()
{int L = 1, R = _n;while(L < R){int mid = (L + R + 1) >> 1;if(judge(mid)) L = mid;else R = mid - 1;}return L;
}int main()
{n = read();build(root[0], 1, n);for(int i = 1; i <= n; ++i) a[i] = b[i] = read();sort(b + 1, b + n + 1);_n = unique(b + 1, b + n + 1) - b - 1;for(int i = 1; i <= n; ++i){a[i] = lower_bound(b + 1, b + _n + 1, a[i]) - b;v[a[i]].push_back(i);}for(int i = 1; i <= _n; ++i){root[i] = root[i - 1];for(int j = 0; j < (int)v[i - 1].size(); ++j)insert(root[i], root[i], 1, n, v[i - 1][j]);}m = read();for(int i = 1, ans = 0; i <= m; ++i){for(int j = 0; j < 4; ++j) q[j] = read();for(int j = 0; j < 4; ++j) q[j] = (q[j] + ans) % n + 1;sort(q, q + 4);ans = b[solve()];write(ans), enter;}return 0;
}

转载于:https://www.cnblogs.com/mrclr/p/10074763.html

[国家集训队]middle相关推荐

  1. P2839 [国家集训队]middle(二分 套 主席树)

    P2839 [国家集训队]middle 有一个长度为nnn的序列,有mmm次询问,每次询问a,b,c,da, b, c, da,b,c,d,为l∈[a,b],r∈[c,d]l \in [a, b], ...

  2. 洛谷P2839 [国家集训队]middle(主席树)

    P2839 [国家集训队]middle 我们可以考虑二分中位数 checkcheckcheck 答案,那么我们对于某个值 midmidmid ,把 [l,r][l,r][l,r] 内的所有小于 mid ...

  3. bzoj 2653 洛谷 P2839 [国家集训队] middle

    2653: middle Time Limit: 20 Sec  Memory Limit: 512 MB Submit: 2381  Solved: 1340 [Submit][Status][Di ...

  4. 可持久化普通线段树 ---- P2839 [国家集训队]middle 可持久化普通线段树 + 二分 求中位数最大值

    题目链接 题目大意: 解题思路: 这个题思路很妙!! 首先我们假设只有一次询问怎么做? 那么我们可以二分出这个最大值midmidmid,然后把大于等于midmidmid设置成111,把小于midmid ...

  5. luogu2839 [国家集训队]middle

    题目链接:洛谷 题目大意:给定一个长度为$n$的序列,每次询问左端点在$[a,b]$,右端点在$[c,d]$的所有子区间的中位数的最大值.(强制在线) 这里的中位数定义为,对于一个长度为$n$的序列排 ...

  6. [bzoj 2653][国家集训队]middle

    传送门 Description 一个长度为\(n\)的序列\(a\),设其排过序之后为\(b\),其中位数定义为\(b[n/2]\),其中\(a,b\)从\(0\)开始标号,除法取下整. 给你一个长度 ...

  7. P2839 [国家集训队]middle

    题面 • 提一下静态区间第k小的nlog2n的做法: 1. 建关于排名的主席树(按排名顺序建树). 2. 二分答案. • 这样做静态区间第k小的虽然有些ZZ,但它的意义在于将线段树   维护的对象改变 ...

  8. P2839 [国家集训队]middle 二分 + 主席树 在值域上建区间

    传送门 文章目录 题意: 思路: 题意: 思路: 我们先解决怎么判断中位数的问题,我们可以二分一个midmidmid,将<mid<mid<mid的值都变成−1-1−1,其他的数都变成 ...

  9. [国家集训队]middle(二分+主席树[中位数思维题])

    文章目录 点击查看 solution code 点击查看 solution 简单口胡一下就跑 考虑二分答案ansansans 区间[x1,x2],x1∈[a,b],x2∈[c,d][x1,x2],x1 ...

最新文章

  1. 程序调试的时候利用Call Stack窗口查看函数调用信息
  2. set_origin_pose算子说明
  3. 微软2008系列 (Orcas + Longhorn Server+SQL2008)将于2008年2月27日发布
  4. java redis 命令_命令界面:使用Java中的动态API处理Redis
  5. 微软中国职位[北京工作]
  6. git 创建webpack项目_webpack项目的搭建及环境构建
  7. 千呼万唤始出来,OpenCV 4.0正式发布!
  8. LeetCode--256. 粉刷房子(动态规划)
  9. 2017.8.14 文本生成器 失败总结
  10. SpringBoot启动流程解析
  11. SQL Server数据库占用过多内存的解决方法
  12. linux远程ssh的端口映射,SSH端口映射
  13. 360安全卫士清理C盘
  14. dell3900机型安装xp64位操作系统
  15. YOLO: 3 步实时目标检测安装运行教程 [你看那条狗,好像一条狗!]
  16. Spring Boot为什么不需要额外安装Tomcat?
  17. Eclipse代码格式化无效解决方案
  18. 这种高逼格的图片效果,居然也可以用PPT制作!
  19. 【评分卡开发】信用评分模型构建流程
  20. 一个高尚的人,一个纯粹的人,一个有道德的人,一个脱离了低级趣味的人,一个有益于人民的人。...

热门文章

  1. CCNA-第十一篇-VTP+STP(上)
  2. 【HDU - 5876】Sparse Graph(补图bfs,STLset)
  3. 【POJ - 2663】Tri Tiling (简单dp)
  4. linux 命令pg,linux下操作PostgreSQL的常用命令
  5. IDEA Tomcat Catalina Log出现乱码
  6. robolectric android studio,Android Studio + Robolectric + AndroidAnnotations 根本框架
  7. 博途中用的是c吗_一只“66 鸭”离奇出走,只为助蛋宝宝C位出道
  8. python websocket服务器https_Socket与WebSocket以及http与https重新总结
  9. 关于阿里云服务器本地访问不了的问题
  10. 手把手maven的功能/安装/使用/idea集成