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

时间限制:15000ms

单点时限:3000ms

内存限制:512MB

描述

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

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

解题方法提示

输入

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

输出

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

样例输入

aab

样例输出

2
1
1

---

---

#include <bits/stdc++.h>
using namespace std;
const int maxn = 2*1e6+10;


int last =1,tot=1;
struct NODE 
{
    int ch[26];
    int len,fa;
    long long  maxlen=0;
    NODE () {memset(ch,0,sizeof ch); len=fa=0;}

}dian[maxn];
long long  cnt[maxn];
int nn=0;
inline void add(int c)
{
   int p = last; int np = last = ++tot; 
   dian[np].maxlen=1; // 这种点包含了一个后缀，所以比他所有的儿子多1，maxlen设为一 
   dian[np].len = dian[p].len+1;
   for( ; p && !dian[p].ch[c]; p=dian[p].fa)dian[p].ch[c] = np;

   if(!p) dian[np].fa = 1,cnt[1]++;
   else 
   {
    //   cout<<"HERE"<<endl;
       int q=dian[p].ch[c];
       if(dian[q].len == dian[p].len+1)dian[np].fa=q,cnt[q]++;
       else
       {
           int nq=++tot;
           dian[nq]=dian[q]; 
           dian[nq].maxlen=0;// 注意这一步，直接复制点的话，要把maxlen清零
           //一开始没注意到这里，wa无数次 
           dian[nq].len=dian[p].len+1;
           dian[q].fa=dian[np].fa=nq;
           cnt[nq]+=2;
           for(;p&&dian[p].ch[c]==q;p=dian[p].fa)dian[p].ch[c]=nq;
       }
   }

}
long long  ans[maxn];

int topsort()
{
      queue<int > q;
      for(int i=1;i<=tot;i++)if(!cnt[i])q.push(i);
      while(!q.empty())
      {
           int t=q.front(); q.pop();
           dian[dian[t].fa].maxlen += dian[t].maxlen;
           if(--cnt[dian[t].fa] == 0)q.push(dian[t].fa);
      }
      // 整个就是向上的递推 
      
      for(int i=1;i<=tot;i++)ans[dian[i].len] = max(ans[dian[i].len],dian[i].maxlen );
      for(int i=nn-1;i>=1;i--)ans[i] = max(ans[i],ans[i+1]);
      for(int i=1;i<=nn;i++)cout<<ans[i]<<endl;
      
}


signed main()
{
    string s;cin>>s; int len=s.length();  nn =len;
    for (int i=0;i<len;i++)add(s[i]-'a');
  //  for(int i=1;i<=tot;i++)cout<<dian[i].maxlen<<" ";
    topsort();
}