【BZOJ】【2946】【POI2000】公共串

后缀数组

---

　　好感动，复习了下后缀数组居然写出来了……（感谢ykz大神）

　　求最长公共子串……WA了一发是因为：【不同字符串之间要用不同的特殊字符隔开】否则就会匹配到相同→_→比如都是aaa结尾，如果用相同特殊字符就会使得最长公共子串变成aaa#这样子……

/**************************************************************
    Problem: 2946
    User: Tunix
    Language: C++
    Result: Accepted
    Time:60 ms
    Memory:4104 kb
****************************************************************/
 
//BZOJ 2946
#include<vector>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)
#define D(i,j,n) for(int i=j;i>=n;--i)
#define pb push_back
using namespace std;
typedef long long LL;
inline int getint(){
    int r=1,v=0; char ch=getchar();
    for(;!isdigit(ch);ch=getchar()) if(ch=='-')r=-1;
    for(; isdigit(ch);ch=getchar()) v=v*10+ch-'0';
    return r*v;
}
const int N=1e5+10,INF=~0u>>2;
/*******************template********************/
int sa[N],rank[N],height[N],belong[N],wa[N],wb[N],c[N],n,m;
int len[6];
char s[N];
bool cmp(int *r,int a,int b,int l){
    return r[a]==r[b] && r[a+l]==r[b+l];
}
void DA(char *s,int *sa,int n,int m){
    int i,j,p,*x=wa,*y=wb;
    rep(i,m) c[i]=0;
    rep(i,n) c[x[i]=s[i]]++;
    F(i,1,m-1) c[i]+=c[i-1];
    D(i,n-1,0) sa[--c[x[i]]]=i;
    for(j=1,p=0;p<n;j<<=1,m=p){
        for(p=0,i=n-j;i<n;i++) y[p++]=i;
        rep(i,n) if (sa[i]>=j) y[p++]=sa[i]-j;
 
        rep(i,m) c[i]=0;
        rep(i,n) c[x[y[i]]]++;
        F(i,1,m-1) c[i]+=c[i-1];
        D(i,n-1,0) sa[--c[x[y[i]]]]=y[i];
        swap(x,y); p=1; x[sa[0]]=0;
        F(i,1,n-1) x[sa[i]]=cmp(y,sa[i-1],sa[i],j) ? p-1 : p++;
    }
}
void calheight(char *s,int *sa,int n){
    int k=0;
    F(i,1,n) rank[sa[i]]=i;
    rep(i,n){
        if (k) k--;
        int j=sa[rank[i]-1];
        while(s[i+k]==s[j+k]) k++;
        height[rank[i]]=k;
    }
}
int main(){
#ifndef ONLINE_JUDGE
    freopen("2946.in","r",stdin);
    freopen("2946.out","w",stdout);
#endif
    n=getint();int l=0;
    F(i,1,n){
        scanf("%s",s+l);
        len[i]=strlen(s)-l;
        l=strlen(s);
        s[l++]='a'-n+i;
    }
    rep(i,l) s[i]=s[i]-'a'+10;
    l--;
    DA(s,sa,l+1,40);
    calheight(s,sa,l);
    int tmp=1;
    rep(i,l){
        if (s[i]<10) {tmp++;continue;}
        belong[rank[i]]=tmp;
    }
 
 
    bool vis[6]={0};
    int L=0,R=2000,mid,cnt,ans=0;
    while(L<=R){
        mid=(L+R)>>1;bool sign=0;
        cnt=0; memset(vis,0,sizeof vis);
        F(i,n,l){
            if (height[i]<mid){
                cnt=1;memset(vis,0,sizeof vis);
                vis[belong[i]]=1; continue;
            }
            if (!vis[belong[i]]) vis[belong[i]]=1,cnt++;
            if (cnt==n) sign=1;
        }
        if (sign) ans=mid,L=mid+1;
        else R=mid-1;
    }
    printf("%d
",ans);
    return 0;
}

2946: [Poi2000]公共串

Time Limit: 3 Sec  Memory Limit: 128 MB
Submit: 160  Solved: 67
[Submit][Status][Discuss]

Description

 

       给出几个由小写字母构成的单词，求它们最长的公共子串的长度。

任务：

l        读入单词

l        计算最长公共子串的长度

l        输出结果

 

Input

 

文件的第一行是整数 n，1<=n<=5，表示单词的数量。接下来n行每行一个单词，只由小写字母组成，单词的长度至少为1，最大为2000。

 

Output

仅一行，一个整数，最长公共子串的长度。

 

Sample Input

3
abcb
bca
acbc

Sample Output

２

HINT

Source

[Submit][Status][Discuss]