传送门

文章目录

  • 题意:
  • 思路:

题意:

思路:

我们先解决怎么判断中位数的问题,我们可以二分一个midmidmid,将<mid<mid<mid的值都变成−1-1−1,其他的数都变成111,那么当全部的和>=0>=0>=0的时候,就说明当前数可以为中位数,且可能会变得更大,让后继续二分判断就好。
以上过程的复杂度是O(qnlogn)O(qnlogn)O(qnlogn)的,显然不能接受,但是可以发现我们的预处理之后的数组是可以递推出来的,不需要每次都算一遍,比如当前midmidmid为111,那么就全都是111,midmidmid为222的时候,只需要将111改成−1-1−1即可,这个过程我们可以用主席树可持久化一下,rootrootroot的下标为中位数的值的位置。
让后我们考虑怎么判断区间选择哪个,首先我们可以发现[b+1,c−1][b+1,c-1][b+1,c−1]这个区间内的数是必选的,所以直接将[b+1,c−1][b+1,c-1][b+1,c−1]区间和加上即可。让后再就是[a,b][a,b][a,b]和[c,d][c,d][c,d],对于[a,b][a,b][a,b]我们只需要取一个最长后缀,[c,d][c,d][c,d]我们只需要取一个最长前缀,让后判断sum+ls+rs>=0sum+ls+rs>=0sum+ls+rs>=0即可。
以下有两种写法。

