DNA sequence HDU - 1560

最新推荐文章于 2024-08-03 22:30:59 发布

原创最新推荐文章于 2024-08-03 22:30:59 发布 · 444 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#迭代搜索 #hdu

搜索专栏收录该内容

33 篇文章

订阅专栏

本文介绍了一种基于深度优先搜索的DNA序列匹配算法，通过枚举不同的搜索深度来寻找最佳匹配方案，并采用剪枝策略提高效率。

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

点击打开链接

算是搜索进阶题吧（发现自己搜索要学的东西还很多啊到现在才开始进阶。。）

枚举搜索深度（可二分）看当前深度下是否可以解决问题

这道题目中还需要剪枝若剩下所需匹配数大于剩余步数则无解提前返回

#include <bits/stdc++.h>
using namespace std;

int len[10];
int n,lim,flag;
char dna[10][10];
char per[4];

void dfs(int step,int* cur);
int judge(int* cur);

int main()
{
    int num[10];
    int t,i;
    scanf("%d",&t);
    per[0]='A',per[1]='T',per[2]='C',per[3]='G';
    while(t--)
    {
        scanf("%d",&n);
        for(i=0;i<n;i++)
        {
            scanf("%s",dna[i]);
            len[i]=strlen(dna[i]);
        }
        lim=1,flag=0;
        while(1)
        {
            memset(num,0,sizeof(num));
            dfs(1,num);
            if(flag==1) break;
            lim++;
        }
        printf("%d\n",lim);
    }
    return 0;
}

void dfs(int step,int* cur)
{
    int num[10];
    int i,j;
    if(judge(cur)>lim-step+1) return;
    if(step==lim+1)
    {
        flag=1;
        return;
    }

    for(i=0;i<4;i++)
    {
        for(j=0;j<n;j++)
        {
            if(cur[j]<len[j]&&dna[j][cur[j]]==per[i])
            {
                num[j]=cur[j]+1;
            }
            else
            {
                num[j]=cur[j];
            }
        }
        dfs(step+1,num);
        if(flag==1) return;
    }

    return;
}

int judge(int* cur)
{
    int sum[4],maxx[4];
    int i,j;
    memset(maxx,0,sizeof(maxx));
    for(i=0;i<n;i++)
    {
        memset(sum,0,sizeof(sum));
        for(j=cur[i];j<len[i];j++)
        {
            if(dna[i][j]=='A')
            {
                sum[0]++;
            }
            else if(dna[i][j]=='T')
            {
                sum[1]++;
            }
            else if(dna[i][j]=='C')
            {
                sum[2]++;
            }
            else
            {
                sum[3]++;
            }
        }
        for(j=0;j<4;j++)
        {
            maxx[j]=max(maxx[j],sum[j]);
        }
    }
    return maxx[0]+maxx[1]+maxx[2]+maxx[3];
}