【IDA*】hdu 1560

本文深入剖析了一道ACM竞赛中的高难度问题,通过使用A*算法和深度优先搜索(DFS)结合剪枝技巧来解决一个关于DNA序列匹配的问题。通过对算法的详细解释和源代码分析,展示了如何有效地解决问题并找到最优解。

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

http://acm.hdu.edu.cn/showproblem.php?pid=1560

分析:这已经是神一般的剪枝了,膜拜……


#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;


const int NM=10;
char str[NM][NM],cs[4]={'A','C','T','G'};
int v[NM],len[NM],sta[NM],dep,n;
bool flag;


int compare()  //A*的估价函数
{
	int i,t,mmax;
	mmax=0;
	for(i=0;i<n;i++){
		t=len[i]-v[i];
		if(t>mmax) mmax=t;
	}
	return mmax;
}


bool DFS(int t)
{
	int i,j,leave;
	int sta[NM];  //
	leave=compare();
	if(leave==0) return true;
	if(t+leave>dep) return false;  //剩下需要匹配的单词个数+当前已匹配的单词个数
	for(i=0;i<n;i++) sta[i]=v[i];
	for(j=0;j<4;j++){
		flag=false;
		for(i=0;i<n;i++){
			if(str[i][v[i]]==cs[j])  //匹配一个字符串的第v[i]个字符
			{
				flag=true;
				v[i]++;
			}
		}
		if(!flag) continue;  //没有一个匹配成功,说明该字符无效
		if(DFS(t+1)) return true;
		for(i=0;i<n;i++) v[i]=sta[i];  //
	}
	return false;
}


int main()
{
	int T,i;
	scanf("%d",&T);
	while(T--){
		scanf("%d",&n);
		dep=0;
		for(i=0;i<n;i++){
			scanf("%s",str[i]);
			len[i]=strlen(str[i]);
			if(dep<len[i]) dep=len[i];
		}
		memset(v,0,sizeof(v));
		while(!DFS(0)) //为了找出最少的所需添加的字母数
			dep++;
		printf("%d\n",dep);
	}
	return 0;
}



评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值