[SPOJ220]PHRASES(后缀数组+二分)

题目:

我是超链接

题意:

给定 n 个字符串,求在每个字符串中至少出现两次且不重叠的最长子串。

题解:

1A
和上道题基本一样啊,只不过不用倒过来了,二分,对于【至少出现两次】这个条件,我们以前也做过,在后缀分组内每组字符串给一个最大值最小值,差只要大于mid就好了

代码:

#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],*x,*y;
char st[N];bool vis[N];
struct hh{int minn,maxx;}zu[20];
void build_sa()
{
    m=300;x=X;y=Y;
    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;
    for (z=0;z<len;z++) zu[z].maxx=0,zu[z].minn=N;
    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;
        zu[id[sa[i]]].maxx=max(zu[id[sa[i]]].maxx,sa[i]);
        zu[id[sa[i-1]]].maxx=max(zu[id[sa[i-1]]].maxx,sa[i-1]);
        zu[id[sa[i]]].minn=min(zu[id[sa[i]]].minn,sa[i]);
        zu[id[sa[i-1]]].minn=min(zu[id[sa[i-1]]].minn,sa[i-1]);
      }
      else 
      {
        bool fff=1;
        for (z=0;z<len;z++) 
          if (!vis[z] || zu[z].maxx-zu[z].minn<mid){fff=0;break;}
        if (fff) return 1; 
        for (z=0;z<len;z++) zu[z].maxx=0,zu[z].minn=N;
        memset(vis,0,sizeof(vis));
      }
    bool fff=1;
    for (z=0;z<len;z++) 
      if (!vis[z] || zu[z].maxx-zu[z].minn<mid){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;fp=0;
    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++);
    }
    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);     
    }
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值