16.题目:编号3904 DNA序列修正

文章讨论了如何通过使用字符数组、map数据结构以及不同策略(如互换和替换)来计算两条DNA链进行最小操作次数,以便使它们匹配。作者展示了多种方法,包括基于数组的计数和map的直接比较,以找到最优解。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目:

 ###本题主要考察模拟

#include<bits/stdc++.h>
using namespace std;
const int S=1010;
#define N 4
char a[S],b[S];
int num[N][N];
int main(){
  ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
  memset(num,0,sizeof(num));
  int n;cin>>n;
  for(int i=0;i<n;++i)cin>>a[i];
  for(int i=0;i<n;++i)cin>>b[i];
  for(int i=0;i<n;++i){
    if(a[i]=='A'){
      if(b[i]=='A')num[0][3]++;
      if(b[i]=='C')num[0][2]++;
      if(b[i]=='G')num[0][1]++;
    }
    if(a[i]=='C'){
      if(b[i]=='A')num[1][3]++;
      if(b[i]=='C')num[1][2]++;
      if(b[i]=='T')num[1][0]++;
    }
    if(a[i]=='G'){
      if(b[i]=='A')num[2][3]++;
      if(b[i]=='G')num[2][1]++;
      if(b[i]=='T')num[2][0]++;
    }
    if(a[i]=='T'){
      if(b[i]=='C')num[3][2]++;
      if(b[i]=='G')num[3][1]++;
      if(b[i]=='T')num[3][0]++;
    }
  }
  int ans=0;
  for(int i=0;i<N;++i){
    for(int j=0;j<i;++j){
      ans+=max(num[i][j],num[j][i]);
    }
  }
  cout<<ans<<'\n';
  return 0;
}

整体思路:

  1. 我采用了开两个字符类型数组和一个int类型的数组来分别存储DNA的两条链的字母和每种情况,第一条链入a数组,第二条链入b数组
  2. 要最小操作次数那么互换操作一定是优于替换操作的,根据规律,这样关于主对角线的两组是可以进行互换的,主对角线上的是已经配好的不需要操作的,我们将每个位置上的个数用num数组都做预处理,将之记录下来
  • 我们接着想这样可以互换的两个位置,将所以情况,两个位置的数取最小值就是我们的互换的操作次数(互换)
  • 那么,替换次数一定就是剩下没有配对的多出来的次数,用所有最小值再加上多出来的就能得到答案了(替换)

我们换个角度来想,其实我们只需要对主对角两侧能匹配的取两个位置的最大值就是替换+互换的次数了 ,如此我们借鉴9 9乘法表的循环用max函数对角标互换的两个位置取最大值都相加就得到最小操作次数了


再来看下这种做法:(用map计数,采用了计数和等于3的方法巧妙的将不相干的两条链有了联系中间的判断与上面写法的最后的循环有异曲同工之处,通过算数层面的判断,实现非交换即替换的操作,对交换的数据实时做出修改

#include<iostream>
#include<map>

using namespace std;

map<char, int> mp{
        {'A', 0},
        {'C', 1},
        {'G', 2},
        {'T', 3}
};

int main() {
    int n;
    cin >> n;
    string a, b;
    cin >> a >> b;
    int cnt = 0;
    for (int i = 0; i < n; i++) {
        if (mp[a[i]] + mp[b[i]] != 3) {
            for (int j = i + 1; j < n; j++) {
                if (mp[a[i]] + mp[b[j]] == 3 && mp[a[j]] + mp[b[i]] == 3) {
                    swap(b[i], b[j]);
                    break;
                }
            }
            cnt++;
        }
    }
    cout << cnt << endl;

    return 0;
}


再看这份:(用map的两个char类型,以函数的形式,在字母层面进行了判断,也是做出了非交换即替换的判断)

#include <bits/stdc++.h>
using namespace std;
int find(string &s1,string &s2,int i,char a,char b) {
    while (i < s1.size()) {
        if (s2[i] == a && s1[i] == b) return i;
        i++;
    }
    return -1;
}
int main() {
    ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
    int n,cnt = 0; cin >> n;
    string s1,s2; cin >> s1 >> s2;
    map<char,char>mp = {{'A','T'},{'C','G'},{'T','A'},{'G','C'}};
    for (int i = 0; i < n; i++) {
        s1[i] = mp[s1[i]];
    }
    for (int i = 0; i < n; i++) {
        if (s1[i] != s2[i]) {
            int x = find(s1,s2,i,s1[i],s2[i]);
            if (x != -1) swap(s2[i],s2[x]);
            else s2[i] = s1[i];
            cnt++; 
        }
    }
    cout << cnt << '\n';
    
    return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值