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，那么这个位置对应的G\C\T都要加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] = '\0';cout << ans << endl << sum << endl;} return 0;}