【JSOI2008】火星人

题面

https://www.luogu.org/problem/P4036

题解

#include<cstdio>
#include<iostream>
#include<cstring>
#define ri register int
#define N 105000
#define uLL unsigned long long
using namespace std;

const uLL p=107;
uLL pp[N];
char s[N],op[3],c[5];
int n,a,b,ptt;

struct Splay{
  int ch[N][2],f[N],siz[N];
  uLL sum[N],v[N];
  int tot;
  void clear() {tot=1; siz[1]=1;}  // &&&***
  bool opt(int x) {return ch[f[x]][1]==x;} 
  int new_node(int x,int ff) {
    ++tot;
    sum[tot]=v[tot]=x; // &&&***
    siz[tot]=1;
    f[tot]=ff;
    return tot;
  }
  void pushup(int x) {
    if (!x) return; // &&&***
    siz[x]=siz[ch[x][0]]+siz[ch[x][1]]+1;
    sum[x]=sum[ch[x][0]]+v[x]*pp[siz[ch[x][0]]]+sum[ch[x][1]]*pp[siz[ch[x][0]]+1];
  }
  void rotate(int x) {
    int y=f[x],z=f[y],s=opt(x),w=ch[x][1^s];
    f[x]=z; ch[z][opt(y)]=x;
    ch[y][s]=w; if (w) f[w]=y;
    f[y]=x; ch[x][1^s]=y;
    pushup(y); pushup(x);
  }
  void splay(int at,int to) {
    while (1) {
      if (f[at]==to) return;
      if (f[f[at]]==to) rotate(at);
      else {
        if (opt(at)==opt(f[at])) rotate(f[at]),rotate(at);
        else rotate(at),rotate(at);
      }
    }
  }
  int findbyrank(int x){
    int les=0,now=1;
    while (1) {
      if (les+siz[ch[now][0]]==x) return now;
      else if (les+siz[ch[now][0]]>x) now=ch[now][0]; else les+=siz[ch[now][0]]+1,now=ch[now][1];
    }
  }
  uLL getval(int l,int r) {
    int x=findbyrank(l-1);
    splay(x,1);
    if (r==ptt) return sum[ch[x][1]];
    else {
      int y=findbyrank(r+1);
      splay(y,x);
      return sum[ch[y][0]]; // &&&***
    }
  }
} t;

int main() {
  //freopen("tt.in","r",stdin);
  pp[0]=1;
  for (ri i=1;i<N;i++) pp[i]=pp[i-1]*p;
  scanf("%s",s+1);
  int L=strlen(s+1);
  t.clear();
  for (ri i=1;i<=L;i++) {
    int x=t.findbyrank(i-1);
    t.splay(x,1);
    t.ch[x][1]=t.new_node(s[i]-'a',x);
    t.pushup(x);
  }
  scanf("%d",&n);
  ptt=L;
  while (n--) {
    scanf("%s",op);
    if (op[0]=='Q') {
      scanf("%d %d",&a,&b);
      int Lb=1,Rb=min(ptt-a+1,ptt-b+1),ans=0;
      while (Lb<=Rb) {
        int mid=(Lb+Rb)/2;
        uLL t1=t.getval(a,a+mid-1);
        uLL t2=t.getval(b,b+mid-1);
        //cout<<t1<<" "<<t2<<endl;
        if (t1==t2) ans=mid,Lb=mid+1; else Rb=mid-1;
      }
      printf("%d
",ans);
    }
    else if (op[0]=='I') {
      scanf("%d %s",&a,c);
      int x=t.findbyrank(a);
      t.splay(x,1);
      if (a==ptt) {
        t.ch[x][1]=t.new_node(c[0]-'a',x);
        t.pushup(x);
      }
      else {
        int y=t.findbyrank(a+1);
        t.splay(y,x);
        t.ch[y][0]=t.new_node(c[0]-'a',y); // &&&***
        t.pushup(y); t.pushup(x);  // &&&***
      }
      ptt++;
    }
    else if (op[0]=='R') {
      scanf("%d %s",&a,c);
      int x=t.findbyrank(a);
      t.splay(x,1);
      t.v[x]=c[0]-'a';
      t.pushup(x);
    }
  }
  return 0;
}