EX—KMP

一、基本概念:
   **问题定义:**给定两个字符串 S 和 T(长度分别为 n 和 m),下标从 0 开始,
定义extend[i]等于S[i]…S[n-1]与 T 的最长相同前缀的长度,求出所有的extend[i]。

   来看下next[i]与extend[i]的定义:
   next[i]: T[i]…T[m - 1]与 T 的最长相同前缀长度;
   extend[i]: S[i]…S[n - 1]与 T 的最长相同前缀长度。
   求解next[i]的过程就是 T 自己和自己的一个匹配过程。

#include <iostream>
#include <string>
using namespace std;

/* 求解 T 中 next[],注释参考 GetExtend() */
void GetNext(string& T, int& m, int next[])
{
    int a=0,p=0;
    next[0]=m;

    for (int i=1;i<m;i++)
    {
        if (i>=p||i+next[i-a]>=p)
        {
            if (i>=p)
                p=i;

            while (p<m&&T[p]==T[p - i])
                p++;

            next[i]=p-i;
            a=i;
        }
        else
            next[i]=next[i-a];
    }
}

/* 求解 extend[] */
void GetExtend(string& S, int& n, string& T, int& m, int extend[], int next[])
{
    int a=0,p=0;
    GetNext(T, m, next);

    for (int i=0;i<n;i++)
    {
        if (i>=p||i+next[i-a]>=p) // i >= p 的作用:举个典型例子,S 和 T 无一字符相同
        {
            if (i>=p)
                p=i;

            while (p<n&&p-i<m&&S[p]==T[p-i])
                p++;

            extend[i]=p-i;
            a=i;
        }
        else
            extend[i]=next[i-a];
    }
}

int main()
{
    int next[100];
    int extend[100];
    string S, T;
    int n, m;
    
    while (cin >> S >> T)
    {
        n = S.size();
        m = T.size();
        GetExtend(S, n, T, m, extend, next);

        // 打印 next
        cout << "next:   ";
        for (int i = 0; i < m; i++)
            cout << next[i] << " ";
 
        // 打印 extend
        cout << "\nextend: ";
        for (int i = 0; i < n; i++)
            cout << extend[i] << " ";

        cout << endl << endl;
    }
    return 0;
}

二、例题: CODEVS——T 1404 字符串匹配
题目描述 Description
给你两个串A,B,可以得到从A的任意位开始的子串和B匹配的长度。
给定K个询问,对于每个询问给定一个x,求出匹配长度恰为x的位置有多少个。
N,M,K<=200000

输入描述 Input Description
第一行三个数 N,M,K,表示A的长度、B的长度和询问数。
第二行为串A。
第三行为串B。
接下来K行,每行1个数X。

输出描述 Output Description
对于每个询问输出一个数。

样例输入 Sample Input
6 2 2
aabcde
ab
0
2

样例输出 Sample Output
4
1

数据范围及提示 Data Size & Hint
各个测试点1s

#include<iostream>
#include<cstring>
#include<string>
#include<cmath>
#include<algorithm>
#include<cstdio>
#include<vector>
#define ll long long
using namespace std;
const int maxn=200010;
char str1[maxn],str2[maxn];
int nt[maxn];
int ex[maxn];
int b[maxn]={0};

int main(void)
{
    int n,m,k;
    int x;
    scanf("%d%d%d",&n,&m,&k);
    scanf("%s%s",str1,str2);
    int a=0,p=0;
    nt[0]=m;
    for(int i=1;i<m;i++)
    {
        if(i>=p||i+nt[i-a]>=p)
        {
            if(i>=p)
                p=i;

            while(p<m&&str2[p]==str2[p-i])
                p++;

            nt[i]=p-i;
            a=i;
        }
        else
            nt[i]=nt[i-a];
    }

    a=0,p=0;
    for(int i=0;i<n;i++)
    {
        if(i>=p||i+nt[i-a]>=p)
        {
            if(i>=p)
                p=i;

            while(p<n&&p-i<m&&str1[p]==str2[p-i])
                p++;

            ex[i]=p-i;
            a=i;
        }
        else
            ex[i]=nt[i-a];

        b[ex[i]]++;
    }

    int g;
    for(int i=1;i<=k;i++)
    {
        scanf("%d",&g);
        printf("%d\n",b[g]);
    }
    return 0;

}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值