51nod 1600 Simple KMP

又被机房神犇肉丝哥哥和glory踩爆了

首先这个答案的输出方式有点套路，当前的答案=上一个答案+每一个后缀的f值=上一个答案+上一次算的每个后缀的f值+当前每个后缀的深度

这个题意给了个根深度为-1有点诡异，考虑它的现实意义是这个后缀在前面出现了几次，这些后缀的深度和就是前面有多少子串和后缀是能匹配的

考虑用SAM把fail树搞出来，那么对于加入一个后缀，就是这个叶子到根的路径的点+1，询问也是问当前点到根的和

拿个LCT维护一下跑路

离线树剖好像更好写 肉丝早就跑路了

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
const int _=1e2;
const int maxn=1e5+_;
const LL mod=1e9+;

struct node
{
    int f,son[];
    LL s,c;//子树内子串总出现数 子树内子串数
    LL tot,num;//管理子串数，节点出现次数
    LL lazy;
}tr[*maxn];
void update(int now)
{
    int lc=tr[now].son[],rc=tr[now].son[];
    tr[now].c=tr[lc].c+tr[rc].c+tr[now].tot;
    tr[now].s=(tr[lc].s+tr[rc].s+tr[now].tot*tr[now].num+(tr[now].c-tr[now].tot)*tr[now].lazy)%mod;
}
void pushdown(int now)
{
    int lc=tr[now].son[],rc=tr[now].son[];
    if(lc!=)tr[lc].num+=tr[now].lazy,tr[lc].lazy+=tr[now].lazy,update(lc);
    if(rc!=)tr[rc].num+=tr[now].lazy,tr[rc].lazy+=tr[now].lazy,update(rc);
    tr[now].lazy=;
}
void rotate(int x,int w)
{
    int f=tr[x].f,ff=tr[f].f;
    int R,r;

    R=f,r=tr[x].son[w];
    tr[R].son[-w]=r;
    if(r!=)tr[r].f=R;

    R=ff,r=x;
         if(tr[ff].son[]==f)tr[R].son[]=r;
    else if(tr[ff].son[]==f)tr[R].son[]=r;
    tr[r].f=R;

    R=x,r=f;
    tr[R].son[w]=r;
    tr[r].f=R;

    update(f);
    update(x);
}
bool isroot(int x,int rt)
{
    int f=tr[x].f;
    if(f==rt|| tr[f].son[]!=x&&tr[f].son[]!=x )return true;
    else return false;
}
int tt,tmp[*maxn];
void splay(int x,int rt)
{
    int i=x; tt=;
    while(!isroot(i,rt))
        tmp[++tt]=i,i=tr[i].f;
    tmp[++tt]=i;
    while(tt>)
    {
        if(tr[tmp[tt]].lazy!=)pushdown(tmp[tt]);
        tt--;
    }

    while(!isroot(x,rt))
    {
        int f=tr[x].f,ff=tr[f].f;
        if(isroot(f,rt))
        {
            if(tr[f].son[]==x)rotate(x,);
            else rotate(x,);
        }
        else
        {
                 if(tr[ff].son[]==f&&tr[f].son[]==x)rotate(f,),rotate(x,);
            else if(tr[ff].son[]==f&&tr[f].son[]==x)rotate(f,),rotate(x,);
            else if(tr[ff].son[]==f&&tr[f].son[]==x)rotate(x,),rotate(x,);
            else if(tr[ff].son[]==f&&tr[f].son[]==x)rotate(x,),rotate(x,);
        }
    }
}
void access(int x)
{
    int y=;
    while(x!=)
    {
        splay(x,);
        tr[x].son[]=y;
        if(y!=)tr[y].f=x;
        update(x);
        y=x;x=tr[y].f;
    }
}
void gotop(int x){access(x),splay(x,);}
void Link(int x,int y){gotop(x),tr[x].f=y;}
void Cut(int x)
{
    gotop(x);
    int lc=tr[x].son[];
    tr[lc].f=;tr[x].son[]=;
    update(x);
}
void add(int x)
{
    gotop(x);
    tr[x].num++;tr[x].lazy++;
    update(x);
}
LL getsum(int x){gotop(x);return tr[x].s;}
void change(int x,LL tot)
{
    gotop(x);
    tr[x].tot=tot;
    update(x);
}

struct SAM
{
    int w[],dep,fail;
}ch[*maxn];int cnt,last;
int gettot(int x){return ch[x].dep-ch[ch[x].fail].dep;}
LL insert(int k,int x)
{
    int now=++cnt,pre=last; LL ret=;
    ch[now].dep=k;
    while(pre!=&&ch[pre].w[x]==)ch[pre].w[x]=now,pre=ch[pre].fail;
    if(pre==)
    {
        ch[now].fail=;
        change(now,gettot(now));
        Link(now,ch[now].fail);
        add(now);
    }
    else
    {
        int nxt=ch[pre].w[x];
        if(ch[nxt].dep==ch[pre].dep+)
        {
            ch[now].fail=nxt;
            change(now,gettot(now));
            Link(now,ch[now].fail);
            ret=getsum(now);
            add(now);
        }
        else
        {
            int nnxt=++cnt; ch[nnxt]=ch[nxt];
            ch[nnxt].dep=ch[pre].dep+;
            ch[nxt].fail=nnxt;

            change(nxt,gettot(nxt));
            tr[nnxt].num=tr[nxt].num;change(nnxt,gettot(nnxt));
            Cut(nxt);
            Link(nnxt,ch[nnxt].fail);
            Link(nxt,ch[nxt].fail);

            //.......先处理nnxt和nxt.......

            ch[now].fail=nnxt;
            change(now,gettot(now));
            Link(now,ch[now].fail);
            ret=getsum(now);
            add(now);

            //........再搞nnxt和now.......

            while(pre!=&&ch[pre].w[x]==nxt)ch[pre].w[x]=nnxt,pre=ch[pre].fail;
        }
    }
    last=now;
    return ret;
}

char ss[maxn];
int main()
{
//    freopen("a.in","r",stdin);
//    freopen("a.out","w",stdout);
    int n;
    scanf("%d%s",&n,ss+);

    LL ans=,sum=;
    cnt=last=;
    for(int i=;i<=n;i++)
    {
        sum=(sum+insert(i,ss[i]-'a'+))%mod;
        ans=(ans+sum)%mod;
        printf("%lld\n",ans);
    }

    return ;
}