HDU Advanced Fruits

Longest Common Subsequence : LCS

dfs 回溯找到公共子序列的位置

题意：

给出两个水果的名字，现在进行转基因造了一个新的品种，请你给它命名，命名方式是按照原来两个字符串中字母出现的位置，如果某一个字母两个字符串都有，就只需记录一次，尽可能的是字符串短。

思路：

很明显是为了找出最长的公共子序列，这个很简单就能找得到，但是即使找到了怎么输出呢？

$$dp[i][j]$$

只是记录路径，而无法输出。这个时候就会发现dfs的巨大魅力，因为是可以回溯！寻着dp 的路径找到那个字母是公共的。记录，然后输出即可。

#include <iostream>
#include <cstdio>
#include <cstring>

using namespace std;

const int maxn = 1005;

int dp[maxn][maxn];
int pre[maxn][maxn];
char s1[maxn],s2[maxn],s[maxn];
int same1[maxn],same2[maxn],cnt;

void LCS(int len1,int len2)
{
    memset(dp,0,sizeof(dp));
    for(int i = 1;i <= len1; i++) {
        for(int j = 1;j <= len2; j++) {
            if(s1[i-1] == s2[j-1]) {
                dp[i][j] = dp[i-1][j-1] + 1;
                pre[i][j] = 0;
            }
            else if(dp[i-1][j] > dp[i][j-1]) {
                dp[i][j] = dp[i-1][j];
                pre[i][j] = 1;
            }
            else {
                dp[i][j] = dp[i][j-1];
                pre[i][j] = 2;
            }
        }
    }
}

void dfs(int i,int j)
{
    if(i == 0 || j == 0) return ;
    else if(pre[i][j] == 1) dfs(i-1,j);
    else if(pre[i][j] == 2) dfs(i,j-1);
    else {
        dfs(i-1,j-1);
        same1[++cnt] = i - 1;
        same2[cnt] = j - 1;
    }
}

int main(int argc, char const *argv[])
{
    //freopen("in.txt","r",stdin);

    while(scanf("%s%s",s1,s2) != EOF) {
        int len1 = strlen(s1);
        int len2 = strlen(s2);
        memset(pre,0,sizeof(pre));
        LCS(len1,len2);
        cnt  = 0;
        dfs(len1,len2);
        int i = 0,j = 0;
        int ans = 0;
        for(int k = 1;k <= cnt; k++) {
            for(;i < same1[k]; i++)
                s[ans++] = s1[i];
            for(;j < same2[k]; j++)
                s[ans++] = s2[j];
            s[ans++] = s1[i];
            i++,j++;
        }
        while(s1[i] != '\0')
            s[ans++] = s1[i++];
        while(s2[j] != '\0')
            s[ans++] = s2[j++];
        s[ans] = '\0';
        printf("%s\n", s);
    }
    return 0;
}