hdu 4300 Clairewd’s message 扩展kmp

S数组 母串

T数组 模式串

A数组 模式串自己的每一个后缀与其本身的最长公共前缀 T[i...i+k]== T[0,,,k]   则 A[i]= k+1;

B数组 母串的每一个后缀与模式串的最长公共前缀     S[i...i+k]== T[0...k]  则  B[i]= k+1;

如果 strlen(T) == B[i]  则模式串一定在母串中出现过，i即起出现的位置，所以扩展kmp为kmp算法的一个扩展

 

hdu 4300 Clairewd’s message

http://acm.hdu.edu.cn/showproblem.php?pid=4300

 

题意为：给定一个翻译表，即第i个字母用哪个字母表示

再给一个串，里面前面为密文，后面为明文，密文一定是完整的，但明文不完整或可能没有

求这个完整的前面密文后面明文的串

 

 

把给定的串全部按表翻译为密文，这样例如样例
qwertabcde

翻译之后为：

 jvtkzqwert

可发现翻译后的后半部分与原传相同，把原串作为T串，翻译后的为S串，求B

题目要求求最短长度， 从B的一半开始， 当找到第一个满足 B[i]== len- i条件的，i即为第一个明文的位置

 

#include <iostream>#include <cstdio>#include <cstring>using namespace std;char mi[30], ming[30]; //mi 下标密文，译为明文， ming  下标明文，译为密文char S[100010], T[100010];int len, A[100010], B[100010];void get_A(int m){   //求A数组，m= strlen（T）    int i, j, k, len, L;    j= 0;    while( T[0+j] == T[1+j])        j= j+1;    A[0]= 0;    A[1]= j;    k= 1;    for( i=2; i<m; i++){        len= k+ A[k] -1;        L= A[i-k];        if( L<len -i+1)            A[i]= L;        else{            j= max(0, len-i+1);            while( T[i+j]== T[0+j])                j= j+1;            A[i]= j;            k= i;        }    }}void get_B(int m, int n){ //求B数组 m= strlen（T），n= strlen（S）    int i, j, k, len, L;    j= 0;    while( T[0+j]== S[0+j] && j<m && j<n ){        j= j+1;    }    B[0]= j;    k= 0;    for( i=1; i<n; i++){        len= k+B[k]-1;        L= A[i-k];        if( L< len-i+1)            B[i]= L;        else{            j= max(0, len-i+1);            while( S[i+j]== T[0+j])                j= j+1;            B[i]= j;            k= i;        }    }}int main(){  // freopen("1.txt", "r", stdin);    int i, j, k;    int TT;    scanf("%d", &TT);    while( TT--){        scanf("%s", ming);        scanf("%s", T);        for( i=0; i<26; i++){            mi[ming[i]-'a']= i+'a';        }        len= strlen(T);        for( i=0; i<len; i++){            S[i]= ming[T[i]-'a'] ;        }        S[len] ='\0';     //注意，S数组为手动赋值得，最后一位加'\0'        get_A(len);        get_B(len, len);/*        for( i=0; i<len; i++)            cout<<T[i]<<"  ";        cout<<endl;        for( i=0; i<len; i++)            cout<<S[i]<<"  ";        cout<<endl;        for( i=0; i<len; i++)            cout<<A[i]<<"  ";        cout<<endl;        for( i=0; i<len; i++)            cout<<B[i]<<"  ";        cout<<endl;*/        // 题目要求求最短长度， 从串的一半开始，        //当找到第一个满足 B[i]== len- i条件的，i即为第一个明文的位置        for( i= (len+1)/2; i<len; i++){            if( B[i]== len- i) break;        }        for( j=0; j<i; j++){            T[j+i]= mi[ T[j]-'a'];        }        T[j+i]='\0';        printf("%s\n", T);    }    return 0;}