DNA sequence HDU – 1560 (IDA*模板)

题目链接:https://cn.vjudge.net/problem/HDU-1560
题目大意:给定数个DNA序列,构造出一个最短的DNA序列,使得所有给定的DNA序列都是它的子序列
解题思路:基于迭代加深的A*算法
代码

#include <iostream>
#include <bits/stdc++.h>
using namespace std;
int n;
int num[11],point[11];
char c[11][11];
char str[]={'A','C','G','T'};
int f;

int judge()
{
    int l=0;
    for(int i=0;i<n;i++)
    {
        l=max(l,num[i]-point[i]);该序列的长度减掉已构造的长度即为该序列距离构造完成所需的最小长度
    }
    return l;
}

void dfs(int s)
{
    if(judge()==0)//如果还需要构造的长度等于0,则构造成功
    {
        f=1;
        return;
    }
    if(judge()>s)//如果距离构造成功还需的长度大于所剩步骤数s,则一定无法完成
    {
        return;
    }
    int ff=0;
    int bf[11];
    for(int i=0;i<n;i++)
    {
        bf[i]=point[i];//将point数组备份
    }
    for(int i=0;i<4;i++)
    {
        ff=0;
        for(int j=0;j<n;j++)//是否可以构造
        {
            if(str[i]==c[j][point[j]])
            {
                point[j]++;
                ff=1;
            }
        }
        if(ff==1)//可以则dfs
        {
            dfs(s-1);
            if(f==1)
                return;
            else//无法成功恢复point数组
            {
                for(int j=0;j<n;j++)
                {
                    point[j]=bf[j];
                }
            }
        }
    }
}

int main()
{
    ios::sync_with_stdio(false);
    int t;
    cin>>t;
    while(t--)
    {
        memset(point,0,sizeof(point));//第i个序列构造完成的长度
        cin>>n;
        int l=0;
        for(int i=0;i<n;i++)
        {
            cin>>c[i];
            num[i]=strlen(c[i]);//记录每个序列的长度
            l=max(l,num[i]);//记录序列中的最长长度
        }
        f=0;
        while(1)
        {
            dfs(l);//从序列中的最长长度开始入手
            if(f==1)
            {
                cout<<l<<endl;
                break;
            }
            l++;//无法构成则将长度+1
        }
    }
    return 0;
}

发表评论

电子邮件地址不会被公开。 必填项已用*标注