AcWing 状态机DP相关问题 1052. 设计密码 KMP状态机模型,有点难想到

from typing import List


def getKmpNext(data: List) -> List[int]:
    next = [0 for _ in range(len(data))]
    next[0] = -1
    j, k = 0, -1

    while j < len(data) - 1:
        if k == -1 or data[j] == data[k]:
            j, k = j + 1, k + 1
            next[j] = k
        else:
            k = next[k]

    return next


# 预处理先把状态转移的边算出来
N = int(input())
edge = {}
T = input()
next = getKmpNext(T)
for j in range(len(T)):
    for ch in 'abcdefghijklmnopqrstuvwxyz':
        if ch == T[j]:
            edge[(j, ch)] = j + 1
        else:
            jj = j
            while T[jj] != ch:
                jj = next[jj]
                if jj == -1:
                    break
            edge[(j, ch)] = jj + 1

M = len(T)

# dp(i, j)表示的是原字符串和模式串的KMP匹配位置为i和j这样一种匹配状态是多少个可能的j的变化序列构造出来的
# 只要有一个j的变化序列能够构造出(i, j)这种状态,就存在一种长度为i的不包含子串T
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
KMP算法是一种用于字符串匹配的高效算法,它可以在O(m+n)的时间复杂度内解决DNA匹配问题。下面是KMP算法的详细设计: 1. 预处理模式串(pattern): - 初始化两个指针i和j,分别指向模式串的第一个字符和第二个字符。 - 创建一个辅助数组next,用于存储每个字符对应的最长公共前后缀的长度。 - 将next初始化为-1,next初始化为0。 - 进入循环,直到i等于模式串的长度-1: - 如果j等于-1或者当前字符和模式串的第j个字符相等,则将i和j分别增加1。 - 如果当前字符和模式串的第j个字符相等,则将next[i+1]赋值为j。 - 否则,将next[i+1]赋值为next[j]。 - 否则,将j赋值为next[j]。 2. 匹配文本串(text): - 初始化两个指针i和j,分别指向文本串的第一个字符和模式串的第一个字符。 - 进入循环,直到i等于文本串的长度或者j等于模式串的长度: - 如果j等于-1或者当前字符和模式串的第j个字符相等,则将i和j分别增加1。 - 否则,将j赋值为next[j]。 3. 返回匹配结果: - 如果j等于模式串的长度,则表示匹配成功,返回匹配的起始位置为i-j。 - 否则,表示匹配失败,返回-1。 下面是一个使用KMP算法求解DNA匹配问题的示例代码: ```python def kmp_search(text, pattern): n = len(text) m = len(pattern) next = [-1] * m # 预处理模式串 i = 0 j = -1 while i < m - 1: if j == -1 or pattern[i] == pattern[j]: i += 1 j += 1 next[i] = j else: j = next[j] # 匹配文本串 i = 0 j = 0 while i < n and j < m: if j == -1 or text[i] == pattern[j]: i += 1 j += 1 else: j = next[j] # 返回匹配结果 if j == m: return i - j else: return -1 text = "ACGTACGTACGTACGT" pattern = "CGTACG" result = kmp_search(text, pattern) print("匹配的起始位置为:", result) # 输出:3 ```

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值