题目链接实现10S 题目链接实现3S(本文代码cf800ms) 后缀数组论文
题目大意是给你n个字符串S1~Sn,然后要你找出一个S1中的最短子串t,而且t不是S2~Sn的子串。
- 首先我们需要字典序最小,那么就可以想到把全部才串连在一起跑一下后缀数组。
- 然后遍历找到起点在S1的suffix(i),易知与suffix(i)前缀最长的是suffix(i-1)与suffix(i+1)。
- 然后查询一下相互之间的最长公共前缀长度加上suffix(i)的起点是否超出S1的长度。
- 若不超出则尝试更新最短答案。
#include<algorithm>
#include<iostream>
#include<math.h>
#include<set>
#include<map>
#include<cstring>
#include<string>
#define qcin; ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
using namespace std;
typedef long long ll;
const int maxn = 300010;
typedef int arr[maxn];
struct ST{
int mm[maxn],dp[maxn][22],k;
const int& (*fun[2])(const int&,const int&)={max,min};
void init(int n,int b[],int fl){
mm[0]=-1;
k=fl;
for(int i=1;i<=n;++i){
mm[i]=mm[i-1]+((i&(i-1))==0?1:0);
dp[i][0]=b[i];
}
for(int j=1;j<=mm[n];j++){
for(int i=1;i+(1<<j)-1<=n;++i){
dp[i][j]=fun[k](dp[i][j-1],dp[i+(1<<(j-1))][j-1]);
}
}
}
inline int qry(int l,int r){
int kk=mm[r-l+1];
return fun[k](dp[l][kk],dp[r-(1<<kk)+1][kk]);
}
};
struct SA{
//h[i]=suffix(sa[i-1])和 suffix(sa[i])的最长公共前缀的长度
//sa[i]表示字典序第i小的后缀起点
//rk[i]表示起点i的后缀的字典序排位
int sa[maxn],rk[maxn],h[maxn];
ST st;
void init(string s,int n){
static arr c,v,q;
int m=127;
for(int i=0;i<m;++i)c[i]=0;
for(int i=0;i<n;++i)++c[s[i]];
for(int i=1;i<m;++i)c[i]+=c[i-1];
for(int i=n-1;~i;--i)sa[--c[s[i]]]=i;
m=rk[sa[0]]=0;
for(int i=1;i<n;++i)rk[sa[i]]=s[sa[i]]!=s[sa[i-1]]?++m:m;
for(int j=1;;j<<=1){
if(++m==n)break;
for(int i=0;i<j;++i)v[i]=n-j+i;
for(int i=0,k=j;i<n;++i)if(sa[i]>=j)v[k++]=sa[i]-j;
for(int i=0;i<n;++i)q[i]=rk[v[i]];
for(int i=0;i<m;++i)c[i]=0;
for(int i=0;i<n;++i)++c[q[i]];
for(int i=1;i<m;++i)c[i]+=c[i-1];
for(int i=n-1;~i;--i)sa[--c[q[i]]]=v[i],v[i]=rk[i];
m=rk[sa[0]]=0;
for(int i=1;i<n;++i)rk[sa[i]]=v[sa[i]]!=v[sa[i-1]]||v[sa[i]+j]!=v[sa[i-1]+j]?++m:m;
}
int k=0;
for(int i=0;i<n-1;++i){
if(k)--k;
if(!rk[i])continue;
while(s[i+k]==s[sa[rk[i]-1]+k])++k;
h[rk[i]]=k;
}
st.init(n-1,h,1);
}
}sa;
int t,n,m,len;
char sp='*';
string s,tmp;
int pre[maxn],nxt[maxn];
void sol(int x){
memset(pre,-1,(n+1)*sizeof(int));
memset(nxt,-1,(n+1)*sizeof(int));
int k=-1,ans=-1,fi=-1;
for(int i=0;i<n;++i){
if(sa.sa[i]>=len)k=i;
else pre[i]=k;
}
k=-1;
for(int i=n-1;~i;--i){
if(sa.sa[i]>=len)k=i;
else nxt[i]=k;
}
for(int i=0;i<n;++i){
int mx=0;
if(sa.sa[i]<len){
if(~pre[i]){
mx=max(mx,sa.st.qry(pre[i]+1,i));
}
if(~nxt[i]){
mx=max(mx,sa.st.qry(i+1,nxt[i]));
}
if((ans>mx+1||ans==-1)&&sa.sa[i]+mx+1<=len){
ans=mx+1;
fi=sa.sa[i];
}
}
}
if(~ans){
cout<<"Case #"<<x<<": ";
for(int i=fi;i<fi+ans;++i){
cout<<s[i];
}
cout<<endl;
return;
}
cout<<"Case #"<<x<<": Impossible"<<endl;
}
int main(){
qcin;
cin>>t;
for(int tt=1;tt<=t;++tt){
cin>>n;s="";s.clear();
for(int i=0;i<n;++i){
cin>>tmp;
if(i)s+=sp;
else len=tmp.length();
s+=tmp;
}
n=s.length();
sa.init(s,s.length());
sol(tt);
}
}