动态规划（DP）中LCS问题合集（有待更新中~

     现在遇到的LCS问题，大概分为下列：

一、求朴素最长公共子序列问题

（推荐题目：hdu 1159 Common Subsequence：http://acm.hdu.edu.cn/showproblem.php?pid=1159 ）

参考代码：

//hdu 1159
#include <iostream>
#include <cstring>
#include <cstdio>
#include <string>
#include <algorithm>
#define MAXN 1010
using namespace std;
char a[MAXN];
char b[MAXN];
int dp[MAXN][MAXN];
int main()
{
    while(~scanf("%s%s",a,b))
    {
        memset(dp,0,sizeof(dp));
        int len1=strlen(a);
        int len2=strlen(b);
        for(int i=1;i<=len1;i++)
        {
            for(int j=1;j<=len2;j++)
            {
                if(a[i-1]==b[j-1])
                    dp[i][j]=dp[i-1][j-1]+1;
                else
                    dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
            }
        }
        printf("%d\n",dp[len1][len2]);
    }
    return 0;
}

二、求最长公共子序列（LCS）输出任意一个符合的子串

（还没找到题，不过推荐个变型题：hdu Advanced Fruits：http://acm.hdu.edu.cn/showproblem.php?pid=1503）

//hdu 1503
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#define MAXN 1010
using namespace std;
char str[MAXN],s[MAXN];

int vis[MAXN][MAXN];
int dp[MAXN][MAXN];
void Lcs(int len1,int len2)
{

    memset(dp,0,sizeof(dp));
    for(int i = 0; i<=len1; i++)
        vis[i][0] = 1;
    for(int i = 0; i<=len2; i++)
        vis[0][i] = -1;
    for(int i=1; i<=len1; i++)
    {

        for(int j = 1; j<=len2; j++)
        {
            if(str[i-1]==s[j-1])
            {
                dp[i][j] = dp[i-1][j-1]+1;
                vis[i][j] = 0;
            }
            else if(dp[i-1][j]>=dp[i][j-1])
            {
                dp[i][j] = dp[i-1][j];
                vis[i][j] = 1;
            }
            else
            {
                dp[i][j] = dp[i][j-1];
                vis[i][j] = -1;
            }
        }
    }
}

void PrintLCS(int i,int j)
{
    if(i==0&&j==0)
        return ;
    else if(vis[i][j]==0)
    {
        PrintLCS(i-1,j-1);
        printf("%c",str[i-1]);
    }
    else if(vis[i][j]==1)
    {
        PrintLCS(i-1,j);
        printf("%c",str[i-1]);//将此条与下面61行注释掉就是输出任意解了
    }
    else
    {
        PrintLCS(i,j-1);
        printf("%c",s[j-1]);  //这两条是输出s1和s2的非公共字符
    }
}


int main()
{
    while(~scanf("%s%s",str,s))
    {
        int len1=strlen(str);
        int len2=strlen(s);
        Lcs(len1,len2);
        PrintLCS(len1,len2);
        printf("\n");
    }


    return 0;
}

三、求公共子序列，并且输出全部符合要求的公共子序列

（推荐题目：poj 1934 Trip：http://poj.org/problem?id=1934）

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <set>
#define MAXN 1010
using namespace std;

