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