POJ1159 回文串DP与MTE，滚动数组

题意：

输入一串字符串，插入一些字符使这个字符串变成回文串，求插入字符个数的最小值

要点：

可以这么看，设字符串为C1,C2,C3……Cn，如果C1==Cn，那么只要再考虑C2是否等于Cn-1即可，如果不相等，有下列两种选择：

在串后加C1'==C1，即C1,C2……Cn,C1'，那么只要继续求C2……Cn即可。

或在串前加Cn'==Cn,即Cn',C1,C2……Cn，那么只要求C1,C2……Cn-1即可。

接下来用记忆化搜索，但提交时说内存超出了，因为递归确实本身很容易爆内存，所以将一些变量改为short型勉强过了。

15033788

Seasonal

1159

Accepted

49300K

1438MS

C++

647B

2015-12-27 13:51:04

最后这内存差一点就爆了，肯定有更好的方法，再研究一下

代码如下：

#include<stdio.h>
#include<string.h>
#define maxn 5000
char a[maxn];
int n;
short dp[maxn][maxn];//会超内存所以用short

short solve(int left, int right)
{
	if (dp[left][right] >= 0)
		return dp[left][right];//判断是否被计算过,计算过下面就不用继续下去了
	if (left >= right)//临界条件，如果中间只有1个是大于，2个是等于
		return dp[left][right] = 0;
	if (a[left] == a[right])
		return dp[left][right] = solve(left + 1, right - 1);
	else
	{
		short x = solve(left + 1, right);//在后面加的情况
		short y = solve(left, right - 1);//在前面加的情况
		return dp[left][right] = (x < y )? (x + 1) : (y + 1);
	}
}

int main()
{
	scanf("%d", &n);
	scanf("%s", a);
	memset(dp, -1, sizeof(dp));
	printf("%d\n", solve(0,n-1));
	return 0;
}

网上大神还有更好的解法

首先，直接输入字符串和它的逆序串求最长公共子序列长度len。n-len即为所求。(n为原串长度)，我验证了一下这是正确的，比如AB3AB，它的逆序串是BA3BA，公共子序列长度是3,字符串长度5-3=2，是正确的，最后可以变成BAB3BAB。证明是这样的，比如AB3AB这个序列，它可以有最大的回文串A3A或B3B，不需要连续。这个回文串就算整个数列逆序也不会改变，因为是最大的，所以可以理解输入字符串和它的逆序串的最长公共子序列就是这个最大的回文串，比如AB3AB对应的最大公共子序列也是A3A或B3B，那么长度自然相等。

再者，为了减少空间复杂度，主要是减少数组的区间，有了一个新概念——滚动数组

这个博客讲的很好：点击打开链接

对于现在这个问题，d[i][j]只依赖于d[i-1][j],d[i][j-1]，所以dp数组可以直接开成dp[2][maxn]，一下子从5010*5010减少到了2*5010，减少了很多内存，滚动数组其实就是循环利用一些空间，从而达到优化空间的目的。

注意滚动数组打印方案十分麻烦，只能输出最后一个值，当然这里直接dp[n%2][j]输出没问题，网上一般都定义一个MAX比较一下，可能怕出错，其实我觉得无所谓。

15034289

Seasonal

1159

Accepted

208K

813MS

C++

633B

2015-12-27 17:17:04

现在无论时间还是空间都好多了，但好像定义一个max稍微快了20几ms，可能从dp所在的栈读取反而慢？

代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn = 5005;

int main()
{
	int n,i,j;
	char a[maxn],b[maxn];
	int dp[2][maxn];//滚动数组
	while (scanf("%d", &n) != EOF)
	{
		scanf("%s", a + 1);//注意这里用%s，否则会读入\n,这种字符串输入以后最好都用%s
		memset(dp, 0, sizeof(dp));
		for (i = 1; i <= n; i++)
			b[i] = a[n +1- i];//取逆序数组
		for (i = 1; i <= n; i++)
			for (j = 1; j <= n; j++)//标准LCS
			{
				if (a[i] == b[j])
					dp[i % 2][j] = dp[(i - 1) % 2][j - 1] + 1;//滚动数组应用，其实只要改牵涉到dp的第一个下标就行
				else
					dp[i % 2][j] = max(dp[i % 2][j - 1], dp[(i - 1) % 2][j]);		
			}
		printf("%d\n", n - dp[n%2][n]);//输出最后一个值
	}
	return 0;
}