[2016ecfinal][Gym-101194F]Mr. Panda and Fantastic Beasts(后缀数组模版题)

7 篇文章 0 订阅
1 篇文章 0 订阅

题目链接实现10S  题目链接实现3S(本文代码cf800ms) 后缀数组论文

题目大意是给你n个字符串S1~Sn,然后要你找出一个S1中的最短子串t,而且t不是S2~Sn的子串。

  • 首先我们需要字典序最小,那么就可以想到把全部才串连在一起跑一下后缀数组。
  • 然后遍历找到起点在S1的suffix(i),易知与suffix(i)前缀最长的是suffix(i-1)与suffix(i+1)。
  • 然后查询一下相互之间的最长公共前缀长度加上suffix(i)的起点是否超出S1的长度。
  • 若不超出则尝试更新最短答案。
#include<algorithm>
#include<iostream>
#include<math.h>
#include<set>
#include<map>
#include<cstring>
#include<string>
#define qcin; ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
using namespace std;

typedef long long ll;
const int maxn = 300010;
typedef int arr[maxn];

struct ST{
    int mm[maxn],dp[maxn][22],k;
    const int& (*fun[2])(const int&,const int&)={max,min};
    void init(int n,int b[],int fl){
        mm[0]=-1;
        k=fl;
        for(int i=1;i<=n;++i){
            mm[i]=mm[i-1]+((i&(i-1))==0?1:0);
            dp[i][0]=b[i];
        }
        for(int j=1;j<=mm[n];j++){
            for(int i=1;i+(1<<j)-1<=n;++i){
                dp[i][j]=fun[k](dp[i][j-1],dp[i+(1<<(j-1))][j-1]);
            }
        }
    }
    inline int qry(int l,int r){
        int kk=mm[r-l+1];
        return fun[k](dp[l][kk],dp[r-(1<<kk)+1][kk]);
    }
};

struct SA{
    //h[i]=suffix(sa[i-1])和 suffix(sa[i])的最长公共前缀的长度
    //sa[i]表示字典序第i小的后缀起点
    //rk[i]表示起点i的后缀的字典序排位
    int sa[maxn],rk[maxn],h[maxn];
    ST st;
    void init(string s,int n){
        static arr c,v,q;
        int m=127;
        for(int i=0;i<m;++i)c[i]=0;
        for(int i=0;i<n;++i)++c[s[i]];
        for(int i=1;i<m;++i)c[i]+=c[i-1];
        for(int i=n-1;~i;--i)sa[--c[s[i]]]=i;
        m=rk[sa[0]]=0;
        for(int i=1;i<n;++i)rk[sa[i]]=s[sa[i]]!=s[sa[i-1]]?++m:m;
        for(int j=1;;j<<=1){
            if(++m==n)break;
            for(int i=0;i<j;++i)v[i]=n-j+i;
            for(int i=0,k=j;i<n;++i)if(sa[i]>=j)v[k++]=sa[i]-j;
            for(int i=0;i<n;++i)q[i]=rk[v[i]];
            for(int i=0;i<m;++i)c[i]=0;
            for(int i=0;i<n;++i)++c[q[i]];
            for(int i=1;i<m;++i)c[i]+=c[i-1];
            for(int i=n-1;~i;--i)sa[--c[q[i]]]=v[i],v[i]=rk[i];
            m=rk[sa[0]]=0;
            for(int i=1;i<n;++i)rk[sa[i]]=v[sa[i]]!=v[sa[i-1]]||v[sa[i]+j]!=v[sa[i-1]+j]?++m:m;
        }
        int k=0;
        for(int i=0;i<n-1;++i){
            if(k)--k;
            if(!rk[i])continue;
            while(s[i+k]==s[sa[rk[i]-1]+k])++k;
            h[rk[i]]=k;
        }
        st.init(n-1,h,1);
    }
}sa;

int t,n,m,len;
char sp='*';
string s,tmp;
int pre[maxn],nxt[maxn];
void sol(int x){
    memset(pre,-1,(n+1)*sizeof(int));
    memset(nxt,-1,(n+1)*sizeof(int));
    int k=-1,ans=-1,fi=-1;
    for(int i=0;i<n;++i){
        if(sa.sa[i]>=len)k=i;
        else pre[i]=k;
    }
    k=-1;
    for(int i=n-1;~i;--i){
        if(sa.sa[i]>=len)k=i;
        else nxt[i]=k;
    }
    for(int i=0;i<n;++i){
        int mx=0;
        if(sa.sa[i]<len){
            if(~pre[i]){
                mx=max(mx,sa.st.qry(pre[i]+1,i));
            }
            if(~nxt[i]){
                mx=max(mx,sa.st.qry(i+1,nxt[i]));
            }
            if((ans>mx+1||ans==-1)&&sa.sa[i]+mx+1<=len){
                ans=mx+1;
                fi=sa.sa[i];
            }
        }
    }
    if(~ans){
        cout<<"Case #"<<x<<": ";
        for(int i=fi;i<fi+ans;++i){
            cout<<s[i];
        }
        cout<<endl;
        return;
    }
    cout<<"Case #"<<x<<": Impossible"<<endl;
}

int main(){
    qcin;
    cin>>t;
    for(int tt=1;tt<=t;++tt){
        cin>>n;s="";s.clear();
        for(int i=0;i<n;++i){
            cin>>tmp;
            if(i)s+=sp;
            else len=tmp.length();
            s+=tmp;
        }
        n=s.length();
        sa.init(s,s.length());
        sol(tt);
    }
}

 

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值