bzoj 1014 [JSOI2008]火星人prefix

填一波坑。。當初懶得挑這題。。(還是我太弱。。寫的時候腦子不好使)。

這題很裸。。但我的錯誤太傻了。。

首先不要相信什麼長度一定小於100000,不開二倍會死的。。。

然後每次增加點的時候要把全長+1.。。。(我就錯這了。。調了1h)。

之後就是很裸的啦。

維護hash值,正常的splay,對於增加點,就是把前驅轉上去,後繼轉到前驅的右子節點,新的點插入左子節點。

修改就是正常splay上去修改就好。

查詢的時候利用二分,對於每一個答案進行驗證就好了。

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#define seed 13131
#define N 100005
using namespace std;
typedef unsigned long long ull;
struct node
{
    int son[2],siz,fa;
    ull has;
    char val;
}tr[2*N];
ull cj[2*N];
char ss[2*N],ins[10];
int root,tot,n,m;
void update(int x)
{
    tr[x].siz=tr[tr[x].son[0]].siz+tr[tr[x].son[1]].siz+1;
    tr[x].has=tr[tr[x].son[0]].has*cj[tr[tr[x].son[1]].siz+1]+cj[tr[tr[x].son[1]].siz]*tr[x].val+tr[tr[x].son[1]].has;
}
void rotate(int x)
{
    int y=tr[x].fa,z=tr[y].fa;
    int typ=(x==tr[y].son[1]);
    tr[y].son[typ]=tr[x].son[typ^1],tr[tr[x].son[typ^1]].fa=y;
    tr[x].son[typ^1]=y,tr[y].fa=x;
    tr[x].fa=z,tr[z].son[tr[z].son[1]==y]=x;
    update(y);update(x);
}
int findkth(int x,int k)
{
    int t=tr[tr[x].son[0]].siz;
    if(k<=t)return findkth(tr[x].son[0],k);
    else if(k==t+1)return x;
    else return findkth(tr[x].son[1],k-t-1);
}
void splay(int x,int goal)
{
    for(int y;(y=tr[x].fa)!=goal;rotate(x))
    {
        if(tr[y].fa!=goal)rotate((x==tr[y].son[0])==(y==tr[tr[y].fa].son[0])?y:x);
    }
    if(!goal)root=x;
}
int biuld(int l,int r,int f)
{
    if(l>r)return 0;
    int mid=(l+r)/2;
    int rt=++tot;
    tr[rt].val=ss[mid];tr[rt].siz=1;tr[rt].fa=f;//has會自己更新進去
 	tr[rt].son[0]=biuld(l,mid-1,rt);
 	tr[rt].son[1]=biuld(mid+1,r,rt);
 	update(rt);
 	return rt;
}
void insert(int po,char s)//x
{
    int x=findkth(root,po+1);
    int y=findkth(root,po+2);
    splay(x,0);splay(y,x);
    int rt=++tot;n++;
    tr[y].son[0]=rt;
    tr[rt].val=s;tr[rt].siz=1;tr[rt].fa=y;/*tr[rt].has=s*/
    update(rt);update(y);update(x);
}
void change(int k,char v)
{
    int x=findkth(root,k+1);
    splay(x,0);
    tr[x].val=v;
    update(x);
}
int LCQ(int xx,int yy)//x->x,x+mid+2
{
    int l=0,r=tot-2,mid;
    while(l<=r)
    {
        mid=(l+r)/2;
        if(yy+mid-1>n)
        {
            r=mid-1;continue;
        }
        int x=findkth(root,xx),y=findkth(root,xx+mid+1);
        splay(x,0);splay(y,x);
        ull tt=tr[tr[y].son[0]].has;
        x=findkth(root,yy),y=findkth(root,yy+mid+1);
        splay(x,0);splay(y,x);
        if(tt==tr[tr[y].son[0]].has)l=mid+1;
        else r=mid-1;
    }
    return r;
}
int main()
{
    scanf("%s",ss+2);n=strlen(ss+2);scanf("%d",&m);
    cj[0]=1;
    for(int i=1;i<=2*N-2;i++)
    {
        cj[i]=cj[i-1]*seed;
    }
    root=biuld(1,n+2,0);
    int t1,t2;char t3[2];
    for(int i=1;i<=m;i++)
    {
        scanf("%s",ins);
        if(ins[0]=='Q')
        {
            scanf("%d%d",&t1,&t2);
            if(t1>t2)swap(t1,t2);
        //	if(t1==t2)printf("%d\n",n-t1+1);
            printf("%d\n",LCQ(t1,t2));
        }else if(ins[0]=='R')
        {
            scanf("%d%s",&t1,t3);change(t1,t3[0]);
        }else 
        {
            scanf("%d%s",&t1,t3);insert(t1,t3[0]);
        }
    }
}

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章