【poj3461】kmp裸题

poj3461

直接看样例:T组数据,每组包括两个字符串AB,求A串在B串所有匹配位置的总数。

KMP裸题……那天刚学KMP于是拿来练手。

【想看代码的直接拉下去】

KMP是一个字符串匹配的算法,可以在O(n)的时间内求出一个串在另一个串的所有匹配位置。
暴力打法是在检查的时候,如果失配则往后移一位再进行匹配…KMP算法认为这样浪费了一些信息。因为若在第k位失配,我们已经知道[0,k-1]位的信息,却只往后移一位进行匹配,显然是很浪费的。
比如匹配abcba时在第四位失配,那么我们已经知道那几位是abcb,然后拿a和b来匹配,这个匹配显然是错误的。

所以说,我们可以用一个nxt数组表示第I位失配后跳到nxt[i]位。
而nxt数组的求法则是自己和自己匹配。

其实KMP主要就是背代码,基本上没法拓展了…我背的大白上的。

下面是代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;


int nxt[1000010];

void getnxt(char *a)
{
    memset(nxt,0,sizeof(nxt));
    int l=strlen(a);
    nxt[1]=nxt[0]=0;
    for(int i=1;i<l;i++)
    {
        int j=nxt[i];
        while(j&&a[i]!=a[j]) j=nxt[j];
        nxt[i+1]= a[i]==a[j]? j+1 :0;
    }
}

int kmp(char *a,char *b)
{
    int n=strlen(a),m=strlen(b);
    getnxt(b);
    int ans=0;
    for(int i=0,j=0;i<n;i++)
    {
        while(j&&a[i]!=b[j]) j=nxt[j];
        if(a[i]==b[j]) j++;
        if(j==m) ans++;
    }
    return ans;
}

char a[1000010],b[1000010];
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%s%s",a,b);
        printf("%d\n",kmp(b,a));
    }
    return 0;
}


2016.2.24

今天又来敲了一遍…这东西敲的少啊QAQ

还有就是有个新发现,strlen()好像不是O(1)的…别问我怎么知道的,T了两次…

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

const int SZ = 1000010;

int nxt[SZ];

void getnxt(char P[])
{
    memset(nxt,0,sizeof(nxt));
    nxt[0] = nxt[1] = 0;
    int n = strlen(P);
    for(int i = 1;i < n;i ++)
    {
        int j = nxt[i];
        while(j && P[i] != P[j]) j = nxt[j];
        nxt[i + 1] = P[i] == P[j] ? j + 1 : 0;
    }
}

int ask(char P[],char T[])
{
    getnxt(P);
    int ans = 0;
    int m = strlen(P),n = strlen(T);
    for(int i = 0,j = 0;i < n;i ++)
    {
        while(j && P[j] != T[i]) j = nxt[j];
        if(P[j] == T[i]) j ++;
        if(j == m) ans ++;
    }
    return ans;
}

char s1[SZ],s2[SZ];

int main()
{
    int T;
    scanf("%d",&T);
    while(T --)
    {
        scanf("%s%s",s1,s2);
        printf("%d\n",ask(s1,s2));
    }
    return 0;
}

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值