POJ 2250 Compromise（最长公共子序列）

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int dp[505][505],num1,num2;
char s[505][505],s1[505][505],s2[505][505];
void lcs(int a,int b)
{
if(a==0||b==0) return;
if(s[a][b]==0)
{
lcs(a-1,b-1);printf("%s ",s1[a]);
}
else if(s[a][b]==-1)
{
lcs(a-1,b);
}
else
{
lcs(a,b-1);
}
}
int main()
{
int i,j,k;
while(scanf("%s",s1[1])!=EOF)
{
num1=2,num2=1;
memset(s,0,sizeof(s));
memset(dp,0,sizeof(dp));
while(scanf("%s",s1[num1])&&s1[num1][0]!='#') num1++;
while(scanf("%s",s2[num2])&&s2[num2][0]!='#') num2++;
for(i=1;i<num1;i++)
{
for(j=1;j<num2;j++)
{
if(!strcmp(s1[i],s2[j]))
{
dp[i][j]=dp[i-1][j-1]+1;s[i][j]=0;
}
else
{
if(dp[i][j-1]>dp[i-1][j])
{
dp[i][j]=dp[i][j-1];s[i][j]=1;
}
else
{
dp[i][j]=dp[i-1][j];s[i][j]=-1;
}
}

}
}
lcs(num1-1,num2-1);
printf("\n");
}
return 0;
}