ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

2019牛客多校第七场E Find the median 离散化+线段树维护区间段

2019-08-27 21:50:22  阅读:248  来源: 互联网

标签:int 第七场 sum median 多校 ++ zz 操作 define


Find the median

题意

刚开始集合为空,有n次操作,每次操作往集合里面插入[L[i],R[i]]的值,问每次操作后中位数是多少

分析

由于n比较大,并且数可以达到1e9,我们无法通过权值线段树来进行操作,那么怎么办呢?题目中还有什么性质?插入的值是一段一段的,那么我们是不是能从这些段中入手?维护这些段,怎么维护呢,如果[1,2][2,5]这两段有一个点重合那该怎么办,此时我们可以用一个常规操作把r加1进行分段,什么意思呢例如上例我们分成[1,2][2,3][3,6]那么我们对[1,2]进行操作,那么只要对第一段进行操作即可,如果我们对[2,5]进行操作,那么只要对第二段和第三段进行操作即可,这样就把区间分开来了,相当于把公共点抽出来了放到了两个段中,这样就不会导致本来只想对一段进行操作,结果相邻的段也被操作了的尴尬情形。起始这种操作在线段树中不少见,刷模板题的时候就有过,目的也相同,只是补完才想起来?
还有因为分段的问题,段树会需要n2的点,也就是需要开n8的空间,RE了1W年,我傻了。

#include<bits/stdc++.h>
#define pb push_back
#define F first
#define S second
#define pii pair<int,int>
#define mkp make_pair
using namespace std;
typedef long long ll;
const int maxn=8e5+4;
ll sum[maxn<<2],lazy[maxn<<2];
ll x[maxn],y[maxn],c[maxn*3];
int id1[maxn],id2[maxn];
ll a1,a2,b1,b2,c1,c2,m1,m2;
int sz;
inline void build(int o,int l,int r){
    lazy[o]=sum[o]=0;
    if(l==r)return ;
    int mid=l+r>>1;
    build(o<<1,l,mid);
    build(o<<1|1,mid+1,r);
}
inline void push_down(int o,int l,int r){
    if(lazy[o]){
        int mid=l+r>>1;
        sum[o<<1]+=lazy[o]*(c[mid+1]-c[l]);
        sum[o<<1|1]+=lazy[o]*(c[r+1]-c[mid+1]);
        lazy[o<<1]+=lazy[o];
        lazy[o<<1|1]+=lazy[o];
        lazy[o]=0;
    }
}
inline void update(int o,int l,int r,int x,int y)
{
    if(x<=l&&y>=r){
        sum[o]+=(c[r+1]-c[l]);
        lazy[o]++;
    }
    else {
        int mid=l+r>>1;
        push_down(o,l,r);
        if(mid>=x)update(o<<1,l,mid,x,y);
        if(mid<y)update(o<<1|1,mid+1,r,x,y);
        sum[o]=sum[o<<1|1]+sum[o<<1];
    }
}
inline ll query(int o,int l,int r,ll v){
    if(l==r){
        int cishu=sum[o]/(c[r+1]-c[l]);
        return c[l]+(v-1)/cishu;
    }
    else {
        int mid=l+r>>1;
        push_down(o,l,r);
        if(sum[o<<1]>=v)return query(o<<1,l,mid,v);
        else return query(o<<1|1,mid+1,r,v-sum[o<<1]);
    }
 
}
 
int main(){
    int n;
    scanf("%d",&n);
    scanf("%lld%lld%lld%lld%lld%lld",&x[1],&x[2],&a1,&b1,&c1,&m1);
    scanf("%lld%lld%lld%lld%lld%lld",&y[1],&y[2],&a2,&b2,&c2,&m2);
     
    for(int i=3;i<=n;i++){
        x[i]=(1ll*a1*x[i-1]+1ll*b1*x[i-2]+c1)%m1;
        y[i]=(1ll*a2*y[i-1]+1ll*b2*y[i-2]+c2)%m2;
    }
    int zz=0;
    for(int i=1;i<=n;i++){
        x[i]++,y[i]++;
        if(x[i]>y[i])swap(x[i],y[i]);
        c[++zz]=x[i];
        c[++zz]=y[i]+1;
    }
    sort(c+1,c+1+zz);
    sz=unique(c+1,c+1+zz)-(c+1);
    for(int i=1;i<=n;i++){
        id1[i]=lower_bound(c+1,c+1+sz,x[i])-c;
        id2[i]=lower_bound(c+1,c+1+sz,y[i]+1)-c;
    }
    sz--;
    build(1,1,sz);
    ll cnt=0;
    for(int i=1;i<=n;i++){
        update(1,1,sz,id1[i],id2[i]-1);
        cnt+=(y[i]-x[i]+1);
        printf("%lld\n",query(1,1,sz,(cnt+1)/2));
    }
    return 0;
}

标签:int,第七场,sum,median,多校,++,zz,操作,define
来源: https://www.cnblogs.com/ttttttttrx/p/11421061.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有