求最小Hamming距离的DNA序列

提示：两个等长字符串的Hamming距离，等于两个字符串相同的位置中，字符不同的个数。（如ACGT和GCGA的Hamming距离为2）

问题：

TATGATAC

TAAGCTAC

AAAGATCC

TGAGATAC

TAAGATGT

求一个与上面所有DNA序列 Hamming距离最小的DNA序列

解：TAAGATAC

思路：找到每一列中出现最多的字符，并把一个字符记录下来，当每一列出现最多的字符都记录下来后，这个记录就是答案。

方法1（小白）：

#include<iostream>
using namespace std;
#define max 256
int n, m;
char DNA[10][10];char maxchar(int index) { //找到 列 中出现最多的字符int i,tmp=0, maxs;int W[max] = { 0 };char t;for (i = 0; i < m; i++) { // 统计列中各个字符的数量t = DNA[i][index];  W[t]++;}if (tmp < W['T']) {tmp = W['T'];maxs = 'T';}if (tmp < W['C']) {tmp = W['C'];maxs = 'C';}if (tmp < W['A']) {tmp = W['A'];maxs = 'A';}if (tmp < W['G']) {tmp = W['G'];maxs = 'G';}return maxs; //返回出现最多次的字符
}int main() {cin >> n >> m;int i, j;char str[10];for (i = 0; i < m; i++)for (j = 0; j < n; j++)cin >> DNA[i][j];for (i = 0; i < n; i++) {cout << maxchar(i); }return 0;
}

方法2（大神版）：

#include<iostream>
#include<string>
using namespace std;
#include<algorithm>
#include<assert.h>
#include<vector>struct ChCnt {int cnt;//与c相同的字符的个数char c;void init(char ch = 'A') {c = ch;cnt = 0;}ChCnt() {init();}bool operator<(const ChCnt& cc2) const {  // sort排序用到的函数return cnt > cc2.cnt || (cnt == cc2.cnt && c < cc2.c);}
};int main() {int T = 1;int n, m;cin >> T;string line;vector<string> seqs; char IDX[256] = { 0 }; // 定好 字符 对应的 数组索引，方便找到 字符 对应的 结构体IDX['A'] = 0; IDX['C'] = 1; IDX['G'] = 2; IDX['T'] = 3;while (T--) {seqs.clear();cin >> m >> n;for (int i = 0; i < m; i++) {cin >> line;assert(line.size() == n); // 判断属否输入 n 个字符seqs.push_back(line);// 加入到 seqs 容器}string ansStr;int ans = 0;vector<ChCnt> ccs(4); // 定义4个ChCnt结构体，分别存放 A C G Tfor (int i = 0; i < n; i++) {// 每一列的开始都要对 结构体数组 进行初始化ccs[0].init('A');ccs[1].init('C');ccs[2].init('G');ccs[3].init('T');for (int j = 0; j < m; j++)ccs[IDX[seqs[j][i]]].cnt++;  // IDX在这里的作用是，能在结构体数组中 找到 字符 对应的 结构体sort(ccs.begin(), ccs.end());ansStr += ccs.front().c;ans += (m - ccs.front().cnt); // 加上每列的 Hamming距离 }cout << ansStr << endl << ans << endl;}return 0;}

求最小Hamming距离的DNA序列相关推荐

21行代码AC——习题3-7 DNA序列(UVa-1368)_解题报告
励志用尽量少的代码做高效表达. 题目(提交)链接→UVa-1368 思路: DNA序列:按列遍历,记录每一列出现次数最多(若同样多,则字典序最小)的字母,录入s串累加. 距离:重新遍历,录入出现次数比 ...
3-7 DNA序列（DNA Consensus String, ACM/ICPC Seoul 2006, UVa1368）
输入m个长度均为n的DNA序列,求一个DNA序列,到所有序列的总Hamming距离尽量小. 两个等长字符串的Hamming距离等于字符不同的位置个数,例如,ACGT和GCGA的 Hamming距离为 ...
DNA序列 UVa1368
原题链接:https://cn.vjudge.net/problem/UVA-1368 大概题意: 输入m个长度均为n的DNA序列,求一个DNA序列,到所有序列的总Hamming距离尽量小.两个等长字 ...
中南大学上机：PIPI的DNA序列
问题 B: PIPI的DNA序列时间限制: 1 Sec 内存限制: 128 MB 题目描述 PIPI有m个长度均为n的DNA序列,求一个DNA序列,使其到所有序列的总HaMM距离尽量小. HaM ...
信息论与编码_信息与编码系列（五）不可靠信道——决策法则、Hamming距离与信息论基本定理...
目录序决策法则(decision rules) Hamming距离信息论基本定理--Shannon 定理 Shannon定理的逆命题序这个算是正本书中关于信息论的最后内容了,后面就是关于编码 ...
dna序列分类数学建模matlab,数学建模DNA序列分类模型(终稿).doc
您所在位置:网站首页 > 海量文档 &nbsp>&nbsp高等教育&nbsp>&nbsp生物学数学建模DNA序列分类模型(终稿).doc32页本文 ...
基于主成分分析与Fisher判别的DNA序列分类
基于主成分分析与Fisher判别的DNA序列分类 [摘要]: 特征提取:对于DNA序列,首先将其分为编码DNA特征与非编码DNA特征.对于非编码DNA,统计A,T,C,G各个碱基以及A+T在其中出现的 ...
Python访问街区10个点，并俩俩绘制一条线，得到5条线，求最短的距离和？
Python访问街区10个点,并俩俩绘制一条线,得到5条线,求最短的距离和? 1. 效果图 2. 源码参考上一篇博客介绍了Python访问街区所有节点最短路径问题,并结合matplotlib可视化 ...
Bridge Across Islands POJ - 3608 旋转卡壳求凸包最近距离
\(\color{#0066ff}{题目描述}\) 几千年前,有一个小王国位于太平洋的中部.王国的领土由两个分离的岛屿组成.由于洋流的冲击,两个岛屿的形状都变成了凸多边形.王国的国王想建立一座桥来连接 ...

求最小Hamming距离的DNA序列

求最小Hamming距离的DNA序列相关推荐

最新文章

热门文章