UVa 10405 & POJ 1458 Longest Common Subsequence

　　求最长公共子序列LCS，用动态规划求解。

　　UVa的字符串可能含有空格，开始用scanf("%s", s);就WA了一次...那就用gets吧，怪不得要一行放一个字符串呢。

　　（本来想用fgets的，可是又放弃了，形式麻烦、代码长是一小方面，另一方面fgets把' '字符也读入，还要做额外的处理...虽然gets有传说中的缓冲区溢出漏洞，不过多加注意一下就好啦，个人认为代码还没大到要用那些工程性的东西的时候）

#include <cstdio>
#include <cstring>
#include <algorithm>
#define MAXN 1000+10
using namespace std;

char s1[MAXN], s2[MAXN];
int c[MAXN][MAXN];

int main()
{
#ifdef LOCAL
    freopen("in", "r", stdin);
#endif
    while (gets(s1) && gets(s2))
    {
        int m = strlen(s1);
        int n = strlen(s2);
        memset(c, 0, sizeof(c));
        for (int i = 1; i <= m; i++)
            for (int j = 1; j <= n; j++)
            {
                if (s1[i-1] == s2[j-1])   c[i][j] = c[i-1][j-1] + 1;
                else c[i][j] = max(c[i][j-1], c[i-1][j]);
            }
        printf("%d
", c[m][n]);
    }
    return 0;
}

　　在POJ中字符串用空格分割，字符串中不含空格，并且两个字符串在同一行上，用scanf("%s", s)替换get(s)就可以了

　　下面是打印LCS的代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#define MAXN 1000+10
using namespace std;

char s1[MAXN], s2[MAXN];
int c[MAXN][MAXN];

void print_LCS(int i, int j)
{
    if (i == 0 || j == 0)   return;
    if (s1[i-1] == s2[j-1])
    {
        print_LCS(i-1, j-1);
        printf("%c", s1[i-1]);
    }
    else if (c[i-1][j] == c[i][j])   print_LCS(i-1, j);
    else if (c[i][j-1] == c[i][j])   print_LCS(i, j-1);
}

int main()
{
#ifdef LOCAL
    freopen("in", "r", stdin);
#endif
    while (scanf("%s%s", s1, s2) != EOF)
    {
        int m = strlen(s1);
        int n = strlen(s2);
        memset(c, 0, sizeof(c));
        for (int i = 1; i <= m; i++)
            for (int j = 1; j <= n; j++)
            {
                if (s1[i-1] == s2[j-1])   c[i][j] = c[i-1][j-1] + 1;
                else  c[i][j] = max(c[i][j-1], c[i-1][j]);
            }
        printf("%d
", c[m][n]);
        // print the LCS
        print_LCS(m, n);
        printf("
");
    }
    return 0;
}