[POJ1226]Substrings(后缀数组+二分)

题目:

我是超链接

题解:

求长度一般二分答案

我们发现如果中间加入一个&的话,前面后面可能构成关于ta的长度链,那么我们每次添加不一样的字符,这就要求我们把s设为一个数数数数组

只需要先将每个字符串都反过来写一遍,中间用一个互不相同的且没有出现在字符串中的字符隔开,再将 n 个字符串全部连起来,中间也是用一个互不相同的且没有出现在字符串中的字符隔开,求后缀数组。然后二分答案,再将后缀分组。判断的时候,要看是否有一组后缀在每个原来的字符串或反转后的字符串中出现。

关于看看是不是每一个字符串都出现过有一个比较好的方法:设置数组id[i]表示i这个位置属于几组,然后用bool数组判断就好了

代码:

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int N=100005;
int c[N],x[N],n,sa[N],y[N],rak[N],height[N],len,m,fp,id[N],s[N];
char st[105];bool vis[105];
void build_sa()
{
    m=300;
    for (int i=0;i<m;i++) c[i]=0;
    for (int i=0;i<n;i++) c[x[i]=s[i]]++;
    for (int i=1;i<m;i++) c[i]+=c[i-1];
    for (int i=n-1;i>=0;i--) sa[--c[x[i]]]=i;

    for (int k=1;k<=n;k<<=1)
    {
        int p=0;
        for (int i=n-k;i<n;i++) y[p++]=i;
        for (int i=0;i<n;i++) if (sa[i]>=k) y[p++]=sa[i]-k;

        for (int i=0;i<m;i++) c[i]=0;
        for (int i=0;i<n;i++) c[x[y[i]]]++;
        for (int i=1;i<m;i++) c[i]+=c[i-1];
        for (int i=n-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i];

        swap(x,y);
        p=1; x[sa[0]]=0;
        for (int i=1;i<n;i++) x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&((sa[i]+k>=n?-1:y[sa[i]+k])==(sa[i-1]+k>=n?-1:y[sa[i-1]+k]))?p-1:p++;
        if (p>n) break;
        m=p;
    }
}
void build_lcp()
{
    for (int i=0;i<n;i++) rak[sa[i]]=i;
    height[0]=0;
    int k=0;
    for (int i=0;i<n;i++)
    {
        if (rak[i]==0) continue;
        if (k) k--;
        int j=sa[rak[i]-1];
        while (i+k<n && j+k<n && s[i+k]==s[j+k]) k++;
        height[rak[i]]=k;
    }
}
bool check(int mid)
{
    int z;
    memset(vis,0,sizeof(vis));
    for (int i=1;i<n;i++)
      if (height[i]>=mid) 
      {
        vis[id[sa[i]]]=1;
        vis[id[sa[i-1]]]=1;
      }
      else 
      {
        bool fff=1;
        for (z=0;z<len;z++) 
          if (!vis[z]){fff=0;break;}
        if (fff) return 1; 
        memset(vis,0,sizeof(vis));
      }
    bool fff=1;
    for (z=0;z<len;z++) 
      if (!vis[z]){fff=0;break;}
    if (fff) return 1; 
    return 0;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
    scanf("%d",&len);
    if (len==1) {scanf("%s",st);printf("%d\n",strlen(st));continue;}
    n=-1;
    for (int i=0;i<len;i++)
    {
        scanf("%s",st);
        int ll=strlen(st);
        for (int j=0;j<ll;j++) s[++n]=st[j],id[n]=i;
        s[++n]='#'+(fp++);
        for (int j=ll-1;j>=0;j--) s[++n]=st[j],id[n]=i;
        s[++n]='#'+(fp++);
    }
    build_sa();
    build_lcp();
    int l=1,r=n,ans=0;
    while (l<=r)
    {
        int mid=(l+r)>>1;
        if (check(mid)) ans=mid,l=mid+1;
        else r=mid-1;
    }
    printf("%d\n",ans);     
    }
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值