hdu 1080 Human Gene Functions

31 篇文章 0 订阅

Problem

acm.hdu.edu.cn/showproblem.php?pid=1080

Meaning

给出一个二维表 similarity[][],表示对应核苷酸配对时的相似度值,横杠 '-' 表示用空格代替一个核苷酸。

给出两个DNA序列 a[] 和 b[],可以任意添加空格,但不能空格对空格。求两串的最大相似度。

Analysis

dp[i][j]:第一串前 i 个和第二串前 j 个配出的最大相似度

在考虑 dp[i][j] 时,a[i] 可能和 b[1] ~ b[j] 的任意一个配,或者一个都不配;b[j] 也可以和 a[1] ~ a[i] 的任何一个配,或一个都不配。

当 a[i] 和 b[k] 配对时,b[k+1] ~ b[j] 的就只能和空格配。

状态转移:dp[i][j] = max { dp[i][0] + dp[0][j] , dp[i-1][s-1] + sim[ a[i] ][ b[s] ] + sum , dp[t-1][j-1] + sim[ a[t] ][ b[j] ] + sum }

其中 sum 表示 b[s+1] ~ b[j] 都和空格配,或 a[t+1] ~ a[i] 都和空格配的总相似度。

注意 a[i] 要和 b[1] ~ b[j] 都匹配一次,b[j] 也要和 a[1] ~ a[i] 都匹配一次。

dp[i][0] 表示第一串前 i 个都和空格配的总相似度;dp[0][j] 类似。

为了好用那个表,可以对DNA串进行格式化,把 A 变成 0,C 变成 1…,空格是 4。

Source code

#include <stdio.h>
#define N 100
#define BIG 12345
#define NUC 4 /* nucleotides */

const int c[5][5] = {
	{5, -1, -2, -1, -3},
	{-1, 5, -3, -2, -4},
	{-2, -3, 5, -2, -2},
	{-1, -2, -2, 5, -1},
	{-3, -4, -2, -1, -BIG}
};
char a[N+3], b[N+3];
int aa[N+1], bb[N+1], dp[N+1][N+1];

void format(char *s, int *d, int l)
{
	for( ; l; --l)
		if(s[l] == 'A')
			d[l] = 0;
		else if(s[l] == 'C')
			d[l] = 1;
		else if(s[l] == 'G')
			d[l] = 2;
		else if(s[l] == 'T')
			d[l] = 3;
}

int max(int x, int y)
{
	return x > y ? x : y;
}

int main()
{
	int t;
	scanf("%d", &t);
	while(t--)
	{
		int la, lb, i, j;
		scanf("%d %s %d %s", &la, a+1, &lb, b+1);
		format(a, aa, la);
		format(b, bb, lb);

		dp[0][0] = 0;
		for(i=1; i<=la; ++i)
			dp[i][0] = dp[i-1][0] + c[aa[i]][NUC];
		for(j=1; j<=lb; ++j)
			dp[0][j] = dp[0][j-1] + c[NUC][bb[j]];

		for(i=1; i<=la; ++i)
			for(j=1; j<=lb; ++j)
			{
				int cnt = 0, k;
				// a[]前i个全和空格配,b[]前j个也全和空格配
				dp[i][j] = dp[i][0] + dp[0][j];
				// a[i]和b[1] ~ b[j]配
				for(k=j; k; --k)
				{
					dp[i][j] = max(dp[i][j],
						dp[i-1][k-1] + c[aa[i]][bb[k]] + cnt);
					cnt += c[NUC][bb[k]];
				}
				// b[j]和a[1] ~ a[i]配
				cnt = 0;
				for(k=i; k; --k)
				{
					dp[i][j] = max(dp[i][j],
						dp[k-1][j-1] + c[aa[k]][bb[j]] + cnt);
					cnt += c[aa[k]][NUC];
				}
			}
		printf("%d\n", dp[la][lb]);
	}
	return 0;
}

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

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值