SPOJ_NSUBSTR

题目意思是给你一个字符串,f[x]是长度为x的子串中,出现个数最多的那个串的出现次数。

给出原串,依次输出f[1],f[2],……。

后缀自动机。对于某一个状态,right[]值的大小就是出现的次数,而且是对于长为step[]的子串的出现次数。

因为小于step值的串在前面已经加了,在pre指针线上面的状态要把当前数量累加上去即可。

不过注意这里要先拓扑排序,这样才能保证更新的正确性。恩,很实用很不错的一种拓扑排序方式。新技能get。

召唤代码君:

#include <iostream>
#include <cstdio>
#define maxn 500500
using namespace std;

int next[maxn][26],pre[maxn],step[maxn],dp[maxn],g[maxn],Q[maxn],cnt[maxn];
int N,last,p,q,np,nq;
char s[maxn];

void insert(int x,int m)
{
	p=last,np=++N,step[np]=m,last=np,g[np]++;
	while (p!=-1 && next[p][x]==0)
        next[p][x]=np,p=pre[p];
	if (p==-1) return ;
	q=next[p][x];
	if (step[q]==step[p]+1)
        { pre[np]=q; return; }
	nq=++N,step[nq]=step[p]+1,pre[nq]=pre[q];
	for (int i=0; i<26; i++)
        next[nq][i]=next[q][i];
      pre[np]=pre[q]=nq;
	while (p!=-1 && next[p][x]==q)
        next[p][x]=nq,p=pre[p];
}

int main()
{
	pre[0]=-1;
	scanf("%s",s);
	for (int i=0; s[i]; i++) insert(s[i]-‘a‘,i+1);
	for (int i=1; i<=N; i++) cnt[step[i]]++;
	for (int i=1; i<=N; i++) cnt[i]+=cnt[i-1];
	for (int i=1; i<=N; i++) Q[cnt[step[i]]--]=i;
	for (int i=N; i>=1; i--) dp[step[Q[i]]]=max(dp[step[Q[i]]],g[Q[i]]),g[pre[Q[i]]]+=g[Q[i]];
	for (int i=1; s[i-1]; i++) printf("%d\n",dp[i]);
	return 0;
}

  

SPOJ_NSUBSTR,布布扣,bubuko.com

时间: 2024-11-15 21:52:58

SPOJ_NSUBSTR的相关文章

后缀自动机/后缀树

只是笔记罢了,不要看 关于DAWG: 见紫书P390 把后缀自动机上所有节点都设为接受态就形成DAWG,可以接受一个字符串的所有子串. 一个子串的end-set是它在原串w中出现位置(从1开始编号)的右端点集合. 在DAWG中,end-set相同的子串属于同一个状态. 原因没原因,这应该算定义吧? 任意两个节点的end-set要么不相交,要么是包含关系. 原因:在DAWG上走一步,当前end-set的变化是将原end-set中各个元素+1(要去掉超出字符串长度的元素),然后拆分成1个或多个新en