51nod1600 Simple KMP

Posted zh-comld

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了51nod1600 Simple KMP相关的知识,希望对你有一定的参考价值。

题目描述

对于一个字符串|S|,我们定义fail[i],表示最大的x使得S[1..x]=S[i-x+1..i],满足(x<i)
显然对于一个字符串,如果我们将每个0<=i<=|S|看成一个结点,除了i=0以外i向fail[i]连边,这是一颗树的形状,根是0
我们定义这棵树是G(S),设f(S)是G(S)中除了0号点以外所有点的深度之和,其中0号点的深度为-1
定义key(S)等于S的所有非空子串S‘的f(S‘)之和
给定一个字符串S,现在你要实现以下几种操作:
1.在S最后面加一个字符
2.询问key(S)

题解

遇到这种对所有子串统计的问题,考虑差分答案数组。

我们将答案数组二次差分之后,只需要算每个字符在所有以它结尾的子串中的贡献即可。

考虑这个\(border\)树的深度有什么意义。

观察或手玩即可发现,对于某个串来说,末尾字符位置的深度就是这个串前缀等于后缀的串的数,就是每一对前缀等于后缀都会对这个位置有1的贡献。

所以对于i位置,考虑前\(i-1\)个位置的字符串,我们现在只需要求所有以这个点结尾的串在前\(i-1\)的字符串的出现次数的和。

用链剖或\(LCT\)维护即可。

代码

#include<bits/stdc++.h>
#define N 200002
using namespace std;
typedef long long ll;
const int mod=1e9+7;
int la[N<<2],num[N],tot,head[N],fa[N],deep[N],son[N],dfn[N],_tag[N],top[N],n;
char s[N];
ll tr[N<<2],size[N<<2],ans[N];
inline ll rd(){
    ll x=0;char c=getchar();bool f=0;
    while(!isdigit(c)){if(c=='-')f=1;c=getchar();}
    while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
    return f?-x:x;
}
struct edge{int n,to;}e[N];
inline void MOD(ll &x){x=x>=mod?x-mod:x;}
struct SAM_t{
    int ch[N][26],fa[N],l[N],cnt,last;
    SAM_t(){cnt=last=1;}
    inline void ins(int x,int id){
        int p=last,np=++cnt;l[np]=l[p]+1;last=np;num[id]=cnt;
        for(;p&&!ch[p][x];p=fa[p])ch[p][x]=np;
        if(!p)fa[np]=1;
        else{
            int q=ch[p][x];
            if(l[p]+1==l[q])fa[np]=q;
            else{
                int nq=++cnt;l[nq]=l[p]+1;
                memcpy(ch[nq],ch[q],sizeof(ch[nq]));
                fa[nq]=fa[q];fa[q]=fa[np]=nq;
                for(;ch[p][x]==q;p=fa[p])ch[p][x]=nq; 
            }
        } 
    }
}sam;
inline void pushdown(int cnt){
    la[cnt<<1]+=la[cnt];
    la[cnt<<1|1]+=la[cnt];  
    MOD(tr[cnt<<1]+=size[cnt<<1]*la[cnt]%mod);
    MOD(tr[cnt<<1|1]+=size[cnt<<1|1]*la[cnt]%mod);
    la[cnt]=0;
}
ll upd(int cnt,int l,int r,int L,int R){;
    if(l>=L&&r<=R){
        ll x=tr[cnt];
        MOD(tr[cnt]+=size[cnt]);
        la[cnt]++;
        return x;
    }
    int mid=(l+r)>>1;
    if(la[cnt])pushdown(cnt);ll ans=0;
    if(mid>=L)MOD(ans+=upd(cnt<<1,l,mid,L,R));
    if(mid<R)MOD(ans+=upd(cnt<<1|1,mid+1,r,L,R));
    MOD(tr[cnt]=tr[cnt<<1]+tr[cnt<<1|1]);
    return ans;
}
void build(int cnt,int l,int r){
    if(l==r){
        size[cnt]=sam.l[_tag[l]]-sam.l[sam.fa[_tag[l]]];
        return;
    }
    int mid=(l+r)>>1;
    build(cnt<<1,l,mid);build(cnt<<1|1,mid+1,r);
    MOD(size[cnt]=size[cnt<<1]+size[cnt<<1|1]);
}
inline void add(int u,int v){e[++tot].n=head[u];e[tot].to=v;head[u]=tot;}
void dfs1(int u){
    size[u]=1;
    for(int i=head[u];i;i=e[i].n)if(e[i].to!=fa[u]){
        int v=e[i].to;fa[v]=u;deep[v]=deep[u]+1;
        dfs1(v);
        size[u]+=size[v];
        if(size[v]>size[son[u]])son[u]=v;
    }
}
void dfs2(int u){
    dfn[u]=++dfn[0];_tag[dfn[0]]=u;
    if(!top[u])top[u]=u;
    if(son[u])top[son[u]]=top[u],dfs2(son[u]); 
    for(int i=head[u];i;i=e[i].n)if(e[i].to!=fa[u]&&e[i].to!=son[u])dfs2(e[i].to);
}
ll work(int x){
    ll ans=0;
    while(top[x]!=top[1]){
        MOD(ans+=upd(1,1,sam.cnt,dfn[top[x]],dfn[x]));
        x=fa[top[x]]; 
    }
    MOD(ans+=upd(1,1,sam.cnt,dfn[1],dfn[x]));
    return ans;
}
int main(){
    n=rd();
    scanf("%s",s+1);
    for(int i=1;i<=n;++i)sam.ins(s[i]-'a',i);
    for(int i=1;i<=sam.cnt;++i)if(sam.fa[i])add(sam.fa[i],i);
    dfs1(1);dfs2(1);
    build(1,1,sam.cnt);
    for(int i=1;i<=n;++i){
        int x=work(num[i]);
        MOD(ans[i]=ans[i-1]+x);
    } 
    for(int i=1;i<=n;++i){
        MOD(ans[i]+=ans[i-1]);
        printf("%lld\n",ans[i]);
    }
    return 0;
}

以上是关于51nod1600 Simple KMP的主要内容,如果未能解决你的问题,请参考以下文章

51nod1006simple KMP

51nod 1600 Simplr KMP(后缀自动机+维护树上的数据结构)

[51nod 1129] 字符串最大值(kmp)

51nod 1554 KMP?????????

51Nod 1554 欧姆诺姆和项链 (KMP)

hdu6153 A Secret CCPC网络赛 51nod 1277 KMP