zoukankan      html  css  js  c++  java
  • UVA

    //法一
    //这题算是初等贪心的思路了,要求总harmming距离最小,通过使每一位上选取的字母,都能使该位对总harming距离的贡献最小,按照定义,也就是所有DNA在该位置上,ATCG字母究竟哪个最多,若多个最多,就选字典序小的...思路确实挺巧妙,值得琢磨
    //顺带提一下,不知道为什么,看了别人的代码,发现有人将变量作为数组维度,即在main函数内部,char DNA[m][n+1];时,居然能AC,不是说数组的维度只能是常量吗?


    #include <bits/stdc++.h>
    using namespace std;
    const int N = 1005;
    const int M = 50;
    char DNA[M][N];
    int num[4];
    
    int main()
    {
    	int t, m, n; //m是数目,n是长度
    	cin >> t;
    	while (t--)
    	{
    		int sum = 0;
    		cin >> m >> n;
    //		char DNA[m][n+1];
    		
    		for (int i = 0; i < m; i++) cin >> DNA[i];
    		for (int i = 0; i < n; i++) 
    		{
    			memset(num, 0, sizeof(num));
    			
    			for (int j = 0; j < m; j++)
    			switch(DNA[j][i])
    			{
    				case 'A' :num[0]++;break;
    				case 'C' :num[1]++;break;
    				case 'G' :num[2]++;break;
    				case 'T' :num[3]++;break;
    			}
    			
    			int temp = num[0], ans = 0;
    			for (int j = 1; j < 4; j++)
    			{
    				if (num[j] > temp)
    				{
    					ans = j; temp = num[j]; 
    				}
    			}
    			sum += m - temp;
    			
    			switch(ans)
    			{
    				case 0: cout << "A";break;
    				case 1: cout << "C";break;
    				case 2: cout << "G";break;
    				case 3: cout << "T";break;
    			}  
    		}
    		cout << endl << sum << endl;
    		
    	} 
    	return 0;
    }
    /*
      法二:
      
      用结构体真的能让代码简单很多很多...
     这个思路有个比较巧妙的地方,在于:将Harming距离所要求的,与对应字符不同的字符个数,转换为与对应字符相同的个数,然后再用总数(DNA个数)减去相同的个数
     
     其实不用这个思路也可以,但是会比较繁琐,例如我假设最终求得的那个DNA的,某个位置是A,那么所有m条DNA,对于每一条而言,如果对应的这个位置不是A,那么这个位置对应的GCT都要加1,繁琐且易错...而如果转变思路,用m减去相同的,需要更新的就只有该位置的A (这种思路下,记录的就是,当前位置下,ATCG分别会和几条DNA的该位置,对总harming距离作出贡献,这就真的是取最小值了)
    */

    #include <bits/stdc++.h>
    using namespace std;
    const int N = 1005;
    char check[N];
    char ans[N];
    struct DNA
    {
    	int a, c, g, t;
    	void clearall()
    	{
    		a = c = g = t = 0;
    	}
    	int getmax(char& ch)//此处很多地方要用偏序的不等号,因为有多个答案时,要输出字典序最小的...一开始忘了=号,结果答案错
    	{
    		ch = ( max(a, c) >= max(g, t) ? (a >= c? 'A':'C'):(g >= t? 'G' : 'T') );
    		return max( max(a, c), max(g, t) );
    	}
    }dna[N];
    
    
    int main()
    {
    	int t, m, n; //m是数目,n是长度
    	cin >> t;
    	while (t--)
    	{
    		int sum = 0;
    		cin >> m >> n;
    		for (int i = 0; i < n; i++)
    			dna[i].clearall();//一定是按照位置清空的,也就是DNA长度,而不是DNA数
    			
    		for (int i = 0; i < m; i++)
    		{
    			cin >> check;
    			for (int j = 0; j < n; j++)
    			switch(check[j])
    			{
    				case'A': dna[j].a++;break;
    				case'G': dna[j].g++;break;
    				case'C': dna[j].c++;break;
    				case'T': dna[j].t++;break;
    			}
    		}
    		
    		for (int i = 0; i < n; i++)
    		{
    			sum += m - dna[i].getmax(ans[i]);
    		}
    		ans[n] = '';
    		
    		cout << ans << endl << sum << endl;
    		
    	} 
    	return 0;
    }





  • 相关阅读:
    离线计算框架 MapReduce
    Hadoop概述
    Linux之rpm/yum
    Linux之磁盘分区
    利用CMD合并多个VOB文件
    android中的simple_list_item
    jquery中的跨域-jsonp格式
    安卓代码中设置ImageView属性
    Android中常用的Adapter的种类和用法
    C#导入excel文件到oracle的方法
  • 原文地址:https://www.cnblogs.com/mofushaohua/p/7789506.html
Copyright © 2011-2022 走看看