习题 3-7 DNA序列

版权声明:本人菜鸟一只,如文章有错误或您有高见,请不吝赐教 https://blog.csdn.net/qq_41138935/article/details/82629063
#include<iostream>
#include<cstring>
using namespace std;
int main(){
	char a[55][1005];
	int i,j,k,sum=0;
	int m,n;
	cin>>m>>n;
	
	string str;
	for(i=0;i<m;i++){
		cin>>str;
		for(j=0;j<n;j++){
			a[i][j]=str[j];
		}
	}
	cout<<"\n\n\n";
	for(j=0;j<n;j++){
		int t,aa,g,c;
		aa=c=g=t=0;
		for(i=0;i<m;i++){
			char b=a[i][j];
			switch(b){
				case 'A':aa++;break;
				case 'C':c++;break;
				case 'G':g++;break;
				case 'T':t++;break;
			}
		}
		sum+=aa+c+g+t;
		if(aa>=c&&aa>=g&&aa>=t){
			cout<<"A";
			sum-=aa;
		}else if(c>aa&&c>=g&&c>=t){
			cout<<"C";
			sum-=c;
		}else if(g>aa&&g>c&&g>=t){
			cout<<"G";
			sum-=g;
		}else if(t>aa&&t>c&&t>g){
			cout<<"T";
			sum-=t;
		}
	}
	cout<<"\n\n"<<sum;	
	return 0;
} 

猜你喜欢

转载自blog.csdn.net/qq_41138935/article/details/82629063