(动态规划)最长回文子序列、回文子序列个数

(动态规划)最长回文子序列、回文子序列个数

主要内容:

1、什么是回文?

2、字符子串和字符子序列的区别

3、最长回文子序列的思路和代码

4、回文子序列个数的思路和代码

1、什么是回文palindrome?

回文指的是正读和反读都一样的字符串,如aba,abba等

2、字符子串和字符子序列的区别

字符字串指的是字符串中连续的n个字符;如palindrome中,pa,alind,drome等都属于它的字串

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

3、最长回文子序列

要求:

给定字符串,求它的最长回文子序列长度。回文子序列反转字符顺序后仍然与原序列相同。例如字符串abcdfcba中,最长回文子序列长度为7,abcdcba或abcfcba。

思路:

动态规划思想

对于任意字符串,如果头尾字符相同,那么字符串的最长子序列等于去掉首尾的字符串的最长子序列加上首尾;如果首尾字符不同,则最长子序列等于去掉头的字符串的最长子序列和去掉尾的字符串的最长子序列的较大者。

因此动态规划的状态转移方程为:

设字符串为str,长度为n,p[i][j]表示第i到第j个字符间的子序列的个数(i<=j),则:

状态初始条件:dp[i][i]=1 (i=0:n-1)

状态转移方程:dp[i][j]=dp[i+1][j-1] + 2  if(str[i]==str[j])

                   dp[i][j]=max(dp[i+1][j],dp[i][j-1])  if (str[i]!=str[j])

代码:

以下代码中的两层循环变量i,j的顺序可以改变,但必须满足i<=j的条件。

计算dp[i][j]时需要计算dp[i+1][*]或dp[*][j-1],因此i应该从大到小,即递减;j应该从小到大,即递增。

#include <iostream>
#include <vector>

usingnamespace std;

int longestPalindromeSubSequence1(string str){
    int n=str.length();
    vector<vector<int> > dp(n,vector<int>(n));

    for(int j=0;j<n;j++){
        dp[j][j]=1;
        for(int i=j-1;i>=0;i--){
            if(str[i]==str[j])
                dp[i][j]=dp[i+1][j-1]+2;
            else
                dp[i][j]=max(dp[i+1][j],dp[i][j-1]);
        }
    }
    return dp[0][n-1];
}

int longestPalindromeSubSequence2(string str){
    int n=str.length();
    vector<vector<int> > dp(n,vector<int>(n));

    for(int i=n-1;i>=0;i--){
        dp[i][i]=1;
        for(int j=i+1;j<n;j++){
            if(str[i]==str[j])
                dp[i][j]=dp[i+1][j-1]+2;
            else
                dp[i][j]=max(dp[i+1][j],dp[i][j-1]);
        }
    }
    return dp[0][n-1];
}

int main()
{
    string s;
    int length;
    while(cin>>s){
        length=longestPalindromeSubSequence2(s);
        cout<<length<<endl;
    }
    return0;
}

4、回文子序列个数

要求:

给定字符串,求它的回文子序列个数。回文子序列反转字符顺序后仍然与原序列相同。例如字符串aba中,回文子序列为"a", "a", "aa", "b", "aba",共5个。内容相同位置不同的子序列算不同的子序列。

思路:

动态规划思想

对于任意字符串,如果头尾字符不相等,则字符串的回文子序列个数就等于去掉头的字符串的回文子序列个数+去掉尾的字符串的回文子序列个数-去掉头尾的字符串的回文子序列个数;如果头尾字符相等,那么除了上述的子序列个数之外,还要加上首尾相等时新增的子序列个数,1+去掉头尾的字符串的回文子序列个数,1指的是加上头尾组成的回文子序列,如aa,bb等。

因此动态规划的状态转移方程为:

设字符串为str,长度为n,p[i][j]表示第i到第j个字符间的最长子序列的长度(i<=j),则:

状态初始条件:dp[i][i]=1 (i=0:n-1)

状态转移方程:dp[i][j]=dp[i+1][j] + dp[i][j-1] - dp[i+1][j-1]  if(str[i]!=str[j])

                   dp[i][j]=dp[i+1][j] + dp[i][j-1] - dp[i+1][j-1]+dp[i+1][j-1]+1=dp[i+1][j] + dp[i][j-1]+1  if (str[i]==str[j])

代码:

#include <iostream>
#include <vector>
usingnamespace std;

