【BZOJ1901】【 ZOJ2112】- Dynamic Rankings(整体二分)

传送门

其实和前面那道K-th Number几乎一模一样吧

只需要在中途添加修改操作就可以了

注意由于整体二分的特殊性

修改操作要先将原来的贡献给排除掉在添加这一次操作的贡献

还有数组大小

具体看代码就可以了

#include<algorithm>
#include<iostream>
#include<cstdio>
#define ll long long
#define int long long
using namespace std;
inline int read(){
    char ch=getchar();
    int res=0,f=1;
    while(!isdigit(ch)){if(ch=='-')f=-f;ch=getchar();}
    while(isdigit(ch))res=(res<<3)+(res<<1)+(ch^48),ch=getchar();
    return res*f;
}
#define lc (u<<1)
#define rc  ((u<<1)|1)
#define mid ((l+r)>>1)
const int N=200005;
struct ask{
    int l,r,k,pos,op;
}q[N],q1[N],q2[N];
int n,m,tot,cnt,a[N],ans[N],tr[N<<2],add[N];
inline void pushup(int u){
    tr[u]=tr[lc]+tr[rc];
}
inline void pushdown(int u,int l,int r){
    tr[lc]+=(mid-l+1)*add[u],tr[rc]+=(r-mid)*add[u],add[lc]+=add[u],add[rc]+=add[u],add[u]=0;
}
inline void update(int u,int l,int r,int st,int des,int k){
    if(st<=l&&r<=des){
        tr[u]+=(r-l+1)*k,add[u]+=k;return;
    }
    pushdown(u,l,r);
    if(st<=mid)update(lc,l,mid,st,des,k);
    if(mid<des)update(rc,mid+1,r,st,des,k);
    pushup(u);
}
inline ll query(int u,int l,int r,int st,int des){
    if(st<=l&&r<=des)return tr[u];
    pushdown(u,l,r);
    ll ans=0;
    if(st<=mid)ans+=query(lc,l,mid,st,des);
    if(des>mid)ans+=query(rc,mid+1,r,st,des);
    pushup(u);
    return ans;
}
#undef mid
inline void solve(int l,int r,int st,int des){
	if(l>r||st>des)return;
    if(l==r){
        for(int i=st;i<=des;i++)if(q[i].op)ans[q[i].pos]=l;
        return;
    }
    int mid=(l+r)>>1,cnt1=0,cnt2=0;
    for(int i=st;i<=des;i++){
        if(q[i].op==1){
            int tmp=query(1,1,n,q[i].l,q[i].r);
            if(q[i].k<=tmp)q1[++cnt1]=q[i];
            else q[i].k-=tmp,q2[++cnt2]=q[i];
        }
        else{
            if(q[i].k<=mid)update(1,1,n,q[i].l,q[i].r,1),q1[++cnt1]=q[i];
            else q2[++cnt2]=q[i];
        }
    }
    for(int i=1;i<=cnt1;i++)if(!q1[i].op)update(1,1,n,q1[i].l,q1[i].r,-1);
    for(int i=1;i<=cnt1;i++)q[st+i-1]=q1[i];
    for(int i=1;i<=cnt2;i++)q[st+cnt1+i-1]=q2[i];
    solve(l,mid,st,st+cnt1-1),solve(mid+1,r,st+cnt1,des);
}
signed main(){
//	freopen("lx.txt","r",stdin);
//	freopen("my.out","w",stdout);
    n=read(),m=read();
    for(int i=1;i<=m;i++){
        int op=read(),l=read(),r=read(),k=read();
        if(op==1)q[i]=(ask){l,r,-k,0,0};
        else q[i]=(ask){l,r,k,++tot,1};
    }
    solve(-n,n,1,m);
    for(int i=1;i<=tot;i++)cout<<-ans[i]<<'\n';
}

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值