zoukankan      html  css  js  c++  java
  • hihocoder1415 重复旋律3

    思路:
    扫一遍height 判一下即可

    //By SiriusRen
    #include <cstdio>
    #include <cstring>
    #include <algorithm>
    using namespace std;
    #define N 200050
    int n,cntA[N],cntB[N],A[N],B[N],tsa[N],sa[N],rk[N],ht[N],ans;
    char s[N],b[N];
    void SA(){
        for(int i=1;i<=n;i++)cntA[s[i]]++;
        for(int i=1;i<=n;i++)cntA[i]+=cntA[i-1];
        for(int i=n;i;i--)sa[cntA[s[i]]--]=i;
        rk[sa[1]]=1;
        for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(s[sa[i]]!=s[sa[i-1]]);
        for(int l=1;rk[sa[n]]<n;l<<=1){
            memset(cntA,0,sizeof(cntA));
            memset(cntB,0,sizeof(cntB));
            for(int i=1;i<=n;i++)
                cntA[A[i]=rk[i]]++,
                cntB[B[i]=i+l<=n?rk[i+l]:0]++;
            for(int i=1;i<=n;i++)cntA[i]+=cntA[i-1],cntB[i]+=cntB[i-1];
            for(int i=n;i;i--)tsa[cntB[B[i]]--]=i;
            for(int i=n;i;i--)sa[cntA[A[tsa[i]]]--]=tsa[i];
            rk[sa[1]]=1;
            for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(A[sa[i]]!=A[sa[i-1]]||B[sa[i]]!=B[sa[i-1]]);
        }
        for(int i=1,j=0;i<=n;i++){
            j=j?j-1:0;
            while(s[i+j]==s[sa[rk[i]-1]+j])j++;
            ht[rk[i]]=j;
        }
    }
    int main(){
        scanf("%s%s",s+1,b+1);
        int r=strlen(s+1),m=strlen(b+1);
        for(int i=1;i<=m;i++)s[r+i+1]=b[i];
        s[r+1]='#';n=r+m+1,SA();
        for(int i=1;i<=n;i++){
            int maxx=max(sa[i],sa[i-1]),minn=min(sa[i],sa[i-1]);
            if(maxx>r&&minn<=r)ans=max(ans,ht[i]);
        }
        printf("%d
    ",ans);
    }
  • 相关阅读:
    python学习笔记(2)
    数据分析工具pandas
    python学习笔记(1)
    python学习笔记(3)
    python学习笔记(2):科学计算及数据可视化入门
    python学习笔记(1):python基础
    js setTimeout 和 setInterval 区别
    C#根据URL生成签名
    jquery.validate.js客户端验证
    redis学习(一)
  • 原文地址:https://www.cnblogs.com/SiriusRen/p/6532070.html
Copyright © 2011-2022 走看看