AcWing 1083. Windy数(数位dp)

这篇博客介绍了如何使用动态规划解决在给定区间[A, B]内计算Windy数(相邻数字差至少为2的正整数)的数量问题。首先初始化状态转移矩阵,接着通过数位DP方法分别计算A-1和B的Windy数,最后通过前缀和计算区间内的Windy数。代码中展示了详细的动态规划实现过程。

在这里插入图片描述
在这里插入图片描述

题意:

Windy 数定义:不含前导零相邻两个数字之差至少为 2正整数

输入两个整数:A 、B

问:在 A 和 B 之间,包括 A 和 B总共有多少个 Windy 数?

思路:

和前两题的思路一样,我们分别对 A - 1、B 两个数进行数位dp操作,求出 dp(B)[0, B] 区间内Windy数的数量),dp(A-1)[0, A-1] 区间内Windy数的数量)。

之后利用前缀和的思想,区间 [A, B]Windy数的数量即为:dp(B) - dp(A-1)

对于一个数 N,我们将其每一位 从高位到低位依次取出放入 num 数组,假设 共有 n

假设我们 当前枚举到第 i,且 i 位上的数字为 x,那么 对于答案中第 i 位数字 j 来说,有两类:

  • 1:0 ~ x - 1 (当然,如果第 i 位是最高位,这里应该改成:1 ~ x - 1,因为题目要求没有前导 0)
    我们使用变量 last,只需 记录上一位的数字,然后枚举 j,如果 abs(j - last) >= 2累加答案 resres += f[i+1][j]。(f[i][j]表示 一共有 i 位,且最高位数字为 j满足windy数定义 的数个数)

  • 2:x

  • 此时显然不可以随意累加 res(因为我们对其后位数字的情况无法判断),使last = x,之后枚举后面的数位

注意一定要有的步骤:上述做完之后,由于上面的答案都是 n 位的,对于数位个数低于 n 的,再累加到答案 res 中就行

在这里插入图片描述

对于f[i][j]的动规分析:

由上方分析过程我们知道我们需要有一种能够直接算出“左分支”答案的方法,考虑使用dp进行f[i][j]的预处理

  • f[i][j]状态表示f[i][j] 表示 一共有 i 位,且最高位数字为 j满足windy数定义 的数的个数。

  • 状态转移: 因为i 位是 j 已经确定,考虑第 i-1 位,设i-1 位数字为 k,那么根据windy数定义只要满足 abs(k-j) >= 2 就可以进行转移。

  • 状态转移方程在这里插入图片描述

代码:

#include<bits/stdc++.h>

using namespace std;

const int N = 15;
typedef vector<int> vi;
#define pb push_back
#define pp pop_back
int f[N][N];

void init()
{
    for(int i=0; i<=9; ++i) f[1][i] = 1;
    for(int i=2; i<N; ++i)
    {
        for(int j=0; j<=9; ++j)
        {
            for(int k=0; k<=9; ++k)
            {
                if(abs(j-k) >= 2) f[i][j] += f[i-1][k];//根据状态转移,满足就加上
            }
        }
    }
}

int dp(int n)
{
    if(!n) return 0;//特判边界 规定 0 不算Windy数,且 0 不在使用区间内
    vi num; while(n) num.pb(n%10), n /= 10;
    int res = 0;//记录答案
    int last = -2;//记录前面某种信息,由于Windy数只和前面相邻数有关,所以last只需记录上一个数即可,由于首位数字可以任意填充,此处只需将last赋值为一个和0~9每个数相比绝对值之差都 >=2的数即可
    
    for(int i=num.size()-1; i>=0; --i)
    {
        int x = num[i];
        for(int j=(i==num.size()-1); j<x; ++j)//枚举左分支
        {
            if(abs(j - last) >= 2) res += f[i+1][j];//从 i 到 0 一共 i+1 位数
        }
        
        if(abs(x - last) < 2) break;//不满足条件,直接break
        
        last = x;//右分支
        if(!i) ++res;
    }
    //<num.size()位的要特殊处理有前导 0 的数
    for(int i=1; i<num.size(); ++i)
    {
        for(int j=1; j<=9; ++j)//i 位数,最高位一定不能为 0,从 1 开始枚举
        {
            res += f[i][j];
        }
    }
    
    return res;
}

int main()
{
    init();
    int l, r;
    cin>>l>>r;
    cout<<dp(r) - dp(l-1)<<'\n';
    
    return 0;
}
### 数位DP算法的学习资源 数位DP是一种用于解决特定范围内的字统计问题的动态规划方法。它通过逐位分析字特性来优化计算效率,适用于许多复杂场景中的计问题。 以下是针对数位DP算法推荐的相关视频教程和学习资料: #### 推荐视频教程 1. **灵神的数位DP通用模板讲解** 灵神提供了详细的数位DP教学视频,涵盖了从基础到高级的应用案例,并附带了多种编程语言实现的代码示例[^1]。该课程适合初学者快速掌握核心概念并应用于实际解题。 2. **据结构与算法一年刷题特训营 - 数位DP模块** 特训营中的第20部分专门讨论了不同类型的动态规划问题,其中章节 `20.7 数位DP` 提供了深入浅出的教学内容以及配套练习题解析[^4]。此资源非常适合希望系统化提升技能的学习者。 #### 推荐学习资料 1. **LeetCode官方文档及相关题目列表** LeetCode上的难度较高的题目如“统计特殊整(No. 2376),“至少有1位重复的字”(No. 1012)等均涉及数位DP的思想[^1]。这些题目不仅帮助理解理论知识,还能锻炼实战能力。 2. **OI竞赛经典题目集合** OI社区贡献了许多经典的数位DP训练题,例如【Luogu1066】2^k 进制、【Bzoj1026】Windy等问题[^3]。这些问题覆盖了各种变体情况下的应用技巧,有助于全面巩固知识点。 3. **今日做题家 - 面试算法题教程系列** 此项目包含了丰富的算法主题介绍及其应用场景说明,虽然是面向面试准备设计的内容,但对于想要深入了解某一领域技术细节的人群同样具有很高的参考价值[^2]。 ```python # 示例代码片段展示如何初始化状态转移方程 (基于 Python 实现) def digit_dp(n, condition_func): s = str(n) @cache def f(i, is_limit, state): if i == len(s): return int(condition_func(state)) res = 0 up = int(s[i]) if is_limit else 9 for d in range(up + 1): next_state = update_state(d, state) res += f(i+1, is_limit and d==up, next_state) return res return f(0,True,initial_state()) ``` 以上即为关于数位DP算法的部分优质视频教程及学习材料汇总信息。建议根据个人需求选取合适的方式展开进一步探索实践。
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值