LA3602DNA序列
生活随笔
收集整理的這篇文章主要介紹了
LA3602DNA序列
小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.
題意:
? ? ?給你一個(gè)一些DNA序列(只有ACGT)然后讓你構(gòu)造一個(gè)序列,使得所有的序列到他的Hamming距離最小,所有的序列包括構(gòu)造的序列長度都是N,Hamming表示兩個(gè)序列的不同字符位置個(gè)數(shù),比如ACCT AACA 第二個(gè)和第四個(gè)位置不同,所以距離是2.
思路:
? ? ?簡單題目,我們一位一位構(gòu)建,每一位都是獨(dú)立的,我們選擇出現(xiàn)次數(shù)最多的字母(如果最多的不止一個(gè),那么我們選擇ASCII小的)作為答案,其他的和加到總的距離里面去。
#include<stdio.h>
#include<string.h>
char map[55][1100];
char Ans[1100];
int main ()
{
? ? int t ,n ,m ,i ,j ,sum;
? ? int A ,T ,G ,C;
? ? scanf("%d" ,&t);
? ? while(t--)
? ? {
? ? ? ?scanf("%d %d" ,&n ,&m);
? ? ? ?for(i = 1 ;i <= n ;i ++)
? ? ? ?scanf("%s" ,&map[i]);
? ? ? ?sum = 0;
? ? ? ?for(i = 1 ;i <= m ;i ++)
? ? ? ?{
? ? ? ? ? A = G = C = T = 0;
? ? ? ? ? for(j = 1 ;j <= n ;j ++)
? ? ? ? ? if(map[j][i-1] == 'A') A ++;
? ? ? ? ? else if(map[j][i-1] == 'G') G ++;
? ? ? ? ? else if(map[j][i-1] == 'C') C ++;
? ? ? ? ? else T ++;
? ? ? ? ? ?int max = A;
? ? ? ? ? ?char ch = 'A';
? ? ? ? ? ?if(max < C) max = C ,ch = 'C';
? ? ? ? ? ?if(max < G) max = G ,ch = 'G';
? ? ? ? ? ?if(max < T) max = T ,ch = 'T';
? ? ? ? ? ?Ans[i-1] = ch;
? ? ? ? ? ?sum = sum + A + G + C + T - max;
? ? }
? ? Ans[m] = '\0';
? ? puts(Ans);
? ? printf("%d\n" ,sum);
??
? ?}
? ?return 0;
}
? ? ? ?
《新程序員》:云原生和全面數(shù)字化實(shí)踐50位技術(shù)專家共同創(chuàng)作,文字、視頻、音頻交互閱讀
? ? ?給你一個(gè)一些DNA序列(只有ACGT)然后讓你構(gòu)造一個(gè)序列,使得所有的序列到他的Hamming距離最小,所有的序列包括構(gòu)造的序列長度都是N,Hamming表示兩個(gè)序列的不同字符位置個(gè)數(shù),比如ACCT AACA 第二個(gè)和第四個(gè)位置不同,所以距離是2.
思路:
? ? ?簡單題目,我們一位一位構(gòu)建,每一位都是獨(dú)立的,我們選擇出現(xiàn)次數(shù)最多的字母(如果最多的不止一個(gè),那么我們選擇ASCII小的)作為答案,其他的和加到總的距離里面去。
#include<stdio.h>
#include<string.h>
char map[55][1100];
char Ans[1100];
int main ()
{
? ? int t ,n ,m ,i ,j ,sum;
? ? int A ,T ,G ,C;
? ? scanf("%d" ,&t);
? ? while(t--)
? ? {
? ? ? ?scanf("%d %d" ,&n ,&m);
? ? ? ?for(i = 1 ;i <= n ;i ++)
? ? ? ?scanf("%s" ,&map[i]);
? ? ? ?sum = 0;
? ? ? ?for(i = 1 ;i <= m ;i ++)
? ? ? ?{
? ? ? ? ? A = G = C = T = 0;
? ? ? ? ? for(j = 1 ;j <= n ;j ++)
? ? ? ? ? if(map[j][i-1] == 'A') A ++;
? ? ? ? ? else if(map[j][i-1] == 'G') G ++;
? ? ? ? ? else if(map[j][i-1] == 'C') C ++;
? ? ? ? ? else T ++;
? ? ? ? ? ?int max = A;
? ? ? ? ? ?char ch = 'A';
? ? ? ? ? ?if(max < C) max = C ,ch = 'C';
? ? ? ? ? ?if(max < G) max = G ,ch = 'G';
? ? ? ? ? ?if(max < T) max = T ,ch = 'T';
? ? ? ? ? ?Ans[i-1] = ch;
? ? ? ? ? ?sum = sum + A + G + C + T - max;
? ? }
? ? Ans[m] = '\0';
? ? puts(Ans);
? ? printf("%d\n" ,sum);
??
? ?}
? ?return 0;
}
? ? ? ?
《新程序員》:云原生和全面數(shù)字化實(shí)踐50位技術(shù)專家共同創(chuàng)作,文字、視頻、音頻交互閱讀
總結(jié)
以上是生活随笔為你收集整理的LA3602DNA序列的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: LA3415保守的老师
- 下一篇: LA3635派