作业9

1.问题
[描述算法问题，首选形式化方式（数学语言），其次才是非形式化方式（日常语言）]
给定两个字符串，求解这两个字符串的最长公共子序列（Longest Common Sequence）。
比如字符串1：BDCABA；字符串2：ABCBDAB则这两个字符串的最长公共子序列长度为4，最长公共子序列是：BCBA。
2.解析
S1的最后一个元素与s2的最后一个元素相同，这说明该元素一定位于公共子序列中。因此，现在只需要找：
LCS(s1n-1，s2m-1);
如果s1n != s2m，那么产生了两个子问题：
LCS(s1n-1，s2m) 和 LCS(s1n，s2m-1)。
解出上面两个子问题，得到的公共子序列后，比较长短，长的就是 LCS（s1,s2）。
对子问题重复以上步骤，就可以得到答案。
这张图对应了三种情况。C[i,j]记录的是子序列长度。
3.设计
[核心伪代码]
void LCS(char *s1, char *s2, int len1, int len2, int c[][MAXLEN], int dir[][MAXLEN]){
int i, j;
for(i = 0; i <= len1; i++) //将第一列设置0
c[i][0] = 0;
for(j = 1; j <= len2; j++) //第一行设置为0
c[0][j] = 0;
for(i = 1; i<= len1; i++){ //用1,0,-1三个整形作为上下斜上的方向标记
for(j = 1; j <= len2; j++){
if(s1[i-1] == s2[j-1]){ //如果最后一个字符相同，那么在剩下的里找
c[i][j] = c[i-1][j-1] + 1;
dir[i][j] = 1;
}
else if(c[i-1][j] >= c[i][j-1]){
c[i][j] = c[i-1][j];
dir[i][j] = 0;
}
else{
c[i][j] = c[i][j-1];
dir[i][j] = -1;
}
}
}
}
4.分析
[算法复杂度推导]
5.源码
#include<stdio.h>
#include<string.h>
const int maxn=1000+10;
char x[maxn],y[maxn];
int dp[maxn][maxn];
int b[maxn][maxn]; //b[i][j]=0为删除两个，b[i][j]=1为删除x，b[i][j]=2为删除y
char lcs[maxn];
int n,m; //n为字符串x的长度，m为字符串y的长度
int main(){
scanf("%d",&n);
scanf("%s",x+1);
scanf("%d",&m);
scanf("%s",y+1);
//初始化
for(int i=0;i<=m;++i)dp[0][i]=0;
for(int i=0;i<=n;++i)dp[i][0]=0;

/*----------以下为求最长公共子序列长度----------*/ 
for(int i=1;i<=n;++i){
	for(int j=1;j<=m;++j){
		if(x[i]==y[j]){
			dp[i][j]=dp[i-1][j-1]+1;
		}
		else{
			if(dp[i-1][j]>dp[i][j-1]){
				b[i][j]=1;
				dp[i][j]=dp[i-1][j];
			}
			else{
				b[i][j]=2;
				dp[i][j]=dp[i][j-1];
			}
		}
	}
}
/*----------以上为求最长公共子序列长度----------*/ 

/*----------以下为求最长公共子序列----------*/
int i=n,j=m;
while(i!=0&&j!=0){     
	if(b[i][j]==0){
		lcs[dp[i][j]]=x[i];
		i--,j--;
	}
	else if(b[i][j]==1){
		i--;
	}
	else{
		j--;
	}
}
/*----------以上为求最长公共子序列----------*/ 

printf("最长公共子序列长度为：%d\n",dp[n][m]);
if(dp[n][m]!=0){
	printf("最长公共子序列为：");
	for(int i=1;i<=dp[n][m];++i){
		printf("%c",lcs[i]);
	}
}

}