【搜索进阶】HDU 1560 DNA sequence

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1560
IDA*: 容易看出最少步数肯定是大于等于最大字符长度的,所以就构造 h() 为所有剩余长度的最大值
#include<cstdio>
#include<cstring>
#include<bitset>
using namespace std;
#define MAX(a,b) (a>b?a:b)
char code[]="ATGC";
int pos[10],n;
int h(){
    int ans=0;
    for(int i=0;i<n;i++){
        ans=MAX(pos[i],ans);
    }
    return ans;
} 
int deep;
char maze[10][8];
bool dfs(int d){
    int th=h();
    if(th==0) return true;
    if(d+th>deep) return false;
    for(int i=0;i<4;i++){
        bool flag=false;
        int tpos[10];
        for(int j=0;j<n;j++){
            tpos[j]=pos[j];
            if(pos[j]==0) continue;
            if(maze[j][pos[j]-1]==code[i]){
                pos[j]--;
                flag=true;
            }
        } 
        if(flag){
            if(dfs(d+1)) return true;
        }
        for(int j=0;j<n;j++){
            pos[j]=tpos[j];
        }
    }
    return false;
} 
int main(){
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%d",&n);
        for(int i=0;i<n;i++){
            scanf("%s",maze[i]);
            pos[i]=strlen(maze[i]);
        }
        deep=0;
        while(true){
            if(dfs(0)) break;
            deep++; 
        }
        printf("%d\n",deep);
    }
    return 0;
}

单广: 开int进行hash,每组cas的hash值不同,这样就省去了每次初始化的消耗
#include "stdio.h"
#include "string"
#define M 1679616+100000
struct H{
    int st;
    int len;
}q[M];
int hash[M],cas =1;
int n , end ,len;
char str[8][8];
char ku[] ="ACGT";
int bfs() {
    int head,tail,i,st;
    int pos[8];
    q[0].len =0;
    q[0].st =0;
    head = tail =0;
    while(head <= tail) {
        st = q[head].st;
        for(i =0 ; i < n ; i ++) {
            pos[i] = st%len;
            st/=len;
        }
        for(int j =0 ; j <4 ; j ++) {
            st =0;
            for(i = n-1; i >=0 ;i --) {
                st = st * len + pos[i] + (str[i][pos[i]] == ku[j]);
            }
            if(hash[st] == cas) { continue; }
            if(st == end) { return q[head].len +1; }
            hash[st] = cas;
            tail ++;
            q[tail].len = q[head].len +1;
            q[tail].st = st;
        }
        head ++;
    }
    return -1;
}
int main() {
    int T , i ;
    scanf("%d",&T);
    while(T --) {
        scanf("%d",&n);
        len =0;
        for(i =0 ; i < n ; i ++) {
            scanf("%s",str[i]);
            if(strlen(str[i]) > len) {
                len = strlen(str[i]);
            }
        }
        len ++;
        end =0;
        for(i = n-1 ; i >=0 ; i --) {
            end = end * len + strlen(str[i]);
        }
        printf("%d\n",bfs());
        cas ++;
    }
    return 0;
}





  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值