【UVA】 1368 --- DNA Consensus String

题意：
输入m个长度均为n的DNA序列， 求一个DNA序列， 到所有序列的总Hamming距离尽量
小。 两个等长字符串的Hamming距离等于字符不同的位置个数， 例如， ACGT和GCGA的
Hamming距离为2（ 左数第1, 4个字符不同） 。
输入整数m和n（ 4≤m≤50, 4≤n≤1000） ， 以及m个长度为n的DNA序列（ 只包含字母
A， C， G， T） ， 输出到m个序列的Hamming距离和最小的DNA序列和对应的距离。 如有多
解， 要求为字典序最小的解。 例如， 对于下面5个DNA序列， 最优解为TAAGATAC。
输入
TATGATAC
TAAGCTAC
AAAGATCC
TGAGATAC
TAAGATGT
输出
TAAGATAC
7

思路：
依次找出每个序列中第i个字符中出现次数最多的那个字符，将其放在dest数组中第i个位置。
将i从0遍历到m，得出目标序列。
然后依次比较每个序列，将每个序列对应的Hamming相加；

注意：程序中字符顺序最好为ACGT，这样在有多解的情况下不用判断字典序最小的序列是哪个。

AC代码：

#include 
#include 
#include 
int Input(int m, int n); //输入二维数组
int Output(int m, int n);//输出目标序列
using namespace std;
char arr[51][1001], dest[1001]; //arr存储输入的序列，dest存储目标序列
int main()
{
	int t;
	cin >> t;
	while (t--)
	{
		int m, n, num = 0;
		cin >> m >> n;
		getchar();
		memset(dest, '\0', sizeof(dest));
		Input(m, n); //输入二维数组
		Output(m, n);//输出目标序列
		for (int i = 0; i < m; i++)
		{
			for (int j = 0; j < n; j++)
			{
				if (arr[i][j] != dest[j])
				{
					num++;
				}
			}
		}
		cout << dest << endl;
		cout << num << endl;
	}
	return 0;
}

int Input(int m, int n) //输入二维数组
{
	for (int i = 0; i < m; i++)
	{
		for (int j = 0; j < n; j++)
		{
			scanf("%c", &arr[i][j]);
		}
		getchar();
	}
	return 0;
}

int Output(int m,int n) //输出目标序列
{
	for (int j = 0; j < n; j++)
	{
		int count[5]{ 0 }; //统计出现次数
		for (int i = 0; i < m; i++)
		{
			switch (arr[i][j])
			{
			case 'A':
				count[0]++;
				break;
			case 'C':
				count[1]++;
				break;
			case 'G':
				count[2]++;
				break;
			case 'T':
				count[3]++;
				break;
			}
		}
		count[4] = count[0]; //用count[4]标记出现次数最多的次数
		int max = 0;
		for (int i = 1; i < 4; i++)//通过遍历找到每列AGCT中出现次数最多的字符
		{
			if (count[4] < count[i])
			{

				max = i;
				count[4] = count[i];
			}
		}
		switch (max) //将每列中出现最多的字符放入dest数组相应位置
		{
		case 0:
			dest[j] = 'A';
			break;
		case 1:
			dest[j] = 'C';
			break;
		case 2:
			dest[j] = 'G';
			break;
		case 3:
			dest[j] = 'T';
			break;
		}
	}
	return 0;
}