题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3555

题目大意:求 \([1,n]\) 范围内有多少数包含“49”。

解题思路:

这个问题我们可以分两种解法来考虑:第一种是求不包含“49”的数的数量,用后减一下;另一种就是直接求包含“49”的数的数量。

解法1:求多少数不包含“49”

这种方法我们先通过数位DP求出 \([0,n]\) 区间范围内有多少数不包含“49”(假设数量为 \(x\) ),然后可以得到答案为 \(n+1-x\) 。

我们可以设计一个函数 dfs(int pos, int stat, bool limit) 来返回区间 \([0,n]\) 范围内有多少数不包含“49”,其中:

  • pos 表示当前所处数位;
  • stat 表示前一位的数是不是‘4’(如果前一位是‘4’当前位是‘9’则凑成“49”);
  • limit 用于标记当前是否受限制(true:受限制;false:不受限制)。

实现代码如下:

#include <bits/stdc++.h>
using namespace std;
long long f[66][2], n;
int T, a[66];
void init() {
memset(f, -1, sizeof(f));
}
long long dfs(int pos, int stat, bool limit) {
if (pos < 0) return 1;
if (!limit && f[pos][stat] != -1) return f[pos][stat];
int up = limit ? a[pos] : 9;
long long tmp = 0;
for (int i = 0; i <= up; i ++) {
if (stat && i == 9) continue;
tmp += dfs(pos-1, i==4, limit && i == up);
}
if (!limit) f[pos][stat] = tmp;
return tmp;
}
long long get_num(long long x) {
int pos = 0;
while (x) {
a[pos++] = x % 10;
x /= 10;
}
return dfs(pos-1, 0, true);
}
int main() {
init();
scanf("%d", &T);
while (T --) {
scanf("%lld", &n);
printf("%lld\n", n+1-get_num(n));
}
return 0;
}

解法2:求多少数包含“49”

相对解法1是间接的方式进行求解。我们现在这种方式则是 直接 求解 \([0,n]\) 范围内有多少个数包含“49”。

我们可以设计一个函数 dfs(int pos, int stat, bool limit) 来返回区间 \([0,n]\) 范围内有多少数不包含“49”,其中:

  • pos 表示当前所处数位;
  • stat 表示前一位的数是不是‘4’(如果前一位是‘4’当前位是‘9’则凑成“49”);
  • limit 用于标记当前是否受限制(true:受限制;false:不受限制)。

虽然是一样的,但是在处理的时候统计结果的方式却不一样。

这里特别需要注意的是限制条件和非限制条件下不同的处理。

假设我们现在要找的是区间 \([0,n]\) 范围内有多少数包含“49” 遍历到 pos 位并且 pos+1 位是‘4’, pos 位是‘9’,则我们可以得知:

  1. 如果当前处于限制状态下,则之后有 \(n \% 10^{pos} + 1\) 个数满足条件;
  2. 如果当前处于非限制状态下,则之后有 \(10^{pos}\) 个数满足条件。

实现代码如下:

#include <bits/stdc++.h>
using namespace std;
long long f[66][2], n;
int T, a[66];
void init() {
memset(f, -1, sizeof(f));
}
long long pow10(int a) {
if (a == 0) return 1;
if (a == 1) return 10;
long long t = pow10(a/2);
return t * t * (a%2 ? 10 : 1);
}
long long dfs(int pos, int stat, bool limit) {
if (pos < 0) return 0; // 返回0说明没找到包含49的
if (!limit && f[pos][stat] != -1) return f[pos][stat];
int up = limit ? a[pos] : 9;
long long tmp = 0;
for (int i = 0; i <= up; i ++) {
if (stat && i == 9) { // 说明接下来pos位的10的pos次方个数都满足条件
// tmp += pow10(pos); // 直接这么写是错的,此时也要考虑边界条件限制
// 修正如下
if (limit) tmp += n % pow10(pos) + 1;
else tmp += pow10(pos);
}
else tmp += dfs(pos-1, i==4, limit && i == up);
}
if (!limit) f[pos][stat] = tmp;
return tmp;
}
long long get_num(long long x) {
int pos = 0;
while (x) {
a[pos++] = x % 10;
x /= 10;
}
return dfs(pos-1, 0, true);
}
int main() {
init();
scanf("%d", &T);
while (T --) {
scanf("%lld", &n);
printf("%lld\n", get_num(n));
}
return 0;
}

