uva1368

题目大意:
给m个长度为n的DNA序列,求一个DNA序列到所有序列的总汉明距离最小,汉明距离为不同字符的位置的个数。

思路:
用贪心的方法,刚开始还以为是在给定的DNA中选出一个,可是并不是,是要求一个DNA序列可以使得它到给定的DNA序列的距离最小。每次都求出一列中最多相同的字母。

代码:

#include <iostream>
using namespace std;
#include <cstring>
#include <stdio.h>

char grid[1005];
struct node {
    int a,c,g,t;
}N[1005];
int main() {
    int n,m;
    int T;
    scanf("%d",&T);
    while(T--) {
        scanf("%d %d",&n,&m);
        for(int i = 0 ; i < m; i++) 
            N[i].a = N[i].c = N[i].g = N[i].t = 0;
        for(int i = 0; i < n; i++) {
            scanf("%s",grid);
            for(int j = 0; j < m; j++) {
                if(grid[j] == 'A') N[j].a++;
                else if(grid[j] == 'C') N[j].c++;
                else if(grid[j] == 'G') N[j].g++;
                else if(grid[j] == 'T') N[j].t++;
            }
        }
        char s[1005];
        int maxx,ans = 0;
        for(int i = 0; i < m; i++) {
            s[i] = 'A';
            maxx = N[i].a;
            if(N[i].c > maxx) {
                s[i] = 'C';
                maxx = N[i].c;
            }
            if(N[i].g > maxx) {
                s[i] = 'G';
                maxx =N[i].g;
            }
            if(N[i].t > maxx) {
                s[i] = 'T';
                maxx = N[i].t;
            }
            ans += n - maxx;
        }
        s[m] = '\0';
        printf("%s\n%d\n",s,ans);
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值