SPOJ8222

传送门

输入一个字符串，要求输出长度为x的子串在字符串中的最多的出现位置。

后缀自动机，由于后缀自动机的每一个节点s的right集大小都表示子串最长为step s的数量，

（而且长度在min(s)到max(s)之间的字符串出现次数必定相同，所以每次只需要更新step s，最后从大到小顺次cmax即可）

所以只需要按topsort序求出每个节点的权值并更新Max[step[s]]即可。

最后需要注意的一点是，在构造的时候，np节点的权值赋为1，nq节点权值赋为0，不要求成了当前子树的size。

代码简单易懂。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<string>
#include<cstring>
#include<cmath>
#include<ctime>
#include<algorithm>
#include<iomanip>
#include<set>
#include<map>
#include<queue>
using namespace std;
#define mem1(i,j) memset(i,j,sizeof(i))
#define mem2(i,j) memcpy(i,j,sizeof(i))
#define LL long long
#define up(i,j,n) for(int i=(j);i<=(n);i++)
#define FILE "dealing"
#define poi vec
#define eps 1e-10
#define db double 
const int maxn=500010,inf=1000000000,mod=1000000007;
int read(){
	int x=0,f=1,ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<1)+(x<<3)+ch-'0',ch=getchar();}
	return f*x;
}
bool cmax(int& a,int b){return a<b?a=b,true:false;}
bool cmin(int& a,int b){return a>b?a=b,true:false;}
int que[maxn],Max[maxn];
char s[maxn];
namespace SAM{
	int cnt=1,now=1;
	int pre[maxn],c[maxn][26],len[maxn],ct[maxn],rank[maxn];
	void expend(int x){
		int np,p,q,nq;
		p=now;np=++cnt;len[np]=len[p]+1;now=np;que[np]=1;
		while(p&&!c[p][x])c[p][x]=np,p=pre[p];
		if(!p)pre[np]=1;
		else {
			q=c[p][x];
			if(len[q]==len[p]+1)pre[np]=q;
			else {
				len[nq=++cnt]=len[p]+1;que[nq]=0;
				mem2(c[nq],c[q]);
				pre[nq]=pre[q];
				pre[q]=pre[np]=nq;
				while(p&&c[p][x]==q)c[p][x]=nq,p=pre[p];
			}
		}
	}
	void getsort(){
		up(i,1,cnt)ct[len[i]]++;
		up(i,1,cnt)ct[i]+=ct[i-1];
		for(int i=cnt;i>=1;i--)rank[ct[len[i]]--]=i;//从后往前即为topsort序
	}
	void build(char* s){
		cnt=1,now=1;int n=strlen(s+1);
		up(i,1,n)expend(s[i]-'a');
	}
};
int main(){
	freopen(FILE".in","r",stdin);
	freopen(FILE".out","w",stdout);
	scanf("%s",s+1);int n=strlen(s+1);
	SAM::build(s);
	SAM::getsort();
	int* t=SAM::rank;
	for(int i=SAM::cnt;i>=1;i--){
		que[SAM::pre[t[i]]]+=que[t[i]];
		cmax(Max[SAM::len[t[i]]],que[t[i]]);
	}
	for(int i=n;i>=1;i--)
		cmax(Max[i],Max[i+1]);//这一步貌似可省，没试过，但感觉没有这个不行啊，哪个大牛指点一下？
	up(i,1,n)printf("%d
",Max[i]);
	return 0;
}