int NumOfPalindromeSubSequence(string str){
    int len=str.length();
    vector<vector<int> > dp(len,vector<int>(len));

    for(int j=0;j<len;j++){
        dp[j][j]=1;
        for(int i=j-1;i>=0;i--){
            dp[i][j]=dp[i+1][j]+dp[i][j-1]-dp[i+1][j-1];
            if(str[i]==str[j])
                dp[i][j]+=1+dp[i+1][j-1];
        }
    }
    return dp[0][len-1];
}

int main()
{
    string str;
    int num;
    while(cin>>str){
        num=NumOfPalindromeSubSequence(str);
        cout<<num<<endl;
    }
    return0;
}

标签:  动态规划回文子序列
id="BAIDUDSPCPRO_WRAP_pu100007" src="http://entry.baidu.com/rp/home?psid=pu100007&pswidth=auto&psheight=260&rsi0=723&rsi1=260&di=pu100007&title=(%E5%8A%A8%E6%80%81%E8%A7%84%E5%88%92)%E6%9C%80%E9%95%BF%E5%9B%9E%E6%96%87%E5%AD%90%E5%BA%8F%E5%88%97%E3%80%81%E5%9B%9E%E6%96%87%E5%AD%90%E5%BA%8F%E5%88%97%E4%B8%AA%E6%95%B0%20-%20AndyJee%20-%20%E5%8D%9A%E5%AE%A2%E5%9B%AD&mobile=true&t=1477361511667&query=%E5%9B%9E%E6%96%87%E5%AD%90%E5%BA%8F%E5%88%97&dsid=&url=http%3A%2F%2Fwww.cnblogs.com%2FAndyJee%2Fp%2F4465696.html&sType=list&tcid=1c0tfr8&ptcid=&newTrans=true&qwtj=www_normal_2_0_10_title&qwpd=undefined&isSf=false&qwsid=-1&ifr=infr%3A1_cross%3A0_drs%3A3_pcs%3A740x523_pss%3A723x5355_cfv%3A0_cpl%3A1_chi%3A1_cce%3A1_cec%3AUTF8_tlm%3A1477361511_ecd%3A1_adw%3A723x260&ltu=http%3A%2F%2Ftranscoder.baidu.com%2Ffrom%3D844b%2Fbd_page_type%3D1%2Fssid%3D0%2Fuid%3D0%2Fpu%3Dusm%25400%252Csz%25401320_2001%252Cta%2540iphone_1_9.3_3_601%2Fbaiduid%3D7FEF7AC862B6F8E6E78D9AA437555D2F%2Fw%3D0_10_%2Ft%3Diphone%2Fl%3D3%2Ftc%3Fref%3Dwww_iphone%26lid%3D18324862759200805339%26order%3D2%26fm%3Dalop%26tj%3Dwww_normal_2_0_10_title%26vit%3Dosres%26m%3D8%26srd%3D1%26cltj%3Dcloud_title%26asres%3D1%26title%3D%2528%25E5%258A%25A8%25E6%2580%2581%25E8%25A7%2584%25E5%2588%2592%2529%25E6%259C%2580%25E9%2595%25BF%25E5%259B%259E%25E6%2596%2587%25E5%25AD%2590%25E5%25BA%258F%25E5%2588%2597%25E5%259B%259E%25E6%2596%2587%25E5%25AD%2590%25E5%25BA%258F%25E5%2588%2597%25E4%25B8%25AA%25E6%2595%25B0-..._%25E5%258D%259A%25E5%25AE%25A2%25E5%259B%25AD%26dict%3D30%26w_qd%3DIlPT2AEptyoA_yk5tv9brQOxAjNSeIsnviG%26sec%3D16167%26di%3D2c05724f9e192cf6%26bdenc%3D1%26tch%3D124.0.294.265.0.0%26nsrc%3DIlPT2AEptyoA_yixCFOxXnANedT62v3IEQGG_ytK1DK6mlrte4viZQRAQz8rOZGTJoCb9n00sqcJxH_e_W5k6so4g43%26eqid%3Dfe4efd9717f2d00010000001580e2989%26wd%3D%26clk_info%3D%257B%2522srcid%2522%253A%2522159&ref=" height="260" scrolling="no" frameborder="0" allowtransparency="true" style="margin: 0px; padding: 0px; width: 1px; min-width: 100%; height: 295px;">
提示: 原网站已由百度转码,以便在移动设备上查看
  • 1
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值