习题3-7 UVA1368 DNA Consensus String（31行AC代码）

紫书刷题进行中，题解系列点这里

习题3-7 UVA1368 DNA Consensus String（31行AC代码）

题目大意

给定若干个DNA序列，仅由ACGT四个字符构成，求解Hamming距离最小的字符串并输出错误数。

Hamming距离：字符不同的位置个数总和

原串1：ATCG
原串2：CTAG
共识串：ATAG
错误数：2 （第一个位置和第三个位置）

若有多个合法字符，则选择字典序最小的，再如上例，第一和第三个位置中AC个数相同，均选择字典序小者A

思路分析

按行读入字符串s，遍历s，按列统计ACTG出现次数，遍历完毕后，按列找出最大值MAX，输出相应字符，m-MAX作为总错误数的贡献，需累加。

可用哈希表完成0-3到ACGT之间双射，这里按字典序排列，保证存在多解时找到的都是字典序最小串

AC代码（C++11，简单模拟，打表）

#include<bits/stdc++.h>
using namespace std;
map<char, int> mp{{'A',0}, {'C',1}, {'G',2}, {'T',3}}; // 便于循环处理
char itc[4] = {'A', 'C', 'G', 'T'}; // [0,4)->[A,C,G,T]
int T, m, n, num[1010][4] = {0}; // ATGC出现次数
string s;
int main() {
    scanf("%d", &T);
    while(T --) {
        scanf("%d %d", &m, &n);
        fill(num[0], num[0] + 4040, 0); // 初始化二维数组
        for (int i = 0; i < m; i ++) { // 统计每一列中每个子符出现次数
            cin >>s;
            for (int j = 0; j < n; j ++) num[j][mp[s[j]]] ++;
        }
        int errNum = 0, idx = 0, MAX = -1;
        for (int i = 0; i < n; i ++) { // 遍历每一列
            idx = 0, MAX = -1; // 初始化
            for (int j = 0; j < 4; j ++) { // ACGT
                if (MAX < num[i][j]) { // 找最大值
                    idx = j;
                    MAX = num[i][j];
                }
            }
            errNum += (m - MAX); // 不同的个数
            printf("%c", itc[idx]);
        }
        printf("\n%d\n", errNum);
    }
    return 0;
}