Luogu P2463 [SDOI2008]Sandy的卡片

题意：求差分串的最长公共子串。

(SA) 二分答案 (k) ，即求是否存在一段连续的 (ht[]) ，且其中包括了所有串的后缀。

#include<iostream>
#include<cstdio>
#include<cstring>
#define R register int
using namespace std;
namespace Luitaryi {
inline int g() { R x=0,f=1;
  register char s; while(!isdigit(s=getchar())) f=s=='-'?-1:f;
  do x=x*10+(s^48); while(isdigit(s=getchar())); return x*f;
} const int N=1010,M=110010,L=1865;
int T,n,k,len[N],a[M],s[M],sa[M],rk[M],ht[M],c[M],x[M],y[M],p[M];
bool vis[M];
inline void pre() {
  R m=L+1000;
  for(R i=1;i<=n;++i) ++c[x[i]=s[i]];
  for(R i=1;i<=m;++i) c[i]+=c[i-1];
  for(R i=n;i;--i) sa[c[x[i]]--]=i;
  for(R t=1,top=0;top<n;m=top,t<<=1) {
    top=0;
    for(R i=n-t+1;i<=n;++i) y[++top]=i;
    for(R i=1;i<=n;++i) if(sa[i]>t) y[++top]=sa[i]-t;
    memset(c,0,m+1<<2);
    for(R i=1;i<=n;++i) ++c[x[i]];
    for(R i=1;i<=m;++i) c[i]+=c[i-1];
    for(R i=n;i;--i) sa[c[x[y[i]]]--]=y[i];
    swap(x,y),x[sa[1]]=top=1;
    for(R i=2;i<=n;++i) 
      x[sa[i]]=(y[sa[i]]==y[sa[i-1]]&&y[sa[i]+t]==y[sa[i-1]+t])?top:++top;
  } for(R i=1;i<=n;++i) rk[sa[i]]=i;
  R k=0;
  for(R i=1,j;i<=n;++i) {
    if(k) --k; j=sa[rk[i]-1];
    while(s[i+k]==s[j+k]) ++k;
    ht[rk[i]]=k; 
  }
}
inline bool calc(int l,int r) {
  if(r-l+1<T) return false;
  R ret=0;
  memset(vis,0,sizeof vis);
  for(R i=l;i<=r;++i) {
    if(!vis[p[sa[i]]]) 
      vis[p[sa[i]]]=true,++ret;
  } return ret==T;
}
inline bool ck(int x) {
	R l=1,r=1;
  while(l<=n) {
    while(l<=n && ht[l]<x) ++l;
    if(l>n) break;
    r=l;
    while(r<n && ht[r+1]>=x) ++r;
    if(calc(l-1,r)) return true;
    l=r+1,r=l;
  } return false;
}
inline void main() {
	T=g(); 
  for(R i=1,l;i<=T;++i) { l=g();
    for(R j=1;j<=l;++j) a[j]=g();
    for(R j=2;j<=l;++j)
      s[++n]=a[j]-a[j-1],p[n]=i;
    s[++n]=L+i,p[n]=i;
  } pre();
  R l=0,r=n,md;
  while(l<r) {
    md=(l+r+1)>>1;
    if(ck(md)) l=md;
    else r=md-1;
  } printf("%d
",l+1);
}
} signed main() {Luitaryi::main(); return 0;}

---

(SAM) 板子？

---

2020.01.09