KMP算法专题题解

已于 2023-03-28 17:10:10 修改
阅读量157 收藏
点赞数
文章标签： 算法 c++
于 2023-03-22 21:24:01 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_73674034/article/details/129694728
版权

A P3375 【模板】KMP字符串匹配

#include <iostream>
#include<map>
#include<set>
#include<vector>
#include<algorithm>
#include<string.h>
using namespace std;
char b[1000005], a[1000005]; int n, m, x, nex[1000005];
int main()
{
    while (cin >> b) {
        if (b[0] == '#') break;
        cin >> a;
        m = strlen(b);
        n = strlen(a);
        nex[1] = 0;
        int N = 0;
        for (int i = 2, j = 0; i <= n; i++) {
            while (j > 0 && a[i - 1] != a[j]) j = nex[j];
            if (a[i - 1] == a[j]) j++;
            nex[i] = j;
        }
        for (int i = 1, j = 0; i <= m; i++) {
            while (j > 0 && (j == n || b[i - 1] != a[j])) j = nex[j];
            if (b[i - 1] == a[j]) j++;
            if (j == n) {
                //开始位置
                cout << i-n+1 << endl;
            }
        }
        for (int i = 1; i <= n; i++) cout << nex[i]<<' ';
    }
    return 0;
}

B Period(HDU1358)

题目要求查找输入字符串所有存在循环元的前缀并返回循环次数(相当于求最小循环元)
设最小循环元长度为k,字符串长度为l,则k,l满足关系式l%(l-k)=0,循环次数为l/(l-k).值得注意的是(l-k)就对应着前后缀匹配上的长度,每次查看长度是否符合关系式即可

#include <iostream>
#include<map>
#include<set>
#include<vector>
#include<algorithm>
#include<string.h>
using namespace std;
char a[1000005]; int n,m,nex[1000005];
int main()
{
    while (cin >> n) {
        if (!n) break;
        m++;
        cin >> a;
        cout << "Test case #" << m << endl;
        for (int i = 2, j = 0; i <= n; i++) {
            while (j > 0 && a[i - 1] != a[j]) j = nex[j];
            if (a[i - 1] == a[j]) j++;
            //出现循环元
            if (j && i % (i - j) == 0) {
                cout << i << ' ' << i / (i - j)<<endl;
            }
        }
        cout << endl;
    }
    return 0;
}

C P4391 [BOI2009]Radio Transmission 无线传输

不需要知道s2,我们可以直接把s1的开头看成s2的开头,设s2’为s2的子串,s1满足:
s1=(s2*k)s2’(k=0,1,2,…)
设s2长为x,s1长为n,则s2=n-next[n]

#include <iostream>
#include<map>
#include<set>
#include<vector>
#include<algorithm>
#include<string.h>
using namespace std;
char a[1000005]; int n,N=1,nex[1000005];
int main()
{
    cin >>n>> a;
    for (int i = 2, j = 0; i <= n; i++) {
        while (j > 0 && a[i - 1] != a[j]) j = nex[j];
        if (a[i - 1] == a[j]) j++;
        nex[i] = j;
    }
    cout << n-nex[n];
    return 0;
}

D P1470 [USACO2.3]最长前缀 Longest Prefix

由kmp得出P中每个元素在序列s中出现的位置,再依靠dp求得能拼凑出的最长前缀

#include <iostream>
#include<map>
#include<set>
#include<vector>
#include<algorithm>
#include<string.h>
using namespace std;
char a[205][15];string x,b=""; int N, n, m, nex[1000005], dp[200005] = {1};
vector<int> t[200005];
void kmp(char a[]) {
    n = strlen(a);
    for (int i = 2, j = 0; i <= n; i++) {
        while (j > 0 && a[i - 1] != a[j]) j = nex[j];
        if (a[i - 1] == a[j]) j++;
        nex[i] = j;
    }
    for (int i = 1, j = 0; i <= m; i++) {
        while (j > 0 && (j == n || b[i - 1] != a[j])) j = nex[j];
        if (b[i - 1] == a[j]) j++;
        if (j == n) t[i].push_back(n);
    }
}
int main()
{
    int i = 0;
    while (cin >> a[i]) 
        if (a[i++][0] == '.')break;
    //注意b可能不只有一行
    while (cin >> x) b += x;
    m = b.size();
    for (int j = 0; j < i - 1; j++) kmp(a[j]);
    for (int i = 1; i <= m; i++) {
        for (int j = 0; j < t[i].size(); j++)
            if (dp[i - t[i][j]]) dp[i] = 1;
        if (dp[i]) N = max(N, i);
    }
    cout << N;
    return 0;
}

E Theme Section (hdu-4763)

S=EAEBE,子字符串E长度所有候选项为(next[m],next[next[m]],…),依次检验候选项是否会存在于字符串S的中间部分,检验成功则返回此时的候选项大小

