UVa 1368 DNA Consensus String

题

1368 DNA Consensus String

要弄明白的

• Hamming distance

  就是不相同的个数。

解

• 一列一列的看，选重复最多的那个。

• 统计字符出现次数？

  计数器

• 把多个变量存入一个字符串

  直接创建一个大空间的字符串，for循环赋值。

• 字典最小

  问题的本质是，一列中如果最多出现次数的字母 有多个，那么取最小的那个就能保证字典最小。在判断的时候，将ACGT从大到小排。

总结

• buf未清空的问题

在程序中我们这样做：

#define maxn 1001
int main()
{
    char buf[max]
}

创建一个比题目给定空间大的maxn，保证空间够用。然后进行如下操作：

for()
{
    int buf[maxn]; //
    for()
    {
        //对buf赋值
    }
}

外层for循环进行第二次的时间，buf第一次被赋的值 不会被清除。所以容易产生问题。使用的时候一定要保证正确的buf大小。

例如第一次给buf赋值为1,2,3,4,5

buf = {1, 2, 3, 4, 5,...}

第二次赋值为1，2，3，4，那么此时的buf为：

buf = {1, 2, 3, 4, 5, ....}

5依然存在。如果第二次用错了buf的大小，那么容易产生错误的结果。所以这个预先声明大空间的方法使用简单，但也空间出错。我在做题的过程中写的getMore函数，使用了strlen来得到它的长，就会出现上面的问题，这里应该直接指定s的大小，不然如果出来多余的字符就会有问题。

mycode

#include <stdio.h>
#include <string.h>
#define maxm 60
#define maxn 1001

//获取出现最多字符的位置和出现个数
int getMore(char *s, int slen, char &resCh, int &num)
{
    int location[4]; //A,C,G,T出现个数
    char res[4] = {'A', 'C', 'G', 'T'}; //方便返回字符，从小到大排保证字典最小
    memset(location, 0, sizeof(location));
    //int len = strlen(s);
    for (int i = 0; i < slen; i++)
    {
        if (s[i] == 'A')
        {
            location[0]++;
        }
        if (s[i] == 'C')
        {
            location[1]++;
        }
        if (s[i] == 'G')
        {
            location[2]++;
        }
        if (s[i] == 'T')
        {
            location[3]++;
        }
    }
    int locMore = 0;
    for (int i = 0; i < 3; i++)
    {
        if (location[locMore] < location[i + 1])
        {
            locMore = i + 1;
        }
    }
    resCh = res[locMore];
    num = location[locMore];
    return locMore;
}
int main()
{
    int T;
    scanf("%d", &T);
    for (int i = 0; i < T; i++)
    {
        int m, n;
        scanf("%d%d", &m, &n);
        char s[maxm][maxn]; //输入的字符串
        char res[maxn]; //返回结果字符串
        for (int im = 0; im < m; im++)
        {
            scanf("%s", s[im]);
        }

        int hamming = 0;
        for (int in = 0; in < n; in++)
        {
            char buf[maxm];
            for (int im = 0; im < m; im++)
            {
                //把多个变量存入一个字符串
                buf[im] = s[im][in];
            }
            char locMore;
            int num=0;

            getMore(buf, m, locMore, num); //获取出现次数最多的字符并得到出现次数
            hamming += m - num;
            res[in] = locMore;

        }
        for (int j = 0; j < n; j++)
        {
            putchar(res[j]);
        }
        printf("\n%d\n", hamming);
    }
    return 0;
}