HDU1560 DNA sequence (IDA*)

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

SaltyFishWei

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

阅读量417

点赞数

分类专栏：搜索文章标签： IDAstar

本文链接：https://blog.csdn.net/cqbzwja/article/details/47087639

版权

搜索专栏收录该内容

21 篇文章 0 订阅

订阅专栏

题目大意：给你n(n<=8)个字符串，每个字符串的长度不超过5，现在让你找出一个最短的母串，使得这n个字符串为这个母串的子串。

思路：搜索。很容易看出来最少步数大于等于最长的字符长度，因此估值函数h() = 同一字符串中所有剩余没匹配字母长度的最大值

#include<cstdio>
#include<cstring>
#define Max(a,b) a>b?a:b
using namespace std;
int n;
struct T
{
	char s[10];
	int len;
}a[10];
int pos[10];//字符串i现在匹配到的位置
char letter[5] = {'A','G','C','T'};
int maxn;
int depth;
int hash[30];

int get_h()//估值函数
{
	int w[4] = {0};//所有剩余长度的最大值
	for(int i = 1; i <= n; i++)
	{
		int t[4] = {0};
		for(int j = pos[i]; j < a[i].len; j++)
		{
			t[ hash[ a[i].s[j] - 'A' ] ]++;
		}
		for(int j = 0; j < 4; j++) w[j] = Max(w[j],t[j]);
	}
	return w[0] + w[1] + w[2] + w[3];
}

bool IDAstar(int cur)
{
	if(cur + get_h() > depth) return false;//当前进行的深度加上至少还要匹配的深度，如果它比限定的深度还要大，那么我们就可以剪枝了
	if(get_h() == 0) return true;
	
	int temp[10];
	bool flag = 0;
	for(int i = 0; i < 4; i++)//枚举字母
	{
		memcpy(temp,pos,sizeof pos);
		for(int j = 1; j <= n; j++)
		{
			if(letter[i] == a[j].s[pos[j]])//当前一个字符可以匹配，因此该字符串向后移动一位
			{
				flag = 1;
				pos[j]++;
			}
		}
		if(flag)
		{
			if(IDAstar(cur+1)) 
				return true;//迭代下一层
			memcpy(pos,temp,sizeof temp);//复原
		}
	}
	return false;
}

int main()
{
	int T;
	scanf("%d",&T);
	hash[0] = 0; hash['G'-'A'] = 1; hash['C'-'A'] = 2; hash['T'-'A'] = 3;
	while(T--)
	{
		scanf("%d",&n);
		maxn = 0;
		for(int i = 1; i <= n; i++)
		{
			scanf("%s",a[i].s);
			a[i].len = strlen(a[i].s);
			if(a[i].len > maxn)
				maxn = a[i].len;
			pos[i] = 0;
		}
		for(depth = maxn; ; depth++)//枚举匹配后长度
		{
			if(IDAstar(0)) 
			{
				printf("%d\n",depth);
				break;
			}
		}
	}
}