初探数位dp 介绍了数位类统计的基础知识。以下列出其中的基础点:

基本问题

统计在区间[l, r]中满足条件的数的个数

思路

1. [l, r] 将问题转换为 在[0, r]中满足条件的个数 -[0, l)满足条件的个数

2. 求解 区间[0, n]满足条件的个数

  • 性质:一个小于n的数m,一定是从高位到低位在某一位小于n的对应位的一个数
  数值 百位 十位 个位  
n 358 3 5 8  
m 350 3 5 0 0<8
m 349 3 4 9 4<5
  • 遍历所有小于n的数:从高位到低位枚举第一次小于n的数位,之后数位的值就不受限制了
  数值 百位 十位 个位
n 358 3 5 8
m   3 4 0...9
m     2 0...9 0...9
  • 遍历可以构成一棵树,类似深度优先遍历。第一个图是小于35的数,第二个图是一般情况,i+1的数位比i的大

  • 递归计算区间[0,n]中满足条件的数
数位 i i-1 ...
数位上的值 digit digit' 0...可以取的最大值 ...

F[i, digit, st]表示长度为i,第i位的位为digit,且满足条件st的数的个数。F[i, digit, st] = F[i, digit, st] + F[i-1, digit', st']

  • 动态规划:函数的参数表明了一个状态
 void dp(int 当前数的长度len, int 当前数位的值 digit, bool 当前数前面数位都是0 begin_zero, bool 下一个数位的值有限制 islimit)
{
if(len == )//当前数的长度为1
{ ...//直接讨论返回
return;
}
//遍历下一位的所有取值情况
int maxValue = (islimit?n[len]:);//下一位的最大取值
for(int i = ; i <= maxValue; i++)
{
bool _begin_zero = (i == && begin_zero);
bool _islimit = (i == maxValue);
dp(len - , i, _begin_zero, _islimit);
}
}

Tips:

  • 假设区间[0, n]中n的位数为 len(n),初始要从len(n)+1开始。这样对第len(n)位讨论的时候就可以和其它位一样。
  • 计算过程中,只有当数位的值没有限制时才把结果存入F[i, digit, st]中,否则会缺失一些结果
  • 前面的数位都为0不全为0需要作为两种状态存储

源码 http://hihocoder.com/problemset/problem/1033

 #include<iostream>
using namespace std; #define ll long long int //需要用long long 来进行存放
const int mod = ;//结果的最大值,用于取模
struct node{
ll s, n;//s: 各数位交错和 n:满足条件的个数
};
node dp[][][];//[长度][以哪个数作为开头][该数的数位交错和]
int bits[];//存放一个数的各位,低到高存放在0到20
ll base[];//每一数位的基准 比如十位数=base[2]=10 //len数位长度, dig是首个数字, begin_zero表示从最高位到当前位是否全部为0, limit表示下一位枚举是否有限制(bit[len-2]或9), sum是要求的数字和
node dfs(int len, int dig, bool begin_zero, bool limit, int sum){
node t;//存放长度为 len的结果
t.s = , t.n = ;
//超过边界值
if (len <= || len >= || dig < || dig > || sum < - || sum >= )
return t;
//返回已有的DP结果,即记忆化搜索
if (!limit && dp[len][dig + (begin_zero ? : )][sum + ].n != -)
return dp[len][dig + (begin_zero ? : )][sum + ];
//长度只有一位,就不需要枚举下一位了,直接讨论返回即可
if (len == ){
if (dig != sum)
return t;
t.n = , t.s = sum;
return t;
}
//开始枚举下一位的数字
int end = limit ? bits[len - ] : ;//下一位数字的最大值
int newsum = dig - sum;
node tmp;
for (int j = ; j < end + ; j++)
{
if (begin_zero){//前面都是0,接下来的就由当前位决定 j是否为0
tmp = dfs(len - , j, j == , limit && (j == end), sum);
}
else{//前面不是全为0
tmp = dfs(len - , j, false, limit && (j == end), newsum);
}
//将tmp的值累加到t上
t.n += tmp.n;//满足条件的个数
//计算满足条件的长度为 len 的data[len]所有数的和, tmp 是长度为 i-1的data[len-1]所有数的和
//每一个: data[len] = dig * base[len] + data[len-1] 且共有n个
//t.s = t.s + tmp.n * (dig * base[len]) + tmp.s 增加 %mod 就得到下面的结果
t.s = ((t.s + tmp.s) % mod + ((tmp.n * dig) % mod * base[len]) % mod) % mod;
}
//当长度为len,且以dig开头的,数位和为sum 的所有结果都计算完成,才将其进行存储
if (!limit)
dp[len][dig + (begin_zero ? : )][sum + ] = t;
//dig + (begin_zero ? 0 : 10) 用来区分两种状态,(1. 前导都为0;2. 前导包含其它数)
return t;
} int solve(ll n, int s){
if (n <= )
return ;
int l = ;
for (int i = ; i < ; i++)
bits[i] = ;
//将n的每一位从低到高放到 bits[0]到bits[l] 中
while (n){
bits[l++] = n % ;
n /= ;
}
//从l+1开始,比n的长度大1,并且第l+1位数置为0
return dfs(l + , , true, true, s).s;
} int main(){
ll l, r, s;
node t;
t.n = -;
t.s = ;
for (int i = ; i < ; i++)//长度
for (int j = ; j < ; j++)//第i位的取值
for (int k = ; k < ; k++)// i到0 所有位上的交错和 + 200
dp[i][j][k] = t;
base[] = ;
for (int i = ; i < ; i++)//base 是作为基准数 比如两位数需要 *10, 三位数需要*100
base[i] = base[i - ] * % mod;
cin >> l >> r >> s;
cout << (solve(r, s) - solve(l - , s) + mod) % mod << endl;
return ;
}

