数位dp思想

一、基础篇——介绍

 

1、概念：

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

     所谓数位dp，字面意思就是在数位上进行dp。

     数位的含义：一个数有个位、十位、百位、千位......数的每一位就是数位啦！

     之所以要引入数位的概念完全就是为了dp。数位dp的实质就是  换一种暴力枚举的方式，使得新的枚举方式满足dp的性质，然      后记忆化就可以了。

     数位DP其实是很灵活的，所以一定不要奢求一篇文章就会遍所有数位DP的题，这一篇只能是讲清楚一种情况，其他情况遇到      再总结，在不断总结中慢慢体会这个思想，以后说不定就能达到一看到题目就能灵活运用的水平。（其实DP都是这样……）

      ex:

      例1.求  a~b  中不包含49的数的个数 . 0 < a、b < 2*10^9

      注意到n的数据范围非常大，暴力求解是不可能的，考虑dp,如果直接记录下数字，数组会开不起，该怎么办呢？

       此时就要用到数位dp.

 

2、用处：

   一般应用于：

 （1）、  求出在给定区间  [A,B]   内，符合条件  P(i)   的数 i 的个数.

  （2）、条件 P(i)  一般与 数 的大小无关，而与  数的组成  有关.

3、这样，我们就要考虑一些特殊的记录方法来做这道题.一般来说，要保存给定数的每个位置的数.然后要记录的状态为当前操作数的位数，剩下的就是根据题目的需要来记录.可以发现，数位dp的题做法一般都差不多，只是定义状态的不同罢了.

 

二、实战篇——例题分析

1、

例：

这一篇要说的数位DP是一道最简单的数位DP：题目链接

题目大意：多组数据，每次给定区间  [n,m]  ，求在n到m中没有  “62“ （连续）或 “4“ 的数的个数。

　　　　　如62315包含62，88914包含4，这两个数都是不合法的。  

                  0 < n<=m < 1000000 

（1）、两种不同的枚举：对于一个求区间[le,ri]满足条件数的个数，最简单的暴力如下：

               

1. for(int i=le;i<=ri;i++)

2. if(right(i)) ans++;

  然而这样枚举不方便记忆化，或者说根本无状态可言。

 （2）新的枚举——————数位 dp

A:

思路分析：

试想：我们如果能有一个函数count(int x)，可以返回[0,x]之间符合题意的数的个数。

那么是不是直接输出count(m)-count(n-1)就是答案？

好，那么下面我们的关注点就在于怎么做出这个函数。我们需要一个数组。（dp原本就是空间换时间）

我们设一个数组 dp[ i ][ j ]，表示 i 位数，最高位是j 的符合题意的个数。

比如dp[ 1 ][ 2 ]=1; （位数是 1,高位是2 的只能是 2本身）

dp[ 1 ][ 4 ]=0;  （ 位数是 0,高位是 4 的不符合题意，个数 是 0）

