最长公共子序列（LCS）

用dp[i][j]表示长度为i的序列和j的序列的最长公共子序列，接下来，状态转移很不好想，我理解了3个小时才搞懂，假若a1[i]==a2[j]也就是说序列a1的最后一个数和序列a2的最后一个数相等， 这种情况下dp[i][j]=dp[i-1][j-1]+1,如果序列a1的最后一个数不等于序列a2的最后一个数，那么最长公共序列肯定不包括这两个序列里面其中一个最后一个数，这时候就看可以扔掉谁了，假若不包含a1的最后一个数，那么就是dp[i][j]=dp[i-1][j]+1,若不包含a2的最后一个数，那么dp[i][j]=dp[i][j-1].
接下来的输出也是有一定技巧的，这个还只是输出其中一个解的版本，还有输出全部解的版本，以后学习…

#include<bits/stdc++.h>
#define maxn 1000
using namespace std;

int a1[maxn];
int a2[maxn];
int c[maxn];  //打印子序列 
int dp[maxn][maxn];
int main(){

int n1,n2;
memset(dp,0,sizeof(dp));
scanf("%d",&n1);
for(int i=0;i<n1;i++)
scanf("%d",&a1[i]);

scanf("%d",&n2);
for(int i=0;i<n2;i++)
scanf("%d",&a2[i]);


int i=0,j=0;

for(int i=0;i<n1;i++)
for(int j=0;j<n2;j++){
	if(a1[i]==a2[j])
	dp[i][j]=dp[i-1][j-1]+1;
	else if(dp[i-1][j]>dp[i][j-1])
	dp[i][j]=dp[i-1][j];
	else
	dp[i][j]=dp[i][j-1];
} 


i=n1-1,j=n2-1;
int z=0;
while(i!=0 && j!=0){
	if(a1[i]==a2[j]){
	c[z++]=a1[i];		
	i--,j--;
	}
	else if(dp[i][j-1]>dp[i-1][j])	
	j--;
	else if(dp[i][j-1]<=dp[i-1][j])
	i--;


}

printf("最长公共子序列为%d",dp[n1-1][n2-1]);
printf("输出：\n");

for(int i=z-1;i>=0;i--)
printf("%d ",c[i]);	

}