传送门:POJ - 2774  (最长相同子串)

题意:(在vj上看到了一个很有意思的描述)这个版本的

L学长喜欢上Z学妹,他发现他们的相似度很高,担心会不会就那么巧合,所以就想办法弄到了Z学妹的基因,然后也把自己的基因一起拿去比对,希望找出基因中完全一样的一段的最大长度,来判断要不要去德国骨科。

 

题解:把两个字符串拼一起,中间加一个分隔的字符(字符串里不会出现的字符)。为什么要加分隔字符呢,因为可能会出现这种结果,前一个字符串的尾部和后一个的首部连一起形成的字符串,和后一个剩下的某部分正好一样,导致判断难度加大。加了分隔符就相当于直接剪断了他们相等的可能。然后求出height数组,遍历 height数组 ,先判断height值是否比当前的最大值大,然后判断此时的sa[i-1]是否小于,sa[i]是否大于第一个字符串的长度,如果都满足就可以将最大值更新。

最近迷上了二分的做法,然后发现二分多此一举了,,,

 

  1 #include<cstdio>
  2 #include<algorithm>
  3 #include<queue>
  4 #include<iostream>
  5 #include<cmath>
  6 #include<cstring>
  7 using namespace std;
  8 
  9 const int INF=0x3f3f3f3f;
 10 const int maxn = 1e6+10;
 11 int wa[maxn],wb[maxn],wsf[maxn],wv[maxn],sa[maxn];
 12 int rnk[maxn],height[maxn];
 13 char s[maxn],t[maxn];
 14 int r[maxn],n,ans=0,len1,len2;
 15 
 16 //sa:字典序中排第i位的起始位置在str中第sa[i]  sa[1~n]为有效值
 17 
 18 //rnk:就是str第i个位置的后缀是在字典序排第几 rnk[0~n-1]为有效值
 19 
 20 //height:字典序排i和i-1的后缀的最长公共前缀  height[2~n]为有效值,第二个到最后一个
 21 
 22 int cmp(int *r,int a,int b,int k)
 23 {
 24     return r[a]==r[b]&&r[a+k]==r[b+k];
 25 }
 26 
 27 void getsa(int *r,int *sa,int n,int m)//n为添加0后的总长
 28 {
 29     int i,j,p,*x=wa,*y=wb,*t;
 30     for(i=0; i<m; i++)  wsf[i]=0;
 31     for(i=0; i<=n; i++)  wsf[x[i]=r[i]]++;
 32     for(i=1; i<m; i++)  wsf[i]+=wsf[i-1];
 33     for(i=n; i>=0; i--)  sa[--wsf[x[i]]]=i;
 34     p=1;
 35     j=1;
 36     for(; p<=n; j*=2,m=p){
 37         for(p=0,i=n+1-j; i<=n; i++)  y[p++]=i;
 38         for(i=0; i<=n; i++)  if(sa[i]>=j)  y[p++]=sa[i]-j;
 39         for(i=0; i<=n; i++)  wv[i]=x[y[i]];
 40         for(i=0; i<m; i++)  wsf[i]=0;
 41         for(i=0; i<=n; i++)  wsf[wv[i]]++;
 42         for(i=1; i<m; i++)  wsf[i]+=wsf[i-1];
 43         for(i=n; i>=0; i--)  sa[--wsf[wv[i]]]=y[i];
 44         swap(x,y);
 45         x[sa[0]]=0;
 46         for(p=1,i=1; i<=n; i++)
 47             x[sa[i]]=cmp(y,sa[i-1],sa[i],j)? p-1:p++;
 48     }
 49 }
 50 
 51 void getheight(int *r,int n)//n为添加0后的总长
 52 {
 53     int i,j,k=0;
 54     for(i=1; i<=n; i++)  rnk[sa[i]]=i;
 55     for(i=0; i<n; i++){
 56         if(k)
 57             k--;
 58         else
 59             k=0;
 60         j=sa[rnk[i]-1];
 61         while(r[i+k]==r[j+k])
 62             k++;
 63         height[rnk[i]]=k;
 64     }
 65 }
 66 
 67 bool check(int k)
 68 {
 69     bool flag=0;
 70     int mx=-INF,mn=INF;
 71     for(int i=2;i<=n;i++){
 72         if(height[i]>=k){
 73             mn=min(mn,min(sa[i],sa[i-1]));
 74             mx=max(mx,max(sa[i],sa[i-1]));
 75             if(mn<len1&&mx>len1) return true;
 76         }
 77         else mx=-INF,mn=INF;
 78     }
 79     return false;
 80 }
 81 
 82 int main()
 83 {
 84     ios::sync_with_stdio(false);
 85     cin.tie(0);
 86     cout.tie(0);
 87     cin>>s;
 88     cin>>t;
 89     n=0;
 90     len1=strlen(s);
 91     len2=strlen(t);
 92     for(int i=0;i<len1;i++)
 93         r[n++]=s[i]-'a'+1;
 94     r[n++]=0;
 95     for(int i=0;i<len2;i++){
 96         r[n++]=t[i]-'a'+1;
 97     }
 98     r[n]=0;
 99     getsa(r,sa,n,150);
100     getheight(r,n);
101     int ans=0;
102     for(int i=2;i<=n;i++){
103         if(height[i]>ans){
104             if(sa[i-1]>len1&sa[i]<len1) ans=height[i];
105             if(sa[i]>len1&&sa[i-1]<len1) ans=height[i];
106         }
107     }
108     cout<<ans<<endl;
109     return 0;
110 }