bzoj 1014 [JSOI2008]火星人prefix（splay+hash）

【题目链接】

http://www.lydsy.com/JudgeOnline/problem.php?id=1014

【题意】

给定一个字符串，要求提供修改一个字符，插入一个字符，查询两个后缀LCP的功能。

【思路】

　 splay维护字符串的哈希值。因为要提供区间，splay采用先查找后调整至根的写法。

一个结点的hash值为:

ch[0]->h * X^(ch[1]->s+1)+v * X^(ch[1]->s)+ch[1]->h

    对于一个询问每次二分长度，提取区间后比较hash值即可。

需要注意的是splay要提前在区间的左右两边各加上一个节点，不然会调用到null。

ull自然溢出相当于模2^64。

【代码】

 #include<set>
 #include<cmath>
 #include<queue>
 #include<vector>
 #include<cstdio>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 #define trav(u,i) for(int i=front[u];i;i=e[i].nxt)
 #define FOR(a,b,c) for(int a=(b);a<=(c);a++)
 using namespace std;

 typedef long long ll;
 typedef unsigned long long ull;
 const int N = 5e5+;
 const int X = ;

 ll read() {
     char c=getchar();
     ll f=,x=;
     while(!isdigit(c)) {
         if(c=='-') f=-; c=getchar();
     }
     while(isdigit(c))
         x=x*+c-'',c=getchar();
     return x*f;
 }

 char s[N];
 int n,q;
 ull powx[N];

 struct Node* null;
 struct Node {
     int s,v; ull h;
     Node* ch[];
     int cmp(int k) {
         if(k==ch[]->s+) return -;
         return k<=ch[]->s? :;
     }
     void init(int x) {
         v=h=x; s=;
         ch[]=ch[]=null;
     }
     void maintain() {
         s=ch[]->s+ch[]->s+;
         h=ch[]->h*powx[ch[]->s+]+v*powx[ch[]->s]+ch[]->h;
     }
 } *root,nodepool[N]; int nodesz=;

 void rot(Node* &o,int d) {
     Node* k=o->ch[d^]; o->ch[d^]=k->ch[d]; k->ch[d]=o;
     o->maintain(),k->maintain(); o=k;
 }
 void splay(Node*& o,int k) {
     int d=o->cmp(k);
     if(d==) k-=o->ch[]->s+;
     if(d!=-) {
         Node* p=o->ch[d];
         int d2=p->cmp(k),k2=(d2==?k:k-p->ch[]->s-);
         if(d2!=-) {
             splay(p->ch[d2],k2);
             if(d==d2) rot(o,d^); else rot(o->ch[d],d);
         }
         rot(o,d^);
     }
 }
 //return range (l,r]
 //加过点后 s[l,r]=range(l,r+1)
 Node*& range(int l,int r) {
     splay(root,l);
     splay(root->ch[],r-l+);
     return root->ch[]->ch[];
 }

 Node* build(int l,int r)
 {
     if(r<l) return null;
     int mid=l+r>>;
     Node* o=&nodepool[++nodesz];
     o->init(s[mid]-'a'+);
     o->ch[]=build(l,mid-);
     o->ch[]=build(mid+,r);
     o->maintain();
     return o;
 }
 void insert(int p,int v)
 {
     splay(root,p+);
     Node* o=&nodepool[++nodesz];
     o->init(v);
     o->ch[]=root->ch[]; o->ch[]=null;
     o->maintain();
     root->ch[]=o; root->maintain();
 }
 void change(int p,int v)
 {
     splay(root,p);
     root->v=v;
     root->maintain();
 }

 int main()
 {
     //freopen("in.in","r",stdin);
     //freopen("out.out","w",stdout);
     null=new Node();
     scanf("%s",s+);
     int n=strlen(s+);
     s[]='z'+; s[++n]='z'+; s[n+]='\0';
     scanf("%d",&q);
     powx[]=;
     FOR(i,,n+q) powx[i]=powx[i-]*X;
     root=build(,n);
     while(q--) {
         char op[],val[];
         int x,y;
         scanf("%s%d",op,&x);
         if(op[]=='R') {
             scanf("%s",val);
             change(x+,val[]-'a'+);
         } else
         if(op[]=='I') {
             scanf("%s",val);
             insert(x+,val[]-'a'+);
         } else {
             scanf("%d",&y);
             int len=root->s,L=,R=;
             R=min(len-y-,len-x-);
             while(L<R) {
                 int M=L+(R-L+)/;
                 ull H=range(x,x+M)->h;
                 H-=range(y,y+M)->h;
                 if(!H) L=M; else R=M-;
             }
             printf("%d\n",L);
         }

     }
     return ;
 }