DNA序列重复模式检测代码故障求助：仅匹配起始位置模式

核心问题 #

你的代码中pattern_count函数的return count语句缩进错误，被放在了for循环内部，导致循环仅执行一次（仅检查序列起始位置）就直接返回结果，因此无法检测后续位置的模式匹配。

修正后的代码 #

将return count移到for循环外部，确保遍历完所有可能的匹配位置后再返回计数。同时优化代码结构，把函数定义移到循环外避免重复定义，并新增匹配位置的输出功能：

def pattern_count(text, pattern):
    count = 0
    positions = []  # 存储匹配的起始位置（索引从0开始）
    pattern_len = len(pattern)
    for i in range(len(text) - pattern_len + 1):
        if text[i: i + pattern_len] == pattern:
            count += 1
            positions.append(i)
    return count, positions

# 目标核苷酸序列
text = 'AGACGCCTGGGAACTGCGGCCGCGGGCTCGCGCTCCTCGCCAGGCCCTGCCGCCGGGCTGCCATCCTTGCCCTGCCATGTCTCGCCGGAAGCCTGCGTCGGGCGGCCTCGCTGCCTCCAGCTCAGCCCCTGCGAGGCAAGCGGTTTTGAGCCGATTCTTCCAGTCTACGGGAAGCCTGAAATCCACCTCCTCCTCCACAGGTGCAGCCGACCAGGTGGACCCTGGCGCTgcagcggctgcagcggccgcagcggccgcagcgCCCCCAGCGCCCCCAGCTCCCGCCTTCCCGCCCCAGCTGCCGCCGCACATA'

while True:
    print("Input Pattern:")
    pattern = input().strip()
    if not pattern:
        print("请输入有效的核苷酸模式")
        continue
    match_count, match_positions = pattern_count(text, pattern)
    print(f"模式重复次数: {match_count}")
    if match_count > 0:
        print(f"匹配起始位置: {match_positions}")

额外优化建议 #

• 大小写兼容处理：原序列包含大小写混合的核苷酸，若需忽略大小写匹配，可在比较前统一转换为同一大小写，例如：text[i: i + pattern_len].upper() == pattern.upper()
• 边界校验：增加了空输入的判断，避免无效计算
• 效率优化：提前计算模式长度，避免循环内重复调用len(pattern)

内容的提问来源于stack exchange，提问作者ClarkThark