[dp] hdu 5282 Senior's String

最新推荐文章于 2026-05-17 06:34:59 发布

原创最新推荐文章于 2026-05-17 06:34:59 发布 · 536 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

动态规划同时被 2 个专栏收录

59 篇文章

订阅专栏

字符串

25 篇文章

订阅专栏

本文详细介绍了如何使用动态规划算法解决最长公共子序列问题，并在此基础上计算给定字符串集合中特定长度子序列在另一个字符串中出现的次数。通过实例演示了从动态规划到最终答案计算的全过程。

题意：

学姐姐非常喜欢字符串，所以学弟送给了她两个字符串作为礼物。

两个字符串分别为X，Y。她非常开心，但在开心之余她还想考考学弟。

她定义L为X与Y的最长公共子序列的长度（子序列在字符串内不一定连续，一个长度为L的字符串有2L个子序列，包括空子序列）。

现在学姐姐取出了X的所有长度为L的子序列，并要求学弟回答在这些子序列中，有多少个是Y的子序列。

因为答案可能很大，所以学弟只需要回答最终答案模109+7。

思路：

首先我们用O(n2)的动态规划算法处理出dp数组，dp[i][j]表示X串的前i个字符和Y串的前j个字符的最长公共子序列的长度，在这个基础上我们再进行一个动态规划。用f[i][j]表示在X串的前i个字符中，有多少个长度为dp[i][j]的子序列在Y的前j个字符中也出现了。转移：若dp[i−1][j]==dp[i][j]，则f[i][j]+=f[i−1][j]，表示i这个字符不选；再考虑选i这个字符，找到Y串前j个字符中最靠后的与X[i]匹配的字符的位置，设为p，若dp[i−1][p−1]+1==dp[i][j]，则f[i][j]+=f[i−1][p−1]。最终的答案即为f[n][m]。复杂度O(n2)。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<iostream>
#include<stack>
using namespace std;
#define ll __int64
int dp[1234][1234],wz[1234][1234],f[1234][1234];
char x[1555],y[1555];
int mod=1000000007;
int main()
{

    int t;
    cin>>t;
    while(t--)
    {
        scanf("%s%s",x+1,y+1);
        memset(dp,0,sizeof(dp));
        int n,m;
        n=strlen(x+1);
        m=strlen(y+1);
        for(int i=1; x[i]; i++)
        {
            for(int j=1; y[j]; j++)
            {
                if(x[i]==y[j]) dp[i][j]=dp[i-1][j-1]+1;
                else dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
            }
        }
        memset(wz,0,sizeof(wz));
        for(int i=1; i<=m; i++)
        {
            for(int j=0; j<26; j++) wz[i][j]=wz[i-1][j];
            wz[i][y[i]-'a']=i;
        }
        memset(f,0,sizeof(f));
        for(int i=0; i<=n; i++)
        {
            for(int j=0; j<=m; j++)
            {
                if(dp[i][j]==0)  //空串
                {
                    f[i][j]=1;
                    continue;
                }
                if(dp[i][j]==dp[i-1][j]) f[i][j]+=f[i-1][j];
                if(f[i][j]>=mod) f[i][j]%=mod;
                int p=wz[j][x[i]-'a'];
                if(p==0) continue;
                if(dp[i][j]==dp[i-1][p-1]+1) f[i][j]+=f[i-1][p-1];
                if(f[i][j]>=mod) f[i][j]%=mod;
            }
        }
        printf("%d\n",f[n][m]);
    }
    return 0;
}