char s1[MAXN],s2[MAXN],tmp[MAXN];
int len1,len2;
int dp[MAXN][MAXN];
int last1[MAXN][30],last2[MAXN][30];
set<string> ans;
void LCS()
{
    memset(dp,0,sizeof(dp));
    int i,j;
    for(i=1;s1[i]!='\0';i++)
    {
        for(j=1;s2[j]!='\0';j++)
        {
            if(s1[i]==s2[j])
                dp[i][j]=dp[i-1][j-1]+1;
            else
                dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
        }
    }
    len1=i-1;
    len2=j-1;
}
//将字母表中字母按是否在s1,s2中出现，若出现，将其最高位存贮在辅助数组last1和last2中
void buildLast(char s[],int len,int last[][30])
{
   for(int i=0;i<26;i++)
    {
        char c='a'+i;
        for(int j=1;j<=len;j++)
        {
            int k;
            for(k=j;k>=1;k--)
            {
                if(c==s[k])
                {
                    last[j][i]=k;
                    break;
                }
            }
            if(k==0)
                last[j][i] = 0;
        }
    }
}
void SaveLCS(int p1,int p2,int len)
{
    if(len<=0)
    {
        ans.insert(&tmp[1]);
        return;
    }
    if(p1>0&&p2>0)
    {
        for(int i=0;i<26;i++)
        {
            int st1=last1[p1][i];   //同样是按字母查询，当dp[st1][st2] == len符合是，加入该字母。
            int st2=last2[p2][i];
            if(dp[st1][st2] == len)
            {
                tmp[len] = 'a' + i;
                SaveLCS(st1-1, st2-1, len- 1);
            }
        }
    }
}


int main()
{
    scanf("%s%s",s1+1,s2+1);
    LCS();
    buildLast(s1,len1,last1);
    buildLast(s2,len2,last2);
    
    int cnt=dp[len1][len2];
    tmp[cnt+1] = '\0';
    SaveLCS(len1,len2,cnt);
    set<string>::iterator p=ans.begin();
    for( ;p!= ans.end();p++)
        printf("%s\n", p->c_str());

    return 0;
}

四、通过LCS来求回文串（同一串反转后求LCS）

（推荐题目：hdu 1513 Palindrome http://acm.hdu.edu.cn/showproblem.php?pid=1513

此题字符串长度有5000，开这么大的二维数组会爆内存的，推荐使用滚动数组的dp[2][MAXN]的写法。

#include <iostream>
#include <cstring>
#include <cstdio>
#define MAXN 5010
using namespace std;
char s1[MAXN],s2[MAXN];
int dp[2][MAXN],len;

void LCS()
{
    memset(dp,0,sizeof(dp));
    for(int i=1;i<=len;i++)
    {
        for(int j=1;j<=len;j++)
        {
            int x = i%2;
            int y = 1-x;
            if(s1[i-1]==s2[j-1])
                dp[x][j] = dp[y][j-1]+1;
            else
                dp[x][j] = max(dp[y][j],dp[x][j-1]);
        }
    }
}
int main()
{

    while(~scanf("%d",&len))
    {
        scanf("%s",s1);
        for(int i=0,j=len-1;i<len;i++,j--)
            s2[i]=s1[j];
        s2[len]='\0';
        LCS();
        printf("%d\n",len-dp[len%2][len]);
    }


    return 0;
}

四、字符文章中求最多单词的LCS（和求字符其实是一样的）

（题目：poj 2250 Compromise : http://poj.org/problem?id=2250）

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

char s1[35][105],s2[35][105];
char ans[35][105];
int dp[105][105],vis[105][105];
int len1,len2,cnt;

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

void PrintLCS(int i,int j)
{
    if(i==0 && j==0)
        return;
    if(vis[i][j]==0)
    {
        PrintLCS(i-1,j-1);
        strcpy(ans[cnt++],s1[i-1]);
    }
    else if(vis[i][j]==1)
        PrintLCS(i-1,j);
    else
        PrintLCS(i,j-1);
}

int main()
{
    int i;
    while(~scanf("%s",s1[0]))
    {
        len1=1,len2=1;
        while(strcmp(s1[len1-1],"#"))
            scanf("%s",s1[len1++]);
        len1--;
        scanf("%s",s2[0]);
        while(strcmp(s2[len2-1],"#"))
            scanf("%s",s2[len2++]);

        LCS();

        cnt = 0;
        PrintLCS(len1,len2);
        for(i=0;i<cnt;i++)
        {
            printf("%s",ans[i]);
            if(i!=cnt-1) printf(" ");
            else printf("\n");
        }
    }
    return 0;
}