问题描述
在生物学中,DNA 序列的相似性常被用来研究物种间的亲缘关系。现在我们有两条 DNA 序列,每条序列由 A、C、G、T 四种字符组成,长度相同。但是现在我们记录的 DNA 序列存在错误,为了严格满足 DNA 序列的碱基互补配对即 A - T 和 C - G,我们需要依据第一条 DNA 序列对第二条 DNA 序列进行以下操作:
-
选择第二条 DNA 序列的任意两个位置,交换他们的字符。
-
选择第二条 DNA 序列任意一个位置,将其字符替换为 A、C、G、T 中的任何一个。
需要注意的是:每个位置上的碱基只能被操作一次!
你的任务是通过最小的操作次数,使第二条 DNA 序列和第一条 DNA 序列互补。并且已知初始两条 DNA 序列长度均为 N。
输入格式
第一行包含一个整数 N(0<N≤1000)表示 DNA 序列的长度。
接下来的两行,每行包含一个长度为 N 的字符串,表示两条 DNA 序列。
输出格式
输出一个整数,表示让第二条 DNA 序列和第一条 DNA 序列互补所需的最小操作次数。
#include "iostream"
#include "vector"
#include "string"
using namespace std;
const int N = 1e3 + 10;
string str1, str2;
bool isMatch(char ch1, char ch2) {
return (ch1 == 'A' && ch2 == 'T') || (ch1 == 'C' && ch2 == 'G') ||
(ch1 == 'T' && ch2 == 'A') || (ch1 == 'G' && ch2 == 'C');
}
int findMinSwap() {
vector<int> unMatch;
int ans = 0;
// 记录未匹配的位置
for (int i = 0; i < str1.length(); i++) {
if (!isMatch(str1[i], str2[i]))
unMatch.push_back(i);
}
/*
* 字符串1的i位置与字符串2的j位置匹配
* 字符串1的j位置与字符串2的j位置匹配
* 进行交换
* */
for (int i = 0; i < unMatch.size(); i++) {
for (int j = i + 1; j < unMatch.size(); j++) {
if (isMatch(str1[unMatch[i]], str2[unMatch[j]]) &&
isMatch(str1[unMatch[j]], str2[unMatch[i]])) {
swap(str2[unMatch[i]], str2[unMatch[j]]);
ans++;
unMatch.erase(unMatch.begin() + j); // 注意这里一定要先释放j位置 如果先释放i位置 那么j位置会向前移动一个位置 导致错误!!!
unMatch.erase(unMatch.begin() + i);
if(i)i--;
break;
}
}
}
// ans 加上没有交换的位置只能改变碱基
ans += unMatch.size();
return ans;
}
int main() {
int n;
cin >> n;
cin >> str1 >> str2;
cout << findMinSwap() << endl;
}