题意:给定一些DNA序列,求一个最短序列能够包含所有序列。

思路:记录第i个序列已经被匹配的长度p[i],以及第i序列的原始长度len[i]。则有两个剪枝:

剪枝1:直接取最长待匹配长度。1900ms

int h = 0;
for(int i = 0; i < n; ++i) {
    h = max(len[i] - p[i], h); //最大待匹配长度
}

剪枝二:统计每个序列里面四种序列值,并求得每种序列值的最长长度。将四种序列值加起来就是最长待匹配长度。180ms

int cal() { //至少还需要匹配的长度
	memset(cost, 0, sizeof(cost));
	memset(tp, 0, sizeof(tp));
	for(int i = 0; i < n; ++i) {
		for(int j = p[i]; j < len[i]; ++j) tp[ha[str[i][j]]]++;
		for(int j = 0; j < 4; ++j) {
			cost[j] = max(cost[j], tp[j]);
			tp[j] = 0;
		}
	}
	int h = 0;
	for(int i = 0; i < 4; ++i) h += cost[i];
	return h;
}

剪枝二优于剪枝一。

AC代码:180ms

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
using namespace std;
map<char, int>ha;
const int maxn = 10;
char str[maxn][maxn];
int len[maxn], p[maxn];
int n, maxd, pd;
int cost[4], tp[4];
char ch[] = {'A', 'G', 'C', 'T'};

int cal() { //至少还需要匹配的长度
	memset(cost, 0, sizeof(cost));
	memset(tp, 0, sizeof(tp));
	for(int i = 0; i < n; ++i) {
		for(int j = p[i]; j < len[i]; ++j) tp[ha[str[i][j]]]++;
		for(int j = 0; j < 4; ++j) {
			cost[j] = max(cost[j], tp[j]);
			tp[j] = 0;
		}
	}
	int h = 0;
	for(int i = 0; i < 4; ++i) h += cost[i];
	return h;
}

int dfs(int cnt) {
	int h = cal();
	if(h == 0) {
		printf("%d\n", maxd);
		return 1;
	}
	if(h + cnt > maxd) {
		pd = min(h + cnt, pd);
		return 0;
	}
	int old[maxn];
	memcpy(old, p, sizeof(old));
	for(int i = 0; i < 4; ++i) {
		char c = ch[i];
		int flag = 0;
		for(int j = 0; j < n; ++j) {
			if(p[j] < len[j] && str[j][p[j]] == c) {
				flag = 1;
				++p[j];
			}
		}
		if(flag && dfs(cnt + 1)) return 1;
		memcpy(p, old, sizeof(old));
	}
	return 0;
}

int main() {
	for(int i = 0; i < 4; ++i) ha[ch[i]] = i;
	int T;
	scanf("%d", &T);
	while(T--) {
		maxd = 0;
		memset(p, 0, sizeof(p));
		scanf("%d", &n);
		for(int i = 0; i < n; ++i) {
			scanf("%s", str[i]);
			len[i] = strlen(str[i]);
			maxd = max(maxd, len[i]);
		}
		while(1) {
			pd = 1 << 30;
			if(dfs(0)) break;
			maxd = pd;
		}
	}
	return 0;
}

如有不当之处欢迎指出!

