hdu 1560 DNA sequence(IDA*)

原创 2015年11月18日 20:12:42

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1560

DNA sequence

Time Limit: 15000/5000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 1505    Accepted Submission(s): 730


Problem Description
The twenty-first century is a biology-technology developing century. We know that a gene is made of DNA. The nucleotide bases from which DNA is built are A(adenine), C(cytosine), G(guanine), and T(thymine). Finding the longest common subsequence between DNA/Protein sequences is one of the basic problems in modern computational molecular biology. But this problem is a little different. Given several DNA sequences, you are asked to make a shortest sequence from them so that each of the given sequence is the subsequence of it.

For example, given "ACGT","ATGC","CGTT" and "CAGT", you can make a sequence in the following way. It is the shortest but may be not the only one.

 

Input
The first line is the test case number t. Then t test cases follow. In each case, the first line is an integer n ( 1<=n<=8 ) represents number of the DNA sequences. The following k lines contain the k sequences, one per line. Assuming that the length of any sequence is between 1 and 5.
 

Output
For each test case, print a line containing the length of the shortest sequence that can be made from these sequences.
 

Sample Input
1 4 ACGT ATGC CGTT CAGT
 

Sample Output
8
 

Author
LL
 

Source
 

Recommend
LL   |   We have carefully selected several similar problems for you:  1667 1043 1813 1226 1401 

题目大意:给n个序列,找到一个包含所有给出序列的最短长度并输出。
解题思路:采用gei_h()得到当前状态下最长的未匹配的长度。在进行深度搜索。每个串的长度不超过5,最多只有8个序列,所以IDA不超过40次。

详见代码。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>

using namespace std;

int n;
char ch[10][10];
int len[10],want;
char dir[10]= {'A','C','G','T'};
int wei[10];//记录第i个序列正在使用第几个位置

int get_h()
{
    int t=0;
    for (int i=1; i<=n; i++)
    {
        t=max(t,len[i]-wei[i]);//得到当前情况下最长的未被匹配的长度
    }
    return t;
}

int IDA(int dep)
{
    if(dep+get_h()>want)//当前长度+估测的长度比我想要的还大的话,就不必继续搜索
    {//cout<<get_h()<<endl;
        return 0;}
    if(dep==want)
        return 1;
    int tmp[10];
    for (int i=0; i<4; i++)
    {
        int flag=0;
        memcpy(tmp,wei,sizeof(wei));//先存一下
        for (int j=1; j<=n; j++)
        {
            if (ch[j][wei[j]]==dir[i])//当前的这一位符合
            {
                flag=1;//标记当前状态可以
                wei[j]++;
            }
        }
        if (flag)
        {
            if(IDA(dep+1))//如果可以就继续搜索
                return 1;
            memcpy(wei,tmp,sizeof(tmp));//还原回来
        }
    }
    return 0;
}

int main()
{
    int T;
    scanf("%d",&T);
    while (T--)
    {
        int Max=0;
        scanf("%d",&n);
        for (int i=1; i<=n; i++)
        {
            scanf("%s",ch[i]);
            len[i]=strlen(ch[i]);
            if (len[i]>Max)
                Max=len[i];
        }
        memset(wei,0,sizeof(wei));
        want=Max;//从最长序列开始查找
        while (1)
        {
            if (IDA(0))
            {
                break;
            }
            want++;
        }
        printf ("%d\n",want);
    }
    return 0;
}




 

版权声明:本文为博主原创文章,未经博主允许不得转载。

hdu 1560 DNA sequence(迭代加深搜索)

题意:找到一个最短的序列,使得输入的所有序列都是这个序列的子序列。输出最短序列的长度 DNA sequence Time Limit: 15000/5000 MS (Java/Ot...
  • t1019256391
  • t1019256391
  • 2013年06月14日 15:28
  • 1259

每天一道算法_8_DNA Sorting

Description One measure of ``unsortedness'' in a sequence is the number of pairs of entries that are...
  • jason0539
  • jason0539
  • 2013年10月30日 17:37
  • 2544

HDU 1003 解题报告

Max Sum Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total ...
  • u013341274
  • u013341274
  • 2014年05月24日 19:09
  • 2586

HDU1560:DNA sequence(IDA星)

Problem Description The twenty-first century is a biology-technology developing century. We know th...
  • libin56842
  • libin56842
  • 2014年04月09日 21:40
  • 3064

HDU1560 DNA sequence (IDA*)

题目大意:给你n(n 思路:搜索。很容易看出来最少步数大于等于最长的字符长度,因此估值函数h() =  同一字符串中所有剩余没匹配字母长度的最大值 #include #includ...
  • cqbzwja
  • cqbzwja
  • 2015年07月27日 17:09
  • 165

HDU1560 DNA sequence(IDA星) (E)

Problem Description The twenty-first century is a biology-technology developing century. We know th...
  • qq_36782366
  • qq_36782366
  • 2017年07月09日 10:02
  • 72

HDU 1560 DNA sequence IDA*搜索

转载请注明出处,谢谢 http://blog.csdn.net/ACM_cxlove?viewmode=contents           by---cxlove 继续IDA*,明显可以想到的一个...
  • ACM_cxlove
  • ACM_cxlove
  • 2012年07月14日 00:15
  • 1591

HDU 1560 DNA sequence(IDA*)

DNA sequence Time Limit: 15000/5000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) ...
  • Yasola
  • Yasola
  • 2017年07月01日 11:14
  • 125

hdu1560 DNA sequence (IDA*)

题意: 给出几个字符串,求一个最短的字符串使得给出的字符串都是该字符串的子序列。 思路: 最多8个字串,长度不超过5,所以所求串的长度不超过40,用IDA*。 每次枚举给主串后加的字符,如果给...
  • u011265346
  • u011265346
  • 2014年12月25日 17:57
  • 422

HDU 1560 DNA sequence(IDA*搜索)

题意:给出n个dna序列,现在要求一个最短的序列,使得给出的dna序列都是它的子序列。 思路:IDA*搜索,一开始用的剪枝估价条件是计算剩余字符串的最大长度,但是总是超时。 后来更改了剪枝条件,估价函...
  • u014664226
  • u014664226
  • 2016年01月31日 21:57
  • 216
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:hdu 1560 DNA sequence(IDA*)
举报原因:
原因补充:

(最多只允许输入30个字)