洛谷P1140 基因匹配 DP真正意義上的一血

2022-04-14 07:43:11 字數 1901 閱讀 5459

大家都知道,基因可以看作乙個鹼基對序列。它包含了44種核苷酸,簡記作a,c,g,ta,c,g,t。生物學家正致力於尋找人類基因的功能,以利用於診斷疾病和發明藥物。

在乙個人類基因工作組的任務中,生物學家研究的是:兩個基因的相似程度。因為這個研究對疾病的**有著非同尋常的作用。

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

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

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

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

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

輸入格式:

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

輸出格式:

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

輸入樣例#1: 

7 agtgatg

5 gttag

輸出樣例#1:

14
其實本蒟蒻以前做過dp題 但是這是第一次自己構造dp方程並且ac 

對於我來說這個題的難點在於判斷是誰插誰= =(聽著好奇怪

對於每乙個f[i][j]

可以由f[i][j-1]轉移而來

可以由f[i-1][j]轉移而來

可以由f[i-1][j-1]轉移而來

現在來討論每一種情況

新來的b[j]要與乙個空鹼基匹配

相當於在a串上插了乙個空鹼基 並與b[j]進行匹配 f[i][j] = f[i][j-1] + trs[b[j][0];

新來的a[i]要與乙個空鹼基匹配

f[i][j] = f[i-1][j] + trs[a[i]][0];

相互匹配 f[i][j] =f[i-1] [j-1] + trs[a[i]][b[j]];

所以構造完成 迴圈順序普通即可(想思考也可以想一下)

那麼邊界?可以想到f[i][0] =f[i-1][0] + trs[a[i]][0];

b也同理

分析完啦 上**:

#includeusing

namespace

std;

int f[105][105

];int a[105],b[105

],len2,len1;

intmain(),,,

,};cin>>len1;

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

cin>>len2;

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

for(int i=1;i<=len1;i++)}f[

0][0] = 0

;

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

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

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

}cout

return0;

}

真的是好不容易啊 嗚嗚嗚嗚嗚我好弱2節課做了個這sb題;

但是感覺dp思維提公升一些 也算沒虧吧

洛谷P1140 相似基因

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

洛谷P1140 相似基因

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

洛谷P1140 相似基因

本題一看就知道是一道動歸,其實和字串距離非常的像,只不過多了題目規定的匹配相似度罷了。int shuzu 6 6 ps 換行效果更佳。int start char c 然後直接在主函式中呼叫即可 for int i 0 i其中s1,s2均為讀入的字串,我們把它們分別逐字元轉化放入a,b陣列。for ...