SPOJ 8222. Substrings

给一个字符串S，令F(x)表示S的长度为x的子串最多在S中出现 了多少次。求F(1)..F(Length(S))
Length(S) <= 250000
我们构造S的SAM，那么对于一个节点s，它的长度范围是 [Min(s),Max(s)]，同时他的出现次数是|Right(s)|。那么我们用
|Right(s)|去更新F(Max(s))的值。 同时最后从大到小依次用F(i)去更新F(i-1)即可。

SAM中root到每个点经过的路径和原串中的子串一一对应。

每个点到接收态的路径数等于这个节点对应的子串在S中出现了多少次。

从最后一个节点往上走，DFS求出每个子串出现了多少次。

因为出现了4次的子串也出现了3,2,1次。。用dp[i]更新dp[i-1]。。。长度大的更新长度小的。。

Substrings

Time Limit: 1000MS

Memory Limit: Unknown

64bit IO Format: %lld & %llu

[Submit]   [Go Back]   [Status]  

Description

You are given a string S which consists of 250000 lowercase latin letters at most. We define F(x) as the maximal number of times that some string with length x appears in S. For example for string 'ababa' F(3) will be 2 because there is a string 'aba' that occurs twice. Your task is to output F(i) for every i so that 1<=i<=|S|.

Input

String S consists of at most 250000 lowercase latin letters.

Output

Output |S| lines. On the i-th line output F(i).

Example

Input:
ababa

Output:
3
2
2
1
1

#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
using namespace std;
#define ll long long
#define prt(k) cout<<#k"="<<k<<" ";

/**Suffix Automaton*/
const int Char=26,N=600022;
struct SAM_Node
{
    SAM_Node *fa,*ch[Char];
    int len,id,pos;
    SAM_Node() {}
    SAM_Node(int _len)
    {
        fa=0; len=_len;
        memset(ch,0,sizeof ch);
    }
}pool[N],*root,*last,*tail;   ///tail
int SAM_size;
#define node SAM_Node
SAM_Node *newnode(int len)
{
    pool[SAM_size]=node(len);
    pool[SAM_size].id=SAM_size;
    return &pool[SAM_size++];
}
node *newnode(node *p)
{
    pool[SAM_size]=*p;
    pool[SAM_size].id=SAM_size;
    return &pool[SAM_size++];
}
void SAM_init()
{
    SAM_size=0;
    root=last=newnode(0);
    pool[0].pos=0;
}
void add(int x,int len)
{
    node *p=last,*np=newnode(p->len+1);
    np->pos=len; last=np;
    for(;p&&!p->ch[x];p=p->fa) p->ch[x]=np;
    if(!p) { np->fa=root;return ; }
    node *q=p->ch[x];
    if(q->len==p->len+1) { np->fa=q;return; }
    node *nq=newnode(q);
    nq->len=p->len+1;
    q->fa=nq; np->fa=nq;
    for(;p&&p->ch[x]==q;p=p->fa) p->ch[x]=nq;
}
void SAM_build(char *s)
{
    SAM_init();
    int len=strlen(s);
    for(int i=0;i<len;i++) add(s[i]-'a',i+1);
}
void Max(int &a,int b) { a=max(a,b); }
void Min(int &a,int b) { a=min(a,b); }
bool vis[N];
int cnt[N];
int dfs(node* p)
{
    if(vis[p->id]) return cnt[p->id];
    vis[p->id]=1;
    cnt[p->id]=0;
    for(int i=0;i<26;i++)
        if(p->ch[i]) cnt[p->id]+=dfs(p->ch[i]);
    return cnt[p->id];
}
void getcnt()
{
    memset(cnt,0,sizeof cnt);
    memset(vis,0,sizeof vis);
    node* p=last;
    vis[p->id]=1;
    for(;p;p=p->fa) dfs(p),cnt[p->id]++;
}
int dp[N];
char str[N];
int main()
{
    while(scanf("%s",str)==1)
    {
        SAM_build(str);
        memset(dp,0,sizeof dp);
        getcnt();
        for(int i=0;i<SAM_size;i++)
            Max(dp[pool[i].len],cnt[pool[i].id]);
        int n=strlen(str);
        for(int i=n;i>=2;i--) Max(dp[i-1],dp[i]);
        for(int i=1;i<=n;i++)
            printf("%d\n",dp[i]);
    }

}