BZOJ 1901 & 整体二分

/*==========================================================================
# Last modified: 2016-02-25 21:09
# Filename: 1901.cpp
# Description: 
==========================================================================*/
#define me AcrossTheSky 
#include <cstdio> 
#include <cmath> 
#include <ctime> 
#include <string> 
#include <cstring> 
#include <cstdlib> 
#include <iostream> 
#include <algorithm> 
  
#include <set> 
#include <map> 
#include <stack> 
#include <queue> 
#include <vector> 
 
#define lowbit(x) (x)&(-x) 
#define FOR(i,a,b) for((i)=(a);(i)<=(b);(i)++) 
#define FORP(i,a,b) for(int i=(a);i<=(b);i++) 
#define FORM(i,a,b) for(int i=(a);i>=(b);i--) 
#define ls(a,b) (((a)+(b)) << 1) 
#define rs(a,b) (((a)+(b)) >> 1) 
#define getlc(a) ch[(a)][0] 
#define getrc(a) ch[(a)][1] 
 
#define maxn 100000 
#define maxm 100000 
#define pi 3.1415926535898 
#define _e 2.718281828459 
#define INF 1070000000 
using namespace std; 
typedef long long ll; 
typedef unsigned long long ull; 
 
template<class T> inline 
void read(T& num) { 
    bool start=false,neg=false; 
    char c; 
    num=0; 
    while((c=getchar())!=EOF) { 
        if(c==‘-‘) start=neg=true; 
        else if(c>=‘0‘ && c<=‘9‘) { 
            start=true; 
            num=num*10+c-‘0‘; 
        } else if(start) break; 
    } 
    if(neg) num=-num; 
} 
/*==================split line==================*/ 
struct Infor{
	int pos,x,id,op,k,l,r,cur; 
	//k,l,r,cur,id专门针对询问,分别表示询问第k大,l,r表示询问的范围,cur表示贡献---->放到程序中更好解释,表示询问的序号
	//pos表示修改或原序列数值在序列中的位置
	//op表示操作类型:1表示在序列中的元素,2表示修改前的元素-----这里又有一个奇技淫巧,我们将所有的信息--->原序列,修改,询问放在一个区间中,那么对于一个修改,到目前为止应该在序列中的元素则为修改后的元素,就在区间中新建一个位置放原来那个数,标记为2,修改后的数标记为1.
	//x表示修改或原序列中的数值
}q[maxn],temp[maxn];
int c[maxn],a[maxn],ans[maxn],tmp[maxn];
bool mark[maxn];
int n,m,cnt=0,num=0;
//树状数组相关==================================================
void add(int x,int t){
	while (x<=n){
		c[x]+=t;
		x+=lowbit(x);
	}
}
int query(int x){
	int ret=0;
	while(x>0){
		ret+=c[x];x-=lowbit(x);
	}
	return ret;
}
//==============================================================
void solve(int l,int r,int L,int R){//l,r表示正在进行二分的区间,L,R表示答案的范围
	if (l>r) return;
	if (L==R) { //如果答案已经确定了,那么这个区间内每一个询问的答案都是L
		FORP(i,l,r) if (q[i].op==3) ans[q[i].id]=L;
		return;
	}
	int mid=rs(L,R);
	FORP(i,l,r){ //对于区间内每个的每个操作
		if (q[i].op==1 && q[i].x<=mid) add(q[i].pos,1); 
		//如果这个元素应该在序列中存在,并且它会对当前的答案产生影响,我们将该位置加上1.
		if (q[i].op==2 && q[i].x<=mid) add(q[i].pos,-1);
		//奇技淫巧
		//如果这个元素的op为2,那么之前一定有一个相等元素其op为1并且其已经对答案造成影响,如今这个元素被修改那么我们相应的要将它对答案的影响减去.
		if (q[i].op==3) tmp[i]=query(q[i].r)-query(q[i].l-1);
		//如果当前操作是询问操作,注意到一个时间性质,在这个询问之后的修改对这个询问没有影响,我们统计这个点询问范围内的贡献.
	}
	FORP(i,l,r){ //重置树状数组
		if (q[i].op==1 && q[i].x<=mid) add(q[i].pos,-1);
		if (q[i].op==2 && q[i].x<=mid) add(q[i].pos,1);
	}
	int tot=0;//统计在将区间二分时靠右操作的个数
	FORP(i,l,r){//划分操作
		if (q[i].op==3){ //这里就能体现cur的作用,如果当前比判断答案小的个数加上已经比判断答案小的个数大于k,那么我们要将判断答案减小,对于所有要这么操作的我们将它划到左边----于是就可以统一将判定答案的范围减小
			if (q[i].cur+tmp[i]>=q[i].k) mark[i]=true,tot++;
			else mark[i]=false,	q[i].cur+=tmp[i];
		}
		else{
			if (q[i].x<=mid) tot++,mark[i]=true; //对答案造成影响的,判定答案减小后仍可能造成影响.所以我们将它划到右边
			else mark[i]=false;
		}
	}
	int la=l,lb=l+tot;
	FORP(i,l,r)
		if (mark[i]) temp[la++]=q[i];
			else temp[lb++]=q[i];
	FORP(i,l,r) q[i]=temp[i]; //复制粘贴的过程
	solve(l,la-1,L,mid);
	solve(la,lb-1,mid+1,R); //整体二分的过程. 非常6
}
int main(){
	//freopen("a.in","r",stdin);
	//freopen("tmp.out","w",stdout);
	read(n); read(m);
	FORP(i,1,n){
		read(a[i]);
		q[++num].x=a[i]; q[num].pos=i; q[num].op=1; q[num].id=0;
	} //读入原序列
	FORP(i,1,m){
		char s[5]; scanf("%s",s);
		if (s[0]==‘Q‘){
			int x,y,z; read(x); read(y); read(z);
			q[++num].l=x; q[num].r=y; q[num].k=z; 
			q[num].id=++cnt; q[num].op=3;
		}//读入询问
		else {
			int x,t;
			read(t); read(x); 
			q[++num].x=a[t]; q[num].pos=t; q[num].id=0; q[num].op=2;//奇技淫巧,增加一个操作
			q[++num].x=x; q[num].pos=t; q[num].id=0; q[num].op=1;
			a[t]=x;
		}//读入修改操作
	}
	solve(1,num,0,INF); //整体二分
	FORP(i,1,cnt) printf("%d\n",ans[i]);//输出
}