HDU 1403 Longest Common Substring(尾碼數組 最長公用子串)

題目串連：http://acm.hdu.edu.cn/showproblem.php?pid=1403另外POJ 的這個題目用相同的代碼也能AC，兩個題目一樣：http://poj.org/problem?id=2774 題意：給定兩個串，要求輸出兩個串的最長公用子串的長度解題思路：尾碼數組，在height數組上做文章其實這個題目和上面兩篇部落格的演算法思路基本上是一致的，首先把兩個字串拼接成一個字串，在拼接之前記錄下拼接的位置，方便後續處理完成之後二分答案，二分的過程中要判斷每個當前要判斷的答案是否滿足題目條件，首先還是和上面兩題目一樣，按照height值分組，這個分組的依據是height值大於m，然後記錄下這個區間的最左邊和最右邊，只要左邊的在分界線左邊，右邊的在分界線右邊，而且左邊的+上m還在分界線的左邊（因為不同串之間的最長公用字串是不可能重疊的），滿足這個條件就滿足題目要求的當前的答案，這個題目也算搞定了其實二分和對height的值進行分組是在做尾碼數組的時候常用的技巧，想用尾碼數組解決問題這兩個技巧就一定要會，否則就是面對大多數題目，有種無從下手的感覺！

#include <iostream>#include <stdio.h>#include <string.h>#include <algorithm>#include <math.h>using namespace std;#define maxn 210000#define ws ws1int wa[maxn],wb[maxn],wv[maxn],ws[maxn];int cmp(int *r,int a,int b,int l){return r[a]==r[b]&&r[a+l]==r[b+l];}void da(char *r,int *sa,int n,int m){int i,j,p,*x=wa,*y=wb,*t;for(i=0;i<m;i++) ws[i]=0;for(i=0;i<n;i++) ws[x[i]=r[i]]++;for(i=1;i<m;i++) ws[i]+=ws[i-1];for(i=n-1;i>=0;i--) sa[--ws[x[i]]]=i;for(j=1,p=1;p<n;j*=2,m=p){for(p=0,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<n;i++) wv[i]=x[y[i]];for(i=0;i<m;i++) ws[i]=0;for(i=0;i<n;i++) ws[wv[i]]++;for(i=1;i<m;i++) ws[i]+=ws[i-1];for(i=n-1;i>=0;i--) sa[--ws[wv[i]]]=y[i];for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++)x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;}return;}int rank[maxn],height[maxn];void calheight(char *r,int *sa,int n){int i,j,k=0;for(i=1;i<=n;i++) rank[sa[i]]=i;//在計算height的時候順便就把rank計算出來了，反正也要用for(i=0;i<n;height[rank[i++]]=k)for(k?k--:0,j=sa[rank[i]-1];r[i+k]==r[j+k];k++);return;}char rec[maxn],rec1[maxn];int sa[maxn];int n;bool check(int m,int len){    int mmax=0,mmin=n;    for(int i=1;i<=n;i++)    {        if(height[i]<m)        {            mmax=sa[i];            mmin=sa[i];        }        else        {            mmax=max(mmax,max(sa[i],sa[i-1]));            mmin=min(mmin,min(sa[i],sa[i-1]));            if(mmax>=len && mmin+m <len)            return true;            //if(mmax-mmin>m) return true;        }    }    return false;}int main(){    int i,j,k;    int left,right;    int mid;    int answer;    int lena,lenb,lent;    while(scanf("%s%s",rec,rec1)!=EOF)    {        lena=strlen(rec);        lenb=strlen(rec1);        lent=lena+lenb;        strcat(rec,rec1);        n=lent;     da(rec,sa,n+1,200);//這裡是n+1 因為看這個函數裡面是 < n 的     calheight(rec,sa,n);//注意這裡面是 n 了因為看函數裡面是 <=n 的，所以這裡要注意     left=0;     right=100000;     while(left <= right)     {         mid=(left+right)>>1;         if(check(mid,lena))         {             answer=mid;             left=mid+1;         }         else         {             right=mid-1;         }     }     printf("%d\n",answer);    }    return 0;}