spoj220 Relevant Phrases of Annihilation

最新推荐文章于 2020-02-07 11:05:28 发布

原创最新推荐文章于 2020-02-07 11:05:28 发布 · 672 阅读

0 ·

CC 4.0 BY-SA版权

SA&&SAM 专栏收录该内容

19 篇文章

订阅专栏

该问题要求找到一组字符串中，每个字符串都至少包含一个长度为L的相同子串，且这些子串在原字符串内的出现位置不能有交叉。采用二分查找方法，针对每个可能的子串长度L，通过维护每个字符串中该子串出现的最远和最近位置，判断是否存在满足条件的子串。

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

题意是给n个串，求出一个最长的子串在每个串中至少出现两次。直接二分答案来判断，因为题目里要求在一个串中，两次出现的位置不能有交叉，所以这要想办法特殊处理下，大体的思路还是循环height数组，按>=当前二分长度来分组，处理一组时，用两个数组分别记录一下本组后缀在各个串中出现的最大位置和最小位置，那么如果某一个串上他的最大位置和最小位置的差大于等于当前长度了，就是两个没有交叉的匹配了。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <string>
#include <cstring>
using namespace std;
typedef long long ll;
const int maxn=220100;
int s[maxn];
int sa[maxn],t[maxn],t2[maxn],c[maxn];
int sa2[maxn];
char s1[maxn],s2[maxn];
int rank[maxn],height[maxn];
int l1,l2;
int len,len1;
int n,m,ans;
int belong[maxn];
inline int idx(char s)
{
   if (s>='a' && s<='z') return s-'a'+11;

}
void getheight(int n)
{
    int i,j,k=0;
    for (i=0; i<=n; i++) rank[sa[i]]=i;

    for (i=0; i<n; i++)
    {
        if (k) k--;
        int j=sa[rank[i]-1];
        while(s[i+k]==s[j+k]) k++;
        height[rank[i]]=k;
    }
}

void build_ss(int m,int n)
{
    n++;
    int i,*x=t,*y=t2;
    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 (i=n-k; i<n; i++) y[p++]=i;
        for (i=0; i<n; i++) if (sa[i]>=k) y[p++]=sa[i]-k;

        for (i=0; i<m; i++) c[i]=0;
        for (i=0; i<n; i++) c[x[y[i]]]++;
        for (i=1; i<m; i++) c[i]+=c[i-1];
        for (i=n-1; i>=0; i--) sa[--c[x[y[i]]]] = y[i];
        swap(x,y);
        p=1;
        x[sa[0]]=0;
        for (i=1; i<n; i++)
        x[sa[i]]=(y[sa[i-1]]==y[sa[i]] && y[sa[i-1]+k]==y[sa[i]+k])? p-1 : p++;
        if (p>=n) break;
        m=p;
    }
}
int d[maxn][22];

void RMQ_init()
{
    for (int i=1; i<=n; i++) d[i][0]=height[i];
     for (int j=1; (1<<j)<=n; j++)
      for (int i=1; (i+(1<<j)-1)<=n; i++)
      d[i][j]=min(d[i][j-1],d[i+(1<<(j-1))][j-1]);
}
int RMQ(int L,int R)
{
    if (L>R) swap(L,R);
    L++;

    int k=0;
    while((1<<(k+1))<=R-L+1) k++;
    return min(d[L][k],d[R-(1<<k)+1][k]);
}
void init()
{
    scanf("%d",&m);
    n=0;
    int ct=1;
    for (int i=0; i<m; i++)
    {
        scanf("%s",s1);
        l1=strlen(s1);
        for (int j=0; j<l1; j++)
        belong[n]=i,s[n++]=idx(s1[j]);

        belong[n]=-1;
        s[n++]=ct++;
    }
    s[n]=0;
    build_ss(40,n);
    getheight(n);
}
bool check(int k)
{
    int mi[11],mx[11];
    memset(mi,0x3f,sizeof mi);
    memset(mx,0,sizeof mx);

    bool ok=false;
    for (int i=1; i<=n; i++)
    {
       if (height[i]>=k)
       {
           mi[belong[sa[i-1]]]=min(mi[belong[sa[i-1]]],sa[i-1]);
           mx[belong[sa[i-1]]]=max(mx[belong[sa[i-1]]],sa[i-1]);
           mi[belong[sa[i]]]=min(mi[belong[sa[i]]],sa[i]);
           mx[belong[sa[i]]]=max(mx[belong[sa[i]]],sa[i]);

           ok=true;
           for (int j=0; j<m; j++)
           if (mx[j]-mi[j]<k) ok=false;
           if (ok) return true;
       }
       else
       {
            memset(mi,0x3f,sizeof mi);
            memset(mx,0,sizeof mx);
            mi[belong[sa[i]]]=sa[i];
            mx[belong[sa[i]]]=sa[i];
       }
    }
    return false;
}
void slove()
{
    int l=1,r=10000;
    ans=0;
    while (l<r)
    {
        int mid=(l+r)>>1;
        if (check(mid)) ans=mid,l=mid+1;
        else r=mid;
    }
    printf("%d\n",ans);
}
int tt;
int main()
{
//    freopen("in.txt","r",stdin);
    scanf("%d",&tt);
    while(tt--)
    {
        init();
        slove();
    }
}