算典03_习题_07_UVA-1368

DNA Consensus String

题意

给出一组DNA序列(即一些字符串),找出与每个DNA序列的差最小的DNA序列
差的意思是序列中位置相同但字符不同的位置的个数
输出这个DNA序列以及最小的差

题解

1.准备

这里有一个小技巧，DNA只有“ATCG”四种，要统计这四种出现的次数，就需要一种对应关系，如让ATCG分别对应0123，那么我就可能用a[0]表示A的出现次数

char N[] = {"ACGT"};int Nid(char c){ return strchr(N, c) - N; } //N[Nid(c)];

我这里还是一如既往地使用常量表来很简便地建立这种对应关系，然后用一个查找函数，使其可以通过字符找到下标，如上面的代码，N[Nid(“A”)]就可以表示A的次数，这样就方便多了

2.实现

有了上一步的准备，代码就很显而易见了，这里我对所有的DNA一起统计某一位，这一位上出现最多的字符就是结果在这一位的字符，而差就用总的DNA数减去其出现的次数即可(每一位都这样处理)

#include <iostream>#include <cstdio>#include <cstring>using namespace std;const int maxn = 1e3 +5;#define met(a,b) memset(a, b, sizeof(a));char N[] = {"ACGT"};int Nid(char c){ return strchr(N, c) - N; } //N[Nid(c)];int m, n, res;char s[maxn][maxn];int ans[maxn][4];char solve(int id){    int Max = 0, ID = 0;    for(int i = 0; i < 4; ++i){        if(Max < ans[id][i]) {            Max = ans[id][i];            ID = i;        }    }    res += n - Max;    return N[ID];}int main(){    #ifdef _LOCAL    freopen("in.txt","r", stdin);    #endif // _LOCAL    int t; cin >> t;    while(t--) {        cin >> n >> m;        for(int i = 0; i < n; ++i) {            scanf("%s", s[i]);        }        met(ans, 0);        for(int i = 0; i < m; ++i){            for(int j = 0; j < n; ++j){                ++ans[i][Nid(s[j][i])];            }        }        res = 0;        for(int i = 0; i < m; ++i){            printf("%c", solve(i));        }        printf("\n%d\n", res);    }    return 0;}