UVa 10405

动态规划的经典问题， 最长公共子序列 .
方程:

这题提交了6次，一直WA，找标程对拍才发现序列可以包含空格...于是用gets，AC.

另外，降维以后通过了样例，但是WA，原因未知 .

这里的 描述 很有启发性，现转载一部分：

本题初看有些类似线形动态规划 ，但经过仔细审题发现用矩阵类动态规划 更好理解， 如下图所显示，我们把两个字符串看作矩阵的横纵轴： 首先读入数据时，把两个字符相同的标记成1，不相同的地方标记成0；

很容易证明，如果能形成公共子序列，则最长子序列一定是从某个位置开始的对角线的长度，所以我们需要做的就是统计对角线的长度，DP方程就为if a[i,j]=1 then d[i,j]:=d[i-1,j-1]+1;

经过循环我们就可以把上表转化成:

因此我们只需记录在循环中的最大值 就是我们要的最长公共子序列的长度 ；

优化：但是在0和1的矩阵中找最长的1对角线序列又要花去一定的时间。通过改进矩阵的生成方式和设置标记变量，可以省去这部分时间。

当字符匹配的时候，我们并不是简单的给相应元素赋上1，而是赋上其左上角元素的值加一。

我们用两个标记变量来标记矩阵中值最大的元素的位置，在矩阵生成的过程中来判断当前生成的元素的值是不是最大的，据此来改变标记变量的值，那么到矩阵完成的时候，

最长匹配子串的位置和长度就已经出来了.这样做速度比较快，但是花的空间太多。我们注意到在改进的矩阵生成方式当中，每生成当前一层，前面的那一层就已经没有用了。

因此我们只需使用一维数组 即可。可以自己写一下。

------转载结束-------

[Code]

#include<iostream> using namespace std; const int MAXN = 1000 + 10; char a[MAXN], b[MAXN]; int f[MAXN][MAXN]; int max(int x, int y){return (x > y) ? x : y;} int main() { int i, j, la, lb; while(gets(a)&& gets(b)) { la = strlen(a); lb = strlen(b); memset(f, 0, sizeof(f)); for (i = 0; i < la; i++) for (j = 0; j < lb; j++) if (a[i] == b[j]) f[i+1][j+1] = f[i][j]+1; else f[i+1][j+1] = max(f[i][j+1], f[i+1][j]); cout << f[la][lb] << endl;; } }