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

最新推荐文章于 2022-03-23 21:46:51 发布
最新推荐文章于 2022-03-23 21:46:51 发布
阅读量148 收藏 1
点赞数
分类专栏： 数据结构 二分 可持久化权值线段树
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_51068403/article/details/115566803
版权

传送门

文章目录

题意：

在这里插入图片描述

思路：

我们先解决怎么判断中位数的问题，我们可以二分一个 m i d mid mid，将 < m i d <mid <mid的值都变成 − 1 -1 −1，其他的数都变成 1 1 1，那么当全部的和 > = 0 >=0 >=0的时候，就说明当前数可以为中位数，且可能会变得更大，让后继续二分判断就好。
以上过程的复杂度是 O ( q n l o g n ) O(qnlogn) O(qnlogn)的，显然不能接受，但是可以发现我们的预处理之后的数组是可以递推出来的，不需要每次都算一遍，比如当前 m i d mid mid为 1 1 1，那么就全都是 1 1 1， m i d mid mid为 2 2 2的时候，只需要将 1 1 1改成 − 1 -1 −1即可，这个过程我们可以用主席树可持久化一下， r o o t root root的下标为中位数的值的位置。
让后我们考虑怎么判断区间选择哪个，首先我们可以发现 [ 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]我们只需要取一个最长前缀，让后判断 s u m + l s + r s > = 0 sum+ls+rs>=0 sum+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
*/
YB Lin
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
2015高中数学二次函数在指定区间上的最值与值域学案无答案新人教A版必修1
08-05
2015高中数学二次函数在指定区间上的最值与值域学案无答案新人教A版必修1
P2839 [国家集训队]middle
weixin_30270561的博客
02-14 97
题面 　　• 提一下静态区间第k小的nlog2n的做法： 　　　　1. 建关于排名的主席树（按排名顺序建树）。 　　　　2. 二分答案。 　　• 这样做静态区间第k小的虽然有些ZZ，但它的意义在于将线段树 　　 维护的对象改变了。 1 #include<bits/stdc++.h> 2 using namespace std; 3 int...
主席树练习6——P2839 [国家集训队]middle
Love_mona的博客
04-03 270
蒟蒻的垂死挣扎这题写得。。。。不想说什么了。首先，对于每次询问，我们可以知道，范围内比中位数小的数是要等于比他大的数的，而这题有重复就意味着小的数与大的数之差要在这个中位数个数的范围内，为了取到最大，我们自然是要它取到极致。我们考虑用二分答案来取到这个值，每次如果大于等于它的数 要大于等于 比他小的数，则这个答案是满足的，就可以继续缩范围，这样我们就得到了一个不错的思路。那么我们如何快速解决che...
洛谷P2839 [国家集训队]middle（主席树）
qq_45863710的博客
09-11 153
P2839 [国家集训队]middle 我们可以考虑二分中位数 checkcheckcheck 答案，那么我们对于某个值 midmidmid ，把 [l,r][l,r][l,r] 内的所有小于 midmidmid 的数赋值为 −1-1−1 ，大于等于 midmidmid 的数赋为 111 ，若 [l,r][l,r][l,r] 的和大于等于 000 就代表这个区间内的中位数至少都是 midmidmid ，否则只可能小于 midmidmid 。答案一定是序列里面的数字，所以我们就考虑对于每个数字建出一颗线段.
P2839 [国家集训队]middle（可持久化线段树+二分）
蒟蒻Steven_Meng的博客
07-11 232
题目 传送门 一个长度为nnn的序列aaa，设其排过序之后为bbb，其中位数定义为b[n/2]b[n/2]b[n/2]，其中a,ba,ba,b从000开始标号,除法取下整。 给你一个长度为nnn的序列sss。 回答QQQ个这样的询问：sss的左端点在[a,b][a,b][a,b]之间,右端点在[c,d][c,d][c,d]之间的子序列中，最大的中位数。 其中a&lt;b&lt;c&...
洛谷 P2839 [国家集训队]middle（可持久化线段树+二分）
C_Dreamy的博客
04-12 171
const int N=2e4+5; int n,m,_; int i,j,k; int a[N]; struct Node { int l,r; int sum,lsum,rsum; void update(int x) { sum = lsum = rsum = x ; } Node (int l=0,int r=0,in...
江苏高二数学复习学案+练习6 函数值域和最值(二) 文 试题.doc
10-26
江苏高二数学复习学案+练习6 函数值域和最值(二) 文 试题.doc
江苏高二数学复习学案+练习5 函数值域和最值(一) 文 试题.doc
10-26
江苏高二数学复习学案+练习5 函数值域和最值(一) 文 试题.doc
关于多值域一一对应，分类后一一对应显示的问题
01-11
关于多值域一一对应，分类后一一对应显示的问题
决策树算法下EMR的测试分类指标
最新发布
05-08
在决策树中，每个非叶节点表示一个特征属性上的测试，每个分支代表这个特征属性在某个值域上的输出，而每个叶节点存放一个类别。 决策树的应用场景非常广泛，包括但不限于以下几个方面： 金融风险评估：用于预测...
P2839 [国家集训队]middle 解题报告
日居月诸的博客
01-21 141
P2839 [国家集训队]middle 解题报告 题目链接 https://www.luogu.com.cn/problem/P2839 题目大意 给一个长度为 nnn 的序列。多次询问，每次给出四个参数 a,b,c,da,b,c,da,b,c,d ，要求找一个子区间 [l,r][l,r][l,r]，满足 a≤l≤b,c≤r≤da\le l\le b,c\le r\le da≤l≤b,c≤r≤d，使得子区间的中位数最大。输出中位数的最大值。 n≤2⋅104n\le2\cdot 10^4n≤2⋅104，询问
P2839 [国家集训队]middle【好题】【主席树+二分】
slass-mot-odet
09-04 199
题目描述 一个长度为n的序列a，设其排过序之后为b，其中位数定义为b[n/2]，其中a,b从0开始标号,除法取下整。 给你一个长度为n的序列s。 回答Q个这样的询问：s的左端点在[a,b]之间,右端点在[c,d]之间的子序列中，最大的中位数。 其中a<b<c<d。位置也从0开始标号。我会使用一些方式强制你在线。 输入格式 第一行序列长度n。 接下来n行按顺序给出a中的数。 接下来...
P2839 畅通工程
weixin_33738982的博客
08-26 78
P2839 畅通工程 题目描述 某省调查城镇交通状况，得到现有城镇道路统计表，表中列出了每条道路直接连通的城镇。省政府“畅通工程”的目标是使全省任何两个城镇间都可以实现交通（但不一定有直接的道路相连，只要互相间接通过道路可达即可）。问最少还需要建设多少条道路？ 输入输出格式 输入格式： 第1行给出两个正整数，分别是城镇数目N ( < 1000 )和道路数目M 随后的M行对...
P2839 [国家集训队]middle 主席树+求中位数技巧
weixin_43824564的博客
07-15 167
首先做这个题 我们要用到一个常用的技巧 二分求中位数 具体步骤就是 先二分一个值 在数组中把大于等于这个值的数都设为1 把小于这个值的数都设为-1 如果数组和>=0 那么 l=mid+1 否则的话 r=mid-1 那么我们如何快速的进行这个操作呢 假如一个询问是 a,b,c,d a<b<c<d 并且对于 区间a~d我们都进行了上述的赋值操作 那么 b+1~c-1之间的值是肯定要算进去的 然后我们需要在 a~b区间内找个最大后缀和...
P3085 [USACO13OPEN]Yin and Yang G 点分治
m0_51068403的博客
02-19 2732
文章目录题意：思路： 传送门 题意： 给你一颗nnn个点的树，每条边为黑色或者白色，问满足以下条件的路径条数：路径上存在一个不是端点的点，使得两端点到该点的两条路径上两种颜色的边数相等。 1≤n≤1000001\le n\le 1000001≤n≤100000 思路： 统计树上路径问题显然需要用到点分治了，这个题维护的信息比较麻烦，想明白了思路还需要考虑如何做才能使代码变得简便好写。 考虑点分治每个步骤需要算的贡献，在选出重心之后，之后的子树中就不会包含重心这个点，所以贡献需要算每颗包含重心的子树内的路径以
The 2021 ICPC Asia Taipei Regional F. What a Colorful Wall 扫描线 + 并查集
m0_51068403的博客
03-23 1156
文章目录题意:思路 传送门 题意: 给你平面nnn个矩形，每个矩形有一种颜色，依次给出矩形以及其的颜色，后面的矩形会覆盖前面的矩形，问最终有多少种颜色。 1≤n≤4000,0≤x1<x2<228,0≤y1<y2<228,1≤c≤n1\le n\le 4000,0\le x_1<x_2< 2^{28},0\le y_1<y_2< 2^{28},1\le c\le n1≤n≤4000,0≤x1​<x2​<228,0≤y1​<y2​<228,
可持久化Splay 学习笔记
m0_51068403的博客
05-12 673
可持久化Splay是怎么回事呢？Splay相信大家都很熟悉，但是可持久化Splay是怎么回事呢，下面就让小编带大家一起了解吧。 　　可持久化Splay，其实就是将Splay持久化一下，大家可能会很惊讶Splay怎么可以持久化呢？但事实就是这样，小编也感到非常惊讶。 　　这就是关于可持久化Splay的事情了，大家有什么想法呢，欢迎在评论区告诉小编一起讨论哦！ ...
XXI Open Cup. Grand Prix of Korea I. Query On A Tree 17 树剖 + 二分 + 树带权重心
m0_51068403的博客
08-23 579
传送门 文章目录题意：思路： 题意： 给你一棵树，每棵树初始权值都为000，现在给你两个操作： (1)(1)(1)将uuu的子树权值全部加111。 (2)(2)(2)将(u,v)(u,v)(u,v)路径上的点权值都加111。 每次输出一个点xxx,满足∑y=1Na[y]∗dis(x,y)\sum_{y=1}^Na[y]*dis(x,y)∑y=1N​a[y]∗dis(x,y)最小，如果多个点相同，输出深度最小的点。 思路： 观察这个式子，发现就是一个带权重心的问题，如果不在重心显然向重心走更优，需要支持修改。
【洛谷】【线段树】P3353 在你窗外闪耀的星星
05-14
对于一个查询操作 $Q(l,r,k)$，需要在值域线段树上找到 $[l,r]$ 对应的区间，然后查询该区间内权值大于 $k$ 的点的个数。这个可以通过线段树的区间查询操作实现。 对于一个修改操作 $C(x,y)$，需要在值域线段树上...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值