[后缀自动机] POJ 1743 Musical Theme

又一道男人八题哦
先差分 然后就是找出现两次的不相交的最大子串
考虑对于每个点 记录它的right集的最值 用差值和maxl的最小值更新答案

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#define cl(x) memset(x,0,sizeof(x))
using namespace std;

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin); if (p1==p2) return EOF; }
  return *p1++;
}

inline void read(int &x){
  char c=nc(),b=1;
  for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
  for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}

const int N=20005;
const int K=87<<1|1;

struct state{
  int len,link,next[K];
}st[N<<1];
int ncnt,last;

int maxv[N<<1],minv[N<<1];

inline void Extend(int c){
  int cur=++ncnt,p; 
  st[cur].len=st[last].len+1; maxv[cur]=minv[cur]=st[cur].len;
  for (p=last;p!=-1 && !st[p].next[c];p=st[p].link)
    st[p].next[c]=cur;
  if (p==-1)
    st[cur].link=0;
  else{
    int q=st[p].next[c];
    if (st[q].len==st[p].len+1)
      st[cur].link=q;
    else{
      int nq=++ncnt; maxv[nq]=-1<<30; minv[nq]=1<<30;
      st[nq].len=st[p].len+1;
      st[nq].link=st[q].link;
      for (int i=0;i<K;i++) st[nq].next[i]=st[q].next[i];
      for (;p!=-1 && st[p].next[c]==q;p=st[p].link)
    st[p].next[c]=nq;
      st[q].link=st[cur].link=nq;
    }
  }
  last=cur;
}

int tmp[N<<1],sum[N<<1];

inline void Sort(){
  for (int i=0;i<=ncnt;i++) sum[st[i].len]++;
  for (int i=1;i<=st[last].len;i++) sum[i]+=sum[i-1];
  for (int i=0;i<=ncnt;i++) tmp[sum[st[i].len]--]=i;
}

int n,a[N];

int main(){
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  while (1){
    read(n); if (!n) break;
    for (int i=1;i<=n;i++) read(a[i]);
    for (int i=1;i<=n;i++) a[i]-=a[i+1];
    for (int i=1;i<n;i++) a[i]+=87;
    st[0].link=-1;
    for (int i=1;i<n;i++) Extend(a[i]);
    Sort();
    for (int i=ncnt+1;i;i--)
      if (tmp[i]){
    maxv[st[tmp[i]].link]=max(maxv[st[tmp[i]].link],maxv[tmp[i]]);
    minv[st[tmp[i]].link]=min(minv[st[tmp[i]].link],minv[tmp[i]]);
      }
    int Ans=0;
    for (int i=1;i<=ncnt;i++){
      int ret=min(maxv[i]-minv[i],st[i].len);
      if (ret>=4) Ans=max(Ans,ret+1);
    }
    printf("%d\n",Ans);
    ncnt=last=0; cl(st); cl(tmp); cl(sum);
  }
  return 0;
}