洛谷1140 相似基因

2021-07-10 20:53:09 字數 1658 閱讀 5280

洛谷1140 相似基因

題目背景

大家都知道,基因可以看作乙個鹼基對序列。它包含了4種核苷酸,簡記作a,c,g,t。生物學家正致力於尋找人類基因的功能,以利用於診斷疾病和發明藥物。 在乙個人類基因工作組的任務中,生物學家研究的是:兩個基因的相似程度。因為這個研究對疾病的**有著非同尋常的作用。

題目描述

兩個基因的相似度的計算方法如下:

對於兩個已知基因,例如agtgatg和gttag,將它們的鹼基互相對應。當然,中間可以加入一些空鹼基-,例如:

這樣,兩個基因之間的相似度就可以用鹼基之間相似度的總和來描述,鹼基之間的相似度如下表所示:

那麼相似度就是:(-3)+5+5+(-2)+(-3)+5+(-3)+5=9。因為兩個基因的對應方法不唯一,例如又有:

相似度為:(-3)+5+5+(-2)+5+(-1)+5=14。規定兩個基因的相似度為所有對應方法中,相似度最大的那個。

輸入輸出格式

輸入格式:

共兩行。每行首先是乙個整數,表示基因的長度;隔乙個空格後是乙個基因序列,序列中只含a,c,g,t四個字母。1<=序列的長度<=100。

輸出格式:

僅一行,即輸入基因的相似度。

輸入輸出樣例

輸入樣例#1:

7 agtgatg

5 gttag

輸出樣例#1:

14

【思路】

線性dp。

可以模擬於求最大公共子串行。設d[i][j]為到a串的i與b串的j為止的最大相似度。則有轉移方程:

d[i][j]=max

需要注意的是初始化d[0][j]和d[i][0],另外字串右移一位讀入方便處理。

【**】

1 #include2 #include3

using

namespace

std;45

const

int maxn = 100+10;6

const

int w[5][5]=,

8 ,

9 ,

10 ,

11

12};

13int h[333

],d[maxn][maxn];

14int

n,m;

15char

a[maxn],b[maxn];

1617 inline int cost(char x,char

y) 20

21int

main()

26for(int j=1;j<=m;j++)

2728

for(int i=1;i<=n;i++)

29for(int j=1;j<=m;j++)

34 printf("%d"

,d[n][m]);

35return0;

36 }

洛谷1140相似基因

兩個基因的相似度的計算方法如下 對於兩個已知基因,例如agtgatg和gttag,將它們的鹼基互相對應。當然,中間可以加入一些空鹼基 例如 這樣,兩個基因之間的相似度就可以用鹼基之間相似度的總和來描述,鹼基之間的相似度如下表所示 那麼相似度就是 3 5 5 2 3 5 3 5 9。因為兩個基因的對應...

洛谷P1140 相似基因

大家都知道,基因可以看作乙個鹼基對序列。它包含了4種核苷酸,簡記作a,c,g,t。生物學家正致力於尋找人類基因的功能,以利用於診斷疾病和發明藥物。在乙個人類基因工作組的任務中,生物學家研究的是 兩個基因的相似程度。因為這個研究對疾病的 有著非同尋常的作用。兩個基因的相似度的計算方法如下 對於兩個已知...

洛谷P1140 相似基因

大家都知道,基因可以看作乙個鹼基對序列。它包含了44種核苷酸,簡記作a,c,g,ta,c,g,t。生物學家正致力於尋找人類基因的功能,以利用於診斷疾病和發明藥物。在乙個人類基因工作組的任務中,生物學家研究的是 兩個基因的相似程度。因為這個研究對疾病的 有著非同尋常的作用。兩個基因的相似度的計算方法如...