【算法竞赛学习笔记】数位DP

最新推荐文章于 2025-03-02 12:33:32 发布
原创 最新推荐文章于 2025-03-02 12:33:32 发布 · 489 阅读 · 2 ·
本内容遵循CC 4.0 BY-SA版权协议
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
标签
#算法 #学习 #动态规划 #c++ #数据结构
本文详细介绍了数位动态规划（DP）的解题思路和常见问题类型，包括数字计数、不降数个数、避免特定数字的数的个数以及满足特定条件的数字平方和等。通过数位DP分析法，将数字分解成一棵树并计算每个部分对答案的贡献，从而解决一系列计数问题。文章提供了多个例题和代码实现，帮助读者理解和掌握数位DP的应用。

title : 数位DP
date : 2022-4-11
tags : ACM,动态规划
author : Linno

数位DP

数位dp是用到数字位数性质来计数的动态规划问题。一般形式比较固定：通常问一个范围内满足某种性质的数的计数问题。我们通过做差可以将问题转化为前n个数中满足条件的数的个数。要解决这个问题，一般步骤分为：①初始化；②计算每个数位对答案的贡献；③累加得到答案。

分析思路

这里强烈推荐yxc的数位dp分析法。将数字n分解为一棵树，然后求每一部分对答案的贡献。在初始化的过程，我们一般会求得前X位数关于数字X的答案。我们定义一个函数f(n)来解决前n个数中满足条件的数的个数。首先特判n是否为0，不为0就按位记录每个数码。我们在计算每位数的贡献时，在[0,num[i]-1]范围枚举数码，这样做是为了保证统计到的是小于n的答案，然后判一下最后一个数是否满足条件即可。

例题

数字计数（模板题）

题意：输出[L,R]范围内每个数码的出现次数。

