简单LCS HDU_1503

学了一下最长公共子串，它是属于dp里面的 dp=max{(i,j-1),(i-1,j),(i-1,j-1)+d}问题，不得不说，规划方向确实厉害，当然这只适用于两个字符串匹配的问题，n个字符串的话，我查了一下，要用后缀数组，这个以后再来弄吧。

比较难的是记录公共字串，在DP的时候记录好搜索的方向，就可以询步骤得到公共串。尤其1503这题呢，除了把公共部分打印出来，还要把非公共部分打印出来。

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
char a[] , b[];
int dp[][];
int rec[][];
int dir[][];
int LCS(int r,int c)
{
    if (dp[r][c]>=) return dp[r][c]; //用DP+记忆化搜索写的LCS
    if (a[r-]==b[c-]){
        dp[r][c]=LCS(r-,c-)+;
        dir[r][c]=;
    }
    else
    {
        int dp1=LCS(r-,c);
        int dp2=LCS(r,c-);
        if (dp1>dp2){
            dp[r][c]=dp1;
            dir[r][c]=; //用dir记录规划的方向
        }
        else
        {
            dp[r][c]=dp2;
            dir[r][c]=;
        }
    }
    return dp[r][c];
}
void output(int r,int c)//已经用dir记录了规划方向，这里用递归将字符串还原出来
{
    if (r==&&c==) return; //一开始这里犯糊涂了，写成或，要知道只有一个为0的话，还要把其他非公共部分打印出来，还得继续下去。
   // cout<<"dir"<<" "<<r<<" "<<c<<" "<<dir[r][c]<<endl;
    if (dir[r][c]==){
        output(r-,c-);
        printf("%c",a[r-]);
    }
    if (dir[r][c]==){ //此时a[r-1]相当于被抛弃了，所以可以打印了，下面同理
        output(r-,c);
        printf("%c",a[r-]);
    }
    if (dir[r][c]==){
        output(r,c-);
        printf("%c",b[c-]);
    }
}
int main()
{
    while (scanf("%s%s",a,b)!=EOF)
    {
        int len1=strlen(a);
        int len2=strlen(b);
        int i,j;
        for (i=;i<=len1;i++){
            for (j=;j<=len2;j++){
                dp[i][j]=-;
                if (i==||j==) dp[i][j]=;
                if (i==) dir[i][j]=;//为了让output把所有非公共部分也打印完，得预先把这些设置好
                if (j==) dir[i][j]=;
            }
        }
        int ans=LCS(len1,len2);
        output(len1,len2);
        //printf("%d",ans);
        putchar('\n');
    }
}