题意：多个字符串的最长公共子串。

解题关键：字符串的任何一个子串都是这个字符串的某个后缀的前缀。求A和B的最长公共子串等价于求A的后缀和B的后缀的最长公共前缀的最大值。

后缀数组的经典例题，连接在一起，二分长度，height数组遍历即可。

注意flag的问题，采用二分小于的方式，可能会出现有最优解但是flag为false的情况，下界需要-1，采用0，而采用等于的话，就不会出现，不过有些题会出现死循环。

还有因为多添加的符号一定不会加入vis数组，所以vis数组只需建立4000即可。

为什么两个字符串不需要二分？而多个字符串需要二分？因为 两个字符串可以直接判定height[i]是属于两个字符串的最大公共子串，而多个必须通过vis数组判定。

法一：

 #include<cstdlib>

 #include<cstring>

 #include<cstdio>

 #include<algorithm>

 #include<iostream>

 #include<cmath>

 using namespace std;

 const int N=;

 int r[N],id[N];

 char tmp[],ans[];

 bool vis[];

 int wa[N],wb[N],wv[N],wc[N],n,m;

 bool cmp(int *r,int a,int b,int l){return r[a]==r[b]&&r[a+l]==r[b+l];}

 void make_sa(int *r,int *sa,int n,int m){

     int i,j,p,*x=wa,*y=wb;

     for(i=;i<m;i++) wc[i]=;

     for(i=;i<n;i++) wc[x[i]=r[i]]++;

     for(i=;i<m;i++) wc[i]+=wc[i-];

     for(i=n-;i>=;i--) sa[--wc[x[i]]]=i;

     for(j=,p=;p<n;j*=,m=p){

         for(p=,i=n-j;i<n;i++) y[p++]=i;

         for(i=;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;

         for(i=;i<n;i++) wv[i]=x[y[i]];

         for(i=;i<m;i++) wc[i]=;

         for(i=;i<n;i++) wc[wv[i]]++;

         for(i=;i<m;i++) wc[i]+=wc[i-];

         for(i=n-;i>=;i--) sa[--wc[wv[i]]]=y[i];

         for(swap(x,y),p=,x[sa[]]=,i=;i<n;i++) x[sa[i]]=cmp(y,sa[i-],sa[i],j)?p-:p++;

     }

     return;

 }

 int rank1[N],height[N],sa[N];

 void make_height(int *r,int *sa,int n){

     int i,j,k=;

     for(i=;i<=n;i++) rank1[sa[i]]=i;

     for(i=;i<n;height[rank1[i++]]=k)

         for(k?k--:,j=sa[rank1[i]-];r[i+k]==r[j+k];k++);

     return;

 }

 bool check(int x){

     memset(vis,,sizeof vis);

     int cnt=;

     for(int i=;i<=n;i++){

         if(height[i]<x){

             memset(vis,,sizeof vis);

             cnt=;

             continue;

         }

         if(!vis[id[sa[i-]]]) vis[id[sa[i-]]]=true,cnt++;

         if(!vis[id[sa[i]]]) vis[id[sa[i]]]=true,cnt++;

         if(cnt==m){

             for(int j=;j<x;j++) ans[j]=r[sa[i]+j]+'a'-;ans[x]=;

             return true;

         }

     }

     return false;

 }

 bool erfen(int l,int r){

     bool flag=false;

     while(l<r){

         int mid=(l+r+)>>;

         if(check(mid)) l=mid,flag=true;

         else r=mid-;

     }

     return flag;

 }

 int main() {

     while(scanf("%d",&m)&&m){

         n=;

         int temp=;

         for(int i=;i<=m;i++){

             scanf("%s",tmp);

             int siz=strlen(tmp);

             for(int j=;j<siz;j++)    id[n]=i,r[n++]=tmp[j]-'a'+;

             id[n]=temp;

             r[n++]=temp++;

         }

         r[n]=;

         make_sa(r,sa,n+,);

         make_height(r,sa,n);

         bool f=erfen(,);//为什么0可以，1不可以

         if(f) printf("%s\n",ans);

         else printf("IDENTITY LOST\n");

     }

 }

法二：

 #include<cstdlib>

 #include<cstring>

 #include<cstdio>

 #include<algorithm>

 #include<iostream>

 #include<cmath>

 using namespace std;

 const int N=;

 int r[N],id[N];

 char tmp[],ans[];

 bool vis[];

 int wa[N],wb[N],wv[N],wc[N],n,m;

 bool cmp(int *r,int a,int b,int l){return r[a]==r[b]&&r[a+l]==r[b+l];}

 void make_sa(int *r,int *sa,int n,int m){

     int i,j,p,*x=wa,*y=wb;

     for(i=;i<m;i++) wc[i]=;

     for(i=;i<n;i++) wc[x[i]=r[i]]++;

     for(i=;i<m;i++) wc[i]+=wc[i-];

     for(i=n-;i>=;i--) sa[--wc[x[i]]]=i;

     for(j=,p=;p<n;j*=,m=p){

         for(p=,i=n-j;i<n;i++) y[p++]=i;

         for(i=;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;

         for(i=;i<n;i++) wv[i]=x[y[i]];

         for(i=;i<m;i++) wc[i]=;

         for(i=;i<n;i++) wc[wv[i]]++;

         for(i=;i<m;i++) wc[i]+=wc[i-];

         for(i=n-;i>=;i--) sa[--wc[wv[i]]]=y[i];

         for(swap(x,y),p=,x[sa[]]=,i=;i<n;i++) x[sa[i]]=cmp(y,sa[i-],sa[i],j)?p-:p++;

     }

     return;

 }

 int rank1[N],height[N],sa[N];

 void make_height(int *r,int *sa,int n){

     int i,j,k=;

     for(i=;i<=n;i++) rank1[sa[i]]=i;

     for(i=;i<n;height[rank1[i++]]=k)

         for(k?k--:,j=sa[rank1[i]-];r[i+k]==r[j+k];k++);

     return;

 }

 bool check(int x){

     memset(vis,,sizeof vis);

     int cnt=;

     for(int i=;i<=n;i++){

         if(height[i]<x){

             memset(vis,,sizeof vis);

             cnt=;

             continue;

         }

         if(!vis[id[sa[i-]]]) vis[id[sa[i-]]]=true,cnt++;

         if(!vis[id[sa[i]]]) vis[id[sa[i]]]=true,cnt++;

         if(cnt==m){

             for(int j=;j<x;j++) ans[j]=r[sa[i]+j]+'a'-;ans[x]=;

             return true;

         }

     }

     return false;

 }

 bool erfen(int l,int r){

     bool flag=false;

     while(l<=r){

         int mid=(l+r+)>>;

         if(check(mid)) l=mid+,flag=true;

         else r=mid-;

     }

     return flag;

 }

 int main() {

     while(scanf("%d",&m)&&m){

         n=;

         int temp=;

         for(int i=;i<=m;i++){

             scanf("%s",tmp);

             int siz=strlen(tmp);

             for(int j=;j<siz;j++)    id[n]=i,r[n++]=tmp[j]-'a'+;

             id[n]=temp;

             r[n++]=temp++;

         }

         r[n]=;

         make_sa(r,sa,n+,);

         make_height(r,sa,n);

         bool f=erfen(,);//为什么0可以，1不可以

         if(f) printf("%s\n",ans);

         else printf("IDENTITY LOST\n");

     }

 }

dc3：依然注意是3倍的问题

 #include<cstdlib>

 #include<cstring>

 #include<cstdio>

 #include<algorithm>

 #include<iostream>

 #include<cmath>

 #define inf 0x3f3f3f3f

 typedef long long ll;

 using namespace std;

 const int N=;

 #define F(x) ((x)/3+((x)%3==1?0:tb))

 #define G(x) ((x)<tb?(x)*3+1:((x)-tb)*3+2)

 int wa[N],wb[N],wv[N],ws1[N];

 int c0(int *r,int a,int b){ return r[a]==r[b]&&r[a+]==r[b+]&&r[a+]==r[b+]; }

 int c12(int k,int *r,int a,int b){

     if(k==) return r[a]<r[b]||r[a]==r[b]&&c12(,r,a+,b+);

     else return r[a]<r[b]||r[a]==r[b]&&wv[a+]<wv[b+];

 }

 void sort(int *r,int *a,int *b,int n,int m){

     int i;

     for(i=;i<n;i++) wv[i]=r[a[i]];

     for(i=;i<m;i++) ws1[i]=;

     for(i=;i<n;i++) ws1[wv[i]]++;

     for(i=;i<m;i++) ws1[i]+=ws1[i-];

     for(i=n-;i>=;i--) b[--ws1[wv[i]]]=a[i];

     return;

 }

 void dc3(int *r,int *sa,int n,int m){

     int i,j,*rn=r+n,*san=sa+n,ta=,tb=(n+)/,tbc=,p;

     r[n]=r[n+]=;

     for(i=;i<n;i++) if(i%!=) wa[tbc++]=i;

     sort(r+,wa,wb,tbc,m);

     sort(r+,wb,wa,tbc,m);

     sort(r,wa,wb,tbc,m);

     for(p=,rn[F(wb[])]=,i=;i<tbc;i++)

     rn[F(wb[i])]=c0(r,wb[i-],wb[i])?p-:p++;

     if(p<tbc) dc3(rn,san,tbc,p);

     else for(i=;i<tbc;i++) san[rn[i]]=i;

     for(i=;i<tbc;i++) if(san[i]<tb) wb[ta++]=san[i]*;

     if(n%==) wb[ta++]=n-;

     sort(r,wb,wa,ta,m);

     for(i=;i<tbc;i++) wv[wb[i]=G(san[i])]=i;

     for(i=,j=,p=;i<ta && j<tbc;p++)

     sa[p]=c12(wb[j]%,r,wa[i],wb[j])?wa[i++]:wb[j++];

     for(;i<ta;p++) sa[p]=wa[i++];

     for(;j<tbc;p++) sa[p]=wb[j++];

     return;

 }

 int rank1[N],height[N],sa[*N];

 void make_height(int *r,int *sa,int n){

     int i,j,k=;

     for(i=;i<=n;i++) rank1[sa[i]]=i;

     for(i=;i<n;height[rank1[i++]]=k)

         for(k?k--:,j=sa[rank1[i]-];r[i+k]==r[j+k];k++);

     return;

 }

 int r[*N],id[N],n,m;

 char tmp[],ans[];

 bool vis[];

 bool check(int x){

     memset(vis,,sizeof vis);

     int cnt=;

     for(int i=;i<=n;i++){

         if(height[i]<x){

             memset(vis,,sizeof vis);

             cnt=;

             continue;

         }

         if(!vis[id[sa[i-]]]) vis[id[sa[i-]]]=true,cnt++;

         if(!vis[id[sa[i]]]) vis[id[sa[i]]]=true,cnt++;

         if(cnt==m){

             for(int j=;j<x;j++) ans[j]=r[sa[i]+j]+'a'-;ans[x]=;

             return true;

         }

     }

     return false;

 }

 bool erfen(int l,int r){

     bool flag=false;

     while(l<r){

         int mid=(l+r+)>>;

         if(check(mid)) l=mid,flag=true;

         else r=mid-;

     }

     return flag;

 }

 int main() {

     while(scanf("%d",&m)&&m){

         n=;

         int temp=;

         for(int i=;i<=m;i++){

             scanf("%s",tmp);

             int siz=strlen(tmp);

             for(int j=;j<siz;j++)    id[n]=i,r[n++]=tmp[j]-'a'+;

             id[n]=temp;

             r[n++]=temp++;

         }

         r[n]=;

         dc3(r,sa,n+,);

         make_height(r,sa,n);

         bool f=erfen(,);//为什么0可以，1不可以

         if(f) printf("%s\n",ans);

         else printf("IDENTITY LOST\n");

     }

 }

[poj3450]Corporate Identity(后缀数组)的更多相关文章

  1. POJ3450 Corporate Identity —— 后缀数组 最长公共子序列

    题目链接:https://vjudge.net/problem/POJ-3450 Corporate Identity Time Limit: 3000MS   Memory Limit: 65536 ...

  2. poj 3518 Corporate Identity 后缀数组->多字符串最长相同连续子串

    题目链接 题意:输入N(2 <= N <= 4000)个长度不超过200的字符串,输出字典序最小的最长公共连续子串; 思路:将所有的字符串中间加上分隔符,注:分隔符只需要和输入的字符不同, ...

  3. POJ-3450 Corporate Identity (KMP+后缀数组)

    Description Beside other services, ACM helps companies to clearly state their “corporate identity”, ...

  4. POJ3450 Corporate Identity 【后缀数组】

    Corporate Identity Time Limit: 3000MS   Memory Limit: 65536K Total Submissions: 7662   Accepted: 264 ...

  5. POJ3450 Corporate Identity

    后缀数组. 解决多个字符串的最长公共子串. 采用对长度的二分,将子串按height分组,每次判断是否在每个字符串中都出现过. 复杂度O(NlogN) By:大奕哥 #include<cstrin ...

  6. [HDU2328]Corporate Identity（后缀数组）

    传送门 求 n 个串的字典序最小的最长公共子串. 和 2 个串的处理方法差不多. 把 n 个串拼接在一起,中间连上一个没有出现过的字符防止匹配过界. 求出 height 数组后二分公共子串长度给后缀数 ...

  7. POJ3080 POJ3450Corporate Identity（广义后缀自动机||后缀数组||KMP）

    Beside other services, ACM helps companies to clearly state their “corporate identity”, which includ ...

  8. hdu2328 Corporate Identity【string库使用】【暴力】【KMP】

    Corporate Identity Time Limit: 9000/3000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Other ...

  9. kuangbin带你飞 后缀数组 题解

    2份模板 DC3 . 空间复杂度O3N 时间复杂度On #define F(x) ((x) / 3 + ((x) % 3 == 1 ? 0 : tb)) #define G(x) ((x) < ...

随机推荐

  1. js jquery 插件

    $(function(){ (function($, document, undefiend){ $.fn.pagination = function(options){ var $this = $( ...

  2. XP，32/64位Win7,32/64位Win8，32/64位Win10系统 【春节版】

    本系统是10月5日最新完整版本的Windows10 安装版镜像,win10正式版,更新了重要补丁,提升应用加载速度,微软和百度今天宣布达成合作,百度成为win10 Edge浏览器中国默认主页和搜索引擎 ...

  3. 【BZOJ4281】[ONTAK2015]Związek Harcerstwa Bajtockiego LCA

    [BZOJ4281][ONTAK2015]Związek Harcerstwa Bajtockiego Description 给定一棵有n个点的无根树,相邻的点之间的距离为1,一开始你位于m点.之后 ...

  4. VS2013 自动添加头部注释 -C#开发

    在团队开发中,头部注释是必不可少的.但在开发每次新建一个类都要复制一个注释模块也很不爽,所以得想个办法让开发工具自动生成我们所需要的模板.....操作方法如下: 方法/步骤 1 找你的vs安装目录, ...

  5. Error524 源站处理超时 Error 524: A timeout occurred

    https://su.baidu.com/helps/index.html#/4/5a61e4b5b34f697f13234a5b Error524 源站处理超时 更新时间:2018-01-19 20 ...

  6. squared-error loss is much more repaidly updated than mean-absolute-deviation when searching for splits

    平方差损失能较绝对值差损失更快地更新

  7. 【题解】CF264B Good Sequences

    [题解]CF264B Good Sequences 具有很明显的无后效性. 考虑\(dp\). 考虑初始条件,显然是\(dp(0)=0\) 考虑转移,显然是\(dp(t)=max(dp[k])+1\) ...

  8. Django的基础操作总结

    1:准备开始 建立一个新的project: django-admin.py startproject XXXXXX(名称) 建立一个新的App:python manage.py startapp XX ...

  9. 不使用 spring-boot-starter-parent 构建 spring boot 应用

    创建 spring-boot 应用通用方法是配置 pom.xml,定义 为 spring-boot-start-parent.如下: <parent> <groupId>org ...

  10. c# XML-Object对象 序列化-反序列化

    using System; using System.Collections.Generic; using System.IO; using System.Linq; using System.Tex ...