经典的lcs问题（输出公共子序列）

这个问题疑惑了很久，只知道求公共子序列长度，不会求子序列。

代码贴出如下：

#include<iostream>
using namespace std;
#include<algorithm>
#include<stack>

int chang[1000][1000]={};
int flag[1000][1000]={};

char * lcs(char *a,char *b)
{
	int len1=strlen(a);
	int len2=strlen(b);
	//for(int i=0;i<=len1;++i)
	//	chang[i][0]=0;
	//for(int i=0;i<=len2;++i)
	//	chang[0][i]=0;
	for(int i=1;i<=len1;++i)
			for(int j=1;j<=len2;++j)
			{
				if(a[i-1]==b[j-1])
				{
					chang[i][j]=chang[i-1][j-1]+1;
					flag[i][j]=0;
				}else if(chang[i][j-1]>=chang[i-1][j]){
					chang[i][j]=chang[i][j-1];
					flag[i][j]=-1;
				}else{
					chang[i][j]=chang[i-1][j];
					flag[i][j]=1;
				}
			}
			int len3=chang[len1][len2];
			char *p=new char[len3+1];
			p[len3--]='\0';
			//stack<char> zifu;
			int i=len1;
			int j=len2;
			while(j>0&&i>0)
			{
				if(flag[i][j]==0){
					//zifu.push(b[j-1]);
					p[len3--]=b[j-1];
					i--,j--;
				}else if(flag[i][j]==-1)
				{
					j--;
				}else if(flag[i][j]==1)
					i--;
			}

			//int len3=zifu.size();
			
			//for(int i=0;i<len3;++i)
			//{
			//	p[i]=zifu.top();
			//	zifu.pop();
			//}

			return p;
}

int main()
{
	char a[]="cafe";
	char b[]="coffee";
	char *p=lcs(a,b);
	cout<<p<<endl;
	return 0;
}