DNA序列（DNA Consensus String, ACM/ICPC seoul 2006, UVa 1368)

题目详情请参考原题↓：

UVa 1368

题目中的定义

注意事项

算法思路：

以下串为例：

TATGATAC

TAAGCTAC

AAAGATCC

TGAGATAC

TAAGATGT

注意，在每一列中，取出现次数最多的字符放入答案字符串，例如第一列中 字符T出现4次为最多。以此类推，8列共取T A A G A T A C。

故上例的最优解为TAAGATAC,

在定义常量数组是，将A C G T按字典序排列，就可以很方便的保存并输出字典序最小的字符，即保存字典序最小的答案串，解决了多解问题。

/*
problem sequence :Dnf序列(DNA Consensus String, ACM/ICPC seoul 2006, UVa1368)
*/
# include <cstdio>
# include <cstring>
const int maxn = 1000+10;
const int maxm = 50+5;
const char c[4] = { 'A','C','G','T' };
int x[4] = { 0 };
int m, n;

void get_maxchar(char (*a)[maxn],char (*d)) {//竖向遍历数组，将出现几率最高的字符赋给d
for (int i = 0; i <n ; i++) {
int max = 0;
for (int j = 0; j < m; j++) {
if (a[j][i] == c[0]) {
x[0]++;
}
else if (a[j][i] == c[1]) {
x[1]++;
}
else if (a[j][i] == c[2]) {
x[2]++;
}
else if (a[j][i] == c[3]) {
x[3]++;
}
for (int k = 0; k < 4; k++) {//获取字符出现最高次数
if (max < x[k])
max = x[k];
}
}
for (int k = 0; k < 4; k++) {

if (max == x[k]) {
d[i] = c[k];
break;
}
}
memset(x, 0, sizeof(x));//清空x数组
}
}
int get_ConErr(char(*a)[maxn], char *d ){
//将数组d与数组s的每行比较，如果有一个position不同，则count+1；便利m次；
int cnt = 0;
for (int i = 0; i < m; i++) {
for (int j = 0; j < n; j++) {
//d与s比较
if (a[i][j] != d[j])
cnt++;

}
}
return cnt;
}
int main() {
//freopen("data1.in", "r", stdin);
//freopen("data1.out", "w", stdout);
int T;
scanf("%d", &T);
while (T--) {

scanf("%d%d", &m, &n);
char s[maxm][maxn];
char d[maxn];
memset(s, 0, sizeof(s));
memset(d, 0, sizeof(d));
//get strings
for (int i = 0; i < m; i++) {
scanf("%s", s[i]);
}
get_maxchar(s, d);
printf("%s\n", d);
printf("%d\n", get_ConErr(s,d));
}
return 0;
}