zoj 3533 Gao the String I(伸展树+字符串hash)

题意：给出一个长度为n的串，有四种操作：

①翻转区间[a,b]②把第p个字符改成c③求以a和b开头的最大公共前缀④求当前串的最大回文子串

思路：维护序列当然伸展树比较方便啦~求LCP可以通过字符串hash，求最大回文比较少，可以直接Manacher搞定。之前做了类似的，所以思路上完全没问题，不过犯了二，debug了好久。。。

代码：

#include<iostream>#include<cstdio>#include<cstring>#include<string>#include<algorithm>#include<map>#include<queue>#include<stack>#include<cmath>#include<vector>#define inf 0x3f3f3f3f#define Inf 0x3FFFFFFFFFFFFFFFLL#define eps 1e-9#define pi acos(-1.0)using namespace std;typedef long long ll;typedef unsigned long long Ull;const int maxn=100000+10;const int rp=233;Ull hash[maxn][2],pw[maxn];int ch[maxn][2],size[maxn],val[maxn],flip[maxn];int pre[maxn],p[maxn<<1],root,tot,cnt;char str[maxn],s[maxn<<1];void NewNode(int &rt,int fa,int v){    rt=++tot;    ch[rt][0]=ch[rt][1]=flip[rt]=0;    pre[rt]=fa;size[rt]=1;    hash[rt][0]=hash[rt][1]=val[rt]=v;}inline void PushUp(int rt){    int ls=ch[rt][0],rs=ch[rt][1];    size[rt]=size[ls]+size[rs]+1;    hash[rt][0]=hash[ls][0]+val[rt]*pw[size[ls]]+hash[rs][0]*pw[(size[ls]+1)];    hash[rt][1]=hash[rs][1]+val[rt]*pw[size[rs]]+hash[ls][1]*pw[(size[rs]+1)];}inline void fson(int rt){    if(!rt) return ;    swap(ch[rt][0],ch[rt][1]);    swap(hash[rt][0],hash[rt][1]);    flip[rt]^=1;}void PushDown(int rt){    if(flip[rt])    {        fson(ch[rt][0]);        fson(ch[rt][1]);        flip[rt]=0;    }}void Rotate(int x,int kind){    int y=pre[x];    PushDown(y);PushDown(x);    ch[y][kind^1]=ch[x][kind];    pre[ch[x][kind]]=y;    if(pre[y]) ch[pre[y]][ch[pre[y]][1]==y]=x;    pre[x]=pre[y];    ch[x][kind]=y;    pre[y]=x;    PushUp(y);}void Splay(int rt,int goal){    PushDown(rt);    while(pre[rt]!=goal)    {        int y=pre[rt];        if(pre[y]==goal)            Rotate(rt,ch[y][0]==rt);        else        {            int kind=ch[pre[y]][0]==y;            if(ch[y][kind]==rt)            {                Rotate(rt,kind^1);                Rotate(rt,kind);            }            else            {                Rotate(y,kind);                Rotate(rt,kind);            }        }    }    PushUp(rt);    if(goal==0) root=rt;}void RotateTo(int k,int goal){    int rt=root;    PushDown(rt);    while(size[ch[rt][0]]!=k)    {        if(size[ch[rt][0]]>k)            rt=ch[rt][0];        else        {            k-=(size[ch[rt][0]]+1);            rt=ch[rt][1];        }        PushDown(rt);    }    Splay(rt,goal);}void Modify(int k,int v){    RotateTo(k-1,0);    RotateTo(k+1,root);    int rt=ch[ch[root][1]][0];    hash[rt][0]=hash[rt][1]=val[rt]=v;    PushUp(ch[root][1]);    PushUp(root);}void Reverse(int L,int R){    if(L>R) swap(L,R);    RotateTo(L-1,0);    RotateTo(R+1,root);    fson(ch[ch[root][1]][0]);    PushUp(ch[root][1]);    PushUp(root);}bool check(int x,int y,int len){    Ull xval,yval;    RotateTo(x-1,0);    if(size[ch[root][1]]-1<len) return 0;    RotateTo(x+len,root);    int rt=ch[ch[root][1]][0];    xval=hash[rt][flip[rt]];    RotateTo(y-1,0);    if(size[ch[root][1]]-1<len) return 0;    RotateTo(y+len,root);    rt=ch[ch[root][1]][0];    yval=hash[rt][flip[rt]];    return xval==yval;}int Lcp(int x,int y){    if(x>y) swap(x,y);    int l=1,r=size[root]-1;    if(!check(x,y,1)) return 0;    while(r-l>1)    {        int m=(l+r)>>1;        if(check(x,y,m)) l=m;        else r=m;    }    return l;}void dfs(int rt){    if(rt)    {        PushDown(rt);        dfs(ch[rt][0]);        str[cnt++]=val[rt];        dfs(ch[rt][1]);    }}int Palindrome(){    cnt=0;    dfs(root);    str[cnt-1]=0;    int n=0;    s[n++]='$';s[n++]='#';    for(int i=1;str[i];++i)    {        s[n++]=str[i];        s[n++]='#';    }    int mx=0,id=0,ans=0;    p[0]=0;    for(int i=1;i<n;++i)    {        p[i]=1;        if(mx>i) p[i]=min(p[id*2-i],mx-i);        while(s[i+p[i]]==s[i-p[i]]) p[i]++;        if(i+p[i]>mx)        {            mx=i+p[i];            id=i;        }        if(p[i]>ans) ans=p[i];    }    return ans-1;}void build(int l,int r,int &rt,int fa){    if(l>r) return ;    int m=(l+r)>>1;    NewNode(rt,fa,str[m]);    build(l,m-1,ch[rt][0],rt);    build(m+1,r,ch[rt][1],rt);    PushUp(rt);}void Init(int n){    ch[0][0]=ch[0][1]=pre[0]=0;    hash[0][0]=hash[0][1]=0;    flip[0]=val[0]=size[0]=0;    root=tot=0;    NewNode(root,0,0);    NewNode(ch[root][1],root,0);    build(0,n-1,ch[ch[root][1]][0],ch[root][1]);    PushUp(ch[root][1]);    PushUp(root);}int main(){    //freopen("in.txt","r",stdin);    //freopen("out.txt","w",stdout);    pw[0]=1;    for(int i=1;i<maxn;++i) pw[i]=pw[i-1]*rp;    int n,m,a,b;    char qs[15];    while(~scanf("%s",str))    {        n=strlen(str);        Init(n);        scanf("%d",&m);        while(m--)        {            scanf("%s",qs);            if(qs[0]=='P') printf("%d\n",Palindrome());            else if(qs[0]=='M')            {                scanf("%d%s",&a,qs);                Modify(a,qs[0]);            }            else            {                scanf("%d%d",&a,&b);                if(qs[0]=='R') Reverse(a,b);                else printf("%d\n",Lcp(a,b));            }        }    }    return 0;}