今天重温了一下LCS,温故而知新,对原先学习的东西有了更深的理解。动态规划,状态划分,根据当前状态即过去的总结来决定下一步该怎么做。求最长公共子序列时候,如果两个序列太长了,那么空间开销是相当大的,因此要使用滚动数组的方式来优化空间复杂度。


#include<iostream>
using namespace std;
int ans[101][101],s[101][101],m,n;
char x[100],y[100];
int LCS()
{
memset(ans,0,sizeof(ans));
m=strlen(x);
n=strlen(y);
int i,j;
for(i=1;i<=m;i++)
{
for(j=1;j<=n;j++)
{
if(x[i-1]==y[j-1])
{
ans[i][j]=ans[i-1][j-1]+1;
s[i-1][j-1]=1;
}
else
if(ans[i-1][j]>=ans[i][j-1])
{
ans[i][j]=ans[i-1][j];
s[i-1][j-1]=2;
}
else
{
ans[i][j]=ans[i][j-1];
s[i-1][j-1]=3;
}
}
}
return ans[m][n];
}
void pr(int i,int j)
{
if(i==-1||j==-1)
return ;
if(s[i][j]==1)
{
pr(i-1,j-1);
printf("%c",y[j]);
}
else
if(s[i][j]==2)
pr(i-1,j);
else
pr(i,j-1);
}
int main()
{
while(scanf("%s%s",x,y)==2)
{
printf("%d\n",LCS());
pr(m-1,n-1);
printf("\n");
}
return 0;
}
using namespace std;
int ans[101][101],s[101][101],m,n;
char x[100],y[100];
int LCS()
{
memset(ans,0,sizeof(ans));
m=strlen(x);
n=strlen(y);
int i,j;
for(i=1;i<=m;i++)
{
for(j=1;j<=n;j++)
{
if(x[i-1]==y[j-1])
{
ans[i][j]=ans[i-1][j-1]+1;
s[i-1][j-1]=1;
}
else
if(ans[i-1][j]>=ans[i][j-1])
{
ans[i][j]=ans[i-1][j];
s[i-1][j-1]=2;
}
else
{
ans[i][j]=ans[i][j-1];
s[i-1][j-1]=3;
}
}
}
return ans[m][n];
}
void pr(int i,int j)
{
if(i==-1||j==-1)
return ;
if(s[i][j]==1)
{
pr(i-1,j-1);
printf("%c",y[j]);
}
else
if(s[i][j]==2)
pr(i-1,j);
else
pr(i,j-1);
}
int main()
{
while(scanf("%s%s",x,y)==2)
{
printf("%d\n",LCS());
pr(m-1,n-1);
printf("\n");
}
return 0;
}