//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid (tr[u].l+tr[u].r>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=100010,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;int n,m;
int x[N];
int root[N],tot;
vector<int>v,pos[N];
struct Query
{int val,id;bool operator < (const Query &W) const{return val<W.val;}
}a[N];
struct Node
{int l,r;int ls,rs,sum;
}tr[N*40];int find(int x)
{return lower_bound(v.begin(),v.end(),x)-v.begin()+1;
}void pushup(int u)
{tr[u].sum=tr[tr[u].l].sum+tr[tr[u].r].sum;tr[u].ls=max(tr[tr[u].l].ls,tr[tr[u].l].sum+tr[tr[u].r].ls);tr[u].rs=max(tr[tr[u].r].rs,tr[tr[u].r].sum+tr[tr[u].l].rs);
}int build(int l,int r)
{int p=++tot;if(l==r){tr[p].ls=tr[p].rs=tr[p].sum=1;return p;}int mid=l+r>>1;tr[p].l=build(l,mid);tr[p].r=build(mid+1,r);pushup(p);return p;
}void insert(int p,int &q,int l,int r,int pos)
{q=++tot; tr[q]=tr[p];if(l==r){tr[q].ls=tr[q].rs=tr[q].sum=-1;return;}int mid=l+r>>1;if(pos<=mid) insert(tr[p].l,tr[q].l,l,mid,pos);else insert(tr[p].r,tr[q].r,mid+1,r,pos);pushup(q);
}int query_sum(int u,int l,int r,int ql,int qr)
{if(!u) return 0;if(l>=ql&&r<=qr) return tr[u].sum;int ans=0,mid=l+r>>1;if(ql<=mid) ans+=query_sum(tr[u].l,l,mid,ql,qr);if(qr>mid) ans+=query_sum(tr[u].r,mid+1,r,ql,qr);return ans;
}Node query_ls(int u,int l,int r,int ql,int qr)
{if(!u) return {0,0,0,0,0};if(l>=ql&&r<=qr) return tr[u];Node ans;int mid=l+r>>1;if(qr<=mid) ans=query_ls(tr[u].l,l,mid,ql,qr);else if(ql>mid) ans=query_ls(tr[u].r,mid+1,r,ql,qr);else{Node a,b;a=query_ls(tr[u].l,l,mid,ql,qr);b=query_ls(tr[u].r,mid+1,r,ql,qr);ans.sum=a.sum+b.sum;ans.ls=max(a.ls,a.sum+b.ls);}return ans;
}Node query_rs(int u,int l,int r,int ql,int qr)
{if(!u) return {0,0,0,0,0};if(l>=ql&&r<=qr) return tr[u];Node ans;int mid=l+r>>1;if(qr<=mid) ans=query_rs(tr[u].l,l,mid,ql,qr);else if(ql>mid) ans=query_rs(tr[u].r,mid+1,r,ql,qr);else{Node a,b;a=query_rs(tr[u].l,l,mid,ql,qr);b=query_rs(tr[u].r,mid+1,r,ql,qr);ans.sum=a.sum+b.sum;ans.rs=max(b.rs,b.sum+a.rs);}return ans;
}int check(int mid)
{int add=0,ls,rs;if(x[2]+1<=x[3]-1) add=query_sum(root[mid],1,n,x[2]+1,x[3]-1);ls=query_rs(root[mid],1,n,x[1],x[2]).rs;rs=query_ls(root[mid],1,n,x[3],x[4]).ls;return (ls+rs+add)>=0;
}int main()
{//  ios::sync_with_stdio(false);
//  cin.tie(0);scanf("%d",&n);for(int i=1;i<=n;i++) scanf("%d",&a[i].val),a[i].id=i;sort(a+1,a+1+n);root[1]=build(1,n);for(int i=2;i<=n;i++) insert(root[i-1],root[i],1,n,a[i-1].id);scanf("%d",&m);int ans=0;while(m--){int aa,b,c,d; scanf("%d%d%d%d",&aa,&b,&c,&d);x[1]=(aa+ans)%n; x[2]=(b+ans)%n;x[3]=(c+ans)%n; x[4]=(d+ans)%n;x[1]++; x[2]++; x[3]++; x[4]++;sort(x+1,x+1+4);int ll=1,rr=n,anss=0;while(ll<=rr){int mid=ll+rr>>1;if(check(mid)) anss=mid,ll=mid+1;else rr=mid-1;}printf("%d\n",ans=a[anss].val);}return 0;
}
/*
5 1 1 1 1 1 1 0 1 2 3
*/
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid (tr[u].l+tr[u].r>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=100010,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;int n,m;
int a[N],x[N];
int root[N],tot;
vector<int>v,pos[N];
struct Node
{int l,r;int ls,rs,sum;
}tr[N*40];int find(int x)
{return lower_bound(v.begin(),v.end(),x)-v.begin()+1;
}void pushup(int u)
{tr[u].sum=tr[tr[u].l].sum+tr[tr[u].r].sum;tr[u].ls=max(tr[tr[u].l].ls,tr[tr[u].l].sum+tr[tr[u].r].ls);tr[u].rs=max(tr[tr[u].r].rs,tr[tr[u].r].sum+tr[tr[u].l].rs);
}int build(int l,int r)
{int p=++tot;if(l==r){tr[p].ls=tr[p].rs=tr[p].sum=1;return p;}int mid=l+r>>1;tr[p].l=build(l,mid);tr[p].r=build(mid+1,r);pushup(p);return p;
}void insert(int p,int &q,int l,int r,int pos,int x)
{q=++tot; tr[q]=tr[p];if(l==r){tr[q].ls=tr[q].rs=tr[q].sum=-1;return;}int mid=l+r>>1;if(pos<=mid) insert(tr[p].l,tr[q].l,l,mid,pos,x);else insert(tr[p].r,tr[q].r,mid+1,r,pos,x);pushup(q);
}int query_sum(int u,int l,int r,int ql,int qr)
{if(!u) return 0;if(l>=ql&&r<=qr) return tr[u].sum;int ans=0,mid=l+r>>1;if(ql<=mid) ans+=query_sum(tr[u].l,l,mid,ql,qr);if(qr>mid) ans+=query_sum(tr[u].r,mid+1,r,ql,qr);return ans;
}Node query_ls(int u,int l,int r,int ql,int qr)
{if(!u) return {0,0,0,0,0};if(l>=ql&&r<=qr) return tr[u];Node ans;int mid=l+r>>1;if(qr<=mid) ans=query_ls(tr[u].l,l,mid,ql,qr);else if(ql>mid) ans=query_ls(tr[u].r,mid+1,r,ql,qr);else{Node a,b;a=query_ls(tr[u].l,l,mid,ql,qr);b=query_ls(tr[u].r,mid+1,r,ql,qr);ans.sum=a.sum+b.sum;ans.ls=max(a.ls,a.sum+b.ls);}return ans;
}Node query_rs(int u,int l,int r,int ql,int qr)
{if(!u) return {0,0,0,0,0};if(l>=ql&&r<=qr) return tr[u];Node ans;int mid=l+r>>1;if(qr<=mid) ans=query_rs(tr[u].l,l,mid,ql,qr);else if(ql>mid) ans=query_rs(tr[u].r,mid+1,r,ql,qr);else{Node a,b;a=query_rs(tr[u].l,l,mid,ql,qr);b=query_rs(tr[u].r,mid+1,r,ql,qr);ans.sum=a.sum+b.sum;ans.rs=max(b.rs,b.sum+a.rs);}return ans;
}int check(int mid)
{int add=0,ls,rs;if(x[2]+1<=x[3]-1) add=query_sum(root[mid],1,n,x[2]+1,x[3]-1);ls=query_rs(root[mid],1,n,x[1],x[2]).rs;rs=query_ls(root[mid],1,n,x[3],x[4]).ls;return (ls+rs+add)>=0;
}int main()
{//  ios::sync_with_stdio(false);
//  cin.tie(0);scanf("%d",&n);for(int i=1;i<=n;i++) scanf("%d",&a[i]),v.pb(a[i]);sort(v.begin(),v.end()); v.erase(unique(v.begin(),v.end()),v.end());for(int i=1;i<=n;i++) pos[find(a[i])].pb(i);root[1]=build(1,n);for(int i=2;i<=(int)v.size();i++){root[i]=root[i-1];for(auto x:pos[i-1])insert(root[i],root[i],1,n,x,-1);}scanf("%d",&m);int ans=0;while(m--){int a,b,c,d; scanf("%d%d%d%d",&a,&b,&c,&d);x[1]=(a+ans)%n; x[2]=(b+ans)%n;x[3]=(c+ans)%n; x[4]=(d+ans)%n;x[1]++; x[2]++; x[3]++; x[4]++;sort(x+1,x+1+4);int ll=1,rr=(int)v.size(),anss=0;while(ll<=rr){int mid=ll+rr>>1;if(check(mid)) anss=mid,ll=mid+1;else rr=mid-1;}printf("%d\n",ans=v[anss-1]);}return 0;
}
/*
5 1 1 1 1 1 1 0 1 2 3
*/

P2839 [国家集训队]middle 二分 + 主席树 在值域上建区间相关推荐

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

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

  2. 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], ...

  3. P2839-[国家集训队]middle【主席树,二分】

    正题 题目链接:https://www.luogu.com.cn/problem/P2839 题目大意 nnn个数字,mmm次询问给出(a,b,c,d)(a,b,c,d)(a,b,c,d)表示左端点在 ...

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

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

  5. P2839 [国家集训队]middle

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

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

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

  7. YBTOJ洛谷P2839:最大中位数(主席树、二分答案)

    遇事不决,二分试试 解析 很好的一道题 真是把主席树玩明白了 一个关于中位数的常用trick: 二分答案mid,把>=mid的看成1,<mid的看成-1,然后看最大子段和是否>=0 ...

  8. 洛谷P4559 [JSOI2018]列队 【70分二分 + 主席树】

    题目链接 洛谷P4559 题解 只会做\(70\)分的\(O(nlog^2n)\) 如果本来就在区间内的人是不用动的,区间右边的人往区间最右的那些空位跑,区间左边的人往区间最左的那些空位跑 找到这些空 ...

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

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

