【集训队互测2015】Robot

题目描述

http://uoj.ac/problem/88

题解

维护两颗线段树，维护最大值和最小值，因为每次只有单点查询，所以可以直接在区间插入线段就可以了。

注意卡常，不要写STL，用链表把同类修改串起来就好了。

代码

%:pragma GCC optimize()
%:pragma GCC optimize()
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<vector>
#include<cstdlib>
#define N 600009
using namespace std;
typedef long long ll;
struct node{ll tim,x;};
int n,m,la[N],pre[N];
ll t[N],pos[N],b[N],top,now[N],k[N],x[N],tag[N];
char s[];
inline ll rd(){
    ll x=;char c=getchar();bool f=;
    while(!isdigit(c)){if(c=='-')f=;c=getchar();}
    while(isdigit(c)){x=(x<<)+(x<<)+(c^);c=getchar();}
    return f?-x:x;
}
inline ll calc(ll l,ll r,ll x,ll k){return x+k*(r-l);}
struct segment1{
  struct seg{
      int l,r;ll x,k;
  }tr[N<<];
  int tot,root;
  void work(int &cnt,int l,int r,ll x,ll k){
      if(!cnt)cnt=++tot;
      int mid=(l+r)>>;
      if(!tr[cnt].x&&!tr[cnt].k){tr[cnt].x=x;tr[cnt].k=k;return;}
      ll x1=calc(b[l],b[mid],x,k),x2=calc(b[l],b[mid],tr[cnt].x,tr[cnt].k);
      if(x1>=x2){
          ll xx=tr[cnt].x,kk=tr[cnt].k;
          tr[cnt].x=x;tr[cnt].k=k;
          x=xx;k=kk;
      }
      if(l==r)return;
      if(x>=tr[cnt].x)work(tr[cnt].l,l,mid,x,k);else work(tr[cnt].r,mid+,r,x+(b[mid+]-b[l])*k,k);
  }
  void upd(int &cnt,int l,int r,int L,int R,ll x,ll k){
      if(!cnt)cnt=++tot;
      if(l>=L&&r<=R){
          work(cnt,l,r,x+(b[l]-b[L])*k,k);
          return;
      }
      int mid=(l+r)>>;
      if(mid>=L)upd(tr[cnt].l,l,mid,L,R,x,k);
      if(mid<R)upd(tr[cnt].r,mid+,r,L,R,x,k);
  }
}T1;
struct segment2{
  struct seg{
      int l,r;ll x,k;
  }tr[N<<];
  int tot,root;
  void work(int &cnt,int l,int r,ll x,ll k){
      if(!cnt)cnt=++tot;
      int mid=(l+r)>>;
      if(!tr[cnt].x&&!tr[cnt].k){tr[cnt].x=x;tr[cnt].k=k;return;}
      ll x1=calc(b[l],b[mid],x,k),x2=calc(b[l],b[mid],tr[cnt].x,tr[cnt].k);
      if(x1<=x2){
          ll xx=tr[cnt].x,kk=tr[cnt].k;
          tr[cnt].x=x;tr[cnt].k=k;
          x=xx;k=kk;
      }
      if(l==r)return;
      if(x<=tr[cnt].x)work(tr[cnt].l,l,mid,x,k);else work(tr[cnt].r,mid+,r,x+(b[mid+]-b[l])*k,k);
  }
  void upd(int &cnt,int l,int r,int L,int R,ll x,ll k){
      if(!cnt)cnt=++tot;
      if(l>=L&&r<=R){
          work(cnt,l,r,x+(b[l]-b[L])*k,k);
          return;
      }
      int mid=(l+r)>>;
      if(mid>=L)upd(tr[cnt].l,l,mid,L,R,x,k);
      if(mid<R)upd(tr[cnt].r,mid+,r,L,R,x,k);
  }
}T2;
ll query(int x){
      int now1=T1.root,now2=T2.root,l=,r=top;ll ans1=-1e15,ans2=1e15;
      while(now1||now2){
          int mid=(l+r)>>;
          if(now1)ans1=max(ans1,calc(b[l],b[x],T1.tr[now1].x,T1.tr[now1].k));
          if(now2)ans2=min(ans2,calc(b[l],b[x],T2.tr[now2].x,T2.tr[now2].k));
        if(x<=mid)now1=T1.tr[now1].l,now2=T2.tr[now2].l,r=mid;
        else now1=T1.tr[now1].r,now2=T2.tr[now2].r,l=mid+;
    }
    ans2=llabs(ans2);
    return max(ans1,ans2);
}
int main(){
    n=rd();m=rd();
    for(int i=;i<=n;++i)pos[i]=rd(),now[i]=;
    for(int i=;i<=m;++i){
        t[i]=rd();scanf("%s",s);b[++top]=t[i];
        if(s[]=='c'){
          k[i]=rd(),x[i]=rd();
          if(!now[k[i]])tag[k[i]]=i;
          pre[i]=now[k[i]];la[now[k[i]]]=i;
          now[k[i]]=i;
        }
    }
    t[m+]=t[m]+;b[++top]=;b[++top]=t[m]+;
    sort(b+,b+top+);
    top=unique(b+,b+top+)-b-;
    for(int i=;i<=m;++i)t[i]=lower_bound(b+,b+top+,t[i])-b;
    for(int i=;i<=n;++i){
        pre[m+]=now[i];la[now[i]]=m+;
        if(!now[i])tag[i]=m+;
    }
    for(int i=;i<=n;++i){
        int l=,r=t[tag[i]];
        T1.upd(T1.root,,top,l,r,pos[i],);
        T2.upd(T2.root,,top,l,r,pos[i],);
    }
    for(int i=;i<=m;++i){
        if(k[i]){
            int l=t[i],r=t[la[i]],la=t[pre[i]];
            pos[k[i]]=pos[k[i]]+(b[l]-b[la])*x[pre[i]];
            T1.upd(T1.root,,top,l,r,pos[k[i]],x[i]);
            T2.upd(T2.root,,top,l,r,pos[k[i]],x[i]);
        }
        else printf("%lld\n",query(t[i]));
    }
    return ;
}