POJ.2774.Long Long Message(後綴數組 倍增)
阿新 • • 發佈:2018-03-14
前綴 中間 題目 std str 後綴 string pri -m
題目鏈接
\(Description\)
求兩個字符串最長公共子串
\(Solution\)
任何一個子串一定是某個後綴的前綴
可以將兩個字符串拼在一起,中間用一個從未出現過的字符隔開,這樣ht[]的最大值就是答案?
不一定,最大的ht[]可能是由同一個字符串得到的,判一下屬於哪個字符串即可
//3772K 516MS
#include <cstdio>
#include <cstring>
#include <algorithm>
const int N=2e5+10;
int n,l,sa[N],rk[N],ht[N],sa2[N],tm[N];
char s[N];
void Get_SA()
{
int *x=rk,*y=sa2,m=30;
for(int i=1; i<=n; ++i) ++tm[x[i]=s[i]-'a'+1];
for(int i=1; i<=m; ++i) tm[i]+=tm[i-1];
for(int i=n; i; --i) sa[tm[x[i]]--]=i;
for(int p=0,k=1; k<n; m=p,p=0,k<<=1)
{
for(int i=n-k+1; i<=n; ++i) y[++p]=i;
for (int i=1; i<=n; ++i) if(sa[i]>k) y[++p]=sa[i]-k;
for(int i=0; i<=m; ++i) tm[i]=0;
for(int i=1; i<=n; ++i) ++tm[x[i]];
for(int i=1; i<=m; ++i) tm[i]+=tm[i-1];
for(int i=n; i; --i) sa[tm[x[y[i]]]--]=y[i];
std::swap(x,y), p=x[sa[1]]=1;
for(int i=2; i<=n; ++i)
x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+k]==y[sa[i]+k]?p:++p;
if(p>=n) break;
}
for(int i=1; i<=n; ++i) rk[sa[i]]=i;
ht[1]=0;
for(int k=0,p,i=1; i<=n; ++i)
{
if(rk[i]==1) continue;
if(k) --k;
p=sa[rk[i]-1];
while(i+k<=n&&p+k<=n&&s[i+k]==s[p+k]) ++k;
ht[rk[i]]=k;
}
}
int main()
{
scanf("%s",s+1), l=strlen(s+1);
s[l+1]='z'+1;
scanf("%s",s+2+l), n=strlen(s+1);
Get_SA();
int res=0;
for(int i=2; i<=n; ++i)
if((sa[i]<=l)^(sa[i-1]<=l)) res=std::max(res,ht[i]);
printf("%d",res);
return 0;
}
POJ.2774.Long Long Message(後綴數組 倍增)