HDU 5412 CRB and Queries (Kth number 整体二分 动态转静态)

题意：

给出长度为n的数列，m次操作

操作分为 1 x y 将x位置的数改成y

2 x y k 查询【x ， y】区间内的第k小数

分析：

有了POJ 2104 那题的静态整体二分的基础后就非常好做了

还是离线处理的方法，将所有数据读入然后再整体二分

对于修改操作，将它分为删除操作和添加操作，仔细想想，这样做并不影响整体二分的结果，因为一次修改操作只对后面的有影响，而后面的影响会随着二分消除掉。

要注意一点，由于修改操作改为了两种操作，所以总空间应该开到 n+2*m

#include<cstring>
#include<cstdio>
#include<iostream>
#include<algorithm>

using namespace std;

#define INF 0x3f3f3f3f
#define maxn 320000

struct node
{
int val,l,r,k;
int index,kind;
} q[maxn],q1[maxn],q2[maxn];

int n,m;
int a[maxn];
int c[maxn];
int ans[maxn];

inline int lowbit(int x)
{
return x&-x;
}

inline void update(int x,int val)
{
for(; x<=n; x+=lowbit(x)) c[x]+=val;
}

inline int query(int x)
{
int sum=0;
for(; x>0; x-=lowbit(x)) sum+=c[x];
return sum;
}

void divide(int s,int t,int l,int r)
{
if(s>t) return ;
if(l==r)
{
for(int i=s; i<=t; i++)
{
if(q[i].kind==2)
{
ans[q[i].index]=l;
}
}
return ;
}
int mid=(l+r)>>1;
int num1=0,num2=0;
int flag1=0,flag2=0;
for(int i=s; i<=t; i++)
{
if(q[i].kind==2)
{
int tmp=query(q[i].r)-query(q[i].l-1);
if(tmp>=q[i].k) q1[num1++]=q[i],flag1=1;
else q[i].k-=tmp,q2[num2++]=q[i],flag2=1;
}
else
{
if(q[i].val<=mid) update(q[i].index,q[i].kind),q1[num1++]=q[i];
else q2[num2++]=q[i];
}
}
for(int i=s; i<=t; i++)
{
if(q[i].kind!=2)
{
if(q[i].val<=mid) update(q[i].index,-q[i].kind);
}
}
for(int i=0; i<num1; i++) q[s+i]=q1[i];
for(int i=0; i<num2; i++) q[s+num1+i]=q2[i];
if(flag1) divide(s,s+num1-1,l,mid);
if(flag2) divide(s+num1,t,mid+1,r);
}

int main()
{
while(scanf("%d",&n)!=EOF)
{
memset(c,0,sizeof c);
int cnt=1;
int mm=0;
for(int i=1; i<=n; i++)
{
scanf("%d",&a[i]);
q[cnt].val=a[i];
q[cnt].index=i;
q[cnt].kind=1;
cnt++;
}
scanf("%d",&m);
int num=0;
for(int i=1; i<=m; i++)
{
int op;
scanf("%d",&op);
if(op==2)
{
int l,r,k;
scanf("%d%d%d",&l,&r,&k);
q[cnt].l=l;
q[cnt].r=r;
q[cnt].k=k;
q[cnt].kind=2;
q[cnt].index=num++;
cnt++;
}
else
{
int x,val;
scanf("%d%d",&x,&val);
q[cnt].val=a[x];
q[cnt].index=x;
q[cnt].kind=-1;
cnt++;
q[cnt].index=x;
q[cnt].val=val;
q[cnt].kind=1;
cnt++;
a[x]=val;
}
}
divide(1,cnt-1,1,INF);
for(int i=0; i<num; i++)
{
printf("%d\n",ans[i]);
}
}
}