SPOJ - NSUBSTR Substrings 后缀自动机

题意：令F(x)为长度为x的子串出现的最多次数，求所有F(x)

方法：构建后缀自动机，并对每个节点len做计数排序，然后从大到小进行拓补更新。注意到小的子串为大的子串所包含。所以有dp[i]=max(dp[i],dp[i+1])。PS:一个节点的深度所代表的串为root到该节点最远路程所经过的点串，即s(pos-len+1,pos)这一串。pos代表其在字符串中对应的相应位置。也就是说，我们可以通过计算s(pos-len+1,pos)这一子串所出现的次数，同时更新其父节点所出现的次数来得到所有s(pos-len+1,pos)所出现的次数。

代码：

#include <stdio.h>
#include <stdlib.h>
#include <vector>
#include <algorithm>
#include <iostream>
#include <math.h>
#include <string.h>
using namespace std;
typedef long long ll;
const int Mod=1e9+7;
const int N=250003;
int n;
struct SAM_Node{
    SAM_Node *fa,*next[26];
    int len;
    int id,pos;
    SAM_Node(){}
    SAM_Node(int _len){
        fa=0;len=_len;
        memset(next,0,sizeof(next));
    }
};
SAM_Node SAM_node[N<<1],*SAM_root,*SAM_last,*b[N<<1];
int SAM_size;
SAM_Node *newSAM_Node(int len){
    SAM_node[SAM_size]=SAM_Node(len);
    SAM_node[SAM_size].id=0;
	SAM_node[SAM_size].pos=SAM_size;
    return &SAM_node[SAM_size++];
}
SAM_Node *newSAM_Node(SAM_Node *p){
    SAM_node[SAM_size]=*p;
	SAM_node[SAM_size].pos=SAM_size;
    return &SAM_node[SAM_size++];
}
void SAM_init(){
    SAM_size=0;
    SAM_root=SAM_last=newSAM_Node(0);
    SAM_node[0].pos=0;
}
void SAM_add(int x,int len){
    SAM_Node *p=SAM_last,*np=newSAM_Node(p->len+1);    
    SAM_last=np;
    for(;p&&!p->next[x];p=p->fa)
		p->next[x]=np;
    if(!p){
        np->fa=SAM_root;
        return;
    }
    SAM_Node *q=p->next[x];
    if(q->len==p->len+1){
        np->fa=q;
        return;
    }
    //puts("3");
    SAM_Node *nq=newSAM_Node(q);
    nq->len=p->len+1;
    q->fa=nq;
    np->fa=nq;
    for(;p&&p->next[x]==q;p=p->fa)p->next[x]=nq;
}
char a[N];
int dp[N];
int rank[N<<1],cnt[N];
void Sort(){
    memset(cnt,0,sizeof(cnt));
    for(int i=0;i<SAM_size;i++)
        cnt[SAM_node[i].len]++;
    for(int i=1;i<=n;i++)
        cnt[i]+=cnt[i-1];
    for(int i=0;i<SAM_size;i++){
        rank[cnt[SAM_node[i].len]-1]=i;
        b[cnt[SAM_node[i].len]-1]=&SAM_node[i];
        cnt[SAM_node[i].len]--;
    }

}
void out(){
	for(int i=0;i<SAM_size;i++)
		printf("%d,",SAM_node[i].id);
	printf("\n");
}
void work(){
    int T;
    while(~scanf("%s",a+1)){
        n=strlen(a+1);
        SAM_init();
        for(int i=1;i<=n;i++)
            SAM_add(a[i]-'a',i);
        memset(dp,0,sizeof(dp));
        Sort();
		//out();
        SAM_Node *Cur=SAM_root;
        for(int i=1;i<=n;i++){
            Cur=Cur->next[a[i]-'a'];
            Cur->id++;
        }
		//out();
        for(int i=SAM_size-1;i>0;i--){
            SAM_Node *temp=&SAM_node[rank[i]];			
            //printf("(%d,%d)\n",temp->id,b[i]->id);
            int tt=temp->len;
            dp[tt]=max(dp[tt],temp->id);
            if(temp->fa)temp->fa->id+=temp->id;
			//printf("(%d,%d)",temp->pos,temp->fa->pos);
			//out();
//            int tt=b[i]->len;
//            dp[tt]=max(dp[tt],b[i]->id);
//            if(b[i]->fa)b[i]->fa->id+=b[i]->id;
        }
		//out();
        for(int i=n-1;i>=1;i--)
            dp[i]=max(dp[i],dp[i+1]);
        for(int i=1;i<=n;i++)
            printf("%d\n",dp[i]);

    }
}
int main(){
    freopen("data_in.txt","r",stdin);
    //freopen("C_out.txt","w",stdout);
    work();
    return 0;
}