【hihoCoder】1033: 交错和的更多相关文章

  1. [hihocoder 1033]交错和 数位dp/记忆化搜索

    #1033 : 交错和 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描写叙述 给定一个数 x,设它十进制展从高位到低位上的数位依次是 a0, a1, ..., an - 1 ...

  2. hihoCoder #1033 : 交错和 (数位Dp)

    题目大意: 给定一个数 x,设它十进制展从高位到低位上的数位依次是 a0, a1, ..., an - 1,定义交错和函数: f(x) = a0 - a1 + a2 - ... + ( - 1)n - ...

  3. hihoCoder 1033 : 交错和 数位dp

    思路:数位dp,dp(i, j, k)表示考虑i位数,每位数可以任意取[0~9],并且这i位数的交错和为j,k=1表示前缀全是0(如000456),k=0表示前缀不为0.注意,前缀是否为0是这道题的一 ...

  4. hihoCoder 1033: 交错和

    (1)题目描述: 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 给定一个数 x,设它十进制展从高位到低位上的数位依次是 a0, a1, ..., an - 1,定义交错 ...

  5. HihoCoder 1033交错和(数位DP第三题)

    (写挂了,有空再补) 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 给定一个数 x,设它十进制展从高位到低位上的数位依次是 a0, a1, ..., an - 1,定义 ...

  6. 【HIHOCODER 1033 】 交错和(数位DP)

    描述 输入 输入数据仅一行包含三个整数,l, r, k(0 ≤ l ≤ r ≤ 1018, |k| ≤ 100). 输出 输出一行一个整数表示结果,考虑到答案可能很大,输出结果模 109 + 7. 提 ...

  7. hihoCoder 1033

    题目链接: http://hihocoder.com/problemset/problem/1033 听说这个题是xiaodao出的~~ 我们要知道dp其实就是一个记忆化搜索的过程,如果某个子结构之前 ...

  8. hiho#1033 : 交错和

    描述 给定一个数 x,设它十进制展从高位到低位上的数位依次是 a0, a1, ..., an - 1,定义交错和函数: f(x) = a0 - a1 + a2 - ... + ( - 1)n - 1a ...

  9. 数位dp/记忆化搜索

    一.引例 #1033 : 交错和 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 给定一个数 x,设它十进制展从高位到低位上的数位依次是 a0, a1, ..., an  ...

随机推荐

  1. 关闭Linux防火墙(iptables) 及 SELinux

    一.关闭防火墙 1.重启后永久性生效: 开启:chkconfig iptables on 关闭:chkconfig iptables off 2.即时生效,重启后失效: 开启:service ipta ...

  2. ACM/ICPC 之 ACM计算机工厂-EK算法(POJ3436)

    题意有点难读懂 //网络流-EK算法-ACM计算机工厂-构图重点 //Time:0Ms Memory:208K #include <iostream> #include<cstrin ...

  3. 【python】确保文件写入结束

    今天遇到了个问题: 我在执行如下代码时发现,文件只写了一半.也就是说,当文件量过大时,下面的代码是不能保证文件被正确写入的. fd = open('test.txt','w') fd.write(&q ...

  4. .NET LINQ 相等运算

    相等运算      如果两个序列的对应元素相等且这两个序列具有相同数量的元素,则视这两个序列相等. 方法 方法名 说明 C# 查询表达式语法 Visual Basic 查询表达式语法 更多信息 Seq ...

  5. 调试asp.net网页时不显示treeview的原因

    在.net中本地调试asp.net网页时,treeview控件显示为文字方式,原因是在http://localhost/下面找不到webctrl_client的路径,解决的方法是把webctrl_cl ...

  6. scrollView的讲解

    今天就讲下UIScrollView的一些事情,这个可以拖动的组件无论在应用还是游戏开发都会经常用到,所以我们就一定要更加熟悉它了.下面我们开始下手咯. (1)初始化 一般的组件初始化都可以alloc和 ...

  7. hive 普通创建表和跟新列操作

    创建表 CREATE TABLE if not exists student ( student_id int, sex int, address String, email String ) 这里需 ...

  8. Note1: windows下Python IDE(Eclipse+Pydev)

    毕设想做一个个性化音乐推荐系统,真的是从零开始学起,通过<推荐系统实践>大致地了解了一些基础知识,目前在啃<集体智慧编程>. 选用的开发环境是 Eclipse+Pydev . ...

  9. css笔记1: html页面的CSS、DIV命名规则

    原地址:http://www.cnblogs.com/rising-fay/archive/2013/02/25/2932592.html CSS命名规则 头:header 内容:content/co ...

  10. Android 笔记 day1