[bzoj1058][ZJOI2007]报表统计 splay+卡常

Description

　　小Q的妈妈是一个出纳，经常需要做一些统计报表的工作。今天是妈妈的生日，小Q希望可以帮妈妈分担一些工

作，作为她的生日礼物之一。经过仔细观察，小Q发现统计一张报表实际上是维护一个可能为负数的整数数列，并

且进行一些查询操作。在最开始的时候，有一个长度为N的整数序列，并且有以下三种操作： INSERT i k 在原数

列的第i个元素后面添加一个新元素k； 如果原数列的第i个元素已经添加了若干元素，则添加在这些元素的最后（

见下面的例子） MIN_GAP 查询相邻两个元素的之间差值（绝对值）的最小值 MIN_SORT_GAP 查询所有元素中最接

近的两个元素的差值（绝对值） 例如一开始的序列为 5 3 1 执行操作INSERT 2 9将得到： 5 3 9 1 此时MIN_GAP

为2，MIN_SORT_GAP为2。 再执行操作INSERT 2 6将得到： 5 3 9 6 1 注意这个时候原序列的第2个元素后面已经

添加了一个9，此时添加的6应加在9的后面。这个时候MIN_GAP为2，MIN_SORT_GAP为1。于是小Q写了一个程序，使

得程序可以自动完成这些操作，但是他发现对于一些大的报表他的程序运行得很慢，你能帮助他改进程序么？

Input

　　第一行包含两个整数N，M，分别表示原数列的长度以及操作的次数。第二行为N个整数，为初始序列。接下来

的M行每行一个操作，即“INSERT i k”，“MIN_GAP”，“MIN_SORT_GAP”中的一种（无多余空格或者空行）。

Output

　　对于每一个“MIN_GAP”和“MIN_SORT_GAP”命令，输出一行答案即可。

Sample Input

3 5
5 3 1
INSERT 2 9
MIN_SORT_GAP
INSERT 2 6
MIN_GAP
MIN_SORT_GAP

Sample Output

2
2
1

HINT

 

N , M ≤500000 对于所有的数据，序列内的整数不超过5*10^8。

 

Source

 

思维难度:NOIP+

代码难度:省选

*卡常难度:NOI

