现在的位置: 首页 > 综合 > 正文

hdu1560 DNA sequence (迭代加深搜索)

2013年01月25日 ⁄ 综合 ⁄ 共 1067字 ⁄ 字号 评论关闭

题意:从n个串中找出一个最短的公共串(也许应该说序列吧,因为不要求连续,即只要保持相对顺序就好)。

分析:一开始自以为是爆搜开了,因为已经知道了可以用迭代加深搜索。所谓迭代加深搜索,就是每次都限制了DFS的深度,若搜不到答案,则加深深度,继续搜索,这样就防止了随着深度不断加深而进行的盲目搜索,而且,对于这种求最短长度之类的题目,只要找到可行解,即是最优解了。所以就这样敲完代码了,敲完之后,悲剧TLE。

少了一步十分重要的剪枝,就是每次DFS的时候,都要判断一下,当前的深度+最少还有加深的深度是否大于限制的长度,若是,则退回。

#include<iostream>
#include<algorithm>
#include<queue>
using namespace std;
char str[10][10];
int n,ans,len,size[10];
char DNA[4]={'A','C','G','T'};
void dfs(int cnt,int len1[])
{
	if(cnt>len)
		return ;
	int max=0;
	for(int i=0;i<n;i++)//找出最长还要加深的深度
	{
		int t=size[i]-len1[i];
		if(t>max)
			max=t;
	}
	if(max==0)
	{
		ans=cnt;
		return ;
	}
	if(cnt+max>len)
		return ;
	for(int i=0;i<4;i++)
	{
		int p[10];	
		int flag=0;
		for(int j=0;j<n;j++)
		{
			if(str[j][len1[j]]==DNA[i])//表示该字母可以往下搜索
			{
				flag=1;
				p[j]=len1[j]+1;
			}
			else p[j]=len1[j];
		}
		if(flag)
			dfs(cnt+1,p);
		if(ans!=-1)
			break;
	}
}
int main()
{
	int cas;
	scanf("%d",&cas);
	while(cas--)
	{
		scanf("%d",&n);
		int max1=0;
		for(int i=0;i<n;i++)
		{
			scanf("%s",str[i]);
			size[i]=strlen(str[i]);
			if(size[i]>max1)//找出最长的串的长度,作为初始时的迭代DFS的限制
				max1=size[i];
		}
		ans=-1;
		len=max1;
		int p[10]={0};//记录当前深度下,每一个串已经匹配过的长度
		while(true)
		{
			dfs(0,p);
			if(ans!=-1)
				break;
			len++;//加深迭代
		}
		printf("%d\n",ans);
	}
}

 

抱歉!评论已关闭.