我们设f(i,j)表示i位数中j数码出现次数，就有以下递推公式
{f(i,j)=0,i=0f(i,j)=10×f(i−1,j)+10i−1,i>0 \begin{cases} f(i,j)=0,i=0\\ f(i,j)=10×f(i-1,j)+10^{i-1},i>0 \end{cases} {f(i,j)=0,i=0f(i,j)=10×f(i−1,j)+10i−1,i>0​
预处理这个出来之后就是常规的数位DP了，处理前导0另外减去10的幂次即可。

#include<bits/stdc++.h>
#define inf 0x3f3f3f3f
#define int long long
using namespace std;
const int N=30;

int pw[N],num[N],cnt[10],dp[N],tmp[10];

void init(){ 
	pw[0]=1;for(int i=1;i<N;i++) pw[i]=pw[i-1]*10;
	dp[0]=0;for(int i=1;i<N;i++) dp[i]=(dp[i-1]*10+pw[i-1]);
}

int solve(int n){
	int len=0,res=0,lst=n;
	while(n){
		num[++len]=n%10;
		n/=10;
	}
	for(int i=len;i>=1;i--){
		for(int j=0;j<=9;j++)	cnt[j]+=dp[i-1]*num[i];
		for(int j=0;j<num[i];j++)	cnt[j]+=pw[i-1];
		lst-=num[i]*pw[i-1];
		cnt[num[i]]+=lst+1;
		cnt[0]-=pw[i-1];
	}
	return res;
}

signed main(){
	ios::sync_with_stdio(0);
	cin.tie(0);cout.tie(0);
	init();
	int L,R;
	cin>>L>>R;
	solve(R);
	for(int i=0;i<=9;i++) tmp[i]=cnt[i];
	memset(cnt,0,sizeof(cnt));
	solve(L-1);
	for(int i=0;i<=9;i++) cout<<tmp[i]-cnt[i]<<" ";
	return 0;
}
数字游戏

题意：求[L,R]的不降数个数。定义不降数为从高位到低位的数字是不下降的。

//#pragma GCC optimize("Ofast", "inline", "-ffast-math")
//#pragma GCC target("avx,sse2,sse3,sse4,mmx")
#include<bits/stdc++.h>
#define inf 0x3f3f3f3f
#define int long long
using namespace std;
const int N=2e5+7;
const int mod=1e9+7;

//int read(){	int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-') f=f*-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}return x*f;}

int a,b,dp[N][40],num[40];

void init(){ //预处理dp[i][j]表示最高位是j，一共有i位的不降数个数 
	for(int i=0;i<=9;i++) dp[1][i]=1;
	for(int i=2;i<35;i++){
		for(int j=0;j<=9;j++){
			for(int k=j;k<=9;k++)
				dp[i][j]+=dp[i-1][k]; 
		}
	}
}

int Solve(int n){
	int len=0,res=0,lst=0;
	if(n==0) return 1;
	while(n){
		num[++len]=n%10;
		n/=10;
	}
	for(int i=len;i>=1;i--){  //从高位向低位枚举 
		for(int j=lst;j<num[i];j++) res+=dp[i][j]; //当前位要填比上一位大的数 
		if(num[i]<lst) break; 
		lst=num[i]; //表示上一位的数 
		if(i==1) res++; //最低位满足的答案 
	}
	return res;
}

signed main(){
	ios::sync_with_stdio(0);
	cin.tie(0);cout.tie(0);
	init();
	while(cin>>a>>b){
		cout<<Solve(b)-Solve(a-1)<<"\n";
	}
	return 0;
}
不要62

题意：求[L,R]中有多少个带有62或者带4的数。

这题数据范围出的巨水，暴力也可以过，数位DP的话多记录前面一个数码，剩下的就是模板了。

/*  这里是暴力的代码
//#pragma GCC optimize("Ofast", "inline", "-ffast-math")
//#pragma GCC target("avx,sse2,sse3,sse4,mmx")
#include<bits/stdc++.h>
using namespace std;
const int N=1e7+7;
const int mod=1e9+7;

int l,r,flag=0,tmp=0,sum[N];

signed main(){
	sum[0]=0;
	for(int i=1;i<=1e7;i++){
		tmp=i,flag=0;
		while(tmp){
			if(tmp%10==4) flag=1;
			else if(tmp>=62&&(tmp-62)%100==0) flag=1; 
			tmp/=10;
		}
		sum[i]=sum[i-1]+flag;
	}
	cin>>l>>r;
	while(l||r){
		cout<<(r-l+1)-(sum[r]-sum[l-1])<<"\n"; 
		cin>>l>>r;
	}
	return 0;
}*/

#include<bits/stdc++.h>
using namespace std;
const int N=10;

int L,R,dp[N][10],num[N];

void init(){
	for(int i=0;i<=9;i++) if(i!=4) dp[1][i]=1;
	for(int i=2;i<N;i++){
		for(int j=0;j<=9;j++){
			if(j==4) continue;
			for(int k=0;k<=9;k++){
				if(k==4||j==6&&k==2) continue;
				dp[i][j]+=dp[i-1][k];
			} 
		} 
	}
} 

int solve(int n){
	int len=0,res=0,lst=0;
	if(!n) return 1;
	while(n){
		num[++len]=n%10;
		n/=10;
	}
	for(int i=len;i>=1;i--){
		for(int j=0;j<num[i];j++){
			if(j==4||lst==6&&j==2) continue;
			res+=dp[i][j];
		}
		if(num[i]==4||lst==6&&num[i]==2) break;
		lst=num[i];
		if(i==1) res++;
	}
	return res; 
}

signed main(){
	init();
	while(cin>>L>>R,L||R){
		cout<<solve(R)-solve(L-1)<<"\n";
	}
}
Windy数

题意：求[L,R]有多少个满足相邻数位差的绝对值大于等于2的数字个数

//#pragma GCC optimize("Ofast", "inline", "-ffast-math")
//#pragma GCC target("avx,sse2,sse3,sse4,mmx")
#include<bits/stdc++.h>
#define inf 0x3f3f3f3f
#define int long long
using namespace std;
const int N=2e5+7;
const int mod=1e9+7;

//int read(){	int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-') f=f*-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}return x*f;}

int L,R,dp[40][15],num[40];

void init(){
	for(int i=0;i<=9;i++) dp[1][i]=1;
	for(int i=2;i<=35;i++){
		for(int j=0;j<=9;j++){
			for(int k=0;k<=j-2;k++) dp[i][j]+=dp[i-1][k];
			for(int k=j+2;k<=9;k++) dp[i][j]+=dp[i-1][k];			
		}
	}
}

int Solve(int n){
	int len=0,res=0,lst=-2;
	if(!n) return 0; 
	while(n){
		num[++len]=n%10;
		n/=10;
	}
	for(int i=len;i>=1;i--){
		if(i==len) for(int j=1;j<num[i];j++) res+=dp[i][j];
		else{
			for(int j=0;j<num[i];j++)
				if(abs(j-lst)>=2) res+=dp[i][j]; 
		}
		if(abs(num[i]-lst)>=2) lst=num[i];
		else break;
		if(i==1) res++;
	}
	for(int i=1;i<len;i++) for(int j=1;j<=9;j++) res+=dp[i][j]; //处理前导零的情况 
	return res;
}

signed main(){
	ios::sync_with_stdio(0);
	cin.tie(0);cout.tie(0);
//  freopen("in.cpp","r",stdin);
//  freopen("out.cpp","w",stdout);
	init();
	cin>>L>>R;
	cout<<Solve(R)-Solve(L-1)<<"\n"; 
	return 0;
}
数字游戏

题意：给定L,R,N，求[L,R]区间中有多少个数每位之和模N等于0。

多记录一维表示取模N的余数，并且记下前面k位的数字和。

#include<bits/stdc++.h>
using namespace std;
const int N=11,M=110;

int dp[N][10][M],L,R,P,num[N];

int mod(int x,int y){
	return (x%y+y)%y;
}

void init(){  //预处理dp[i][j][k]表示i位数最高位为j，取模为k的方案数 
	memset(dp,0,sizeof(dp));
	for(int i=0;i<=9;i++) dp[1][i][i%P]=1;
	for(int i=2;i<N;i++){
		for(int j=0;j<=9;j++){
			for(int k=0;k<P;k++){
				for(int x=0;x<=9;x++){
					dp[i][j][k]+=dp[i-1][x][mod(k-j,P)]; 
				}
			}
		}
	}
}

int solve(int n){
	int len=0,res=0,lst=0;
	if(!n) return 1;
	while(n){
		num[++len]=n%10;
		n/=10;
	}
	for(int i=len;i>=1;i--){
		for(int j=0;j<num[i];j++) res+=dp[i][j][mod(-lst,P)];	
		lst+=num[i];
		if(i==1&&lst%P==0) res++;
	}
	return res;
}

signed main(){
	while(cin>>L>>R>>P){
		init();
		cout<<solve(R)-solve(L-1)<<"\n";	
	}
}
恨7不成妻

题意：给定[L,R]，求区间内满足（如下条件的）数字的平方和。

  1. 整数中某一位是 7 ；
  2. 整数的每一位加起来的和是 7 的整数倍；
  3. 这个整数是 7 的整数倍。

这道题性质比较多，难度也较大。我们可以用记忆化搜索的写法，来处理dp[i][j][k]表示前i为每一位加起来模7为j,整个数模7为k的数字的是否存在，记录个数和平方。

#include<bits/stdc++.h>
//#define int long long
using namespace std;
typedef long long ll;
const int N=25,P=1e9+7;

struct E{
	ll s0,s1,s2;
	E(){s0=-1;s1=s2=0;}
	E(ll s0,ll s1,ll s2):s0(s0),s1(s1),s2(s2){}
}dp[N][10][10];

ll num[N],pw[N];

E _DP(int len,int r1,int r2,bool flag){
	E ans,nxt;
	if(!len) return (r1!=0&&r2!=0)?E(1,0,0):E(0,0,0);
	if(!flag&&dp[len][r1][r2].s0!=-1) return dp[len][r1][r2];
	int mxflag=(flag)?num[len]:9;
	ans.s0=0;
	for(int i=0;i<=mxflag;i++){
		if(i==7) continue;
		E nxt=_DP(len-1,(r1+i)%7,(r2*10+i)%7,flag&&i==mxflag);
		ans.s0+=nxt.s0;
		ans.s0%=P;
		ans.s1+=(nxt.s1+((pw[len]*i)%P)*nxt.s0%P)%P;
		ans.s1%=P;
		ans.s2+=(nxt.s2+((2*pw[len]*i)%P)*nxt.s1)%P;
		ans.s2%=P;
		ans.s2+=((nxt.s0*pw[len])%P*pw[len]%P*i*i%P);
		ans.s2%=P; 
	}
	if(!flag) dp[len][r1][r2]=ans;
	return ans;
}

ll solve(ll n){
	int len=0;
	while(n){
		num[++len]=n%10;
		n/=10;
	}
	return _DP(len,0,0,1).s2;
}

signed main(){
	ll t,l,r,ans;
	pw[1]=1;for(int i=2;i<=20;i++) pw[i]=(pw[i-1]*10)%P;
	cin>>t;
	while(t--){
		cin>>l>>r;
		ans=solve(r)-solve(l-1);
		cout<<(ans%P+P)%P<<"\n";
	}
	return 0;
}

参考文献

https://www.bilibili.com/video/BV1yT4y1u7jW

https://www.luogu.com.cn/blog/writeSTL/solution-p4999

https://blog.csdn.net/hzf0701/article/details/116717851

夜深人静写算法（二十九）- 数位DP
英雄哪里出来
02-08 9360
动态规划进阶，她来了！
数位DP学习整理（数位DP看完这篇你就会了）
热门推荐
pursuit的博客
05-12 4万+
文章目录数位DP数位DP介绍数位DP解法数位DP经典例题例题1：度的数量例题2：计数问题例题3：数字游戏例题4：windy数数位DP总结 数位DP 数位DP介绍 数位DP往往都是这样的题型，给定一个闭区间[l,r][l,r][l,r]，让你求这个区间中满足某种条件的数的总数。而这个区间可能很大，简单的暴力代码如下： int ans=0; for(int i=l;i<=r;i++){ if(check(i))ans++; } 我们发现，若区间长度超过1e81e81e8，我们暴力枚举就会超时了，
关于数位dp的几点总结 | 附经典例题与讲解
一只酷酷光儿的博客
08-12 680
把数位dp的简单题目做了一下,今天做一个总结与整理 其实所谓的数位dp做了几道题后发现就是一个模板类型问题,(指的是记忆化搜索强推状态转移就有点difficult了),我的经验下面几点: 1.首先确定上界与限制,就可以使得我们把当前这个数以下的所有数都进行遍历. 2.在过程中我们确定如何保留满足题意的答案,比如说 不要 62,我们需要记录前一位是否为6,如果前一位为6我们标记一下,如果下一位...
《算法竞赛进阶指南》0x50动态规划 - 数位DP
中学生信竞
06-12 336
《算法竞赛进阶指南》0x50动态规划 - 数位DP
位数（digit），计算一个整数是几位数（C++）
weixin_52038497的博客
03-14 1万+
计算一个整数是几位数（C++） 文章目录计算一个整数是几位数（C++）前言题目：解题思路：代码：运行结果总结 前言 刘汝佳算法《算法竞赛入门经典》第二章2.4.6上机练习第一题 题目： 输入一个不超过10的九次方的正整数，输出它的位数。例如：12735的位数是5.请不要使用任何数学函数，用四则运算和循环语句实现。 解题思路： 判断位数可以从第一位开始数也可以从最后一位开始数。从一个整数中获取其中一位数可以用除法或者取余来得到。而我们要求的就是位数，所以并不知道第一位是第几位，应该除以几去得到它。所以我
【数位dp】[Scoi2014] bzoj3598 方伯伯的商场之旅
每天心塞一点点
03-13 2606
题目点这里 和方伯伯的斗争终于结束了。。。我也是快要死了。。。。（请自动忽视那两道非人哉的题 = =） 之前听学长说这题很水的数位dp ：） 对 水到我都做不起了。  = =题解看了三遍啊！！！还是没研究明白这踏马是什么鬼！！ ……总有种我数位dp白学了的感觉（其实确实也是白学了）  思路……当然不是我想的………… 思路请见省选rank3大神：http
算法竞赛备赛进阶之数位DP训练
Williamtym的博客
01-16 3639
数位DP的思想就是对每一位进行DP，计算时记忆化每一位可以有的状态，其作用是减少运算时间，避免重复计算。数位DP是一种计数用的DP，一般就是要统计一个区间[A,B]内满足一些条件数的个数。以1e9甚至1e18、1e100的问题为例，因为在统计情况下有很多重复的计算，数位DP实现了相同状态只计算一次，从而大幅减少运算时间。数位DP：技巧1：[X, Y] => f(Y) - f(X-1)技巧2：用树进行排列。
dp 就 dp ，数位dp是什么意思 ？
仰望天空，妳我亦是行人.
03-28 2779
🐬初学一门技术时，总有些许的疑惑，别怕，它们是我们学习路上的点点繁星，帮助我们不断成长。🦄文章粗浅，希望对大家有帮助！
【算法竞赛】数位dp（超详细版，新手易懂）
最新发布
m0_74140776的博客
03-02 1169
数位dp是一个在竞赛中非常常用的思想。本文将从一道经典的竞赛题目出发，详解数位dp的思路及相关知识点的用法，非常适合初学者，简单易懂，知识点涵盖全面。
(数位dp) 算法竞赛入门到进阶 书本题集
cuber-lotus的博客
09-28 1933
hdu2089 不要62; hdu3555 Bomb; hdu3652 B-number; hdu6148 Valley Numer; hdu4507 吉哥系列故事——恨7不成妻; 本文中的题目来自罗勇军老师和郭卫斌老师的《算法竞赛入门到进阶》俗称小黑书 数位dp一般有两种实现方式，递推和记忆化dfs递推比较重视推到能力，而记忆化dfs也非常好写且模板话程度高。
【算法竞赛题目 & 题解收集】数位DP
m0_73500785的博客
01-20 960
ACM题目 & 题解 录
4.6.1 蓝桥杯动态规划之数位DP
tang7mj的博客
01-26 1651
数位DP是动态规划的一个重要分支，在处理数字的数位问题时非常有效。熟练掌握数位DP，对于编程竞赛和算法解题是非常有帮助的。希望本文能为你在学习数位DP的道路上提供一些指导和帮助。
数位dp入门题目总结——区间中满足要求的数的个数
启人zhr的博客
02-27 3577
引言 总体策略 不要62 Bomb B-number Balanced Number引言 在算法竞赛中，有一类求出给定区间中符合要求的数的个数问题，这类问题往往区间范围较大，无法通过枚举区间中数再判断条件这种方式来求解，数位dp就是一种解决这种方式的策略。 给出一篇写的很好地文章链接总体策略若区间符合可加减， 求解[l,r][l, r]满足条件的数个数可以通过[0,r]−[0,l−1][0,
数位DP讲解
zhangxian___的博客
07-18 1万+
转载自：http://www.cnblogs.com/itlqs/p/5935308.html 数位DP其实是很灵活的，所以一定不要奢求一篇文章就会遍所有数位DP的题，这一篇只能是讲清楚一种情况，其他情况遇到再总结，在不断总结中慢慢体会这个思想，以后说不定就能达到一看到题目就能灵活运用的水平。（其实DP都是这样……）   这一篇要说的数位DP是一道最简单的数位DP：htt
数位dp。
WQhuanm的博客
04-29 1301
在处理1e9甚至1e18,1e100的问题时，因为在统计情况下有很多重复的计算，数位dp实现了相同状态只计算一次，从而大幅减少运算时间，思想就是对每一位进行dp，计算时记忆化每一位可以有的状态。如我们在统计1234的状态时，可以拆成统计0~10000，0~2000，0~300，0~40数位统计我们用bit数组由低到高存储每一位，bit[1]=4,bit[2]=3,bit[3]=2,bit[4]=1.然后dp从高位到低位进行。
数位dp学习笔记
Tudou_Pika的博客
10-23 507
将一个数字按照数位拆开，接着关注每一位数字。数位dp通常可以用计数问题的技巧如[l,r]=[0,r]−[0,l−1]等。数位中，无论是1~3位、10~12位，还是1000~1002位，性质在大多数题目中都相同。所以通常将一些数据处理在一个通用的数组内，然后利用这个数组加加减减乘乘除除求解答案。求解答案的过程可以采用记忆化搜索，或循环dp，具体视具体题目的实现难度而定。接下来看几道题目。
数位DP 详解
weixin_30302609的博客
09-29 965
序 天堂在左，战士向右 引言 数位DP在竞赛中的出现几率极低，但是如果不会数位DP，一旦考到就只能暴力骗分。 以下是数位DP详解，涉及到的例题有： [HDU2089]不要62 [HDU3652]B-number 概述 首先我们要理清的是，到底数位DP是什么。 事实上，一般数位DP的题目题面描述都会有以下内容： 求出一段区间\([l,r]\)中，满足某一特殊条件的数有多少个 在例题...
数位DP——学习笔记
流水
02-16 256
题目特征 数位DP题目一般都是求[l,r][l,r][l,r]间满足某种情况的数有几个, 一般l,rl,rl,r都比较大 或者直接就是算[1,N][1,N][1,N]的 题目解法 数位DP虽然是DP，但是大部分题都能用一个相同的方法来解决 虽然很多题目能用递推写，但是递推难度较记忆化搜索较高，因此最好选用记忆化搜索来写题 模板 ll dfs(int pos, ... , bool lead, bo...
算法：数位DP
GoldenSpace
10-12 418
数位DP 问题 求区间[L,R][L, R][L,R]中满足条件的数有多少个， 0≤L≤R0 \le L \le R0≤L≤R，该条件与数位有关，比如不包含数字444。 思路 时间复杂度 O(log⁡n)O(\log n)O(logn)。实际上为nnn的位数，而nnn的位数可以表示为lg⁡n\lg nlgn，即log⁡nlog⁡10\frac {\log n}{\log 10}log10logn​...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

RWLinno

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值