模板

<pre name="code" class="cpp"><span style="font-size:18px;">#include<cstdio>
#include<stdlib.h>
#include<string.h>
#include<string>
#include<map>
#include<cmath>
#include<iostream>
#include <queue>
#include <stack>
#include<algorithm>
#include<set>
using namespace std;
#define INF 1e8
#define eps 1e-8
#define LL long long 
#define N 1005
#define mod  1000000007 
char a[N],b[N],str[N];  
int c[N][N];  
//求LCS的长度：  
int lcs_len(char *a,char *b,int c[][N])  
{  
	int m=strlen(a),n=strlen(b),i,j;  
	for(i=0;i<=m;i++)  
		c[i][0]=0;  
	for(j=0;j<=m;j++)  
		c[0][j]=0;  
	for(i=1;i<=m;i++)  
	{  
		for(j=1;j<=n;j++)  
		{  
			if(a[i-1]==b[j-1]) c[i][j]=c[i-1][j-1]+1;  
			else if(c[i-1][j]>=c[i][j-1]) c[i][j]=c[i-1][j];  
			else c[i][j]=c[i][j-1];  
		}  
	}  
	return c[m][n];  
}  
//求LCS的字符串：  
char *build_lcs(char s[],char *a,char *b)  
{  
	int i=strlen(a),j=strlen(b);  
	int k=lcs_len(a,b,c);  
	s[k]='\0';  
	while(k>0)  
	{  
		if(c[i][j]==c[i-1][j]) i--;  
		else if(c[i][j]==c[i][j-1]) j--;  
		else  
		{  
			s[--k]=a[i-1];  
			i--;  
			j--;  
		}  
	}  
	return s;  
}  
int main()  
{  
	while(~scanf("%s%s",a,b))  
	{  
		printf("%d\n",lcs_len(a,b,c));  
		printf("%s\n",build_lcs(str,a,b));  
	}  
	return 0;  
}  </span>

************************************************************************ 设有二维数组 f[i][j] 表示 X 的 i 位和 Y 的 j 位之前的最长公共子序列的长度，则有：f[1][1] = same(1,1)f[i][j] = max{f[i-1][j-1] + same(i,j),f[i-1][j],f[i][j-1]}其中，same(a,b)当 X 的第 a 位与 Y 的第 b 位完全相同时为“1”，否则为“0”。此时，f[i][j]中最大的数便是 X 和 Y 的最长公共子序列的长度，依据该数组回溯，便可找出最长公共子序列。该算法的空间、时间复杂度均为O(n^2)，经过优化后，空间复杂度可为O(n)，时间复杂度为O(nlogn)。HDU1080、PKU1080，ZJU1027：初始化：dp[i][0]=dp[i-1][0]+score[ s1[i-1] ]['-']; dp[0][j]=dp[0][j-1]+score['-'][ s2[j-1] ]; 设dp[i][j]为取s1第i个字符，s2第j个字符时的最大分值则决定dp为最优的情况有三种（score[][]为s1[i]和s2[j]两符号的分数）：1、 s1取第i个字母，s2取“ - ”： dp[i-1][j]+score[ s1[i-1] ]['-'];2、 s1取“ - ”，s2取第j个字母：dp[i][j-1]+score['-'][ s2[j-1] ];3、 s1取第i个字母，s2取第j个字母：dp[i-1][j-1]+score[ s1[i-1] ][ s2[j-1] ]; 即dp[i][j]=max( dp[i-1][j]+score[ s1[i-1] ]['-'],dp[i][j-1]+score['-'][ s2[j-1] ],dp[i-1][j-1]+score[ s1[i-1] ][ s2[j-1] ] );for(i=1;i<=len1;i++) for(j=1;j<=len2;j++) { int temp1=dp[i-1][j]+score[ s1[i-1] ]['-']; int temp2=dp[i][j-1]+score['-'][ s2[j-1] ]; int temp3=dp[i-1][j-1]+score[ s1[i-1] ][ s2[j-1] ]; dp[i][j]=max(temp1,temp2,temp3); }