【UOJ #88】【集训队互测2015】Robot（线段树分治 / 斜率优化 / 李超树）

传送门

对于函数$f_i=k_{i}x+b_i$求最大或最小

显然是李超树板子

考虑对时间线段树分治一下

然后再在线段树内求最大最小值
把值列出来发现是一个斜率优化的式子
分别对最大/小值维护上/下凸壳即可

注意对于最小值，第一个点显然就是在$-\infty$处最小的
所以按斜率优化中的横坐标从大到小即从右往左做

考虑建凸包时要先排序
可以先最开始的时候排一次再插入，这样可以做到$O(nlogn)$

但由于常数太大跑不过两个$log$的李超树
喜提$uojrk$倒一

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define bg begin
#define pii pair<int,int>
#define fi first
#define cs const
#define se second
#define pb push_back
#define y1 shinkle
cs int RLEN=1<<20|1;
inline char gc(){
	static char ibuf[RLEN],*ib,*ob;
	(ib==ob)&&(ob=(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
	return (ib==ob)?EOF:*ib++;
}
inline int read(){
	char ch=gc();
	int res=0;bool f=1;
	while(!isdigit(ch))f^=ch=='-',ch=gc();
	while(isdigit(ch))res=(res+(res<<2)<<1)+(ch^48),ch=gc();
	return f?res:-res;
}
inline ll readll(){
	char ch=gc();
	ll res=0;bool f=1;
	while(!isdigit(ch))f^=ch=='-',ch=gc();
	while(isdigit(ch))res=(res+(res<<2)<<1)+(ch^48),ch=gc();
	return f?res:-res;
}
inline int readstring(char *s){
	int top=0;char ch=gc();
	while(isspace(ch))ch=gc();
	while(!isspace(ch)&&ch!=EOF)s[++top]=ch,ch=gc();
	return top;
}
template<typename tp>inline void chemx(tp &a,tp b){a=max(a,b);}
template<typename tp>inline void chemn(tp &a,tp b){a=min(a,b);}
cs int N=100005,M=600005;
char xxx;
typedef pair<int,ll> pt;
inline double slope1(cs pt &x,cs pt &y){
	if(x.fi==y.fi)return x.se<y.se?5e18:-5e18;
	return 1.0*(y.se-x.se)/(y.fi-x.fi);
}
inline double slope2(cs pt &x,cs pt &y){
	if(x.fi==y.fi)return x.se<y.se?-5e18:5e18;
	return 1.0*(y.se-x.se)/(y.fi-x.fi);
}
struct robot{
	int x,t,pre;
	ll a;
}p[N];
namespace tr{
	vector<pt> tr[M<<2],tl[M<<2];
	pt q[M];
	int pl[M<<2],pr[M<<2];
	#define lc (u<<1)
	#define rc ((u<<1)|1)
	#define mid ((l+r)>>1)
	inline void getc(vector<pt> tr,vector<pt> &ans){
		pt x;int tp=0;
		for(int i=0;i<tr.size();i++){
			x=tr[i];
			while(tp>=2&&slope1(q[tp-1],q[tp])<=slope1(q[tp],x))tp--;
			q[++tp]=x;
		}
		ans.resize(tp);
		for(int i=1;i<=tp;i++)ans[i-1]=q[i];
	}
	inline void getd(vector<pt> tr,vector<pt> &ans){
		pt x;int tp=0;
		for(int i=0;i<tr.size();i++){
			x=tr[i];
			while(tp>=2&&slope2(q[tp-1],q[tp])<=slope2(q[tp],x))tp--;
			q[++tp]=x;
		}
		ans.resize(tp);
		for(int i=1;i<=tp;i++)ans[i-1]=q[i];
	}
	inline ll F(pt x,int t){return (ll)x.fi*t+x.se;}
	void update(int u,int l,int r,int st,int des,pt k){
		if(st<=l&&r<=des){
			tr[u].pb(k);return;
		}
		if(st<=mid)update(lc,l,mid,st,des,k);
		if(mid<des)update(rc,mid+1,r,st,des,k);
	}
	void build(int u,int l,int r){
		vector<pt>tmp=tr[u];
		getc(tmp,tr[u]);reverse(tmp.bg(),tmp.end());
		getd(tmp,tl[u]),pr[u]=pl[u]=0;
		if(l==r)return;
		build(lc,l,mid),build(rc,mid+1,r);
	}
	struct opt{
		int l,r;pt k;
		friend inline bool operator <(cs opt &a,cs opt &b){
			return a.k.fi<b.k.fi;
		}
	}qe[M];int num;
	void update(int st,int des,pt k){
		num++;qe[num].l=st,qe[num].r=des,qe[num].k=k;
	}
	void pushupd(int l,int r){
		sort(qe+1,qe+num+1);
		for(int i=1;i<=num;i++)
		update(1,l,r,qe[i].l,qe[i].r,qe[i].k);
	}
	ll query(int u,int l,int r,int p,int k){
		while(pr[u]+1<tr[u].size()&&F(tr[u][pr[u]],k)<=F(tr[u][pr[u]+1],k)){
			pr[u]++;
		}
		while(pl[u]+1<tl[u].size()&&F(tl[u][pl[u]],k)>=F(tl[u][pl[u]+1],k)){
			pl[u]++;
		}
		ll res=0;
		if(pr[u]<tr[u].size())chemx(res,F(tr[u][pr[u]],k));
		if(pl[u]<tl[u].size())chemx(res,-F(tl[u][pl[u]],k));
		if(l==r)return res;
		if(p<=mid)return max(res,query(lc,l,mid,p,k));
		else return max(res,query(rc,mid+1,r,p,k));
	}
	#undef lc
	#undef rc
	#undef mid
}
int n,m,qn;
char str[15];
pii qr[600005];
char yyy;
signed main(){
	#ifdef Stargazer
	freopen("robot9.in","r",stdin);
	freopen("my.out","w",stdout);
	#endif
	n=read(),m=read();;
	for(int i=1;i<=n;i++)p[i].a=read(),p[i].x=0,p[i].t=0,p[i].pre=0;
	for(int i=1;i<=m;i++){
		int t=read();
		readstring(str);
		if(str[1]=='c'){
			int k=read(),x=read();
			tr::update(p[k].pre+1,i,pt(p[k].x,p[k].a));
			p[k].a=p[k].a+1ll*p[k].x*t-(ll)t*x,p[k].t=t,p[k].x=x,p[k].pre=i;
		}
		else{
			qr[++qn]=pii(i,t);
		}
	}
	for(int i=1;i<=n;i++)if(p[i].pre<m){
		tr::update(p[i].pre+1,m,pt(p[i].x,p[i].a));
	}
	tr::pushupd(1,m);
	tr::build(1,1,m);
	for(int i=1;i<=qn;i++){
		cout<<tr::query(1,1,m,qr[i].fi,qr[i].se)<<'\n';
	}return 0;
}