zoukankan      html  css  js  c++  java
  • POJ 3261 后缀数组+二分

    思路:
    论文题…
    二分+对后缀分组

    这块一开始不用基数排序 会更快的(其实区别不大)

    //By SiriusRen
    #include <cstdio>
    #include <cstring>
    #include <algorithm>
    using namespace std;
    #define N 20050
    int n,k,s[N],cntA[N],cntB[N],A[N],B[N],rk[N],sa[N],tsa[N],ht[N];
    struct Node{int pos,wei;}nd[N];
    bool cmp(Node a,Node b){return a.wei<b.wei;}
    void SA(){
        rk[nd[1].pos]=1;
        for(int i=2;i<=n;i++)rk[nd[i].pos]=rk[nd[i-1].pos]+(nd[i].wei!=nd[i-1].wei);
        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;
        }
    }
    bool check(int x){
        int bgn=0;
        for(int i=1;i<=n;i++){
            if(ht[i]>=x){if(i-bgn==k-1)return 1;}
            else bgn=i;
        }return 0;
    }
    int main(){
        scanf("%d%d",&n,&k);
        for(int i=1;i<=n;i++){
            scanf("%d",&s[i]);
            nd[i].pos=i,nd[i].wei=s[i];
        }sort(nd+1,nd+1+n,cmp);
        SA();
        int l=0,r=n,ans;
        while(l<=r){
            int mid=(l+r)>>1;
            if(check(mid))l=mid+1,ans=mid;
            else r=mid-1;
        }printf("%d
    ",ans);
    }

    这里写图片描述

  • 相关阅读:
    苹果审核被拒,Guideline 1.1.6
    Java基础:类与继承
    Java中的static关键字解析
    Entity Framework 预备知识
    Entity Framework Code First执行SQL语句、视图及存储过程
    java多线程技能
    流与文件
    多线程
    集合
    断言
  • 原文地址:https://www.cnblogs.com/SiriusRen/p/6532079.html
Copyright © 2011-2022 走看看