后缀数组（一）（模板题poj2774）

<pre name="code" class="cpp"><span style="font-size:18px;">#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int maxn=200010;
char ss[maxn];
int sa[maxn],t[maxn],t2[maxn],c[maxn],n,s[maxn];
int Rank[maxn],height[maxn];

void get_sa(int m)
{
   int i,*x=t,*y=t2;
   for(i=0;i<m;i++) c[i]=0;
   for(i=0;i<n;i++) c[x[i]=s[i]]++;
   for(i=1;i<m;i++) c[i]+=c[i-1];
   for(i=n-1;i>=0;i--) sa[--c[x[i]]]=i;
   for(int k=1;k<=n;k<<=1)
   {
      int p=0;
      for(i=n-k;i<n;i++) y[p++]=i;
      for(i=0;i<n;i++) if(sa[i]>=k) y[p++]=sa[i]-k;
      for(i=0;i<m;i++) c[i]=0;
      for(i=0;i<n;i++) c[x[y[i]]]++;
      for(i=1;i<m;i++) c[i]+=c[i-1];
      for(i=n-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i];

      swap(x,y);
      p=1;x[sa[0]]=0;
      for(i=1;i<n;i++)
         x[sa[i]]=y[sa[i-1]]==y[sa[i]] && y[sa[i-1]+k]==y[sa[i]+k]?p-1:p++;
      if(p>=n) break;
      m=p;
   }
}

void get_height()
{
   int i,j,k=0;
   for(i=1;i<n;i++) Rank[sa[i]]=i;
   for(i=0;i<n-1;i++)
   {
      if(k)k--;
      j=sa[Rank[i]-1];
      while(s[i+k]==s[j+k]) k++;
      height[Rank[i]]=k;
   }
}

int main()
{
   while(scanf("%s",ss)!=EOF)
   {
      int len=strlen(ss);
      ss[len]='|';
      scanf("%s",ss+len+1);
      n=strlen(ss);
      //memset(s,-1,sizeof(s));
      for(int i=0;i<n;i++) s[i]=ss[i]-'a';
      s[n++]=0;
      get_sa(30);
      get_height();
      int maxn=0;
      for(int i=2;i<n-1;i++)
      {
         if((sa[i]>=len && sa[i-1]<len && sa[i-1]>=0)||(sa[i]<len && sa[i-1]>=len && sa[i]>=0))
            maxn=max(maxn,height[i]);
      }
      printf("%d\n",maxn);
   }
   return 0;
}</span>