UVa1368 DNA序列(DNA Consensus String)

思路：

　　比较同一列A，C，G，T四个字母出现的个数，找到出现数最大的字母即答案，一样的话优先'A'（按字典顺序进行判断就可以得到字典序最小的解）

AC代码

#include<stdio.h>
#include<string.h>

char str[55][1005];
char ans[1005];

int findMax(int a, int b, int c, int d)
{ 
    int tmp[4] = {a, b, c, d};
    int max = 0;
    for(int i = 0; i < 4; i++)
    {
        if(max < tmp[i])
            max = tmp[i];
    }
    return max;
    
}
int main()
{
    int T, m, n, count;
    scanf("%d", &T);
    while(T--)
    {
        count = 0;
        scanf("%d %d", &m, &n);
        getchar();    //吃掉换行符
        for(int i = 0; i < m; i++)
        {
            for(int j = 0; j < n; j++)
            {
                scanf("%c", &str[i][j]);
            }
            getchar();    //吃掉换行符
        }
        //统计每一列字符出现的字数，找到最大的，若一样优先'A'
        for(int i = 0; i < n; i++)
        {
            int a = 0, c = 0, g = 0, t = 0;
            for(int j = 0; j < m; j++)
            {
                if(str[j][i] == 'A')
                    a++;
                if(str[j][i] == 'C')
                    c++;
                if(str[j][i] == 'G')
                    g++;
                if(str[j][i] == 'T')
                    t++;
            }
            int max = findMax(a, c, g, t);
            if(a == max)
            {
                ans[i] = 'A';
                continue;
            }    
            if(c == max)
            {
                ans[i] = 'C';
                continue;
            }    
            if(g == max)
            {
                ans[i] = 'G';
                continue;
            }
            if(t == max)
            {
                ans[i] = 'T';
                continue;
            }    
        }
        //计算总Hamming
        for(int i = 0; i < m; i++)
        {
            for(int j = 0; j < n; j++)
            {
                if(ans[j] != str[i][j])
                    count++;    
            }    
        } 
        for(int i = 0; i < n; i++)
            printf("%c", ans[i]);
        printf("
");
        printf("%d
", count);
    }    
}

PS：

　　其实没有必要开ans[]数组，每一列确定好字母时直接输出即可。