codeforces-Goodbye2015- New Year and Ancient Prophecy

本文探讨了解决复杂问题时使用动态规划方法的策略,包括寻找最长连续子串、状态转移方程的设计以及如何避免重复计算。通过实例分析,展示了如何通过动态规划有效地解决这类问题,最终实现算法优化。

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

第一次做这种难度的题刚开始做一点思路没有,看完别人的思路做出来以后感觉好简单我居然当时没有思路!
1.先求最长连续子串。因为后面要比较两个字串的大小如果在下面DP的时候在一次次比较两个数的大小妥妥会爆的。
2.DP
(1)dp[i][j]代表最末尾为s[i]s[i+1]…s[j]这个数时的方案数,如果末尾为s[i]s[i+1]…s[j]这个数的方案可以,那这些方案同样适用s[i]s[i+1]…s[j]s[j+1]
因为s[i]s[i+1]…s[j]s[j+1]>s[i]s[i+1]…s[j]
所以dp[i][j] = dp[i][j]+dp[i][j-1] (1<=j<=n) 所以dp到i时先通过这个式子叠加一下
(2)当叠加以后就是真正的dp[i][j]了;
然后开始推以后的状态了
取j从i开始取使它满足j-i+1<=n-j 因为后面那个数比前面那个数大所以位数肯定不能小于!
看lcp[i][j+1]的值如果值大于等于j-i+1 那么s[i]s[i+1]…s[j]与s[j+1]…s[2*j+i+1]相等 那只能选s[j+1]…s[2*j+i+1]s[2*j+i+2]才能比前者大,所以dp[j+1][2*j-i+2] = (dp[j+1][2*j-i+2]+dp[i][j])
如果lcp[i][j+1]的结果小于j-i+1,那么比较s[i]s[i+1]…s[j]和s[j+1]…s[2*j+i+1]第一个不同的数来比较大小如果前者小于后者那么
dp[j+1][2*j-i+1] = (dp[j+1][2*j-i+1]+dp[i][j])
否则dp[j+1][2*j-i+2] = (dp[j+1][2*j-i+2]+dp[i][j])
(3)注意取摸

#include<cstdio>
#include<cmath>
#include<cstring>
#define mod 1000000007
#define MAXN 5005
#define LL long long
int dp[MAXN][MAXN],lcp[MAXN][MAXN],n; 
char s[5005];
int main()
{
   while(scanf("%d",&n)!=EOF)
   {
     memset(dp,0,sizeof(dp));
     memset(lcp,0,sizeof(lcp));
     scanf("%s",s+1);
     //从后往前求连续子串
     for(int i=n;i>=1;i--)
        for(int j=n;j>=1;j--)
           if(s[i]==s[j])lcp[i][j] = lcp[i+1][j+1]+1;
           else lcp[i][j] = 0;
      //dp
     dp[1][1] = 1;
     for(int i=1;i<=n;i++)
        {
           if(s[i]=='0')continue;
           for(int j=i;j<=n;j++)
            dp[i][j] = (dp[i][j]+dp[i][j-1])%mod;
           for(int j=i;j-i+1<=n-j;j++)
            {
                 if(s[j+1]=='0')continue;
                 int t = lcp[i][j+1];
                 if(t>=j-i+1)dp[j+1][2*j-i+2] = (dp[j+1][2*j-i+2]+dp[i][j])%mod;
                 else {
                    if(s[i+t]<s[j+1+t])dp[j+1][2*j-i+1] = (dp[j+1][2*j-i+1]+dp[i][j])%mod;
                    else dp[j+1][2*j-i+2] = (dp[j+1][2*j-i+2]+dp[i][j])%mod;
                      }
            }
        }
        int ans = 0;
        for(int i=1;i<=n;i++)
          ans = (dp[i][n]+ans)%mod;
        printf("%d\n",ans);  
   }
   return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值