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);    }       }