HDU3555 Bomb 题解 数位DP的更多相关文章

  1. [Hdu3555] Bomb(数位DP)

    Description 题意就是找0到N有多少个数中含有49. \(1\leq N \leq2^{63}-1\) Solution 数位DP,与hdu3652类似 \(F[i][state]\)表示位 ...

  2. 【Hdu3555】 Bomb(数位DP)

    Description 题意就是找0到N有多少个数中含有49. \(1\leq N \leq2^{63}-1\) Solution 数位DP,与hdu3652类似 \(F[i][state]\)表示位 ...

  3. HDU 3555 Bomb(数位DP)

    Bomb Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 131072/65536 K (Java/Others) Total Subm ...

  4. HDU 3555 Bomb(数位DP模板啊两种形式)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3555 Problem Description The counter-terrorists found ...

  5. 动态规划晋级——HDU 3555 Bomb【数位DP详解】

    转载请注明出处:http://blog.csdn.net/a1dark 分析:初学数位DP完全搞不懂.很多时候都是自己花大量时间去找规律.记得上次网络赛有道数位DP.硬是找规律给A了.那时候完全不知数 ...

  6. POJ-2282题解&数位DP总结

    一.题意 给定一个区间[a, b](注意输入的时候可能a > b,所以,在数据输入后,要先比较a和b,如果a > b,交换a和b的值),统计这个区间里面,数位上有多少个0.多少个1.--. ...

  7. luogu2657-Windy数题解--数位DP

    题目链接 https://www.luogu.org/problemnew/show/P2657 分析 第一道数位DP题,发现有点意思 DP求\([L,R]\)区间内的XXX个数,很套路地想到前缀和, ...

  8. 洛谷P2602 [ZJOI2010]数字计数 题解 数位DP

    题目链接:https://www.luogu.com.cn/problem/P2602 题目大意: 计算区间 \([L,R]\) 范围内 \(0 \sim 9\) 各出现了多少次? 解题思路: 使用 ...

  9. 洛谷P3413 SAC#1 - 萌数 题解 数位DP

    题目链接:https://www.luogu.com.cn/problem/P3413 题目大意: 定义萌数指:满足"存在长度至少为2的回文子串"的数. 求区间 \([L,R]\) ...

随机推荐

  1. oracle函数 ASCII(x1)

    [功能]:返回字符表达式最左端字符的ASCII 码值. [参数]:x1,字符表达式 [返回]:数值型 [示例] SQL> select ascii('A') A,ascii('a') a,asc ...

  2. Python深入:stevedore简介

    stevedore是用来实现动态加载代码的开源模块.它是在OpenStack中用来加载插件的公共模块.可以独立于OpenStack而安装使用:https://pypi.python.org/pypi/ ...

  3. Linux 正文处理命令及tar命令 利用vi编辑器创建和编辑正文文件

    要点回顾 1) 将用户信息数据库文件和组信息数据库文件纵向合并为一个文件/1.txt(覆盖) cp /etc/passwd . cat ./passwd >1.txt cp /etc/group ...

  4. Eclipsed的SVN插件不能识别之前工作空间的项目

    今天用Eclipse的SVN插件,由于之前项目是用"tortoisesvn"进行checkout的,并且已经在工作空间中使用.后来用来Ecplise的SVN插件,可是发现右键Tea ...

  5. 学习CSS预处理器Less

    1.Sass与Less的区别:Sass与Less的区别 2.Less的使用:Less的使用 3.Less的相关网址:Less.js.Less中文网 Less的HSL函数-lighten():HSL函数 ...

  6. Python--day25--面向对象之多态

    多态(Python天生支持多态) 多态指的是一类事物有多种形态 动物有多种形态:人,狗,猪 import abc class Animal(metaclass=abc.ABCMeta): #同一类事物 ...

  7. Python--day25--接口类

    接口类,起规范作用

  8. Spring Data Jpa一对多单向映射

    /** @author StormMaybin @date 2017-01-17 */ 生命不息,奋斗不止! 一对多映射关系 在JPA中,用@OneToMany来标识一对多的关系.实现一对多的单向关联 ...

  9. centos7搭建集群必知:centos7已经无iptables,只有firewall

    1.防火墙概述 centos7搭建集群,发现没有iptables,需要安装.防火墙为firewalle CentOS7默认的防火墙不是iptables,而是firewalle. CentOS 7.0默 ...

  10. 无限调用函数add(1)(2)(3)......

    无限调用函数,并且累计结果 其实这也算一道面试题吧,笔者曾经被提问过,可惜当时没能答上来...