[ACM_动态规划] Palindrome

http://acm.hust.edu.cn/vjudge/contest/view.action?cid=28415#problem/D

 

题目大意：给一个长为n的字符串，问最少插入几个字符成回文串

解题思路：总长-最长公共（原来的和其倒过来的串）子序列（LCS） 

知识详解——LCS:给出两个子序列A,B，求长度最大的公共子序列（如152687和2356984——>568或268）；不妨设d(i,j)为A,B的LCS,   则最有子结构为::A[i]=B[j]时，d(i,j)=d(i-1,j-1)+1;否则，d(i,j)=max{d(i-1,j),d(i,j-1)};复杂度为O（m*n),也可用滚动数组法优化。

#include<iostream>
#include<cmath>
#include<string>
#include<algorithm>
using namespace std;
int c[5005][5005]; 
int main(){
    int n; 
    for(string a,b;cin>>n>>a;){  
       
       b=a;                      //输入a并获得其颠到序列+对c[][]初始化0
       reverse(b.begin(),b.end());
       memset(c,0,sizeof(c));

       for(int j=0;j<n;j++){     //利用最优子结构求解过程
          for(int i=0; i<n; i++){  
               if(a[i]==b[j])  c[i+1][j+1]=c[i][j]+1;  
               else   c[i+1][j+1] = c[i][j+1]>c[i+1][j] ? c[i][j+1]:c[i+1][j];  
           }
       }  

       cout<<n-c[n][n]<<endl;
    }
    return 0;
}

 

 

注意： 上面的代码虽然实现了求解，但是 由于开了一个c[2005][2005]的数组所以不幸“Memory Limit Exceeded”，通过观察不难发
现
其最优子结构中每次计算只用到两层数据，因此可以考录建立一个c[2][2005]的动态表来节省空间，代码如下：

#include<iostream>
#include<cmath>
#include<string>
#include<algorithm>
using namespace std;
int c[2][5005];//开一个2层的数组，0层存先前数值，1层存当前数值，每次计算完一层要用1层更新0层  
int main(){
    int n; 

    for(string a,b;cin>>n>>a;){  
       
       b=a;
       reverse(b.begin(),b.end());
       memset(c,0,sizeof(c));

       for(int j=0;j<n;j++){ 

          for(int i=0; i<n; i++){                  //用前一层计算下一层  
               if(a[i]==b[j])  c[1][i+1]=c[0][i]+1;  
               else   c[1][i+1] = c[1][i]>c[0][i+1] ? c[1][i]:c[0][i+1];  
          }  
          for(int i=1; i<=n; i++) c[0][i]=c[1][i]; //每次计算完一层要用1层更新0层  
       }  

       cout<<n-c[1][n]<<endl;
    }
    return 0;
}