Ural 1635 Mnemonics and Palindromes(DP)

题目地址:Ural 1635

又是输出路径的DP。。。连着做了好多个了。。

状态转移还是挺简单的。要先预处理出来所有的回文串,tag[i][j]为1表示字符串i--j是一个回文串,否则不是回文串。预处理时要用n^2的方法,即枚举回文串中间,可以分奇数和偶数两种分别求一次。

然后dp转移方程为,若tag[j][i]==1,dp[i]=min(dp[i],dp[j-1]+1);

对于最令人讨厌的路径输出,可以用pre来记录回文串前驱分裂点,然后根据这些分裂点来输出。

代码如下:

#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <stdlib.h>
#include <math.h>
#include <ctype.h>
#include <queue>
#include <map>
#include <set>
#include <algorithm>

using namespace std;
#define LL __int64
const int INF=0x3f3f3f3f;
char s[5000];
int dp[5000], tag[4001][4001], len, pre[4001], a[4001];
void init()
{
    int i, j, k, l, r;
    memset(tag,0,sizeof(tag));
    for(i=0;i<len;i++)
        tag[i][i]=1;
    for(i=1;i<len;i++)
    {
        for(j=1;j<=i;j++)
        {
            l=i-j;r=i+j;
            if(r>=len) break;
            if(s[l]==s[r])
            {
                tag[l][r]=1;
            }
            else break;
        }
    }
    for(i=0;i<len-1;i++)
    {
        if(s[i]==s[i+1])
        {
            tag[i][i+1]=1;
            for(j=1;j<=i;j++)
            {
                l=i-j;r=i+j+1;
                if(r>=len) break;
                if(s[l]==s[r])
                {
                    tag[l][r]=1;
                }
                else
                    break;
            }
        }
    }

}
int main()
{
    int i, j, k, cnt;
    scanf("%s",s);
    len=strlen(s);
    init();
    dp[0]=0;
    memset(pre,-1,sizeof(pre));
    for(i=1;i<=len;i++)
    {
        dp[i]=dp[i-1]+1;
        pre[i]=i-1;
        for(j=1;j<i;j++)
        {
            if(tag[j-1][i-1])
            {
                if(dp[i]>dp[j-1]+1)
                {
                    dp[i]=dp[j-1]+1;
                    pre[i]=j-1;
                }
            }
        }
    }
    printf("%d\n",dp[len]);
    cnt=0;
    for(i=len;i!=-1;i=pre[i])
    {
        a[++cnt]=i;
        //printf("%d ",a[cnt-1]);
    }
    sort(a+1,a+cnt+1);
    /*for(i=1;i<=cnt;i++)
    {
        printf("%d ",a[i]);
    }*/
    for(i=2;i<=cnt;i++)
    {
        for(j=a[i-1];j<a[i];j++)
        {
            printf("%c",s[j]);
        }
        if(i!=cnt)
        printf(" ");
        else
            puts("");
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值