最长公共子串（LCS）

1、先科普下最长公共子序列 & 最长公共子串的区别：

找两个字符串的最长公共子串，这个子串要求在原字符串中是连续的。而最长公共子序列则并不要求连续。
（最长公共子序列详细解释）

2、最长公共子串

其实这是一个序贯决策问题，可以用动态规划来求解。我们采用一个二维矩阵来记录中间的结果。这个二维矩阵怎么构造呢？直接举个例子吧：”bab”和”caba”(当然我们现在一眼就可以看出来最长公共子串是”ba”或”ab”)
　　 b　　a　　b
c　　0　　0　　0
a　　0　　1　　0
b　　1　　0　　1
a　　0　　1　　0
我们看矩阵的斜对角线最长的那个就能找出最长公共子串。
不过在二维矩阵上找最长的由1组成的斜对角线也是件麻烦费时的事，下面改进：当要在矩阵是填1时让它等于其左上角元素加1。
　　 b　　a　　b
c　　0　　0　　0
a　　0　　1　　0
b　　1　　0　　2
a　　0　　2　　0
这样矩阵中的最大元素就是 最长公共子串的长度。

动态转移方程为：
如果xi == yj， 则 c[i][j] = c[i-1][j-1]+1
如果xi ! = yj, 那么c[i][j] = 0
最后求Longest Common Substring的长度等于
max{ c[i][j], 1<=i<=n， 1<=j<=m}

代码：

#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
char s1[1010],s2[1010];
int maxlen[1010][1010];
int main()
{
    while(cin>>s1>>s2)
    {
        int length1=strlen(s1);
        int length2=strlen(s2);
        memset(maxlen,0,sizeof(maxlen));
        int maxl=0,x,y;
        for(int i=1;i<=length1;i++)
        {
            for(int j=1;j<=length2;j++)
            {
                if(s1[i-1]==s2[j-1])maxlen[i][j]=maxlen[i-1][j-1]+1;
                else maxlen[i][j]=0;
                if(maxl<maxlen[i][j])
                {
                    maxl=maxlen[i][j];
                    x=i;y=j;
                }
            }
        }
        char s[1000];//输出最长公共子串
        memset(s,0,sizeof(s2));
        int k=maxl;
        s[k--]='\0';
        for(int i=x-1,j=y-1;i>=0&&j>=0;i--,j--)
        {
            if(s1[i]==s2[j])s[k--]=s1[i];
            else break;
        }
        cout<<maxl<<endl;
        cout<<s<<endl;
        memset(s1,0,sizeof(s1));
        memset(s2,0,sizeof(s2));
    }
    return 0;
}