【BZOJ】1901: Zju2112 Dynamic Rankings-BIT套线段树

传送门:bzoj1901


题解

新姿势:BIT套线段树


代码

sro hwzer orz

#include<iostream>
#include<cstdio>
#include<algorithm>
#define N 2200001
using namespace std;
int n,m,tot,top,sz;
int v[10001],num[20001],hash[20001];
int flag[10001],A[10001],B[10001],K[10001],root[10001];
int sum[N],ls[N],rs[N];
int L[30],R[30],a,b;
int lowbit(int x)
{return x&(-x);}
int find(int x)
{
    int l=1,r=tot,mid;
    while(l<=r)
    {
         int mid=(l+r)>>1;
         if(hash[mid]<x)l=mid+1;
         else r=mid-1;
        } 
    return l;
}
void updata(int last,int l,int r,int &rt,int w,int x)
{
     rt=++sz;
     sum[rt]=sum[last]+x;ls[rt]=ls[last];rs[rt]=rs[last];
     if(l==r)  return;
     int mid=(l+r)>>1;
     if(w<=mid)  updata(ls[last],l,mid,ls[rt],w,x);
     else  updata(rs[last],mid+1,r,rs[rt],w,x);
}
int query(int l,int r,int k)
{
     if(l==r)  return l;
     int i,suml=0,sumr=0;
     for(i=1;i<=a;i++)  suml+=sum[ls[L[i]]];
     for(i=1;i<=b;i++)  sumr+=sum[ls[R[i]]];
     int mid=(l+r)>>1;
     if(sumr-suml>=k) 
     {
         for(i=1;i<=a;i++)  L[i]=ls[L[i]];
         for(i=1;i<=b;i++)  R[i]=ls[R[i]];
         return query(l,mid,k);
     }
     else
     {
         for(i=1;i<=a;i++)  L[i]=rs[L[i]];
         for(i=1;i<=b;i++)  R[i]=rs[R[i]];
         return query(mid+1,r,k-(sumr-suml));
     }
}
int main()
{
	char s[3];
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
	{
	    scanf("%d",&v[i]);
	    num[++top]=v[i];
    }
    for(int i=1;i<=m;i++)
    {
    	scanf("%s",s);
    	scanf("%d%d",&A[i],&B[i]);
    	if(s[0]=='Q'){scanf("%d",&K[i]);flag[i]=1;}
    	else num[++top]=B[i];
    }
    sort(num+1,num+top+1);
    hash[++tot]=num[1];
    for(int i=2;i<=top;i++)
       if(num[i]!=num[i-1])
          hash[++tot]=num[i];
    for(int i=1;i<=n;i++)
    {
    	int t=find(v[i]);
    	for(int j=i;j<=n;j+=lowbit(j))
    	    updata(root[j],1,tot,root[j],t,1);
    }
    for(int i=1;i<=m;i++)
    	if(flag[i])
    	{
    		a=0;b=0;A[i]--;
    		for(int j=A[i];j>0;j-=lowbit(j))
    		   L[++a]=root[j];
 		    for(int j=B[i];j>0;j-=lowbit(j))
 		       R[++b]=root[j];
            printf("%d\n",hash[query(1,tot,K[i])]);
    	}
    	else
        {
             int t=find(v[A[i]]);
             for(int j=A[i];j<=n;j+=lowbit(j))
			     updata(root[j],1,tot,root[j],t,-1);
             v[A[i]]=B[i];
             t=find(B[i]);
             for(int j=A[i];j<=n;j+=lowbit(j))
                 updata(root[j],1,tot,root[j],t,1);
        }
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值