阿里笔试之最长公共子串

题目描述：给定一个query和一个text，均由小写字母组成。要求在text中找出以同样的顺序连续出现在query中的最长连续字母序列的长度。例如，query为“acbac”，text为“acaccbabb”，那么text中的“cba”为最长的连续出现在query中的字母序列，因此，返回结果应该为其长度3。请注意程序效率。

思路：用一个矩阵来记录两个字符串中所有位置的两个字符之间的匹配情况，若是匹配则为1，否则为0。然后求出对角线最长的1序列，其对应的位置就是最长匹配子串的位置.

当字符匹配的时候，不是简单的给相应元素赋上1，而是赋上其左上角元素的值加1。我们用两个标记变量来标记矩阵中值最大的元素的位置，在矩阵生成的过程中来判断当前生成的元素的值是不是最大的，据此来改变标记变量的值，那么到矩阵完成的时候，最长匹配子串的位置和长度就已经出来了。

实例如下：

     a    c    b    a    c

a   1    0    0    1    0   

c   0    2    0    0    2

a   1    0    0    1    0

c   0    2    0    0    2

c   0    1    0    0    1

b   0    0    2    0    0

a   1    0    0    3    0    

b   0    0    1    0    0

b   0    0    1    0    0

（红色为最长公共子串）

#include <cstring>
#include <cstdio>
#define M 1010
 
int LCS(char query[], char text[])
{
    int len_query=strlen(query),len_text=strlen(text);
 
    //数组c记录匹配情况,模拟二维矩阵
    char c[len_text];
    int len, i, j;
    len=0;
 
    for(i=0; i<len_query; i++)
    {
        //不反过来会把之前数组元素冲掉的，因为后面的元素需要根据前面的元素计算
        for(j=len_text-1; j>=0; j--)
        {
            if(query[i] == text[j])
            {
                if(i==0 || j==0)
                    c[j]=1;
                else
                    c[j]=c[j-1]+1; 
            }
            else
                c[j] = 0;
 
            if(c[j] > len)
                len=c[j];
        }
    }
    return len;
}
 
int main()
{
    char str1[M],str2[M];
    printf("请输入字符串query:");
    scanf("%s", str1);
    printf("请输入字符串text:");
    scanf("%s", str2);
    printf("所求长度为:");
 
    printf("%d\n",LCS(str1,str2));
    return 0;
}

以上的代码是求出了最长公共子序列的长度，而没有找出最长公共子序列是什么。

下面贴上求出最长公共子序列的代码，供参考

#include <iostream>
using namespace std;
int c[20][20];
int Get_LCS(const char *X, const char * Y,int &len)
{	
	int i, j;
	if ( NULL == X || NULL == Y)
		return 0;
	int m = strlen(X);
	int n = strlen(Y);
	for (i=1;i<=m;++i)
		c[i][0] = 0;
	for(i=1;i<=n;++i)
		c[0][i] = 0;
	c[0][0] = 0;
	for (i=1;i<=m;++i)
	{
		for (j = 1; j<=n;++j)
		{
			if (X[i-1] == Y[j-1])
				c[i][j] = c[i-1][j-1] +1;
			else
				//c[i][j] = c[i][j-1] >= c[i-1][j] ? c[i][j-1]:c[i-1][j];
				c[i][j] = max(c[i][j-1],c[i-1][j]);
		}
	}
	for (i=0;i<=m;++i)
	{
		for (j=0;j<=n;++j)
		{
			cout<<c[i][j]<<"   ";
		}
		cout<<endl;
	}
	len = c[m][n];	
	char result[20];
	int count=0;
	for (i = m, j=n; i>=1 && j>=1;)
	{
		if (X[i-1] == Y[j-1])
		{
			result[count] = X[i-1];//将结果保存起来，待会儿倒着打印即可
			++count;
			--i;
			--j;
		}
		else
		{
			if (c[i][j-1] >c[i-1][j])//若是>= 则输出BDAB,若是> ，则输出BCBA
				--j;
			else
				--i;
		}
	}
	for(i=count-1;i>=0;--i)
		cout<<result[i]<<" ";
	cout<<endl;
	return 0;
}