SPOJ8222 NSUBSTR

传送门[洛谷]

SAM模板题 可惜我不会

本来想大力上线段树的

然后发现题解的做法非常精妙

f[i]=max(f[i],f[i+1])

因为 很显然长度为i的字符串出现长度一定>=长度为i+1的字符串

[作为子串出现就可以取到=]

那么直接记录最后倒着更新一遍就可以啦

附代码。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define inf 20021225
#define ll long long
#define mxn 250010
using namespace std;

struct node{int ch[26],fa,len,sz;}t[mxn*4];
struct edge{int to,lt;}e[mxn*4];
char ch[mxn];int f[mxn],in[mxn*4],cnt,poi,lt,rt;
void add(int x,int y){e[++cnt].to=y;e[cnt].lt=in[x];in[x]=cnt;}
int id(char c){return c-'a';}
void insert(int c)
{
	int p=lt,np; np=lt=++poi; t[np].sz=1; t[np].len=t[p].len+1;
	for(;p&&!t[p].ch[c];p=t[p].fa)	t[p].ch[c]=np;
	if(!p){t[np].fa=rt;return;}
	int q=t[p].ch[c];
	if(t[q].len==t[p].len+1){t[np].fa=q;return;}
	int nq=++poi; t[nq].len=t[p].len+1;//稳定转移 
	memcpy(t[nq].ch,t[q].ch,sizeof(t[q].ch));
	t[nq].fa=t[q].fa; t[q].fa=t[np].fa=nq;
	for(;p&&t[p].ch[c]==q;p=t[p].fa)	t[p].ch[c]=nq;
}
void build(){for(int i=2;i<=poi;i++)add(t[i].fa,i);}
void dfs(int x)
{
	for(int i=in[x];i;i=e[i].lt){dfs(e[i].to);t[x].sz+=t[e[i].to].sz;}
	f[t[x].len]=max(f[t[x].len],t[x].sz);
}
int main()
{
	scanf("%s",ch+1);int n=strlen(ch+1);
	rt=lt=++poi; 
	for(int i=1;i<=n;i++)	insert(id(ch[i]));
	build();dfs(rt);
	for(int i=n;i;i--)	f[i]=max(f[i],f[i+1]);
	for(int i=1;i<=n;i++)	printf("%d
",f[i]);
	return 0;
}
原文地址:https://www.cnblogs.com/hanyuweining/p/10321919.html