dp[ 2 ][ 0 ]=9;    ( 当 j 等于 4 的时候不符合题意，就是 9个，00,01,02,03,05,06,07,08,09 ）

dp[ 2][ 6 ]=8       (60,61,63,65,66,67,68,69).

我们先不关注这个 dp 有什么用，我们先关注 dp 本身怎么求。

首先

1. if(i==4)

2. dp[1][i]=0;

3. else

4. dp[1][i]=1;

这一步是很显然的，那么根据这个题的数据范围，只需要递推到  dp[ 7 ][ i ] 就够用了。那么稍微理解一下，可以想出递推式：

1. if(j==4)

2. dp[i][j]=0;

3. else if(j==6)

4. {

5. for(int k=0; k<=9; k++)

6. {

7. if(k!=2)

8. dp[i][j]+=dp[i-1][k];

9. }

10. }

11. else

12. {

13. for(int k=0; k<=9; k++)

14. dp[i][j]+=dp[i-1][k];

15. }

上面的式子也是很显然的，如果觉得不显然可以这样想：

i  位数，最高位是  j   的符合条件的数，如果 j  是4，肯定都不符合条件（因为题目不让有4），所以直接是0；

如果  j 不是6，那么它后面随便取，只要符合题意就可以，所以  dp[ i - 1 ][ k ] ,k可以随便取的和；

如果  j 是6，后面只要不是 2 就行，所以是 dp[ i -1 ][ k ],   k除了2都可以，求和即可

 

这里要说明一下，认为00052是长度为5，首位为0的符合条件的数，052是长度为3首位为0符合条件的数。

 

那么现在我们已经得到了  dp 数组，再重申一下它的含义：i位数，最高位是j的数，符合题意的数有多少个。

 

现在我们就要关注怎么利  dp 数组做出上面我们说的那个函数 count(int x) ，它可以求出 [0,x] 中符合题意的数有多少个

那么我们做这样一个函数  int solve(int x)  它可以返回   [0,x)  中符合题意的有多少个。

那么solve(x+1)  实际上与count(x)是等价的 ( solve(x) 求得是不等于 x 的个数，那么 solve(x+1) 求的是小于 x+1 的个数，也就是小于等于 x 的个数）

那么现在问题转化成了：  小于x，符合题意的数有多少个？

很简单，既然小于，从最高位开始比，必定有一位要严格小于x（前面的都相等）。所以我们就枚举哪一位严格小于（前面的都相等）。

假设我们现在把x分成了  a1,a2,…,aL （ai  是每位数）这样一个数组，长度为 L，aL 是最高位。

那么结果实际上就是这样：长度为L，最高位取 [0,aL-1] 的所有的符合题意数的和；再加上长度为L-1，最高位取aL，次高位取[0,a（L-1）-1]的所有符合题意数的和；再加上……；一直到第一位。

B:

特判：

最高位aL如果是4，那么这句话直接就可以终止了，4 以后的都不符合题意

最高位 aL如果是6，次高位 取 2，，也可以终止了 。加上这些条件之后就很严谨了。C: 我们要求[a,b]不包含49的数的个数，可C:求和

想到利用 前缀和来做，具体来说，就是  [a,b]  = solve(b+1)-solve(a)  （等于 b的数 减去小于 a的数，就是 [a,b] 区间的值）

                                    solve(x)  求解的是 不等于 x的数

2、CODE:

1. #include<iostream>

2. #include<cstdio>

3. #include<cstring>

4. #include<string>

5. #include<cmath>

6. #include<iomanip>

7. #include<algorithm>

8. #define max(a,b) (a>b?a:b)

9. #define min(a,b) (a<b?a:b)

10. #define swap(a,b) (a=a+b,b=a-b,a=a-b)

11. #define memset(a,v) memset(a,v,sizeof(a))

12. #define X (sqrt(5)+1)/2.0 //Wythoff

13. #define Pi acos(-1)

14. #define e 2.718281828459045

15. #define eps 1.0e-8

16. using namespace std;

17. typedef long long int LL;

18. typedef pair<int,int>pa;

19. const int MAXL(1e5);

20. const int INF(0x3f3f3f3f);

21. const int mod(1e9+7);

22. int dir[4][2]= {{-1,0},{1,0},{0,1},{0,-1}};

23. LL dp[10][20];

24. int len;

25. void getDp()

26. {

27. memset(dp,0);

28. dp[0][0]=1;

29. for(int i=1;i<=7;i++)

30. {

31. for(int j=0;j<=9;j++)

32. {

33. if(j==4)

34. dp[i][j]=0;

35. else if(j==6)

36. {

37. for(int k=0;k<=9;k++)

38. {

39. if(k!=2)

40. dp[i][j]+=dp[i-1][k];

41. }

42. }

43. else

44. {

45. for(int k=0;k<=9;k++)

46. dp[i][j]+=dp[i-1][k];

47. }

48. }

49. }

50. }

51. int a[20];

52. LL solve(int n) // 求的是 小于本身个数

53. {

54. // len=1;

55. // while(n)

56. // a[len++]=n%10,n/=10;

57. // a[len+1]=0; // 必不可少的 ，否则当下文 a[i+1] 的时候不知道是什么值

58. // LL ans=0;

59. // for(int i=len-1;i>=1;i--)

60. // {

61. // for(int j=0;j<a[i];j++)

62. // {

63. // if(j==4||a[i+1]==6&&j==2)

64. // continue;

65. // else

66. // ans+=dp[i][j];

67. // }

68. // if(a[i]==4)

69. // break;

70. // if(a[i+1]==6&&a[i]==2)

71. // break;

72. // }

73. len=0; —————————— 和注释的都可以使用

74. while(n)

75. a[++len]=n%10,n/=10;

76. a[len+1]=0;

77. LL ans=0;

78. for(int i=len;i>=1;i--)

79. {

80. for(int j=0;j<a[i];j++) // 这里小于是求的小于它本身的个数

81. {

82. if(j==4||a[i+1]==6&&j==2)

83. continue;

84. else

85. ans+=dp[i][j];

86. }

87. if(a[i]==4) // ————————————— 1

88. break;

89. if(a[i+1]==6&&a[i]==2)

90. break;

91. }

92. return ans;

93.  

94. }

95. int main()

96. {

97. getDp();

98. int n,m;

99. while(cin>>m>>n&&(m||n))

100. {

101. LL ans1=solve(m);

102. LL ans2=solve(n+1);

103. cout<<ans2<<endl;

104. cout<<ans1<<endl;

105. cout<<ans2-ans1<<endl;

106. }

107. }

108.  

1——————:  放在这里 break 的原因是
                            当 4 或者是 62 连续开头的时候，以后的都不符合题意了，所以就直接停止循环