#include<cstdio>
#include<iostream>
#define ll long long
using namespace std;
const int Maxn=1000005;
const ll inf=1000000000;
ll n,m,a[Maxn],up[Maxn],low[Maxn];
inline ll mn(ll x,ll y){
	return x<y?x:y;
}
inline ll mx(ll x,ll y){
	return x>y?x:y; 
}
inline ll read(){
	ll x=0,f=1;char c=getchar();
	while(c<'0'||c>'9'){
		if(c=='-')f=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9'){
		x=x*10+c-'0';
		c=getchar();
	} 
	return x*f;
}
namespace mg{
	ll ch[Maxn][2],fa[Maxn],cnt[Maxn],val[Maxn],tcnt=0,root=0;
	void rotate(ll x){
		ll y=fa[x];ll z=fa[y],k=ch[y][1]==x;ll w=ch[x][!k];
		if(z)ch[z][ch[z][1]==y]=x;
		ch[x][!k]=y;ch[y][k]=w;
		if(w)fa[w]=y;fa[y]=x;fa[x]=z;
	}
	void splay(ll x){
		ll y,z;
		while(fa[x]){
			y=fa[x];z=fa[y];
			if(z)rotate((ch[y][0]==x)^(ch[z][0]==y)?x:y);
			rotate(x);
		}
		root=x;
	}
	ll find(ll x){
		ll rt=root;
		while(1){
			if(x==val[rt])return rt;
			if(x<val[rt])rt=ch[rt][0];
			else rt=ch[rt][1];
		}
	}
	void insert(ll x){
		if(root==0){
			tcnt++;
			root=tcnt;
			val[tcnt]=x;
			cnt[tcnt]=1;
			return;
		}
		ll rt=root;
		while(1){
			if(x==val[rt]){
				cnt[rt]++;
				return;
			}
			if(x<val[rt]){
				if(!ch[rt][0]){
					val[++tcnt]=x;
					ch[rt][0]=tcnt;
					fa[tcnt]=rt;
					cnt[tcnt]=1;
					splay(tcnt);
					return;
				}
				rt=ch[rt][0];
			}
			else{
				if(!ch[rt][1]){
					val[++tcnt]=x;
					ch[rt][1]=tcnt;
					fa[tcnt]=rt;
					cnt[tcnt]=1;
					splay(tcnt);
					return;
				}
				rt=ch[rt][1];
			}
		}
	}
	void delet(ll x){
		ll rt=find(x);
		if(cnt[rt]>1){
			cnt[rt]--;
			return;
		}
		while(ch[rt][0]||ch[rt][1]){
			if(ch[rt][0])rotate(ch[rt][0]);
			else rotate(ch[rt][1]);
		}
		ch[fa[rt]][ch[fa[rt]][1]==rt]=0;
		splay(fa[rt]);
		fa[rt]=0;
	}
	ll query(){
		ll rt=root;
		while(ch[rt][0])rt=ch[rt][0];
		return val[rt];
	}
} 
namespace msg{
	ll ch[Maxn][2],fa[Maxn],cnt[Maxn],val[Maxn],tcnt=0,root=0;
	void rotate(ll x){
		ll y=fa[x];ll z=fa[y],k=ch[y][1]==x;ll w=ch[x][!k];
		if(z)ch[z][ch[z][1]==y]=x;ch[x][!k]=y;ch[y][k]=w;
		if(w)fa[w]=y;fa[y]=x;fa[x]=z;
	}
	void splay(ll x){
		ll y,z;
		while(fa[x]){
			y=fa[x];z=fa[y];
			if(z)rotate((ch[y][0]==x)^(ch[z][0]==y)?x:y);
			rotate(x);
		}
		root=x;
	}
	ll find(ll x){
		ll rt=root;
		while(1){
			if(val[rt]==x)return rt;
			if(x<val[rt])rt=ch[rt][0];
			else rt=ch[rt][1];
		}
	}
	void insert(ll x){
		if(root==0){
			tcnt++;
			root=tcnt;
			val[tcnt]=x;
			cnt[tcnt]=1;
			return;
		}
		ll rt=root;
		while(1){
			if(val[rt]==x){
				cnt[rt]++;
				splay(rt);
				return;
			}
			if(x<val[rt]){
				if(!ch[rt][0]){
					tcnt++;
					val[tcnt]=x;
					cnt[tcnt]=1;
					ch[rt][0]=tcnt;
					fa[tcnt]=rt;
					splay(tcnt);
					return;
				}
				rt=ch[rt][0];
			}
			else{
				if(!ch[rt][1]){
					tcnt++;
					val[tcnt]=x;
					cnt[tcnt]=1;
					ch[rt][1]=tcnt;
					fa[tcnt]=rt;
					splay(tcnt);
					return;
				}
				rt=ch[rt][1];
			} 
		}
	}
	ll lower(ll x){
		ll rt=root,ans=-inf;
		while(rt){
			if(x>=val[rt]&&val[rt]>ans)ans=val[rt];
			if(x>val[rt])rt=ch[rt][1];
			else rt=ch[rt][0];
		}
		return ans;
	}
	ll upper(ll x){
		ll rt=root,ans=inf;
		while(rt){
			if(x<=val[rt]&&ans>val[rt])ans=val[rt];
			if(x<val[rt])rt=ch[rt][0];
			else rt=ch[rt][1];
		}
		return ans;
	}
}
inline ll ab(ll x){
	return x>0?x:-x;
}
int main(){
	char inc[25];ll u,v,ans=inf,key;
	n=read();m=read();
	for(int i=1;i<=n;++i){
		a[i]=read();
		low[i]=up[i]=a[i];
		ans=mn(ans,mn(msg::upper(a[i])-a[i],a[i]-msg::lower(a[i])));
		msg::insert(a[i]);
		if(i>1)mg::insert(ab(a[i]-a[i-1])); 
	}
	while(m--){
		scanf("%s\n",inc+1);
		if(inc[1]=='I'){
			u=read();v=read();
			ans=mn(ans,mn(msg::upper(v)-v,v-msg::lower(v)));
			msg::insert(v);
			mg::insert(ab(v-up[u]));if(u!=n)mg::insert(ab(v-low[u+1]));
			mg::delet(ab(low[u+1]-up[u]));
			up[u]=v;
		}
		else{
			if(inc[5]=='S'){
				printf("%lld\n",ans);
			}
			else{
				printf("%lld\n",mg::query());
			}
		}
	}
	return 0;
}