LCS最长公共子序列

Description
两个字符串， 求两个字符串最长公共子序列的长度
字符串的长度不超过5000

Sample Input Copy
ABCBDAB
BACBBD
Sample Output Copy
4
公共子串为BCBD

两个字符串a和b。对比每一个字符串。
如果a[i] == b[j] ，dp[i][j] == 左上角+1
否者 dp[i][j] == 左边或者上边最大的

#include<bits/stdc++.h>
using namespace std;
char a[5009];
char b[5009];
int dp[5009][5009];
int main()
{
     scanf("%s",&a);
     scanf("%s",&b);
     int L = strlen(a);
     int R = strlen(b);
     dp[0][0] = 0;
     for(int i=1;i<=R;i++)
     for(int j=1;j<=L;j++){
        if(a[j-1]==b[i-1])
        dp[i][j] = dp[i-1][j-1]+1;
        else
        dp[i][j] = max(dp[i-1][j],dp[i][j-1]);//左边与上面
     }
         cout<<dp[R][L]<<endl;

	return 0;
}

若找到公共子串则需要记录一下方向，
不过一般会有多个比如样例中：
ABCBDAB
BACBBD
包含BCBB，BCDB，ABBD。
先输入a串，或者先输入b串公共子串会不同，当dp[i-1] [j] == dp[i] [j-1]时候，记录为左方向还是上方向也得到的公共子串也不同。图片是其中一种。
为1时向左上走，为2时向上走，为3时向左走。等于1时为子串其中一个字母。

#include<bits/stdc++.h>
using namespace std;
char a[5009];
char b[5009];
int dp[5009][5009];
int vis[5009][5009];

void LCS(int i,int j)
{
	if(vis[i][j] == 0)
		return;
	if(vis[i][j]==1){
		LCS(i-1,j-1);
		cout<<b[i-1]<<" "; //a[j-1]
	}
	else if(vis[i][j]==2)
		LCS(i-1,j);
	else
		LCS(i,j-1);
}

int main()
{
     memset(dp,0,sizeof dp);
     memset(vis,0,sizeof vis);
     scanf("%s",&a);
     scanf("%s",&b);
     int L = strlen(a);
     int R = strlen(b);
     dp[0][0] = 0;
     for(int i=1;i<=R;i++)
     for(int j=1;j<=L;j++){
        if(a[j-1]==b[i-1]){
        dp[i][j] = dp[i-1][j-1]+1;
        vis[i][j] = 1;
        }
        else if(dp[i-1][j]>dp[i][j-1]){
        dp[i][j] = dp[i-1][j];
        vis[i][j] = 2;
        }
        else {
        dp[i][j] = dp[i][j-1];
        vis[i][j] = 3;
        }
     }
    cout<<dp[R][L]<<endl;
    LCS(R,L);

	return 0;
}