poj 2758 && BZOJ 2258 Checking the Text 文本校对

Description

 

为了给Wind买生日礼物，Jiajia不得不找了一份检查文本的工作。这份工作很无聊：给你一段文本 
要求比对从文本中某两个位置开始能匹配的最大长度是多少。但比无聊更糟糕的是，Jiajia的经理 
还可能往文本里面插入一些字符。 
Jiajia想用一个程序来解决这些繁琐的工作。这个程序的速度要足够快，因为Wind的生日就快要到了 
Jiajia必须赚到足够多的钱，也就是处理足够多的文本。

Input

输入文件第一行是原始文本。 
输入文件第二行是操作数n。此后n行，每行描述一条命令，命令有两种形式： 
I ch p：表示将一个字符ch插入到当前文本的第p个字符之前，如果p大于当前文本长度则表示插入到当前文本末尾； 
Q i j：表示询问当前文本从原始文本的第i个和第j个字符现在所在的位置开始能匹配的字符是多少。 
你可以认为文本初始长度不超过50000，I命令最多200条，Q命令最多20000条。

Output

对于每条Q命令输出一行，为最长匹配长度。

Sample Input

abaab
5
Q 1 2
Q 1 3
I a 2
Q 1 2
Q 1 3

Sample Output

0
1
0
3

 

直接hash暴力重构==1196 ms（bzoj）

treap维护hash==5512 ms（bzoj）

spaly维护hash==TLE

gg……

感谢居神贡献权限号

#include<cstdio>
#include<algorithm>
#include<cstring>
#define MN 100001
using namespace std;
int read_p,read_ca;
inline int read(){
    read_p=;read_ca=getchar();
    while(read_ca<''||read_ca>'') read_ca=getchar();
    while(read_ca>=''&&read_ca<='') read_p=read_p*+read_ca-,read_ca=getchar();
    return read_p;
}
unsigned long long hash[MN],mi[MN];
char c[MN],p[];
int n,m,ne[MN],pos,len,x,y;
inline void build(int x){
    for (int i=x;i<=len;i++) hash[i]=hash[i-]*+c[i];
}
inline int query(int x,int y){
    if (x>y) swap(x,y);
    int l=,r=len-y+,mid;
    while (l<r){mid=(l+r+)>>;if (hash[x+mid-]-hash[x-]*mi[mid]==hash[y+mid-]-hash[y-]*mi[mid]) l=mid;else r=mid-;}
    return l;
}
int main(){
    scanf("%s",c+);
    n=read();
    m=len=strlen(c+);
    mi[]=;
    for (int i=;i<=m+;i++) mi[i]=mi[i-]*;
    for (int i=;i<=len;i++) ne[i]=i;
    build();
    for (int i=;i<=n;i++){
        scanf("%s",p);
        if (p[]=='I'){
            scanf("%s",p);pos=read();
            if (pos>len) pos=len+;
            for (int i=len;i>=pos;i--) c[i+]=c[i];
            c[pos]=p[];len++;
            build(pos);
            for (int i=m;i;i--) if (ne[i]>=pos) ne[i]++;else break;
        }else{
            x=read();y=read();
            printf("%d\n",query(ne[x],ne[y]));
        }
    }
}

大暴力 2856 kb 1196 ms 1350 B

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#define MN 100001
using namespace std;
int read_p,read_ca;
inline int read(){
    read_p=;read_ca=getchar();
    while(read_ca<''||read_ca>'') read_ca=getchar();
    while(read_ca>=''&&read_ca<='') read_p=read_p*+read_ca-,read_ca=getchar();
    return read_p;
}
unsigned long long mi[MN];
char c[MN],p[];
int n,m,ne[MN],pos,len,x,y;
struct tree{
    int l,r,s,f,fa;
    unsigned long long h,c;
    tree(){
        l=r=s=h=c=;
    }
};
struct s_tree{
    int root,size;
    tree t[];
    s_tree(){
        root=size=;
    }
    inline void up(int x){
        t[x].h=(t[t[x].l].h*+t[x].c)*mi[t[t[x].r].s]+t[t[x].r].h;
    }
    inline void rir(int &x){
        int k=t[x].r;
        t[x].r=t[k].l;
        t[k].l=x;
        t[t[x].r].fa=x;
        t[k].fa=t[x].fa;
        t[x].fa=k;
        t[k].s=t[x].s;
        t[x].s=+t[t[x].l].s+t[t[x].r].s;
        t[k].h=t[x].h;up(x);
        x=k;
    }
    inline void lir(int &x){
        int k=t[x].l;
        t[x].l=t[k].r;
        t[k].r=x;
        t[t[x].l].fa=x;
        t[k].fa=t[x].fa;
        t[x].fa=k;
        t[k].s=t[x].s;
        t[x].s=+t[t[x].l].s+t[t[x].r].s;
        t[k].h=t[x].h;up(x);
        x=k;
    }
    inline void insert(int &p,int x,int c,int f){
        if (!p){
            p=++size;
            t[p].c=c;
            t[p].h=c;
            t[p].s=;
            t[p].f=rand();
            t[p].fa=f;
        }else{
            t[p].s++;
            if (t[t[p].l].s+>=x){
                insert(t[p].l,x,c,p);
                if (t[t[p].l].f<t[p].f) lir(p);
            }else{
                insert(t[p].r,x--t[t[p].l].s,c,p);
                if (t[t[p].r].f<t[p].f) rir(p);
            }
        }
        up(p);
    }
    inline int ne(int x){
        int mmh=;
        while (x){
            mmh+=t[t[x].l].s+;
            while (t[x].fa&&t[t[x].fa].l==x) x=t[x].fa;
            x=t[x].fa;
        }
        return mmh;
    }
    inline unsigned long long hash(int p,int x){
        if ((!x)||(!p)) return ;
        if (t[t[p].l].s==x-) return t[t[p].l].h*+t[p].c;else
        if (t[t[p].l].s>=x) return hash(t[p].l,x);else return (t[t[p].l].h*+t[p].c)*mi[x-t[t[p].l].s-]+hash(t[p].r,x--t[t[p].l].s);
    }
}t;
inline int query(int x,int y){
    x=t.ne(x);y=t.ne(y);
    if (x>y) swap(x,y);
    int l=,r=len-y+,mid;
    while (l<r){mid=(l+r+)>>;if (t.hash(t.root,x+mid-)-t.hash(t.root,x-)*mi[mid]==t.hash(t.root,y+mid-)-t.hash(t.root,y-)*mi[mid]) l=mid;else r=mid-;}
    return l;
}
int main(){
    scanf("%s",c+);
    n=read();
    m=len=strlen(c+);
    mi[]=;
    for (int i=;i<=1e5;i++) mi[i]=mi[i-]*;
    for (int i=;i<=len;i++) ne[i]=i;
    for (int i=;i<=len;i++) t.insert(t.root,i,c[i],);
    for (int i=;i<=n;i++){
        scanf("%s",p);
        if (p[]=='I'){
            scanf("%s",p);pos=read();
            if (pos>len) pos=len+;t.insert(t.root,pos,p[],);
            len++;
        }else{
            x=read();y=read();
            printf("%d\n",query(x,y));
        }
    }
}

treap 5604 kb 5512 ms 2929 B

暴力又短又快……有毒……