遇事不决，二分试试

解析

很好的一道题
真是把主席树玩明白了
一个关于中位数的常用trick：

二分答案mid，把>=mid的看成1，<mid的看成-1，然后看最大子段和是否>=0

然而如果对离散化后的每一个值建一棵小白逛公园那样的线段树显然时空双炸
这里主席树就用上了
注意到关于相邻的值域所建的两棵主席树只有个别值从1变成-1
所以就用主席树优化时空即可
时间复杂度mlogn2mlogn^2mlogn2，空间复杂度nlognnlognnlogn

代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N=20050;
inline ll read(){ll x=0,f=1;char c=getchar();while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}while(isdigit(c)){x=x*10+c-'0';c=getchar();}return x*f;
}
int n,m;
int tot,r[N];
struct node{int sum,l,r;
};
struct tree{int ls,rs;node o;
}tr[N*20];
inline node merge(node a,node b){node res;res.sum=a.sum+b.sum;res.l=max(a.l,a.sum+b.l);res.r=max(b.r,b.sum+a.r);return res;
}
#define pushup(k) tr[k].o=merge(tr[tr[k].ls].o,tr[tr[k].rs].o)
inline int copy(int x){tr[++tot]=tr[x];return tot;
}
#define mid ((l+r)>>1)
void build(int &k,int l,int r){k=copy(k);if(l==r){tr[k].o=(node){1,1,1};return;}build(tr[k].ls,l,mid);build(tr[k].rs,mid+1,r);pushup(k);//printf("k=%d (%d %d) sum=%d l=%d r=%d\n",k,l,r,tr[k].o.sum,tr[k].o.l,tr[k].o.r);
}
void upd(int &k,int l,int r,int p,int v){k=copy(k);if(l==r){tr[k].o=(node){v,max(v,0),max(v,0)};return;}if(p<=mid) upd(tr[k].ls,l,mid,p,v);else upd(tr[k].rs,mid+1,r,p,v);pushup(k);return;
}
node ask(int k,int l,int r,int x,int y){//printf("  k=%d (%d %d) x=%d y=%d sum=%d\n",k,l,r,x,y,tr[k].o.sum);if(x<=l&&r<=y) return tr[k].o;else if(y<=mid) return ask(tr[k].ls,l,mid,x,y);else if(x>mid) return ask(tr[k].rs,mid+1,r,x,y);else return merge(ask(tr[k].ls,l,mid,x,y),ask(tr[k].rs,mid+1,r,x,y));
}
int a[N],q[N],num;
vector <int>v[N];
int w[5];
bool check(int k,int a,int b,int c,int d){//printf("check:k=%d x=%d (%d %d) (%d %d)\n",k,q[k],a,b,c,d);node o=ask(r[k],1,n,b,c),u=ask(r[k],1,n,a,b-1),v=ask(r[k],1,n,c+1,d);//printf("res=%d+%d+%d\n\n",o.sum,u.r,v.l);return o.sum+u.r+v.l>=0;
}
int main(){n=read();for(int i=1;i<=n;i++) a[i]=q[i]=read();sort(q+1,q+1+n);num=unique(q+1,q+1+n)-q-1;for(int i=1;i<=n;i++){a[i]=lower_bound(q+1,q+1+num,a[i])-q;v[a[i]].push_back(i);}build(r[1],1,n);for(int i=1;i<num;i++){r[i+1]=r[i];for(int j=0;j<v[i].size();j++){int pl=v[i][j];upd(r[i+1],1,n,pl,-1);}}m=read();int lst=0;for(int i=1;i<=m;i++){for(int j=1;j<=4;j++) w[j]=(read()+lst)%n;sort(w+1,w+1+4);int a=w[1],b=w[2],c=w[3],d=w[4];a++;b++;c++;d++;int st=1,ed=num;while(st<ed){int o=(st+ed+1)>>1;if(check(o,a,b,c,d)) st=o;else ed=o-1;}printf("%d\n",lst=q[st]);}return 0;
}

YBTOJ洛谷P2839：最大中位数（主席树、二分答案）相关推荐

洛谷P1182 数列分段Section II 二分答案
洛谷P1182 数列分段Section II 二分答案题意:将 n 个数分为 m段求一种方案,使这m段中最大的和最小额..可能有点拗口,其实就是说每一种方案,都有对应的每段和的最大值, ...
洛谷 P3302 [SDOI2013]森林主席树+启发式合并
click here~:https://www.luogu.org/problem/P3302 emmm这个题是真的烦一看题发现这题不是count on a tree的升级版么如果一点思路没有的话 ...
BZOJ5343[Ctsc2018]混合果汁——主席树+二分答案
题目链接: CTSC2018混合果汁显然如果美味度高的合法那么美味度低的一定合法,因为美味度低的可选方案包含美味度高的可选方案. 那么我们二分一个美味度作为答案然后考虑如何验证? 选择时显然要贪心的 ...
洛谷 [P1024]一元三次方程求解【二分答案】
题目链接:https://www.luogu.org/problemnew/show/P1024 题目描述有形如:ax3+bx2+cx+d=0 这样的一个一元三次方程.给出该方程中各项的系数(a,b ...
洛谷 P2468 粟粟的书架二分（主席树+前缀和）
传送~:https://www.luogu.org/problem/P2468 看了一下数据也发现是两道题,后边当他是一个序列(n==1)的时候直接主席树二分区间前k大和就行了但是有一个细节我觉得就 ...
洛谷 P3384 【模板】树链剖分-树链剖分(点权)(路径节点更新、路径求和、子树节点更新、子树求和)模板-备注结合一下以前写的题目，懒得写很详细的注释...
P3384 [模板]树链剖分题目描述如题,已知一棵包含N个结点的树(连通且无环),每个节点上包含一个数值,需要支持以下操作: 操作1: 格式: 1 x y z 表示将树从x到y结点最短路径上所有节 ...
bzoj 2653 洛谷 P2839 [国家集训队] middle
2653: middle Time Limit: 20 Sec Memory Limit: 512 MB Submit: 2381 Solved: 1340 [Submit][Status][Di ...
洛谷P4315 月下“毛景树” 题解
洛谷P4315 月下"毛景树" 题解题目链接:P4315 月下"毛景树" 题意:请维护一个数据结构,支持改第 kkk 条边的边权结点 uuu 到 vvv ...
洛谷——P1047 校门外的树 python实现
洛谷--P1047 校门外的树 python实现 l, m = map(int, input().split(' ')) trees = [True]*(l+1) for i in range(m): ...

YBTOJ洛谷P2839：最大中位数（主席树、二分答案）

解析

代码

YBTOJ洛谷P2839：最大中位数（主席树、二分答案）相关推荐

最新文章

热门文章