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]);
        }
    }
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值