【hihocoder1449】后缀自动机三·重复旋律6 后缀自动机

描述

小Hi平时的一大兴趣爱好就是演奏钢琴。我们知道一个音乐旋律被表示为一段数构成的数列。

现在小Hi想知道一部作品中所有长度为K的旋律中出现次数最多的旋律的出现次数。但是K不是固定的，小Hi想知道对于所有的K的答案。

解题方法提示

输入

共一行，包含一个由小写字母构成的字符串S。字符串长度不超过 1000000。

输出

共Length(S)行，每行一个整数，表示答案。

样例输入
aab
样例输出
2
1
1

代码

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
#define ll long long
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
char s[1000005];
int cnt,last,len,S,a[1000005];
int l[2000005],trans[2000005][27];
int fa[2000005],ans[2000005];
int du[2000005];queue<int>q;
void add(int x)
{
    int p=last,np=++cnt,c=s[x]-'a';last=np;l[np]=x;ans[np]=1;
    for (;p&&!trans[p][c];p=fa[p]) trans[p][c]=np;
    if (!p) fa[np]=S,du[S]++;
    else
    {
        int q=trans[p][c];
        if (l[p]+1==l[q]) fa[np]=q,du[q]++;
        else
        {
            int nq=++cnt;l[nq]=l[p]+1;
            memcpy(trans[nq],trans[q],sizeof trans[q]);
            fa[nq]=fa[q];fa[q]=fa[np]=nq;du[nq]=2;
            for (;trans[p][c]==q;p=fa[p]) trans[p][c]=nq;
        }
    }
}
int main()
{
    scanf("%s",s+1);
    len=strlen(s+1);
    last=S=++cnt;
    for (int i=1;i<=len;i++) add(i);
    for (int i=1;i<=cnt;i++) if (!du[i]) q.push(i);
    while (!q.empty())
    {
        int now=q.front();q.pop();
        ans[fa[now]]+=ans[now];
        du[fa[now]]--;if (!du[fa[now]]) q.push(fa[now]);
    }
    for (int i=1;i<=cnt;i++) a[l[i]]=max(a[l[i]],ans[i]);
    for (int i=len-1;i;i--) a[i]=max(a[i],a[i+1]);
    for (int i=1;i<=len;i++) printf("%d\n",a[i]);
    return 0;
}