Good BYe 2015 D New Year and Ancient Prophecy

36 篇文章 0 订阅

传送门:
http://codeforces.com/contest/611/problem/D

题意:把一个数字字符串分割成几个严格无前导0的单调递增的序列,问答案mod(1e9+7)
很明显的dp
以下的个人初始想法完全可以忽略掉~~~
粗略的想法就是:用dp[i]表示以i结尾,且数字以i结束的方法数,那么转移方程就是
如果i到j这段字符串要比以 j结尾的最大字符串大,的话dp[i]就等于dp[j],
这样的复杂度就已经是n平方了,那么现在应该如何快速判断O(1)判断出字符串的大小呢,
我们可以构造这个方法,即记录i位置和j位置开始第一个不同的字符串的位置,那么比较的时候取出这个位置然后比较大小就可以了,完美地实现了O(1)的复杂度
貌似应该开两个数组,一个记录的是以一个字符结尾的最后一个串的最小长度,
那么这个要怎么去更新呢?
就是要去看i到j字符串可否拼接在dpj到到j之后嘛,判断一下就好了
还有就是首先可以通过字符串的长短来判断是不是可以直接拼在后面,我说的那个需要比较长度的情况应该是二者长度相等的情况!!!
另一个则是记录以一个字符结尾的答案数

那么事实是这样显然概括不了所有的状态,即我的那个记录最小长度是完全不可取的!!!!

好吧冥想时间已然结束了,其实我那样子做是不行的,因为状态就已经结束了,能用二维记录状态为什么要用一维呢,岂不是自找麻烦?
用二维会记录的什么清晰

ps:有的时候写完之后真的是要仔细看一眼代码,说不定哪个地方就写错了那么一点,还不好用数据检查出来
题解只是用了更新一个前缀和数组就ok了,传引用的方法逼格的确是高
我觉得最重要的还是要注意边界条件,初值等等的问题,还有就是状态更改的方法,如何用正确的方式去记录状态才能保证解得以更新且不会有遗漏
还有就是记录nxt数组的那个位置,本来就应该倒着往前面去写,注意写这种东西的时候一定要注意用的 是dp的思想,从前往后更新的话会很麻烦,而且啰嗦!!!!!!

每当下标要增减的时候就一定要注意特殊判断一下会不会有情况会发生越界
题解:

#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
const int maxn=5e3+10;
const int inf=0x3f3f3f3f;
int n;
//pref[a][b]=dp[1...a][b]
int pref[maxn][maxn];
int nxt[maxn][maxn];
char s[maxn];int t[maxn];
int dp(int a,int b){
    int x=pref[a][b]-pref[a-1][b];
    if(x<0) x+=mod;
    if(x>=mod) x-=mod;
    return x;
}
int main(){
   scanf("%d%s",&n,s+1);
   for(int i=1;i<=n;i++){
    t[i]=s[i]-'0';
   }
   for(int a=n;a>=1;a--){
    for(int b=n;b>a;b--){
        int &x=nxt[a][b];
        if(t[a]!=t[b]) x=a;
        else if(b==n) x=inf;
        /*why inf  just think about preones come to this situation it should stop,in other words , never ends!!!!
        just beacause of the situation now is illegal,and the condition we judeg below is that whree<c ,so we should put this one to inf*/
        else x=nxt[a+1][b+1];
    }
   }
   for(int i=1;i<=n;i++){
    //dp[1][i]=1 and we starts from 2 next
    for(int j=1;j<=n;j++){
        pref[i][j]=1;
    }
   }
   // c starts from 2
   for(int c=2;c<=n;c++){
    for(int d=c;d<=n;d++){
        int &x=pref[c][d];
        int b=c-1;
        int a=b-(d-c);
        // be warn of this,it is not x+= because x is 1
        x=pref[b][b]-pref[max(a,0)][b];
        if(x>=mod) x-=mod;
        //miss this
        if(x<0)  x+=mod;
        if(a>=1){
            int whree=nxt[a][c];
            //miss whree<c
            if(whree<c&&t[whree]<t[whree+c-a]){
                x+=dp(a,b);
                if(x>=mod) x-=mod;
            }
        }
        if(s[c]=='0') x=0;
        pref[c][d]=x+pref[c-1][d];
        if(pref[c][d]>=mod) pref[c][d]-=mod;
    }
   }
   int ans=0;
   for(int i=1;i<=n;i++){
    ans=(ans+dp(i,n))%mod;
   }
   printf("%d\n",ans);
   return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值