1. 程式人生 > >nyoj 36 最長公共子序列(動態規劃)

nyoj 36 最長公共子序列(動態規劃)

時間限制:3000 ms | 記憶體限制:65535 KB
難度:3
描述
咱們就不拐彎抹角了,如題,需要你做的就是寫一個程式,得出最長公共子序列。
tip:最長公共子序列也稱作最長公共子串(不要求連續),英文縮寫為LCS(Longest Common Subsequence)。其定義是,一個序列 S ,如果分別是兩個或多個已知序列的子序列,且是所有符合此條件序列中最長的,則 S 稱為已知序列的最長公共子序列。
輸入
第一行給出一個整數N(0< N<100)表示待測資料組數
接下來每組資料兩行,分別為待測的兩組字串。每個字串長度不大於1000.
輸出
每組測試資料輸出一個整數,表示最長公共子序列長度。每組結果佔一行。
樣例輸入
2
asdf
adfsd
123abc
abc123abc
樣例輸出
3
6

經典的dp問題
注意,子序列和子串的定義又是不一樣的,子串要求在原字串中是連續的。而最長公共子序列則並不要求連續。
有兩個字串a和b,長度分別是lena和lenb,那麼
動態轉移方程: 如果a[i-1]==b[j-1],則dp[i][j]=dp[i-1][j-1]+1;
否則,dp[i][j]=max(dp[i-1][j],dp[i][j-1]);

#include <iostream>
#include <string.h>
using namespace
std; int dp[1001][1001]; int main() { int N; cin>>N; while(N--) { string a,b; cin>>a>>b; int i,j,lena=a.length(),lenb=b.length(); memset(dp,0,sizeof(dp)); for(i=1; i<=lena; i++) { for(j=1; j<=lenb; j++) { if
(a[i-1]==b[j-1]) dp[i][j]=dp[i-1][j-1]+1; else if(dp[i-1][j]>=dp[i][j-1]) dp[i][j]=dp[i-1][j]; else dp[i][j]=dp[i][j-1]; } } cout<<dp[lena][lenb]<<endl; } return 0; }