c++查找DNA中子序列

DNA子序列的选择和分析是基因组学、分子生物学和生物信息学等领域研究的重要组成部分，对理解基因组结构和功能、进化过程以及疾病机理具有重要意义。以下是简单的函数示例，用于返回给定子序列在DNA中的所有位置。

#include <iostream>
#include <string>
#include <vector>

std::vector<size_t> findSubstringPositions(const std::string& str, const std::string& sub) {
    std::vector<size_t> positions;
    size_t pos = str.find(sub, 0);
    while (pos != std::string::npos) {
        positions.push_back(pos);
        pos = str.find(sub, pos + 1);
    }
    return positions;
}

int main() {
    std::string str = "TGACGCTGACGTCGCTAG";
    std::string sub = "GCT";

    std::vector<size_t> positions = findSubstringPositions(str, sub);

    if (positions.empty()) {
        std::cout << "子序列未找到" << std::endl;
    } else {
        std::cout << "子序列在字符串中的位置：" << std::endl;
        for (size_t pos : positions) {
            std::cout << pos << std::endl;
        }
    }

    return 0;
}