#include <iostream>
#include<map>
#include<set>
#include<vector>
#include<algorithm>
#include<string.h>
using namespace std;
char a[1000005]; int  nex[1000005],m;
bool check(int n) {
    //检验在字符串的中间部分是否存在对应的子字符串
    for (int i = n+1, j = 0; i <= m-n; i++) {
        while (j > 0 && (j == n || a[i - 1] != a[j])) j = nex[j];
        if (a[i - 1] == a[j]) j++;
        if (j == n)  return 1;
    }
    return 0;
}
int main()
{
    std::ios::sync_with_stdio(0);
    int n, t; cin >> t;
    while (t--) {
        cin >> a; m = strlen(a);
        nex[1] = 0;
        for (int i = 2, j = 0; i <= m; i++) {
            while (j > 0 && a[i - 1] != a[j]) j = nex[j];
            if (a[i - 1] == a[j]) j++;
            nex[i] = j;
        }
        int N = nex[m];
        while (N) {
            if (check(N)) break;
            //若匹配失败,再利用nex数组得到剩下的最大"候选项"
            N = nex[N];
        }
        cout << N << endl;
    }
    return 0;
}

F P2375 [NOI2014] 动物园

先求一个前后缀可以重叠的num0数组,再仿照next数组的求解过程,由num0数组得出num数组

#include <iostream>
#include<map>
#include<set>
#include<vector>
#include<algorithm>
#include<string.h>
using namespace std;
char a[1000005]; int n, nex[1000005],num0[1000005], k = 1000000007;
int main()
{
    int t; cin >> t;
    while (t--) {
        cin >> a;
        n = strlen(a);
        nex[1] = 0; num0[1] = 1;
        for (int i = 2, j = 0; i <= n; i++) {
            while (j > 0 && a[i - 1] != a[j]) j = nex[j];
            if (a[i - 1] == a[j]) j++;
            nex[i] = j; num0[i] = num0[j] + 1;
        }
        int N = 1;
        for (int i = 2,j=0; i <= n; i++) {
            while (j > 0 && a[i - 1] != a[j]) j = nex[j];
            if (a[i - 1] == a[j]) j++;
            while (j>i/2) j = nex[j];
            N = (long long)(num0[j] + 1) * N % k;
        }
        cout << N << endl;
    }
    return 0;
}

G P3808【模板】AC 自动机（简单版）

#include <iostream>
#include <string>
#include <cstring>
using namespace std;
const int M = 1e6 + 5;
char a[M], t[M]; int q[M], nex[M], cnt[M], tr[M][28], N = 1; bool vis[M];
//字典树
void trie() {
    int T; cin >> T;
    while (T--) {
        cin >> a;
        int n = strlen(a), p = 0, ch;
        for (int i = 0; i < n; i++) {
            ch = a[i] - 'a';
            if (!tr[p][ch])tr[p][ch] = ++N;
            p = tr[p][ch];
        }
        cnt[p]++;
    }
}
//求"fail"数组
void build() {
    int hh = 0, tt = -1;
    for (int i = 0; i < 26; i++)
        if (tr[0][i]) q[++tt] = tr[0][i];
    while (hh <= tt) {
        int s = q[hh++];
        for (int i = 0; i < 26; i++) {
            int p = tr[s][i];
            if (!p) tr[s][i] = tr[nex[s]][i];
            else {
                nex[p] = tr[nex[s]][i];
                q[++tt] = p;
            }
        }
    }
}
int main()
{
    cin.tie(0)->sync_with_stdio(0);
    cout.tie(0)->sync_with_stdio(0);
    trie();
    build();
    int res = 0; cin >> t;
    for (int i = 0, j = 0; t[i]; i++) {
        int s = t[i] - 'a';
        j = tr[j][s];
        int p = j;
        //若p已出现过,不再累加
        while (p&&!vis[p]) {
            res += cnt[p];
            cnt[p] = 0; vis[p] = 1;
            p = nex[p];
        }
    }
    cout << res;
    return 0;
}

H P3796 【模板】AC 自动机（加强版）

