最长回文子序列

原文链接：https://blog.csdn.net/sugarbliss/article/details/80730929

问题描述
回文序列(Palindromic sequence, Palindrome)是指正向遍历和反向遍历完全相同的序列，例如字符串“AAAAA”显然是一个回文序列，又如字符串“ABC@CBA”也是一个回文序列。现在，我们要在一个（字符）序列中找出最长回文子序列的长度。例如字符序列"BBABCBCAB"，最长回文子序列是“BACBCAB”（可能不唯一），它的长度是7；子序列"BBBBB"和"BBABB"虽然也是回文序列，但却不是最长的，因此不合题意。

字符子串和字符子序列的区别
字符字串指的是字符串中连续的n个字符；如palindrome中，pa，alind，drome等都属于它的字串

而字符子序列指的是字符串中不一定连续但先后顺序一致的n个字符；如palindrome中，plind，lime属于它的子序列，而mod，rope则不是，因为它们与字符串的字符顺序不一致。

分析
对任意字符串，如果头和尾相同，那么它的最长回文子序列一定是去头去尾之后的部分的最长回文子序列加上头和尾。如果头和尾不同，那么它的最长回文子序列是去头的部分的最长回文子序列和去尾的部分的最长回文子序列的较长的那一个。

str[0...n-1]是给定的字符串序列，长度为n，假设f(0,n-1)表示序列str[0...n-1]的最长回文子序列的长度。

1.如果str的最后一个元素和第一个元素是相同的，则有：f(0,n-1)=f(1,n-2)+2；例如字符串序列“AABACACBA”，第一个元素和最后一个元素相同，其中f(1,n-2)表示红色部分的最长回文子序列的长度；

2.如果str的最后一个元素和第一个元素是不相同的，则有：f(0,n-1)=max(f(1,n-1),f(0,n-2))；例如字符串序列“ABACACB”，其中f(1,n-1)表示去掉第一元素的子序列，f(0,n-2)表示去掉最后一个元素的子序列。

 设字符串为s，f(i,j)表示s[i..j]的最长回文子序列。 

状态转移方程如下： 

当i>j时，f(i,j)=0。 

当i=j时，f(i,j)=1。 

当i<j并且s[i]=s[j]时，f(i,j)=f(i+1,j-1)+2。   注：f(i+1,j-1)是相对(i,j)位置的左下方

当i<j并且s[i]≠s[j]时，f(i,j)=max( f(i,j-1), f(i+1,j) )。  注：(i,j-1)是相对(i,j)位置的左侧第一个数，(i+1,j)是相对(i,j)位置的下方第一个数

因此，大多数情况下（排除i>j和i=j的情况），最终(i,j)位置的数值只与它的左下方、左侧和下方有关。

由于f(i,j)依赖i+1，所以循环计算的时候，第一维必须倒过来计算，从s.length()-1到0。 

最后，s的最长回文子序列长度为f(0, s.length()-1)。

注：本程序的填表方向斜向左上，从最后一行最后一个数向左填写。

代码

（时间复杂度O（n^2）），空间复杂度O（n^2））

/**
     * 思路：
     *  设dp[i][j]表示s[i:j]中的最长回文子序列长度，则
     *  (1) s[i] == s[j] ==> dp[i][j] = dp[i+1][j-1] + 2;
     *  (2) s[i] != s[j] ==> dp[i][j] = max(dp[i+1][j], dp[i][j-1])
     *  初始化: dp[i][i] = 1
     *  return dp[0][n-1]
     */
public int longestPalindromeSubseq(String s) {
    if(s == null || s.length() == 0) {
        return 0;
    }
    int n = s.length();
    int[][] dp = new int[n][n];
    for (int i = n - 1; i >= 0; i--) {
        dp[i][i] = 1;
        for (int j = i + 1; j < n; j++) {
            if (s.charAt(i) == s.charAt(j)) {
                dp[i][j] = dp[i + 1][j - 1] + 2;
            } else {
                dp[i][j] = Math.max(dp[i + 1][j], dp[i][j - 1]);
            }
        }
    }
    return dp[0][n - 1];
}