经典的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;
}