hdu4300Clairewd’s message(kmp||exkmp)( 模板 )

题目链接
给你一个明文to密文的变化表,然后给你一行字符串,前面输入的是密文后面接明文,但是明文不一定是完整的,求最短要发送的字符串(密文+明文).
假设我们知道字符串以某位开始的后缀与另一字符串的最长公共前缀长ex[],
那么设s[]为原串,t[]为转换回去的串,ex[i]为s[i]和t的最长公共前缀长度
那么密文和明文的分界点一定是满足i+ex[i]==l(字符串从0计),而最短则是顺序扫一遍找到最小的i,注意i也不能小于L/2。然后就是扩展kmp的模板题了
关于exkmp:
首先这个算法能够求出ex[i],也即主串的任意一个后缀与子串的最长公共前缀,看个ppt什么的学起来会快一些,转移的原理大概是这个样子的:
假设p为最右匹配过的位置,a为匹配到q时的起点,那么有这样的
s[a..p]=t[0…p-a]
s[a…k…p]=t[0…k-a…p-a]
s[k…p]=t[k-a…p-a]
l=nxt[k-a]
if(l+k < p+1)ex[k]=l
if(l+k > = p+1){
cmp( s[p+1] , t[p-k+1]);
update(a);
}

#include <stdio.h>
#include<cstring>
#include<iostream>
using namespace std;
char s[100010];
char t[100010];
int ls,lt;
int nxt[100010];
int ex[100010];
void ekmp(){
    int  i,j,l,p;
    memset(nxt,0,sizeof(nxt));
    memset(ex,0,sizeof(ex));
    ls=strlen(s);
    lt=strlen(t);

    nxt[0]=lt,j=0;
    while(j+1<lt&&t[j]==t[1+j])j++;
    nxt[1]=j;

    int a=1;
    for(int i=2;i<lt;i++){
        p=a+nxt[a]-1;
        l=nxt[i-a];
        if(l+i<p+1)nxt[i]=l;
        else {
            j=max(0,p-i+1);
            while(i+j<lt&&t[i+j]==t[0+j])j++;
            nxt[i]=j;
            a=i;
        }
    }

    j=0;
    while(j<ls&&j<lt&&s[j]==t[j])j++;
    ex[0]=j;

    a=0;
    for(int i=1;i<ls;i++){
        p=a+ex[a]-1;
        l=nxt[i-a];
        if(l+i<p+1)ex[i]=l;
        else {
            j=max(0,p-i+1);
            while(i+j<ls&&j<lt&&s[i+j]==t[j])j++;
            ex[i]=j;
            a=i;
        }
    }
}
char tran[30];
int main(void)
{
    int n;
    scanf("%d",&n);
    while(n--){
        scanf("%s%s",tran,s);
        int l=strlen(s);
        for(int i=0;i<l;i++){
            for(int j=0;j<26;j++){
                if(s[i]==tran[j])t[i]=j+'a';
            }
        }
        ekmp();
        int k;
        for( k=0;k<l;k++){
            if(k+ex[k]==l&&ex[k]<=k)break;
        }
        for(int i=0;i<k;i++)printf("%c",s[i]);
        for(int i=0;i<k;i++)printf("%c",t[i]);
        puts("");
    }
    return 0;
}

当然这题有一个更加巧妙地解法
把还原过后的串的前半部分(上取整)拿去和原串的后半部分匹配找到的第一个能够使得原串匹配结束的地方即可,并将后面的部分输出

#include <stdio.h>
#include<cstdio>
#include<set>
#include<cstring>
#include<iostream>
#include<cstdlib>
using namespace std;
char a[100010];
char b[100010];
int la,lb;
int nxt[100010];
void gnx(){
    memset(nxt,0,sizeof(nxt));
    int i=0,j=-1;
    nxt[0]=-1;
    while(i<lb){
        if(j==-1||b[i]==b[j])nxt[++i]=++j;
        else j=nxt[j];
    }
}
int kmp(){
    int i=0,j=0;
    la=strlen(a);
    while(i<la&&j<lb){
        if(j==-1||a[i]==b[j]){
            i++;
            j++;
        }
        else j=nxt[j];
        if(i==la)return j;
    }
    return 0;
}
char s[30];
int main(void)
{
    int n;
    scanf("%d",&n);
    while(n--){
        scanf("%s%s",s,b);
        lb=strlen(b);
        strcpy(a,b+(lb+1)/2);
        printf("%s",b);
        for(int i=0;i<lb;i++){
            for(int j=0;j<26;j++){
                if(b[i]==s[j]){b[i]='a'+j;break;}
            }
        }
        gnx();
        int kk=kmp();
        for(int i=kk;i<lb-kk;i++){
            printf("%c",b[i]);
        }
        puts("");
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值