• 数位dp(模板+例题)


    文章参考:数位dp之总结

    首先,什么是数位dp?它是干什么的?

    数位dp是一种计数用的dp,一般就是要统计一个区间[le,ri]内满足一些条件数的个数。

    举个栗子:
    加入我们要枚举所有上界不超过231的数,那么我们一般的做法就是:

    1 for(int i=le;i<=ri;i++)  
    2         if(right(i)) ans++;  

    这里的le和ri分别是0和231。
    方法可以是可以,但是有时候枚举量会太大,这里采用新的枚举方法:

    从高位起,每次枚举每一位的数,要求每次枚举都不能超过231,比如第一位你的上界是2,那么你就只能枚举0或1或2。当你第一位枚举1的时候,你第二位就可以随便枚举0~9;但当你第一位枚举的是2,接下来那一位便只能枚举0~3。以此类推,总之就是枚举的数不能超过规定的这个上界。与此同时要注意能否出现前导0,这个看题目,代码中的lead代表了前导0。


    模板代码:

     1 typedef long long ll;  
     2 int a[20];  
     3 ll dp[20][state];//不同题目状态不同  
     4 //有人可能想到把dp状态改一下dp[pos][state][limit]就是分别记录不同limit下的个数
     5 ll dfs(int pos,/*state变量*/,bool lead/*前导零*/,bool limit/*数位上界变量*/)//不是每个题都要判断前导零  
     6 {  
     7     //递归边界,既然是按位枚举,最低位是0,那么pos==-1说明这个数我枚举完了  
     8     if(pos==-1) return 1;/*这里一般返回1,表示你枚举的这个数是合法的,那么这里就需要你在枚举时必须每一位都要满足题目条件,也就是说当前枚举到pos位,一定要保证前面已经枚举的数位是合法的。不过具体题目不同或者写法不同的话不一定要返回1 */  
     9     //第二个就是记忆化(在此前可能不同题目还能有一些剪枝)  
    10     if(!limit && !lead && dp[pos][state]!=-1) return dp[pos][state];  
    11     /*常规写法都是在没有限制的条件记忆化,这里与下面记录状态是对应,具体为什么是有条件的记忆化后面会讲*/  
    12     int up=limit?a[pos]:9;//根据limit判断枚举的上界up;这个的例子前面用213讲过了  
    13     ll ans=0;  
    14     //开始计数  
    15     for(int i=0;i<=up;i++)//枚举,然后把不同情况的个数加到ans就可以了  
    16     {  
    17         if() ...  
    18         else if()...  
    19         ans+=dfs(pos-1,/*状态转移*/,lead && i==0,limit && i==a[pos]) //最后两个变量传参都是这样写的  
    20         /*这里还算比较灵活,不过做几个题就觉得这里也是套路了 
    21         大概就是说,我当前数位枚举的数是i,然后根据题目的约束条件分类讨论 
    22         去计算不同情况下的个数,还有要根据state变量来保证i的合法性,比如题目 
    23         要求数位上不能有62连续出现,那么就是state就是要保存前一位pre,然后分类, 
    24         前一位如果是6那么这意味就不能是2,这里一定要保存枚举的这个数是合法*/  
    25     }  
    26     //计算完,记录状态  
    27     if(!limit && !lead) dp[pos][state]=ans;  
    28     /*这里对应上面的记忆化,在一定条件下时记录,保证一致性,当然如果约束条件不需要考虑lead,这里就是lead就完全不用考虑了*/  
    29     return ans;  
    30 }  
    31 ll solve(ll x)  
    32 {  
    33     int pos=0;  
    34     while(x)//把数位都分解出来  
    35     {  
    36         a[pos++]=x%10;//个人老是喜欢编号为[0,pos),看不惯的就按自己习惯来,反正注意数位边界就行  
    37         x/=10;  
    38     }  
    39     return dfs(pos-1/*从最高位开始枚举*/,/*一系列状态 */,true,true);//刚开始最高位都是有限制并且有前导零的,显然比最高位还要高的一位视为0嘛  
    40 }  
    41 int main()  
    42 {  
    43     ll le,ri;  
    44     while(~scanf("%lld%lld",&le,&ri))  
    45     {  
    46         //初始化dp数组为-1,这里还有更加优美的优化,后面讲  
    47         printf("%lld
    ",solve(ri)-solve(le-1));  
    48     }  
    49 }  

    例题:bzoj1026 Windy数
    题意:
    windy定义了一种windy数。不含前导零且相邻两个数字之差至少为2的正整数被称为windy数。 windy想知道,在A和B之间,包括A和B,总共有多少个windy数?
    输入为两个整数,输出windy数的个数。

    代码如下:

     1 #include<cstdio>
     2 #include<algorithm>
     3 #include<string.h>
     4 #include<cmath>
     5 using namespace std;
     6 typedef long long ll;
     7 ll l, r;
     8 ll dp[25][25][2];
     9 int dig[25];
    10 
    11 ll dfs(int pos, int pre, bool limit, bool lead) {
    12     if (pos < 0)return 1;
    13     if (!limit && dp[pos][pre][lead] != -1)
    14         return dp[pos][pre][lead];
    15     ll ans = 0;
    16     int up = limit ? dig[pos]: 9;
    17     for (int i = 0; i <= up; i++) {
    18         if (lead||abs(i - pre) >= 2) {
    19             ans += dfs(pos - 1, i, limit && (i == up), lead && (i == 0));
    20         }
    21     }
    22     if (!limit)dp[pos][pre][lead] = ans;
    23     return ans;
    24 }
    25 
    26 ll cal(ll n) {
    27     int len = 0;
    28     memset(dp, -1, sizeof(dp));
    29     while (n) {
    30         dig[len++] = n % 10;
    31         n /= 10;
    32     }
    33     return dfs(len - 1, 0, true,1);
    34 }
    35 
    36 int main() {
    37     scanf("%lld%lld", &l, &r);
    38     printf("%lld", cal(r) - cal(l - 1));
    39     return 0;
    40 }
  • 相关阅读:
    textarea宽度、高度自动适应处理方法
    Table嵌套去掉子table的外边框
    发现原来自己挺能给自己找理由开脱的
    Life is not the amount of breath you take.
    在遍历ResultSet的循环中再执行SQL会发生什么(前提:同一个Statement)
    按月查询数据
    Oracle SQL 判断某表是否存在
    在Python程序中执行linux命令
    在Oracle中十分钟内创建一张千万级别的表
    Redis Sentinel结构 及相关文档
  • 原文地址:https://www.cnblogs.com/romaLzhih/p/9489812.html
Copyright © 2020-2023  润新知