51Nod 1006 最长公共子序列Lcs

给出两个字符串A B，求A与B的最长公共子序列（A,B的长度 <= 1000，子序列不要求是连续的）。

比如两个串为：

A：abcicba B：abdkscab

ab是两个串的子序列，abc也是，abca也是，其中abca是这两个字符串最长的子序列。。

输入

第1行：字符串A
第2行：字符串B

输出

输出最长的子序列，如果有多个，随意输出1个。

输入样例

abcicba
abdkscab

输出样例

abca

参考程序

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
int main()
{
	short dp[1005][1005];
	char luJing[1005][1005];
	char a[1005],b[1005];
	memset(dp,0,sizeof(dp));
	memset(luJing,0,sizeof(luJing));
	 gets(a);
	 gets(b);
	int la=strlen(a),lb=strlen(b);
	for(int i=0;i<la;i++)
	{
		for(int j=0;j<lb;j++)
		{
			if(a[i]==b[j])
			{                                                         //'d' 'z'  's'  分别代表来自对角 左边 和上边
				dp[i+1][j+1]=dp[i][j]+1;
				luJing[i][j]='d';
			}
			else
			{
				if(dp[i+1][j]>dp[i][j+1])
				{
					dp[i+1][j+1]=dp[i+1][j];
					luJing[i][j]='z';
				}
				else
				{
					dp[i+1][j+1]=dp[i][j+1];
					luJing[i][j]='s';
				}
			}
		}
	}
	/*puts(a);
	puts(b);
	for(int i=0;i<=la;i++)
	{
		for(int j=0;j<=lb;j++)
		{
			cout<<luJing[i][j]<<" ";
		}
		cout<<endl;
	}*/
	int x=la-1,y=lb-1,ls=0;
	char s[1005];
	memset(s,0,sizeof(s));
	while(x>=0&&y>=0)                        //回溯找路径
	{
		if(luJing[x][y]=='d')
		{
			s[ls++]=a[x];
			//cout<<x<<" "<<y<<endl;
			x--; 
			y--;
		}
		else if(luJing[x][y]=='s')
		{
			x--;
		}
		else
		{
			y--;
		}
	}
	//cout<<s[0]<<endl;
	while(ls--)
		cout<<s[ls];
	cout<<endl;
	return 0;
}