poj 1226 求多串最长公共子串 或 回文子串 KMP&&strstr&&后缀数组

题意：求给出的字符串中的最长公共子串 ，其中子串可以反转，比如rose 和 orchid ， ro 和第一个匹配 ， ro 的回文串or 和orchid 匹配 。最后输出最长子串的长度。

思路：找出长度最小的子串， 枚举其所有子串 ，然后将其反转再保存 ，最后用这两个串去匹配其它的字符串。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
char str[52][505];
char ss[305] , ts1[305] , ts2[305];
int next[205];
void get_next(char *T , int len) {      /*****得到next数组*****/ 
    int j = 0, k = -1;
    next[0] = -1;
    while (j < len) {
        if(T[j] == T[k] || k == -1) {
            next[j+1] = k+1;
            j++;k++;    
        }
        else k = next[k];    
    }    
}
int kmp(char *T , int h , int len2 , int len1) { /*****KMP*****/
    int i , j;
    i = j = 0; 
    while(i <len1 && j <len2) {
        if(T[j]==str[h][i] || j == -1) {
            j++;i++;    
        }else {
        j = next[j];
        }
    } 
    if(j == len2) return 1;
    return 0;   
}
int main() {
    int T , i , j , k , len , n , g , temp;
    scanf("%d",&T);
    while(T--) {
        scanf("%d",&n);
        for(i = 0 ; i < n ; i ++)
        scanf("%s",str[i]);
        len = 120;
        for(i = 0 ; i < n ; i ++) {         /*****找长度最小的字符串*****/ 
            if(strlen(str[i]) <= len) {
                len = strlen(str[i]);
                temp = i;    
            }    
        }
        int tt , ans , cnt , l1 , l2 , lenn , f , flag , r , t;
        ans = 0;
        for(i = 1 ; i <= len ; i ++) {
            for(j = 0 ; j <= len - i; j ++) {
                for(g = 0 ,k = j ; k < j + i ; k ++ , g++) { /*****暴力枚举长度是i的最小字符串子串ts1*****/ 
                    ts1[g] = str[temp][k];            
                }
                ts1[g] = '\0';
               for(tt = 0 ; tt < g ; tt ++) {                 /*****求这个子串的反转子串ts2*****/ 
                    ts2[tt] = ts1[g - tt - 1];    
                }
                ts2[tt] = '\0';
                get_next(ts1 , i);                             /*****得到ts1的next数组*****/ 
                cnt = 0;
                bool vis[100];                             /*****标记数组*****/ 
                memset(vis , false , sizeof(vis));
                for(f = 0 ; f < n; f ++) {
                    lenn = strlen(str[f]);
                    if(kmp(ts1 , f , i , lenn) && !vis[f]) {
                        vis[f] = 1; cnt++;} 
                }
                get_next(ts2 , i);
                for(f = 0 ; f < n ; f ++) {
                    lenn = strlen(str[f]);
                    if(kmp(ts2 , f , i , lenn)&&!vis[f]) {
                        vis[f] = 1; cnt++;} 
                }                                              /*****分别用ts1 、 ts1 去匹配其它串*****/ 
                if(cnt == n) {                                 //满足条件中取最长的子串长度 
                    if(i > ans) ans = i;    
                }
            }    
        }  
        printf("%d\n",ans); 
    }    
}