HDU1560 DNA sequence IDA* + 强力剪枝 [kuangbin带你飞]专题二的更多相关文章

  1. 【算法系列学习三】[kuangbin带你飞]专题二 搜索进阶 之 A-Eight 反向bfs打表和康拓展开

    [kuangbin带你飞]专题二 搜索进阶 之 A-Eight 这是一道经典的八数码问题.首先,简单介绍一下八数码问题: 八数码问题也称为九宫问题.在3×3的棋盘,摆有八个棋子,每个棋子上标有1至8的 ...

  2. HDU - 3567 IDA* + 曼哈顿距离 + 康托 [kuangbin带你飞]专题二

    这题难度颇大啊,TLE一天了,测试数据组数太多了.双向广度优先搜索不能得到字典序最小的,一直WA. 思路:利用IDA*算法,当前状态到达目标状态的可能最小步数就是曼哈顿距离,用于搜索中的剪枝.下次搜索 ...

  3. [kuangbin带你飞]专题二十 斜率DP

            ID Origin Title   20 / 60 Problem A HDU 3507 Print Article   13 / 19 Problem B HDU 2829 Lawr ...

  4. HDU - 3533 bfs [kuangbin带你飞]专题二

    看了好久的样例才看懂. 题意:有一个人要从(0,0)走到(n,m),图中有k个碉堡,每个碉堡可以向某个固定的方向每隔t秒放一次炮,炮弹不能穿越另一个碉堡,会被阻挡.人在移动的过程中不会被炮弹打到,也就 ...

  5. HDU - 1043 A* + 康托 [kuangbin带你飞]专题二

    这题我第一次用的bfs + ELFhash,直接TLE,又换成bfs + 康托还是TLE,5000ms都过不了!!我一直调试,还是TLE,我才发觉应该是方法的问题. 今天早上起床怒学了一波A*算法,因 ...

  6. [kuangbin带你飞]专题二十二 区间DP

            ID Origin Title   17 / 60 Problem A ZOJ 3537 Cake   54 / 105 Problem B LightOJ 1422 Hallowee ...

  7. 【算法系列学习】[kuangbin带你飞]专题二 搜索进阶 D - Escape (BFS)

    Escape 参考:http://blog.csdn.net/libin56842/article/details/41909459 [题意]: 一个人从(0,0)跑到(n,m),只有k点能量,一秒消 ...

  8. HDU - 3001 Travelling 状压dp + 三进制 [kuangbin带你飞]专题二

    终于刷完搜索专题了. 题意:给定n个城市,每个城市参观不能超过两次,两个城市之间有道路通过需要花费X,求通过能所有城市的最小花费. 思路:每个城市有三个状态0,1,2,可用三进制存储所有城市的访问状态 ...

  9. HDU - 3085 双向BFS + 技巧处理 [kuangbin带你飞]专题二

    题意:有两只鬼,一个男孩女孩被困在迷宫中,男孩每秒可以走三步,女孩只能1步,鬼可以两步且可以通过墙.问男孩女孩是否可以在鬼抓住他们之前会合? 注意:每秒开始鬼先移动,然后两人开始移动. 思路:以男孩和 ...

随机推荐

  1. DAY10-万物皆对象-2018-2-2

    许久没有写了,虽然每天都有在学,但是学的东西也少了,后面难度慢慢加大,学习速度也是变慢了.这是许多天积累下来的笔记,从第一次接触对象,到慢慢去了解,现在处于还待深入了解的状态.万物皆对象,那是不是说没 ...

  2. IEEE754 处理数据变换

    public  class IEEE754 { /// <summary> /// 将二进制值转ASCII格式十六进制字符串 /// </summary> /// <pa ...

  3. python字符串问题

    相关知识点: 字符串在Python内部的表示是unicode编码,因此,在做编码转换时,通常需要以unicode作为中间编码,即先将其他编码的字符串解码(decode)成unicode,再从unico ...

  4. 【转】国际GNSS服务组织IGS

    国际GNSS服务组织The International GNSS Service,简称IGS,前身为国际GPS服务组织.IGS提供的高质量数据和产品被用于地球科学研究等多个领域. IGS组织由卫星跟踪 ...

  5. unicode文件处理(如果是ANSI编码就不需要了)

    1.unicode文件的打开必须用rb模式. 3.wchar_t str[100] = { 0 }; 这个占200个字节. 2.宽字符对应的处理 fgetc fgetwc fputc fputwc f ...

  6. 一个 div 手写红绿灯- 分别用css3 和 js 实现

    [要求] 一个div,配合 css3 或者 js 实现红绿灯切换的效果. [思路] 使用 css3,要实现红绿灯颜色的变换必然要用到 animation 动画,通过 keyframes 控制颜色的渐变 ...

  7. 关于HTML的两个实例

    例子一:使用HTML语言编写一个菜单 代码如下: <!DOCTYPE html> <html> <head lang="en"> <met ...

  8. Yii框架里用grid.CGridView调用pager扩展不显示最后一页按钮的解决

    有如下一例,调用zii.widgets.grid.CGridView显示Blog信息,代码如下: $this->widget('zii.widgets.grid.CGridView', arra ...

  9. 2018/2/11 ELK技术栈之ElasticSearch学习笔记二

    终于有时间记录一下最近学习的知识了,其实除了写下的这些还有很多很多,但懒得一一写下了: ElasticSearch添加修改删除原理:ElasticSearch的倒排索引和文档一旦生成就不允许修改(其实 ...

  10. 发放春节福利,ASP.NET Core断点续传

    ASP.NET Core断点续传 在ASP.NET WebAPi写过完整的断点续传文章,目前我对ASP.NET Core仅止于整体上会用,对于原理还未去深入学习,由于有园友想看断点续传在ASP.NET ...