后缀数组 poj. 1743. Musial Theme

最新推荐文章于 2020-04-03 16:26:24 发布

原创最新推荐文章于 2020-04-03 16:26:24 发布 · 292 阅读

0 ·

CC 4.0 BY-SA版权

字符串专栏收录该内容

22 篇文章

订阅专栏

本文介绍了一种使用后缀数组和二分查找的方法来解决寻找给定数字串中是否存在长度大于五的相同子串的问题。该方法首先初始化后缀数组，并通过倍增策略构建后缀数组，接着利用高度数组进行二分查找，最终确定最长相同子串的长度。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题意感觉有点难懂。但就是问你，给你一个数字串，然后问有没有相同子串（不能重合且长度大于五），或者不用相同的，全部加上一个或减去一个数之后相等的就好。
先是二分长度，然后通过height数组来判断是否可行

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include<string>
using namespace std;
const int  maxn =4e4+10;
int n;
struct Suffix
{
    int s[maxn];
    int SA[maxn],ranked[maxn],height[maxn];
    int x[maxn],y[maxn],buc[maxn],len,m;
    void init(int *str)
    {
        //len=strlen(str);
        len=n;
        for(int i=0;i<len;i++) s[i]=(int)str[i];
        m=200;//
    }
    void build()
    {
        for (int i = 0; i < m; i++) buc[i] = 0; // buc 是一个桶
        for (int i = 0; i < len; i++) buc[x[i] = s[i]]++;
        for (int i = 1; i < m; i++)
             buc[i] += buc[i - 1];
        for (int i = len - 1; i >= 0; i--)
            SA[--buc[x[i]]] = i;
        for (int k = 1; k <= len; k <<= 1) { // k 倍增
          int p = 0;
          //对第二关键字排序，y[i]:第i大的第二关键字是谁
          // 后缀 len - k 及之后的所有后缀第二关键字最小。为0
          for (int i = len - 1; i >= len - k; i--)
            y[p++] = i;
          for (int i = 0; i < len; i++)
                if (SA[i] >= k)
                    y[p++] = SA[i] - k;
         //总体来排个序，求出SA
          for (int i = 0; i < m; i++) buc[i] = 0;
          for (int i = 0; i < len; i++)
                buc[x[y[i]]]++;
          for (int i = 1; i < m; i++)
                buc[i] += buc[i - 1];
          for (int i = len - 1; i >= 0; i--)
                SA[--buc[x[y[i]]]] = y[i];
          swap(x, y);
          p = 1; x[SA[0]] = 0;
          // 重新计算每个一元的名次。则x数组里存的是总体的顺序
          for (int i = 1; i < len; i++) {
            if (y[SA[i - 1]] == y[SA[i]] && y[SA[i - 1] + k] == y[SA[i] + k])
              x[SA[i]] = p - 1;
            else x[SA[i]] = p++;
          }
          if (p >= len) break; // 每个后缀的名次已经完全不同，不需要继续倍增
              m = p; // 更新名次的最大值。
        }
    }
    void LCP()
    {
        for(int i=0;i<len;i++)
            ranked[SA[i]]=i;
        int k=0;
        for(int i=0;i<len;i++)
        {
            if(ranked[i]==0) {height[0]=0;continue;}
            if(k) k--;
            int j=SA[ranked[i]-1];
            while(s[i+k]==s[j+k]&&i+k<len&&j+k<len) k++;
            height[ranked[i]]=k;
        }
    }
 /*  int LCS(char *s1)
    {
        int l=strlen(s1);
        s1[l]=' ';
        scanf("%s",s1+l+1);
        init(s1);
        build();
        LCP();
        int ans=0;
        for(int i=1;i<len;i++)
        if((SA[i-1]<l&&SA[i]>l)||(SA[i-1]>l&&SA[i]<l))
                ans=max(ans,height[i]);
        return ans;
    }*/
};
Suffix a;
int ok(int mid)
{
    int i=1;
    while(1){
    for(;i<n;i++)
        if(a.height[i]>=mid)
            break;
    if(i>=n) break;
    int minn=a.SA[i-1];
    int maxx=a.SA[i-1];
    while(i<n&&a.height[i]>=mid)
    {
        minn=min(a.SA[i],minn);
        maxx=max(a.SA[i],maxx);
        i++;
    }
    if(maxx-minn>=mid) return 1;
    }
    return 0;

}
int num[maxn];
int str[maxn];
int main()
{
    while(scanf("%d",&n)!=EOF&&n)
    {
       memset(str,0,sizeof(str));
      memset(num,0,sizeof(num));
        for(int i=0;i<n;i++)
            scanf("%d",&num[i]);
        if(n<10) printf("0\n");
        else
        {
            for(int i=1;i<n;i++)
                str[i-1]=num[i]-num[i-1]+100;
             n--;
         //   str[n-1]=0;
            a.init(str);
            a.build();
            a.LCP();
            int l=4,r=n/2,res=0;
            while(l<=r)
            {
                int mid=l+(r-l)/2;
                if(ok(mid)){
                    res=mid;
                    l=mid+1;
                }
                else r=mid-1;
            }
            if(res>=4)
            printf("%d\n",res+1);
            else printf("0\n");
        }
    }
    return 0;
}