NYOJ 36 最长公共子序列

http://acm.nyist.net/JudgeOnline/problem.php?pid=36

今天突然想到那天文哥给我看的最长公共子序列的问题....然后自己找了这道题来看看..
那天文哥叫我看的时候, 自己就这么一下扫过去.   没什么概念.

今天自己从头按照它的分析方法来分析了下. 

f[i][j] 表示 s1 的 i 个字符和 s2 的 j 个字符内最长的公共子序列

if(s1[i] == s2[j])   f[i][j] = f[i-1][j-1] + 1;

else f[i][j] = max(f[i-1][j], f[i][j-1]);

当不相同的时候，则考虑在 s1 的 i 个字符和 s2 的 j 个字符内最长的公共子序列.

可用  s1 的 i -1 个字符和 s2 的 j 个字符内最长的公共子序列长度 和 s1 的 i 个字符和 s2 的 j -1 个字符内最长的公共子序列长度的最大值.

#include <stdio.h>
#include <string.h>
#define maxn 1001
char s1[maxn], s2[maxn];
int ans[maxn][maxn];
int max(int a, int b)
{
    return a > b ? a : b;
}
int main()
{
    int i, j, t, len1, len2, mark;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%s%s",s1,s2);
        len1 = strlen(s1);
        len2 = strlen(s2);
        for(i = 0; i <= len1; i++)
        for(j = 0; j <= len2; j++)
        ans[i][j] = 0;
        mark = 0;
        for(i = 1; i <= len1; i++)
        {
            for(j = 1; j <= len2; j++)
            {
                if(s1[i-1]==s2[j-1])
                {
                    ans[i][j] = ans[i-1][j-1]+1;
                }
                else
                   ans[i][j] = max(ans[i-1][j], ans[i][j-1]);
                if(mark < ans[i][j])  mark = ans[i][j];   
            }
        }
        printf("%d\n",mark);
    }
    return 0;
}