最长回文子序列 LPS 算法

我初次接触最长回文子序列是因为求最长回文子串，子串是需要连续的，而子序列可以不连续。求LPS用到的是区间dp，dp[i][j],在这里i,j代表从i开始到j的最长回文子序列的长度。

理解：对于每一个回文子序列，看成从短向外扩，如果使回文子序列能够扩大，那么毫无疑问，两端点就是让其长度增加的因素。所以dp[i][j] = dp[i+1][j-1]+2;如果不能使其扩大，那么就考虑当前最长子序列长度为多少，可以看到，当前状态，肯定是由两个端点的前一状态得来，即左端+1，右端不变或者左端不变，右端点-1，即
max(dp[i+1][j],dp[i][j-1]);

附上两道hdu练习：

http://acm.hdu.edu.cn/showproblem.php?pid=4745

code

#include<bits/stdc++.h>
using namespace std;
int a[1100];
int dp[1010][1010];
int main ()
{
    int n;
    while(~scanf("%d",&n)&&n)
    {
        for(int i=0;i<n;i++)
        {
            scanf("%d",&a[i]);
        }
        memset(dp,0,sizeof(dp));
        for(int i=0;i<=n;i++)
            dp[i][i] = 1;

        for(int i=n-1;i>=0;i--)
        {
            for(int j=i+1;j<n;j++)
            {
                if(a[i]==a[j]) dp[i][j] = dp[i+1][j-1]+2;
                else dp[i][j] = max(dp[i+1][j],dp[i][j-1]);
            }
        }
        int ans = 0;
        for(int i=0;i<n;i++)
        {
            ans = max(ans,dp[0][i]+dp[i+1][n-1]);
        }
        printf("%d\n",ans);
    }
    return 0;
}

http://acm.hdu.edu.cn/showproblem.php?pid=4632

code

#include<bits/stdc++.h>
using namespace std;
const int mod = 10007;
string s;
int dp[1010][1010];
int main ()
{
    int t;
    int n;
    int cas=0;
    scanf("%d",&t);
    while(t--)
    {
        cin>>s;
        memset(dp,0,sizeof(dp));
        for(int i=0;i<=s.length();i++)
            dp[i][i] = 1;
        long long ans=0;
        for(int i=s.length()-1;i>=0;i--)
        {
            for(int j=i+1;j<s.length();j++)
            {
                dp[i][j] = (dp[i][j-1]+dp[i+1][j]-dp[i+1][j-1]+mod)%mod;
                if(s[i]==s[j])
                    dp[i][j] = (dp[i][j]+dp[i+1][j-1]+1)%mod;
            }
        }
        printf("Case %d: %d\n",++cas,dp[0][s.length()-1]);
    }
    return 0;
}