基因相关性（C++）简单实用

参考代码:

#include <iostream>

#include <string>

using namespace std;

string compareDNA(const string& dna1, const string& dna2, double threshold) {

    int totalPairs = 0;

    int samePairs = 0;

    for (int i = 0; i < dna1.length(); i++) {

        if (dna1[i] != '-' && dna2[i] != '-') {

            totalPairs++;

            if (dna1[i] == dna2[i]) {

                samePairs++;

            }

        }

    }

    double similarityRatio = static_cast<double>(samePairs) / totalPairs;

    if (similarityRatio >= threshold) {

        return "yes";

    } else {

        return "no";

    }

}

int main() {

    double threshold;

    cin >> threshold;

    string dna1, dna2;

    cin >> dna1 >> dna2;

    string result = compareDNA(dna1, dna2, threshold);

    cout << result << endl;

    return 0;

}