【后缀数组】HDU-1403

/* 最长公共子串的问题
 *
 * 
 * File:   main.cpp
 * Author: lenovo
 *
 * Created on 2011年10月18日, 上午9:01
*/

#include <cstdlib>
#include <cstdio>
#include <cstring>
#define maxn 200200
using namespace std;

/*
 * 求A和B的最长公共子串等价于求A的后缀和B的后缀的最长公共前缀的最大值
 * 用到了height数组(排名相邻的两个后缀的最长公共前缀)
 * 不是height值中的最大值就是答案
 * 要保证suffix(sa[i-1])和suffix(sa[i])不是同一个字符串的两个后缀
 * 这里就用到了check函数
 *
 *
 * 
*/


int wa[maxn] , wb[maxn] , wv[maxn] , ws[maxn];
int rank[maxn] , height[maxn] , sa[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 ;
}

void callheight(char *r , int *sa , int n)   //构造height数组
{
int i , j , k = 0;
for(i = 1 ; i <= n ; i++)
        rank[sa[i]] = i;
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 ;
}

int check(int i , int n , int *sa)  //判断两个后缀是不是属于同一个字符串
{
int c1 , c2;
if(sa[i] <= n)                  //直接看排名的值   如果都是小于第一个字符串的长度的话  显然在同一个子串里面
        c1 = 1;
else
        c1 = 0;
if(sa[i-1] <= n)
        c2 = 1;
else
        c2 = 0;
if(c1 != c2)
return 1;
else
return 0;
}

int main(int argc, char** argv) {
char r[maxn] , b[maxn];
int i , j;
while(scanf("%s%s",r,b) != EOF)
    {
int lenr = strlen(r);
int lenb = strlen(b);
int sum_len = lenr+lenb+1;
        r[lenr] = '$';
for(i = 0 ; i < lenb ; i++)
            r[i+lenr+1] = b[i];
//for(i = 0 ; i < sum_len ; i++)
//    printf("%c",r[i]);
        da(r,sa,sum_len,'{');
        callheight(r,sa,sum_len);
//for(i = 0 ; i < sum_len ; i++)
//    printf("%d ",sa[i]);
        int max = 0;
for(i = 1 ; i < sum_len ; i++)
if(height[i] > max && check(i , lenr , sa))
                max = height[i];
        printf("%d\n",max);
    }
return 0;
}

转载于:https://www.cnblogs.com/zuckerTT/archive/2011/10/18/2216155.html