strstr 函数水的~   注意strncpy枚举子串的方法 

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
char str[120][120];
char s1[120] , s2[120];
char zz1[120] , zz2[120];
int main() {
    int n , T , i , j , k , t , temp;
    scanf("%d",&T);
    while(T--) {
        scanf("%d",&n);
        int len = 120;
        for(i = 0 ; i < n ; i ++) {
            scanf("%s",str[i]);
            int ll = strlen(str[i]);
            if(ll < len) {
                len = ll;
                temp  = i;
            }  
        }
        strcpy(s1 , str[temp]);
        int ans = 0;
        for(i = 1 ; i <= len ; i ++) { //枚举子串长度
            for(j = 0 ; j <= len - i ; j ++) {
                strncpy(zz1, s1+j , i);
                zz1[i] = '\0';      //子串 
                //printf("%s\n",zz1);
                for(k = 0 ; k < i ; k ++)
                zz2[k] = zz1[i - k - 1];
                zz2[k] = '\0';     //反转子串 
                //printf("%s\n",zz2);
                for(t = 0 ; t < n ; t ++) {
                    if(strstr(str[t] , zz1)||strstr(str[t] , zz2)) continue;
                    else break;    
                }
                if(t == n) ans = i>ans?i:ans; 
            }   
        }
        printf("%d\n",ans);     
    }    
}

后缀数组的方法：和求多串公共子串差不多的方法，不过这题要连接的不只给出的字符串，也要把给出的子串的反转串连接起来，并且loc数组的状态要和正串是一样的。

最后二分枚举最短字符串的子串长度，判断即可。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn = 21000;
char str[maxn];
int num[maxn] , loc[maxn];
int sa[maxn] , rank[maxn] , height[maxn];
int wa[maxn] , wb[maxn] , wv[maxn] , wd[maxn];
int n;
int cmp(int *r , int a , int b , int l) {
    return r[a] == r[b] && r[a+l] == r[b+l];    
}
void da(int *r , int n , int m) {
    int i , j , p;
    int *x = wa , *y = wb , *t;
    for(i = 0 ; i < m ; i ++) wd[i] = 0;
    for(i = 0 ; i < n ; i ++) wd[x[i]=r[i]]++;
    for(i = 1 ; i < m ; i ++) wd[i] += wd[i-1];
    for(i = n-1 ; i >= 0 ; i --) sa[--wd[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 ++) wd[i] = 0;
        for(i = 0 ; i < n ; i ++) wd[wv[i]] ++;
        for(i = 1 ; i < m ; i ++) wd[i] += wd[i-1];
        for(i = n - 1 ; i >= 0 ; i --) sa[--wd[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 ++;        
        }   
    }    
}
void calheight(int *r , int n) {
    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 ++);    
    } 
}
int vis[1300];
bool check(int mid , int len) {
    int i , j , cnt;
    cnt = 0;
    memset(vis , 0 , sizeof(vis));
    for(i = 2 ; i <= len ; i ++) {
        //printf("%d\n",height[i]);
        if(height[i] < mid) {
            memset(vis,0,sizeof(vis));
            cnt = 0;
            continue;    
        }    
        if(!vis[loc[sa[i-1]]]) {
            vis[loc[sa[i-1]]] = 1;
            cnt ++;        
        }
        if(!vis[loc[sa[i]]]) {
            vis[loc[sa[i]]] = 1;
            cnt ++;    
        }
        if(cnt == n) return 1;
    }
    return 0;    
}
int main() {
    int T , i , j , g , k , temp;
    scanf("%d",&T);
    while(T--) {
        scanf("%d",&n);
        g = 0;
        temp = 130;
        int ll = 120;
        for(i = 0 ; i < n ; i ++) {
            scanf("%s",str);
            int len = strlen(str);
            if(len < ll) ll = len;
            for(j = 0 ; j < len ; j ++) {
                loc[g] = i;    
                num[g++] = str[j];
            }
            loc[g] = temp;
            num[g++] = temp++;
            for(j = 0 ; j < len ; j ++) {
                loc[g] = i;    
                num[g++] = str[len - j - 1];
            }
            loc[g] = temp;
            num[g++] = temp ++;    
        }   
        num[g] = 0;
        da(num , g+1 , temp);
        calheight(num , g);
        int left , right , mid , ans;
        left = 0;right = ll;
        while(left <= right) {
            mid = (left+right)/2;
            if(check(mid , g)) {
                left = mid + 1;
                ans = mid;        
            } else {
                right = mid - 1;
                //ans = mid;    
            }
        }
        printf("%d\n",ans);
    }       
}