POJ 1159 Palindrome 之解法三 LCS (最长公共子序列)

题目链接

题意：给定一个长度为N（3 <= N <= 5000）的字符串，只用插入字符这种操作，问要使其变成回文字符串，最少需多少步的插入？

思路：找出所给字符串中最大的回文子序列长度L，N - L即为答案。L = 所给字符串和其逆序串的LCS。

LCS 背景：

LCS(Longest Common Subsequences)最长公共子序列用一般的动态规划时间复杂度O(N^2), 但经过优化可以达到O(NlogN),

其中O(NlogN)的算法是最长公共子序列向最长递增子序列退化：

　　设有序列A，B。记序列A中各个元素在B 中的位子(降序排列)，然后按在A中的位置依次列出按后求A的最长递增子序列。

　　例如：有A={a,b,a,c,x}，B={b,a,a,b,c,a}则有a={6,3,2},b={4,1},c={5};x=/;（注意降序排列）

然后按A中次序排出{a(6,3,2),b(4,1),a(6,3,2),c(5),x()}={6,3,2,4,1,6,3,2,5}；对此序列求最长递增子序列即可

#include<stdio.h>
#include <stdlib.h>
#include <iostream>
using namespace std;
#define MAX 5000   
#define max(x,y) ((x) > (y) ? (x) : (y))
short int dp[2][MAX + 5];      
char arr[MAX + 5]; 
int main()
{
	int  i, j, n, x, y;
	scanf("%d", &n);
    scanf("%s", arr + 1);
	滚动数组求LCS
	for( i = 1; i <= n; i ++ ){
		x = i & 1;
		y = ( i + 1 ) & 1;
		for( j = 1; j <= n; j ++ ){
			if( arr[i] == arr[n + 1 - j] )
				dp[x][j] = 1 + dp[y][j - 1];
			else dp[x][j] = max( dp[y][j], dp[x][j - 1] );
		}
	}
	printf( "%d\n", n - dp[x][n] );
	return 0;
}