UVA 1368 - DNA Consensus String(贪心)

最新推荐文章于 2023-01-19 17:29:00 发布

HelloWorld10086

最新推荐文章于 2023-01-19 17:29:00 发布

阅读量687

点赞数

CC 4.0 BY-SA版权

分类专栏：贪心文章标签： uva 1368

本文链接：https://blog.youkuaiyun.com/HelloWorld10086/article/details/43954599

贪心专栏收录该内容

36 篇文章

订阅专栏

本文介绍了一种针对多个DNA序列寻找最优目标序列的方法，通过贪心策略选取每位置上出现频率最高的碱基来构建目标序列，并计算了将所有原始序列转换为目标序列所需的最小变化次数。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题意：

给定n个DNA，求一个目标串，使得每个序列变换成目标所需次数最少，且字典序最小。

思路：

贪心。每一位选存在最多的变化即可。

AC代码

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdlib>
using namespace std;
typedef long long ll;
char dna[55][1010];
const char s[] = "ACGT";
int map[256];
char str[1010];

int m,n;
void init() {
    memset(map,0,sizeof(map));
    map['A'] = 0; map['C'] = 1;
    map['G'] = 2; map['T'] = 3;
}
int main() {
    init();
    int T;
    scanf("%d",&T);
    while(T--) {
        scanf("%d%d",&m,&n);
        for(int i = 0; i < m; i++) {
            scanf("%s",dna[i]);
        }
        int vis[5];
        int cnt = 0;
        for(int j = 0; j < n; j++) {
            memset(vis, 0,sizeof(vis));
            for(int i = 0; i < m; i++) {
                vis[map[dna[i][j]]]++;
            }
            int flag = 0;
            int Max = vis[0];
            for(int i = 0; i < 4; i++) {
                if(Max < vis[i]) {
                    Max = vis[i];
                    flag = i;
                }
            }
            str[j] = s[flag];
            cnt += (m - Max);
        }
        str[n] = '\0';
        printf("%s\n",str);
        printf("%d\n",cnt);
    }
    return 0;
}