題目鏈接：點擊查看

題目大意：給出n個字符串，求出n個字符串中最長的公共子串，如果沒有，輸出IDENTITY LOST

題目分析：可以直接用二分+后綴數組來做，先將n個字符串連接起來，中間用一個不同的字符連接，然后二分枚舉最大長度，轉換成check判斷的方式，來利用height數組將其分為不同的組，若其中有一個組中的sa數組分布在所有的n個字符串中，說明滿足條件，不過這個題目卡后綴數組的時間比較緊，如果二分的范圍是1~inf會超時，所以我們直接從1~len就能勉強卡過了

代碼：

#include<iostream> #include<cstdio> #include<string> #include<ctime> #include<cstring> #include<algorithm> #include<stack> #include<queue> #include<map> #include<set> #include<sstream> using namespace std;typedef long long LL;const int inf=0x3f3f3f3f;const int N=1e6+100;char s1[N],s2[N],ans[N];int sa[N]; //SA數組，表示將S的n個后綴從小到大排序后把排好序的 //的后綴的開頭位置順次放入SA中 int t1[N],t2[N],c[N];int rk[N],height[N],belong[N],len,n;int s[N];bool vis[4010];void build_sa(int s[],int n,int m)//n為添加0后的總長 {int i,j,p,*x=t1,*y=t2;for(i=0;i<m;i++) c[i]=0;for(i=0;i<n;i++) c[x[i]=s[i]]++;for(i=1;i<m;i++) c[i]+=c[i-1];for(i=n-1;i>=0;i--) sa[--c[x[i]]]=i;for(j=1;j<=n;j<<=1) {p=0;for(i=n-j;i<n;i++) y[p++]=i;for(i=0;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;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]+j]==y[sa[i]+j]?p-1:p++;if(p>=n) break;m=p;} }void get_height(int s[],int n)//n為添加0后的總長 {int i,j,k=0;for(i=0;i<=n;i++)rk[sa[i]]=i;for(i=0;i<n;i++) {if(k) k--;j=sa[rk[i]-1];while(s[i+k]==s[j+k]) k++;height[rk[i]]=k;} }void solve(int base=128) {build_sa(s,len+1,base);get_height(s,len); }bool check(int x) {memset(vis,false,sizeof(vis));int cnt=0;for(int i=2;i<=len;i++){if(height[i]<x){memset(vis,false,sizeof(vis));cnt=0;continue;}if(belong[sa[i-1]]!=-1&&!vis[belong[sa[i-1]]]){vis[belong[sa[i-1]]]=true;cnt++;}if(belong[sa[i]]!=-1&&!vis[belong[sa[i]]]){vis[belong[sa[i]]]=true;cnt++;}if(cnt==n){for(int j=0;j<x;j++)ans[j]=s[sa[i]+j]+'a'-1;ans[x]=0;return true;}}return false; }int main() { // freopen("input.txt","r",stdin); // ios::sync_with_stdio(false);while(scanf("%d",&n)!=EOF&&n){len=0;int pos=30;for(int i=1;i<=n;i++){char str[210];scanf("%s",str);for(int j=0;str[j];j++){s[len]=str[j]-'a'+1;belong[len]=i;len++;}s[len]=pos;belong[len]=-1;pos++;len++;}belong[len]=-1;s[len]=0;solve(pos);int l=1,r=len;bool flag=false;while(l<=r){int mid=l+r>>1;if(check(mid)){flag=true;l=mid+1;}elser=mid-1;}if(flag)puts(ans);elseputs("IDENTITY LOST");}return 0; }

?

總結

以上是生活随笔為你收集整理的POJ - 3450 Corporate Identity(二分+后缀数组)的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。