POJ - 2774 Long Long Message (后缀数组/后缀自动机模板题)

后缀数组：

 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #include<vector>
 using namespace std;
 typedef long long ll;
 const int N=2e5+,mod=;
 char s[N];
 int sa[N],buf1[N],buf2[N],c[N],n,n1,n2,rnk[N],ht[N],ST[N][],Log[N];
 void Sort(int* x,int* y,int m) {
     for(int i=; i<m; ++i)c[i]=;
     for(int i=; i<n; ++i)++c[x[i]];
     for(int i=; i<m; ++i)c[i]+=c[i-];
     for(int i=n-; i>=; --i)sa[--c[x[y[i]]]]=y[i];
 }
 void da(char* s,int n,int m=) {
     int *x=buf1,*y=buf2;
     x[n]=y[n]=-;
     for(int i=; i<n; ++i)x[i]=s[i],y[i]=i;
     Sort(x,y,m);
     for(int k=; k<n; k<<=) {
         int p=;
         for(int i=n-k; i<n; ++i)y[p++]=i;
         for(int i=; i<n; ++i)if(sa[i]>=k)y[p++]=sa[i]-k;
         Sort(x,y,m),p=,y[sa[]]=;
         for(int i=; i<n; ++i)y[sa[i]]=x[sa[i-]]==x[sa[i]]&&x[sa[i-]+k]==x[sa[i]+k]?p-:p++;
         if(p==n)break;
         swap(x,y),m=p;
     }
 }
 void getht() {
     for(int i=; i<n; ++i)rnk[sa[i]]=i;
     ht[]=;
     for(int i=,k=; i<n; ++i) {
         if(k)--k;
         if(!rnk[i])continue;
         for(; s[i+k]==s[sa[rnk[i]-]+k]; ++k);
         ht[rnk[i]]=k;
     }
 }
 void initST() {
     for(int i=; i<n; ++i)ST[i][]=ht[i];
     for(int j=; (<<j)<=n; ++j)
         for(int i=; i+(<<j)-<n; ++i)
             ST[i][j]=min(ST[i][j-],ST[i+(<<(j-))][j-]);
 }
 int lcp(int l,int r) {
     if(l==r)return n-sa[l];
     if(l>r)swap(l,r);
     l++;
     int k=Log[r-l+];
     return min(ST[l][k],ST[r-(<<k)+][k]);
 }
 int main() {
     Log[]=-;
     for(int i=; i<N; ++i)Log[i]=Log[i>>]+;
     scanf("%s",s),n1=strlen(s),s[n1]='|',scanf("%s",s+n1+),n=strlen(s);
     da(s,n),getht(),initST();
     vector<int> v1,v2;
     for(int i=; i<n1; ++i)v1.push_back(rnk[i]);
     for(int i=n1+; i<n; ++i)v2.push_back(rnk[i]);
     sort(v1.begin(),v1.end());
     sort(v2.begin(),v2.end());
     int ans=;
     for(int i=,j=; i<v1.size(); ++i) {
         for(; j<v2.size()&&v2[j]<v1[i]; ++j);
         if(j->=)ans=max(ans,lcp(v1[i],v2[j-]));
         if(j<v2.size())ans=max(ans,lcp(v1[i],v2[j]));
     }
     printf("%d\n",ans);
     return ;
 }

后缀自动机：

 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 using namespace std;
 typedef long long ll;
 const int N=2e5+,M=,mod=;
 char s[N];
 int n,fa[N],go[N][M],mxl[N],last,tot;
 int newnode(int l) {int u=++tot; mxl[u]=l,memset(go[u],,sizeof go[u]); return u;}
 void add(int ch) {
     int p=last,np=last=newnode(mxl[p]+);
     for(; ~p&&!go[p][ch]; p=fa[p])go[p][ch]=np;
     if(!~p)fa[np]=;
     else {
         int q=go[p][ch];
         if(mxl[q]==mxl[p]+)fa[np]=q;
         else {
             int nq=newnode(mxl[p]+);
             memcpy(go[nq],go[q],sizeof go[q]);
             fa[nq]=fa[q],fa[q]=fa[np]=nq;
             for(; ~p&&go[p][ch]==q; p=fa[p])go[p][ch]=nq;
         }
     }
 }
 int main() {
     fa[]=-;
     scanf("%s",s),n=strlen(s);
     for(int i=; i<n; ++i)add(s[i]-'a');
     scanf("%s",s),n=strlen(s);
     int ans=;
     for(int i=,u=,l=; i<n; ++i) {
         int ch=s[i]-'a';
         if(!go[u][ch]) {
             for(; u&&!go[u][ch]; u=fa[u]);
             l=mxl[u];
         }
         if(go[u][ch])u=go[u][ch],++l;
         ans=max(ans,l);
     }
     printf("%d\n",ans);
     return ;
 }