标签：200005 int Luogu5546 long POI2000 BZOJ2946 sa Sandy define

题目传送门

算法分析

和Sandy的卡片几乎是同一题我们证明了中国OI的POI性。

具体算法分析见Luogu2463/BZOJ4698 [SDOI2008]Sandy的卡片

要注意是否存在完全一样的字符串，此时的分隔符要注意。

代码实现

#include<bits/stdc++.h>
using namespace std;
#define maxn 200005
#define maxm 200005
#define inf 0x3f3f3f3f
#define LL long long
#define mod 1000000007
#define local
template <typename Tp>
void read(Tp &x){
    x=0;int fh=1;char c=getchar();
    while(c>'9'||c<'0'){if(c=='-'){fh=-1;}c=getchar();}
    while(c>='0'&&c<='9'){x=(x<<1)+(x<<3)+(c&15);c=getchar();}x*=fh;
}
int n,m,s[maxn],tot,id[maxn];
int sa[2][maxn],rk[2][maxn],v[maxn],h[maxn];
int st[30][maxn],lg=20,ans;
int get_lcp(int x,int y){
    x++;
    int kk=log(y-x+1)/log(2);
    return min(st[kk][x],st[kk][y-(1<<kk)+1]);
}
int c[maxn],su;
char ss[maxn];
int main(){
    int l,x,xx;
    read(m);
    for(int i=1;i<=m;i++){
        scanf("\n%s",ss);
        l=strlen(ss);
        for(int j=0;j<l;j++){
            s[++tot]=ss[j];id[tot]=i;
        }
        s[++tot]=i;id[tot]=i;
    }
    n=tot;int p=0,q=1;
    for(int i=1;i<=n;i++)v[s[i]]++;
    for(int i=1;i<=128;i++)v[i]+=v[i-1];
    for(int i=1;i<=n;i++)sa[0][v[s[i]]--]=i;
    for(int i=1;i<=n;i++)rk[0][sa[0][i]]=rk[0][sa[0][i-1]]+(s[sa[0][i]]!=s[sa[0][i-1]]);
    for(int k=1;k<=n;k<<=1,swap(p,q)){
        for(int i=1;i<=n;i++)v[rk[p][sa[p][i]]]=i;
        for(int i=n;i;i--)if(sa[p][i]>k)sa[q][v[rk[p][sa[p][i]-k]]--]=sa[p][i]-k;
        for(int i=n-k+1;i<=n;i++)sa[q][v[rk[p][i]]--]=i;
        for(int i=1;i<=n;i++)rk[q][sa[q][i]]=rk[q][sa[q][i-1]]+(rk[p][sa[q][i]]!=rk[p][sa[q][i-1]]||rk[p][sa[q][i]+k]!=rk[p][sa[q][i-1]+k]);
        if(rk[q][sa[q][n]]==n)break;
    }
    for(int i=1,k=0;i<=n;i++){
        if(rk[q][i]==1)continue;
        if(k)k--;
        int j=sa[q][rk[q][i]-1];
        while(s[i+k]==s[j+k])k++;
        h[rk[q][i]]=k;
    }
    for(int i=1;i<=n;i++)st[0][i]=h[i];
    for(int j=1;j<=lg;j++){
        for(int i=1;i+(1<<j)-1<=n;i++){
            st[j][i]=min(st[j-1][i],st[j-1][i+(1<<(j-1))]);
        }
    }
    int r=0;
    for(int i=1;i<=n;i++){
        while(r<=n&&su<m){
            c[id[sa[q][++r]]]++;
            if(c[id[sa[q][r]]]==1)su++;
        }
        if(su<m)break;
        ans=max(ans,get_lcp(i,r));
        c[id[sa[q][i]]]--;
        if(!c[id[sa[q][i]]])su--;
    }
    printf("%d\n",ans);
    return 0;
}

标签：200005,int,Luogu5546,long,POI2000,BZOJ2946,sa,Sandy,define
来源： https://www.cnblogs.com/ZigZagKmp/p/11518851.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。