1. 程式人生 > >SPOJ8222 NSUBSTR - Substrings

SPOJ8222 NSUBSTR - Substrings

傳送門[洛谷]

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\n",f[i]);
	return 0;
}