1. 程式人生 > >【SPOJ】Substrings(後綴自動機)

【SPOJ】Substrings(後綴自動機)

tps 大小 substring 子串 ges 自動 好處 can printf

【SPOJ】Substrings(後綴自動機)

題面

Vjudge
題意:給定一個長度為\(len\)的串,求出長度為1~len的子串中,出現最多的出現了多少次

題解

出現次數很好處理,就是\(right/endpos\)集合的大小
那麽,直接構建\(SAM\)
求出每個位置的\(right\)集合大小
直接更新每個節點的\(longest\)就行了
最後短的可以由長的更新過來就好

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm> #include<set> #include<map> #include<vector> #include<queue> using namespace std; #define ll long long #define RG register #define MAX 2001000 inline int read() { RG int x=0,t=1;RG char ch=getchar(); while((ch<'0'||ch>'9')&&ch!='-'
)ch=getchar(); if(ch=='-')t=-1,ch=getchar(); while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar(); return x*t; } char ch[MAX]; struct Node { int son[26]; int ff,len; }t[MAX<<1]; int size[MAX]; int tot=1,last=1,c[MAX],a[MAX],ans[MAX]; void extend(int
c) { int p=last,np=++tot;last=np; t[np].len=t[p].len+1; while(p&&!t[p].son[c])t[p].son[c]=np,p=t[p].ff; if(!p)t[np].ff=1; else { int q=t[p].son[c]; if(t[q].len==t[p].len+1)t[np].ff=q; else { int nq=++tot; t[nq]=t[q]; t[nq].len=t[p].len+1; t[q].ff=t[np].ff=nq; while(p&&t[p].son[c]==q)t[p].son[c]=nq,p=t[p].ff; } } size[np]=1; } int main() { scanf("%s",ch+1); int l=strlen(ch+1); for(int i=1;i<=l;++i)extend(ch[i]-97); for(int i=1;i<=tot;++i)c[t[i].len]++; for(int i=1;i<=tot;++i)c[i]+=c[i-1]; for(int i=1;i<=tot;++i)a[c[t[i].len]--]=i; for(int i=tot;i;--i) { int u=a[i]; size[t[u].ff]+=size[u]; ans[t[u].len]=max(ans[t[u].len],size[u]); } for(int i=l-1;i;--i)ans[i]=max(ans[i],ans[i+1]); for(int i=1;i<=l;++i)printf("%d\n",ans[i]); return 0; }

【SPOJ】Substrings(後綴自動機)