hdu 2328 Corporate Identity (后缀数组应用)

来源：互联网发布：ubuntu增加根目录空间编辑：程序博客网时间：2024/06/06 09:40

题意：

给出n个字符串，求字典序最小的最长公共子串。

题解：

二分长度，然后用后缀数组判断是否满足，根据height分组，找出分组中包含n个字符串的后缀的分组。字典序最小也容易处理，直接看代码就好了。先前忘记了n有4000多了，结果桶排序的上限搞错了，一直wa。

#include<iostream>#include<math.h>#include<stdio.h>#include<algorithm>#include<string.h>#include<vector>#include<queue>#include<map>#include<set>using namespace std;#define B(x) (1<<(x))typedef long long ll;const int oo=0x3f3f3f3f;const ll OO=1LL<<61;const int MOD=10007;const int maxn=910000;int Rank[maxn],SA[maxn],height[maxn];int t1[maxn],t2[maxn],t3[maxn],t4[maxn];char str[205];int ss[maxn],which[maxn],n;int res,mi;set<int>g;void Swap(int*& x,int*& y){    int *temp=x;    x=y;    y=temp;}int cmp(int t[],int a,int b,int l){    return t[a]==t[b]&&t[a+l]==t[b+l];}void build_SA(int s[],int len,int up){    int *k1=t1,*k2=t2,*r=t3,*cnt=t4;    for(int i=0;i<up;i++)cnt[i]=0;    for(int i=0;i<len;i++)cnt[k1[i]=s[i]]++;    for(int i=1;i<up;i++)cnt[i]+=cnt[i-1];    for(int i=len-1;i>=0;i--)SA[--cnt[k1[i]]]=i;    for(int d=1,p=1;p<len;d<<=1,up=p){        p=0;        for(int i=len-d;i<len;i++)k2[p++]=i;        for(int i=0;i<len;i++)if(SA[i]>=d)k2[p++]=SA[i]-d;        for(int i=0;i<len;i++)r[i]=k1[k2[i]];        for(int i=0;i<up;i++)cnt[i]=0;        for(int i=0;i<len;i++)cnt[r[i]]++;        for(int i=1;i<up;i++)cnt[i]+=cnt[i-1];        for(int i=len-1;i>=0;i--)SA[--cnt[r[i]]]=k2[i];        Swap(k1,k2);        k1[SA[0]]=0;        p=1;        for(int i=1;i<len;i++){            k1[SA[i]]= cmp(k2,SA[i-1],SA[i],d) ? p-1 : p++;        }    }}void get_height(int s[],int len){    for(int i=1;i<=len;i++)Rank[SA[i]]=i;    for(int i=0,p=0;i<len;i++){        int j=SA[Rank[i]-1];        while(s[i+p]==s[j+p])p++;        height[Rank[i]]=p;        if(p)p--;    }}bool ok(int tag,int len){    int i=1;    mi=oo;    while(1){        while(i<=len&&height[i]<tag)i++;        if(i>len)break;        g.clear();        g.insert(which[SA[i-1]]);        mi=i-1;        while(i<=len&&height[i]>=tag){            g.insert(which[SA[i]]);            i++;        }        if(g.size()>=n)return true;    }    return false;}void solve(int len){    g.clear();    res=oo;    int l=1,r=len,mid,ans=0;    while(l<=r){        mid=(r+l)>>1;        if(ok(mid,len)){            ans=mid;            res=mi;            l=mid+1;        }else{            r=mid-1;        }    }    if(res==oo||ans==0)        printf("IDENTITY LOST\n");    else{        for(int i=SA[res];i<SA[res]+ans;i++)            printf("%c",(char)(ss[i]+'a'-1));        puts("");    }}int main(){    int l,T,tlen;    while(scanf("%d",&n)!=EOF){        if(n==0)break;        l=0;        for(int i=1;i<=n;i++){            scanf("%s",str);            tlen=strlen(str);            for(int j=0;j<tlen;j++){                which[l]=i;                ss[l++]=str[j]-'a'+1;            }            which[l]=i;            ss[l++]=i+30;        }        l--;        ss[l]=0;        build_SA(ss,l+1,4300);        get_height(ss,l);        solve(l);    }    return 0;}/**3abacabababababacbaaaabacababc*/

0 0