中南大学上机：PIPI的DNA序列

问题 B: PIPI的DNA序列

时间限制: 1 Sec 内存限制: 128 MB

题目描述

PIPI有m个长度均为n的DNA序列，求一个DNA序列，使其到所有序列的总HaMM距离尽量小。
HaMM距离: 两个等长字符串的HaMM距离等于字符不同的位置个数。

输入

输入包含多组测试用例。
对于每一组测试用例，第一行输入整数m和n (4<=m<=50,4<=n<=1000)。
接下来m行每行一个长度为n的DNA序列(只包含A、C、G、T)。

输出

对于每一组测试用例，输出让HaMM距离最小的DNA序列和其对应的距离。若有多解，输出字典序最小的解。

样例输入

3
5 8
TATGATAC
TAAGCTAC
AAAGATCC
TGAGATAC
TAAGATGT
4 10
ACGTACGTAC
CCGTACGTAG
GCGTACGTAT
TCGTACGTAA
6 10
ATGTTACCAT
AAGTTACGAT
AACAAAGCAA
AAGTTACCTT
AAGTTACCAA
TACTTACCAA

样例输出

TAAGATAC
7
ACGTACGTAA
6
AAGTTACCAA
12

代码：

#include<stdio.h>
int main(){int a[6];//a1 A a2 C a3 G a4 Tint times;scanf("%d",&times);while(times--){int n,m;scanf("%d%d",&n,&m);char str[51][1002];int i,j; int max = 0;int k;int count = 0;char xin[m];for(i=1;i<=n;i++){scanf("%s",str[i]);}for(i=0;i<m;i++){a[1]=0; a[2]=0;a[3]=0;a[4]=0; max=0;for(j=1;j<=n;j++){if(str[j][i]=='A')a[1]++;else if(str[j][i]=='C')a[2]++;else if(str[j][i]=='G')a[3]++;elsea[4]++;}for(j=1;j<=4;j++){if(a[j]>max){max=a[j];k=j;}}count = count + max;if(k==1)xin[i]='A';else if(k==2)xin[i]='C';else if(k==3)xin[i]='G';elsexin[i]='T';}for(i=0;i<m;i++){printf("%c",xin[i]);}printf("\n");printf("%d\n",n*m-count);}}

中南大学上机：PIPI的DNA序列相关推荐

求最小Hamming距离的DNA序列
提示:两个等长字符串的Hamming距离,等于两个字符串相同的位置中,字符不同的个数.(如ACGT和GCGA的Hamming距离为2) 问题: TATGATAC TAAGCTAC AAAGATCC T ...
UA MATH567 高维统计III 随机矩阵12 整数环上的区间的应用：DNA序列突变点侦测的统计量及假设检验
UA MATH567 高维统计III 随机矩阵12 整数环上的区间的应用:拐点侦测的统计量及假设检验整数环上的区间作为随机变量的下标最大值的概率不等式应用:DNA序列突变点侦测整数环上的区间作 ...
21行代码AC——习题3-7 DNA序列(UVa-1368)_解题报告
励志用尽量少的代码做高效表达. 题目(提交)链接→UVa-1368 思路: DNA序列:按列遍历,记录每一列出现次数最多(若同样多,则字典序最小)的字母,录入s串累加. 距离:重新遍历,录入出现次数比 ...
Leetcode 187.重复的DNA序列
重复的DNA序列所有 DNA 由一系列缩写为 A,C,G 和 T 的核苷酸组成,例如:"ACGAATTCCG".在研究 DNA 时,识别 DNA 中的重复序列有时会对研究非常有帮 ...
A1051. DNA序列
问题描述人类基因组计划的第一阶段于2000年6月26日胜利结束,我国的科研工作者圆满地完成了其中的1%的测序工作. 众所周知,对于任意两个人来说,他们的染色体上的DNA序列大部分是相同的(否则就不是 ...
3-7 DNA序列（DNA Consensus String, ACM/ICPC Seoul 2006, UVa1368）
输入m个长度均为n的DNA序列,求一个DNA序列,到所有序列的总Hamming距离尽量小. 两个等长字符串的Hamming距离等于字符不同的位置个数,例如,ACGT和GCGA的 Hamming距离为 ...
python文件处理，将DNA序列转换为RNA序列
1 #!/usr/bin/python 2 #-*- coding:utf-8 -*- 3 "将DNA序列转换为RNA序列,即将T转换为U即可,利用字符串的replace方法" 5 ...
DNA序列编码中Hairpin的定义和计算
DNA序列编码中Hairpin的定义和计算觉得有用的话,欢迎一起讨论相互学习~ 参考文献 [1] 张凯. DNA计算核酸编码优化及算法设计[D]. 2008. [2] Shin, Soo Yong ...
基因家族进化分析之DNA序列批量获取
使用方法参考: 基于全基因组的基因家族分析(3):SlNRAMP家族基因CDS和Genomic DNA序列获取 Gene序列要获得的ID号不仅仅是seqid,而且还需要在染色体上的位置信息--起始和终 ...

中南大学上机：PIPI的DNA序列

问题 B: PIPI的DNA序列

中南大学上机：PIPI的DNA序列相关推荐

最新文章

热门文章