#include <iostream>
#include <string>
#include <cstring>
using namespace std;
const int M = 1e6 + 5;
char a[155][71], t[M]; int T, q[M], nex[M], cnt[M], tr[M][28], N, vis[M];
void init() {
    memset(cnt, 0, sizeof(cnt));
    memset(vis, 0, sizeof(vis));
    memset(tr, 0, sizeof(tr));
    memset(nex, 0, sizeof(nex));
    N = 1;
}
void trie() {
    for (int i = 1; i <= T; i++) {
        cin >> a[i];
        int n = strlen(a[i]), p = 0, ch;
        for (int j = 0; j < n; j++) {
            ch = a[i][j] - 'a';
            if (!tr[p][ch])tr[p][ch] = ++N;
            p = tr[p][ch];
        }
        //离散化表示各模式串
        vis[p] = i;
    }
}
void build() {
    int hh = 0, tt = -1;
    for (int i = 0; i < 26; i++)
        if (tr[0][i]) q[++tt] = tr[0][i];
    while (hh <= tt) {
        int s = q[hh++];
        for (int i = 0; i < 26; i++) {
            int p = tr[s][i];
            if (!p) tr[s][i] = tr[nex[s]][i];
            else {
                nex[p] = tr[nex[s]][i];
                q[++tt] = p;
            }
        }
    }
}
int main()
{
    cin.tie(0)->sync_with_stdio(0);
    cout.tie(0)->sync_with_stdio(0);
    while (cin >> T) {
        if (!T) break;
        init();
        trie();
        build();
        cin >> t;
        for (int i = 0, j = 0; t[i]; i++) {
            int s = t[i] - 'a';
            j = tr[j][s];
            int p = j;
            while (p) {
            //给每次出现的模式串++
                cnt[vis[p]]++;
                p = nex[p];
            }
        }
        int mmax = 0;
        for (int i = 1; i <= T; i++) mmax = max(cnt[i], mmax);
        cout << mmax << endl;
        for (int i = 1; i <= T; i++)
            if (mmax == cnt[i])
                cout << a[i] << endl;
    }
    return 0;
}

题目链接:

【模板】KMP字符串匹配
Period
Radio Transmission
最长前缀 Longest Prefix
Theme Section
动物园
AC 自动机（简单版）
AC 自动机（加强版）

