NOIP模拟回文子串【动态规划】

最新推荐文章于 2024-03-08 22:54:14 发布

Neo__Z

最新推荐文章于 2024-03-08 22:54:14 发布

阅读量323

点赞数

分类专栏： --------动态规划-------- 文章标签：动态规划 NOIP模拟

本文链接：https://blog.csdn.net/cdsszjj/article/details/78428045

版权

--------动态规划-------- 专栏收录该内容

114 篇文章

订阅专栏

本文介绍了一种求解两个字符串合并后的最长回文子串的算法，通过四维动态规划实现，适用于字符串长度不超过50的情况。文章提供了完整的代码实现。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目大意：

给定两个长度不超过50的大写字母字符串s，t，求他们的所有并字符串中的最长回文子串长度。
如：“CLEVER”和“JAVA”的一个并字符串为“CLJEAVVAER”，其最长回文子串为“EAVVAE”,长度为6。

解题思路：

考试时只想到四维dp就不知道怎么做了……蒟蒻。

设想分别从两边开始取到中间合并。
设dp[l][r][L][R]表示表示s串从左往右取到s[l]，从右往左取到s[r]，t串从左往右取到t[L]，从右往左取到t[R]时取得的最长回文子串的长度，则：
if(s[l]==s[r])dp[l][r][L][R]=dp[l-1][r+1][L][R]+2;
if(s[l]==t[R])dp[l][r][L][R]=dp[l-1][r][L][R+1]+2;
if(t[L]==t[R])dp[l][r][L][R]=dp[l][r][L-1][R+1]+2;
if(t[L]==s[r])dp[l][r][L][R]=dp[l][r+1][L-1][R]+2;
四个值取最大值。
这是回文长度为偶数的情况，奇数的话只用最后在中间随便插入一个字符即可。
注意两个字符串都可能有从一边不取的情况，所以y要从0枚举到sn（tn）+1，而初始时可将s[0],s[sn+1],t[0],t[tn+1]赋为不同的字符避免错误。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#define ll long long
using namespace std;

int getint()
{
    int i=0,f=1;char c;
    for(c=getchar();(c<'0'||c>'9')&&c!='-';c=getchar());
    if(c=='-')f=-1,c=getchar();
    for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
    return i*f;
}

const int N=55;
char s[N],t[N];
int sn,tn,ans,dp[N][N][N][N];

int main()
{
    //freopen("palin.in","r",stdin);
    //freopen("palin.out","w",stdout);
    int l,r,L,R;
    scanf("%s%s",s+1,t+1);
    sn=strlen(s+1),tn=strlen(t+1);
    s[0]='!',s[sn+1]='@',t[0]='#',t[tn+1]='$';
    for(l=0;l<=sn;l++)
        for(r=sn+1;r>l;r--)
            for(L=0;L<=tn;L++)
                for(R=tn+1;R>L;R--)
                {
                    int x=dp[l][r][L][R];
                    if(s[l]==s[r])x=max(x,dp[l-1][r+1][L][R]+2);
                    if(s[l]==t[R])x=max(x,dp[l-1][r][L][R+1]+2);
                    if(t[L]==t[R])x=max(x,dp[l][r][L-1][R+1]+2);
                    if(t[L]==s[r])x=max(x,dp[l][r+1][L-1][R]+2);
                    dp[l][r][L][R]=x;
                }
    for(l=0;l<=sn;l++)
        for(L=0;L<=tn;L++)
        {
            ans=max(ans,dp[l][l+1][L][L+1]);
            ans=max(ans,dp[l][l+2][L][L+1]+1);
            ans=max(ans,dp[l][l+1][L][L+2]+1);
        }
    cout<<ans<<'\n';
    return 0;
}