最新文章

  1. R语言基础练习与入门实践
  2. 列举web开发中,为满足高性能的架构技术实现
  3. 【错误记录】eclipse,android,logcat日志无法打印,真机调试
  4. 无法加载 DLL“oramts.dll”: 找不到指定的模块。 (异常来自 HRESULT:0x8007007E)
  5. IIS的安装和配置全过程
  6. [转]js 取得 Unix时间戳(Unix timestamp)
  7. MAP的get与containskey
  8. C++包扩展_利用 MATLAB Coder 将M代码生成C/C++代码
  9. 火狐浏览器中文乱码怎么办 Firefox中文乱码解决方法
  10. 学生网页设计代码成品 大学生汽车网页代做 红旗汽车网页设计模板下载 HTML汽车网站制作 dreamweaver学生汽车网页设计成品
  11. 计算ndvi值需要的数据_ENVI中计算植被覆盖指数
  12. 应用迁移方案_Kubernetes如何成为迁移旧版应用程序的解决方案
  13. 这首歌,竟然是AI生成的,太狠了...
  14. .Net 内存溢出(System.OutOfMemoryException)的常见情况和处理方式总结
  15. CET UTC GMT CST 各种时区区别
  16. RWEQ模型的土壤风蚀模数估算、其变化归因分析
  17. 非线性激励函数sigmoid,tanh,softplus,Relu
  18. Unity3d培训中Rotation和EularAngles的正确使用方法
  19. gamit中SV不识别问题
  20. 感觉-勇气-缘分 爱的永恒

热门文章

  1. css实战手册第四版 pdf_你真的了解CSS继承吗?看完必跪
  2. 竟然有如何奇葩的如厕方式......
  3. 马斯克投15亿,价值47000美元的比特币,其实是一群数学家搞出来的……
  4. 清华体质优良可降5分录取;窃取密钥者奖百万;阿里投入1亿保护方言;腾讯不正当竞争被罚;这就是今天的大新闻...
  5. 趣读:程序员泪流满面的20个瞬间!
  6. 为什么离开学校后,学习能力直线下降?
  7. python葡萄酒数据集_利用python分析红葡萄酒数据
  8. 红旗河最早设计计算机的目的,论红旗河的利弊及其替代方案
  9. Android 封装handler,android封装工作线程跟Handler工具类
  10. python实战excel_实战python 读写EXCEL表