DNA序列修正

这篇文章详细描述了一个C++程序,用于统计给定字符串中特定字符对(如A-T,C-G等)的出现次数,并计算最大计数值,最后输出结果。
摘要由CSDN通过智能技术生成

#include<bits/stdc++.h>
using namespace std;
const int N=1e3+10;
int a2,c2,g2,t2;
int l1,l2,l3,l4,l11,l22,l33,l44;
int main(){
	int n,j=-1;
	cin>>n;
	int ans=0,sum=0;
	string a,b;
	cin>>a>>b;
	for(auto &i:a){
		j++;
		if(i=='A'&&b[j]=='T'||i=='T'&&b[j]=='A'||i=='C'&&b[j]=='G'||i=='G'&&b[j]=='C')continue;
		if(i==b[j]){
			if(i=='A')a2++;
			if(i=='C')c2++;
			if(i=='G')g2++;
			if(i=='T')t2++;
		}
		if(i=='A'&&b[j]=='G')l1++;
		if(i=='C'&&b[j]=='A')l2++;
		if(i=='G'&&b[j]=='A')l3++;
		if(i=='A'&&b[j]=='C')l4++;
		
		if(i=='C'&&b[j]=='T')l11++;
		if(i=='T'&&b[j]=='G')l22++;
		if(i=='T'&&b[j]=='C')l33++;
		if(i=='G'&&b[j]=='T')l44++;
		
	}
	ans+=max(a2,t2)+max(c2,g2)+max(l1,l11)+max(l2,l22)+max(l3,l33)+max(l4,l44);
	cout<<ans;
	return 0;
}

有疑问可以在下面留言哦

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值