bzoj3796(后缀数组)

题目描述

Mushroom最近看上了一个漂亮妹纸。他选择一种非常经典的手段来表达自己的心意——写情书。考虑到自己的表达能力，Mushroom决定不手写情书。他从网上找到了两篇极佳的情书，打算选择其中共同的部分。另外，Mushroom还有个一个情敌Ertanis，此人也写了封情书给妹子。

Mushroom不希望自己的情书中完整的出现了情敌的情书。（这样抄袭的事情就暴露了）。

Mushroom把两封情书分别用字符串s1和s2来表示，Ertanis的情书用字符串s3来表示，他要截取的部分用字符串w表示。

需满足：

1、w是s1的子串

2、w是s2的子串

3、s3不是w的子串

4、w的长度应尽可能大

所谓子串是指：在字符串中连续的一段。

题解

我真是有看见字符串想SAM综合征。。

要求匹配的串只有两个，可以把它们连起来求一遍SA。

然后KMP求出c在大串中出现的位置。

维护一个数组f[i]表示从i位置开始至多延长多少个位置保证不包含c。

然后扫两遍，一遍算第一个串对第二个串的贡献，再算第二个串对第一个串的贡献，滚动维护一下min(height)就好了。

又TM开小数组了。、

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#define N 200009
using namespace std;
int n,m,sa[N],rnk[N],tong[N],y[N],height[N],n1,fail[N],f[N],ans;
bool b[N],tag[N];
char s[N],s1[N];
inline int rd(){
    int x=0;char c=getchar();bool f=0;
    while(!isdigit(c)){if(c=='-')f=1;c=getchar();}
    while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
    return f?-x:x;
}
inline void qsort(){
    for(int i=0;i<=m;++i)tong[i]=0;
    for(int i=1;i<=n;++i)tong[rnk[i]]++;
    for(int i=1;i<=m;++i)tong[i]+=tong[i-1];
    for(int i=n;i>=1;--i)sa[tong[rnk[y[i]]]--]=y[i];
}
inline void SA(){
    m=200;
    for(int i=1;i<=n;++i)rnk[i]=s[i],y[i]=i;
    qsort();
    for(int w=1,p=0;p<n;m=p,w<<=1){
        p=0;
        for(int i=n-w+1;i<=n;++i)y[++p]=i;
        for(int i=1;i<=n;++i)if(sa[i]>w)y[++p]=sa[i]-w;
        qsort();swap(rnk,y);rnk[sa[1]]=p=1;
        for(int i=2;i<=n;++i)rnk[sa[i]]=((y[sa[i]]==y[sa[i-1]])&&(y[sa[i]+w]==y[sa[i-1]+w]))?p:++p;
    }    
    for(int i=1;i<=n;++i){
        if(rnk[i]==1)continue;
        int j=max(0,height[rnk[i-1]]-1);
        while(s[i+j]==s[sa[rnk[i]-1]+j])j++;
        height[rnk[i]]=j;
    }   
}
inline void check(int pos,int ln){
    int j=0;
    for(int i=pos;i<=pos+ln-1;++i){
        while(j&&s[i]!=s1[j+1])j=fail[j];
        if(s[i]==s1[j+1])j++;
        if(j==n1){cout<<"gan";return;}
    } 
    cout<<"yeah";
}
int main(){
    scanf("%s%s",s+1,s1+1);
    n=strlen(s+1);n1=strlen(s1+1);
    for(int i=1;i<=n;++i)f[i]=n-i+1;s[++n]='#'; 
    for(int i=1;i<=n1;++i)s[++n]=s1[i],b[n]=1,f[n]=n1-i+1;
    SA();
    scanf("%s",s1+1);n1=strlen(s1+1);
    int j=0;
    for(int i=2;i<=n1;++i){
        while(j&&s1[i]!=s1[j+1])j=fail[j];
        if(s1[i]==s1[j+1])j++;fail[i]=j;
    }
    j=0;
    for(int i=1;i<=n;++i){
        while(j&&s[i]!=s1[j+1])j=fail[j];
        if(s[i]==s1[j+1])j++;
        if(j==n1)tag[i-j+1]=1;
    }
    int xian=n+1;
    for(int i=n;i>=1;--i){
        if(tag[i])xian=min(xian,i+n1-1);
        f[i]=min(f[i],xian-i);
    }
    xian=0;
    for(int i=1;i<=n;++i){
        int pos=sa[i];
        if(!b[pos])xian=min(xian,height[i]),ans=max(ans,min(f[pos],xian));
        else xian=2e9;
    }
    xian=0;
    for(int i=1;i<=n;++i){
        int pos=sa[i];
        if(b[pos])xian=min(xian,height[i]),ans=max(ans,min(f[pos],xian));
        else xian=2e9;
    }
    cout<<ans;
    return 0;
}

猜你喜欢