喜欢的对象是坤呢
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
kmp专题
sjs_caomei的博客
08-18 281
kmp即字符串匹配，时间复杂度O(M+N)。 核心为求匹配的字符串的next[]数组（以当前位置结束的最长前缀的位置），代码如下： void getnext() { int i=0,j=-1,len=strlen(str2); a[0]=-1; while(i&lt;len) { if(j==-1||str2[i]==str2[j]) a[++i]=++j; ...
KMP专题
SCUEC__zhazhahui的博客
08-25 325
以下内容适合已经掌握KMP算法的同学观看。 KMP应用一：字符串匹配 这是KMP最常见的应用，我就不详细解释。贴一道模板题， 洛谷3375 AC代码 //字符串匹配模板 #include<stdio.h> #include<cstring> #define maxn 1000010 char s[maxn],t[maxn]; //NEXT是用于字符串匹配的优化版：优化了模式串为aaa...aaa.a这种情况 //NEXT是正常版本，用处更广 int NEXT[maxn],NEXT2
kmp习题大全
AI_xue
02-19 3547
http://acm.hust.edu.cn/vjudge/contest/view.action?cid=70325#overview kuangbin带你飞专题十六kmp&扩展kmp&manacher A. 题意：给T组数据，每组有长度为n和m的母串和模式串。判断模式串是否是母串的子串，如果是输出最先匹配完成的位置，否则输出-1. 做法：直接套用模板。把char改成int。kmp函数中
AcWing 831. KMP字符串 【 KMP算法 】题解
林深不见鹿 的博客
08-05 462
目录1.题目2.提示3.代码 1.题目 给定一个模式串S，以及一个模板串P，所有字符串中只包含大小写英文字母以及阿拉伯数字。 模板串P在模式串S中多次作为子串出现。 求出模板串P在模式串S中所有出现的位置的起始下标。 输入格式 第一行输入整数N，表示字符串P的长度。 第二行输入字符串P。 第三行输入整数M，表示字符串S的长度。 第四行输入字符串S。 输出格式 共一行，输出所有出现位置的起始下标（下标从0开始计数），整数之间用空格隔开。 数据范围 1≤N≤105 1≤M≤106 输入样例： 3 aba 5 a
kmp算法题解
yy254117440的博客
11-11 573
题目是leetcode上的Implement strStr()，题目大意就是在字符串haystack中寻找字符串needle的是否存在，并返回位置。 地址是： https://leetcode.com/problems/implement-strstr/，可以用来训练一下KMP算法的实现。学习KMP算法的过程中，参考了这篇文章： http://blog.csdn.net/yutianzuiji
图解kmp算法
前端学习博客
08-15 1620
但是对于kmp算法，如果遇到了不同字符，其实此次你已经知道了前面比较过的字符是什么，那么就设法利用这个已知信息，不要把"搜索位置"移回已经比较过的位置(注意，这里的比较过的位置指的是文本串中已经进行比较过的位置)，继续把它向后移，这样就提高了效率。可以看到，当比较到 f 时，比较错误，然后找到指针左边(不包括该指针)的最长公共前后缀的长度为2，那么指针就移动到下标为2的地方继续向后比较，这个过程就相当于将前缀的位置移动到后缀上。)，从next数组中可以发现，f 左边的长度为2，所以要移动到下标为2的地方。.
算法/题解-KMP
WjNaG的博客
08-11 402
题目链接： luogu 题目描述： 给出两个字符串s1和s2，其中s2为s1的子串，求出s2在s1中所有出现的位置。然后输出next数组 题解： 这篇题解鸽了好久，因为我一直都没有完全理解看猫片，这东西是真的难。KMP的用处就是在母串里面找子串。 朴素做法是一位一位地匹配，判断子串和母串的字符是否相同，所以时间复杂度是O(mn)。（m、n是母串和字串的长度） KMP就是通过一些神奇的操作跳过一些字...
深入浅出系列之——KMP算法详解【吐血整理】
JA_yichao欢迎你 cnblogs账号：https://www.cnblogs.com/mayichao/
03-13 2965
KMP算法详解——秒杀！
KMP算法例题Oulipo题解
08-17
KMP算法的一个经典题题解，
动态规划之KMP字符匹配算法.md
11-17
动态规划之KMP字符匹配算法.md
AlgorithmProblem:算法题解
04-09
1. KMP算法：快速匹配字符串中的子串，避免了不必要的回溯，提高了效率。 2. Rabin-Karp算法：基于哈希的字符串匹配方法，通过滚动哈希减少比较次数。 五、数学与逻辑 1. 质数判断：快速判断一个整数是否为质数，...
algorithms:各种各种算法题题解
03-06
这个压缩包"algorithms-main"很可能包含了一系列的JavaScript实现的算法题解，通过学习和实践这些题解，开发者可以提升逻辑思维能力，更好地应对实际项目中的挑战。同时，熟悉这些算法也有助于准备像LeetCode、...
华为OD算法题解，并采用模块化代码，形成解题思路
05-01
本压缩包文件包含的是一系列使用Python语言编写的华为OD算法题解，这为我们提供了一个学习和练习算法的好资源。以下是这些题目涉及的一些核心知识点和解题策略： 1. **数据结构**：数据结构是算法的基础，包括数组...
备战秋招60天算法挑战，Day17
溜达虎的博客
08-13 396
给定两个字符串 和 ，编写一个函数来判断 是否是 的字母异位词。注意： 若 和 中每个字符出现的次数都相同，则称 和 互为字母异位词。举个例子： 视频题解 有效的字母异位词 思路来源 思路来源 思路解析 方法一 排序 分别对和进行排序，然后比较排序后的字符串是否相等。 java代码 python代码 复杂度分析 时间复杂度： 快排的时间复杂度是O(nlogn)，其中是字符串的长度。空间复杂度： 快排的空间复杂度是O(logn)。假设， 。 java代码 python代码 复杂度分析 时
算法的学习笔记—删除链表中重复的结点（牛客JZ76）
最新发布
业精于勤，荒于嬉
08-14 1201
本文介绍了一种删除排序链表中重复节点的高效算法。通过递归的方式，算法能够有效地删除所有重复节点，并保留链表中仅出现一次的节点。此方法的实现简单且易于理解，同时在时间和空间复杂度方面表现优异，适用于处理大规模数据链表中的重复节点问题。😁热门专栏推荐想学习vue的可以看看这个java基础合集数据库合集redis合集nginx合集linux合集手写机制微服务组件spring_尘觉springMVCmybits🤔欢迎大家加入我的社区尘觉社区。
【秋招笔试】8.11大疆秋招(第二套)-测开岗
清隆学长的博客
08-12 788
> 🍒 本专栏已收集 70+套笔试题，`笔试真题` 会在第一时间跟新 > > 🍹 近期会进行一波价格调整，没订阅的朋友抓紧冲！！！
【LeetCode每日一题】——623.在二叉树中增加一行
IronmanJay的博客
08-14 1134
给定一个二叉树的根 root 和两个整数 val 和 depth ，在给定的深度 depth 处添加一个值为 val 的节点行。 注意，根节点 root 位于深度 1 。 加法规则如下: 给定整数 depth，对于深度为 depth - 1 的每个非空树节点 cur ，创建两个值为 val 的树节点作为 cur 的左子树根和右子树根。 cur 原来的左子树应该是新的左子树根的左子树。 cur 原来的右子树应该是新的右子树根的右子树。 如果 depth == 1 意味着 depth - 1 根本没有深度
程序算法题解大全：洛谷与CF经典题目详解
"《程序算法题解实用知识库分享》是一份全面的IT资料，专注于算法题目的解析和解答。这份文档涵盖了多种编程竞赛平台（如Codeforces和洛谷）上的经典题目，包括CF1294F、CF1632E1/E2、CF1646D等，涉及到的问题类型...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值