HihoCoder1449 重复旋律6(后缀自动机)

描述

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

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

输入

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

输出

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

样例输入

aab

样例输出

2

1

1

1，很好的利用了拓扑。

2，利用了单调性。

  for(i=tot;i>;i--){

        ans[i]=max(ans[i],ans[i+]);

    }

2017-11-24 代码是舶来品，参考hihocoder，等把后缀自动机消化了自己再写一遍吧。

2017-11-25 这几个题理解得差不多了。

#include<iostream>

#include<cstring>

#include<cstdio>

#include<cstring>

#include<algorithm>

#include<cmath>

#include<string>

#include<queue>

using namespace std;

typedef long long ll;

const int inf=0x3f3f3f3f;

const int N=1e6+;

const int M=1e6+;

int tot,slink[*N],trans[*N][],minlen[*N],maxlen[*N],edpts[*N],n;

char str[*N];

int blue[*N],ind[*N],ans[*N+];

int newstate(int _maxlen,int _minlen,int* _trans,int _slink) {

    maxlen[++tot]=_maxlen;

    minlen[tot]=_minlen;

    slink[tot]=_slink;

    if(_trans)

        for(int i=; i<; i++)

            trans[tot][i]=_trans[i];

    return tot;

}

int add(char ch,int u) {

    int c=ch-'a',v=u;

    int z=newstate(maxlen[u]+,-,NULL,);

    blue[z]=;//绿色

    while(v&&!trans[v][c]) {

        trans[v][c]=z;

        v=slink[v];

    }

    if(!v) {

        minlen[z]=;

        slink[z]=;

        ind[]++;

        return z;

    }

    int x=trans[v][c];

    if(maxlen[v]+==maxlen[x]) {

        slink[z]=x;

        minlen[z]=maxlen[x]+;

        ind[x]++;

        return z;

    }

    int y=newstate(maxlen[v]+,-,trans[x],slink[x]);

    slink[z]=slink[x]=y;

    ind[y]+=;

    minlen[x]=minlen[z]=maxlen[y]+;

    while(v&&trans[v][c]==x) {

        trans[v][c]=y;

        v=slink[v];

    }

    minlen[y]=maxlen[slink[y]]+;

    return z;

}

void count() {

    queue<int> q;

    for( int i=; i <=tot; i++ )if( !ind[i] ) {

            q.push(i);

        }

    while( !q.empty() ) {

        int u = q.front();

        q.pop();

        if(blue[u] ) edpts[u]++;

        edpts[ slink[u]] += edpts[u];

        if( !--ind[slink[u]] ) q.push(slink[u]);

    }

}

int main() {

    int i;

    scanf("%s",str);

    int len=strlen(str),pre=;

    tot=;

    for(i=; i<len; i++) {

        pre=add(str[i],pre);

    }

    count();

    for(i=;i<=tot;i++){

        ans[maxlen[i]] = max(ans[maxlen[i]],edpts[i]);

    }

    for(i=tot;i>;i--){

        ans[i]=max(ans[i],ans[i+]);

    }

    for(i=;i<=len;i++)printf("%d\n",ans[i]);

    return ;

}