poj1159

题目不难，主要是数据大，这样就必须格外注意时间复杂度和空间复杂度。即使是用dp解题，不同的dp思想可能产生不同的时间复杂度和空间复杂度。这题就是一个典型例子。题意要求给定一个字符串，删除最少的字符，使之成为回文串。要求这个最小删除个数。这里给出两种不同的dp思想：

1）令dp[i][j]表示使第i个字符到第j个字符组成的字符串成为回文串必须删除的最小字符个数。则分析如下：

（1）若将第j个字符删除，那么dp[i][j]=dp[i]][j-1]+1;

（2）若不删除，则说明第j个字符为结尾字符，那么我们寻找首字符，可以从第i个字符开始查找，顺序扫描直到第j-1个字符，若扫描过程中发现存在扫描字符与尾字符相等，则将该字符最为首字符，去掉该字符前的字符（第i个字符至该字符之间的所有字符），并且最少删除个数还应加上首字符与结尾字符之间的字符组成的字符串要构成回文串所必须删除的最小字符个数。即dp[i][j]=min{k-i+dp[k+1][j-1]}

这样的算法朴素执行（dp没有优化）需要O(n^3),铁定超时（从这里也可看出即使是动态规划，也有优略之分）

而如果加入set数组，则时间复杂度有所降低但是仍然为O(n*log(n!))，时间复杂度依旧很高

下面是代码：（TLE）

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#define Max 5010
#define Min(a,b) (a)<(b)?(a):(b)
char record[Max];
char help[Max];
int dp[Max][Max];
int set[Max]; 
int n;
int main(){
	scanf("%d",&n);
	getchar();
	scanf("%s",help);
	int i,j,k;
	for(i=1;i<=n;i++){
		record[i]=help[i-1];
		dp[i][i]=0;
	}
	memset(set,0,sizeof(set));
	for(i=1;i<=n;i++) // 求解每个字符之前与其相等的最近的字符下标，保存在数组set中，优化dp
		for(j=i-1;j>=1;j--)
			if(record[j]==record[i]){
				set[i]=j;
				break;
			}
	for(i=2;i<=n;i++){
		for(j=1;j<=n-i+1;j++){
			dp[j][j+i-1]=dp[j][j+i-2]+1;
			/*for(k=j;k<j+i-2;k++){ // 线性扫描
				if(record[k]==record[j+i-1])
					dp[j][j+i-1]=Min(dp[j][j+i-1],k-j+dp[k+1][j+i-2]);
			}*/
			k=set[j+i-1];
			while(j<=k && k<j+i-2){ // 优化查找过程，将线性扫描改为只扫描可行位置，利用set数组
				dp[j][j+i-1]=Min(dp[j][j+i-1],k-j+dp[k+1][j+i-2]);
				k=set[k];
			}
			if(record[j+i-2]==record[j+i-1])
			dp[j][j+i-1]=Min(dp[j][j+i-1],i-2);
		}
	}
	printf("%d\n",dp[1][n]); // 输出结果
	return 0;
}

2）同样：令dp[i][j]表示使第i个字符到第j个字符组成的字符串成为回文串必须删除的最小字符个数。则分析如下：

（1）若第i个字符与第j个字符相等，则dp[i][j]=dp[i+1][j-1];

（2）若不相等，则dp[i][j]=min{dp[i][j-1]+1,dp[i+1][j]+1}

这样的算法思想要简单一些，也更容易实现，而且时间复杂度为0(n^2),但是若将dp设置成int则会MLE，必须设置为short,

当然若采用其它dp思想（好像还可以用lcs解决，则可以用到滚动数组，空间复杂度直线下降。只需要两个数组就可以了）

下面为代码：40660K+1407MS

#include <stdio.h>
#include <stdlib.h>
#define Max 5010
#define Min(a,b) (a)<(b)?(a):(b)
char record[Max];
char help[Max];
short dp[Max][Max]; // 若为int,则MLE
int n;
int main(){
	scanf("%d",&n);
	getchar();
	scanf("%s",help);
	int i,j;
	for(i=1;i<=n;i++){
		record[i]=help[i-1];
		dp[i][i]=0;  
	}
	for(i=1;i<=n-1;i++) // 初始化长度为2的dp
		if(record[i]==record[i+1])
			dp[i][i+1]=0;
		else
			dp[i][i+1]=1;
	for(i=3;i<=n;i++) // 求解dp
		for(j=1;j<=n-i+1;j++) 
			if(record[j]==record[j+i-1])
				dp[j][j+i-1]=dp[j+1][j+i-2];
			else
				dp[j][j+i-1]=Min(dp[j][j+i-2]+1,dp[j+1][j+i-1]+1);
	printf("%d\n",dp[1][n]);
	return 0;
}