「一本通」数位dp学习笔记

最新推荐文章于 2020-11-05 07:35:33 发布
最新推荐文章于 2020-11-05 07:35:33 发布
阅读量458 收藏 3
点赞数
分类专栏： loj DP 数位dp 一本通提高篇 文章标签： loj 一本通 数位dp
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_36038511/article/details/82794261
版权
DP 同时被 3 个专栏收录
20 篇文章 0 订阅
订阅专栏
loj
8 篇文章 0 订阅
订阅专栏
一本通提高篇
6 篇文章 1 订阅
订阅专栏
废话：

从数位dp只会windy数变成什么都不会 从只会写递推变成只会写记忆化搜索…

总结：

大概就是解决一些对数位有要求而且上下限特别的大…一般有两种实现方法，递推（dp，比较好理解一般都是先学这一种）/记忆化搜索（暴力，方便，容易写，sb方法）

板子：

https://www.luogu.org/blog/mak2333/shuwei

loj#10163. 「一本通 5.3 例 1」Amount of Degrees

https://loj.ac/problem/10163

几个次幂相加，等于转到b进制时那一位为1。转成多叉树的形式对每一棵含有答案的子树乱搞。注意对应的次幂
分析&论文：https://wenku.baidu.com/view/d2414ffe04a1b0717fd5dda8.html

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int f[32][32],a[32];//f:有i层（有i位）有j个1的个数 （整数次幂相加 对应的b进制数位上应为1）
int k,b;
int solve(int x)
{
	int len=0,kk=k;
	while (x!=0)
	{
		a[len++]=x%b; //注意对应次幂 
		x/=b;
	}
	int ans=0; len--;
	for (int i=len;i>=0;i--)
	{
		if (a[i]==1)
		{
			ans+=f[i][kk];
			kk--;
		}
		else if (a[i]>1)
		{
			ans+=f[i+1][kk];
			break;
		}
		if (kk<0) return ans;
	}
	if (kk==0) ans++; //x刚好也是k个1
	return ans;
}
int main()
{
	for (int i=0;i<31;i++)
	{
		f[i][0]=1;
		for (int j=1;j<=i;j++) f[i][j]=f[i-1][j]+f[i-1][j-1];
	}
	int x,y;
	scanf("%d%d%d%d",&x,&y,&k,&b);
	printf("%d\n",solve(y)-solve(x-1));
	return 0;
}
loj#10164. 「一本通 5.3 例 2」数字游戏

https://loj.ac/problem/10164

递推做法：
//递推 
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int a[32],f[32][32];//f[i位][开始的数字]（递推用开始的数  记忆化搜索用最后一位数？）
int solve(int x)
{
	int len=0; memset(a,0,sizeof(a));
	while (x)
	{
		a[++len]=x%10;
		x/=10;
	}
	int ans=0;
	for (int i=len;i;i--)
	{
		if (a[i+1]>a[i]) break;
		for (int j=a[i+1];j<a[i];j++)
		{
			ans+=f[i+1-1][j];
		}
		if (i==1) ans++;//本身也是 
	}
	return ans;
}
int main()
{
	for (int i=1;i<=31;i++) f[1][i]=1;
	for (int i=2;i<=31;i++)
		for (int j=0;j<=9;j++)
			for (int k=j;k<=9;k++)
				f[i][j]+=f[i-1][k];
	int x,y;
	while (scanf("%d%d",&x,&y)!=EOF)
	{
		printf("%d\n",solve(y)-solve(x-1));
	}
	return 0;
}
记忆化搜索：
//记忆化搜索
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int f[32][32],a[32];//f[第i高位][填的数字]   
int dfs(int len,int st,int lim)
{
	if (len==0) return 1;
	if (lim==0&&f[len][st]!=-1) return f[len][st];
	int ans=0,ed; if (lim==1) ed=a[len]; else ed=9;
	for (int i=st;i<=ed;i++)
	{
		if (lim==1&&i==ed) ans+=dfs(len-1,i,1);
		else ans+=dfs(len-1,i,0);
	}
	return ans;
}
int solve(int x)
{
	memset(f,-1,sizeof(f));
	int len=0;
	while (x)
	{
		a[++len]=x%10;
		x/=10;
	}
	return dfs(len,0,1);
}
int main()
{
	int x,y;
	while (scanf("%d%d",&x,&y)!=EOF)
	{
		printf("%d\n",solve(y)-solve(x-1));
	}
	return 0;
}
loj#10165. 「一本通 5.3 例 3」Windy 数

https://loj.ac/problem/10165

递推做法：
//递推
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
int f[32][32],a[32]; //f[位数i][第i位(最高位)的数字] 
void init()
{
	for (int i=0;i<=9;i++) f[1][i]=1;
	for (int i=2;i<=31;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 solve(int x)
{
	memset(a,0,sizeof(a));
	int len=0;
	while (x)
	{
		a[++len]=x%10;
		x/=10;
	}
	int ans=0;
	for (int i=1;i<len;i++)//不足len位的部分 
		for (int j=1;j<=9;j++) 
			ans+=f[i][j];
	for (int i=1;i<a[len];i++) //第len位不足a[len]的部分
	{
		ans+=f[len][i];
	}
	for (int i=len-1;i;i--) 
	{
		for (int j=0;j<a[i];j++) if (abs(j-a[i+1])>=2) ans+=f[i][j]; //跟前一位比较 
		if (abs(a[i+1]-a[i])<2) break;
		if (i==1) ans++;
	}
	return ans;
}
int main()
{
	init();
	int a,b;
	scanf("%d%d",&a,&b);
	printf("%d\n",solve(b)-solve(a-1));
	return 0;
}
记忆化搜索做法：
//记忆化搜索
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int f[10][10],a[10];//f[第i高位][填的数字]
int dfs(int len,int st,int ze,int lim)//ze是否有前导零  lim:是否达到上限 
{
	if (len==0) return 1;
	if (lim==0&&f[len][st]!=-1) return f[len][st];
	int ans=0,ed; 
	if (lim==1) ed=a[len]; else ed=9; 
	for (int i=0;i<=ed;i++) 
	{ 
		if (ze==1)
		{
			int zz=0;
			if (i==0) zz=1; 
			if (lim==1&&i==ed) ans+=dfs(len-1,i,zz,1);
			else ans+=dfs(len-1,i,zz,0);
		}
		else if (abs(st-i)>=2)
		{
			if (lim==1&&i==ed) ans+=dfs(len-1,i,0,1);
			else ans+=dfs(len-1,i,0,0);
		}
	}
	if (lim==0&&st!=0) f[len][st]=ans;
	return ans;
}
int solve(int x)
{
	memset(a,0,sizeof(a));
	int len=0;
	while (x)
	{
		a[++len]=x%10;
		x/=10;
	}
	return dfs(len,0,1,1);
}
int main()
{
	memset(f,-1,sizeof(f));
	int a,b;
	scanf("%d%d",&a,&b);
	printf("%d\n",solve(b)-solve(a-1));
	return 0;
}
loj#10166. 「一本通 5.3 练习 1」数字游戏

https://loj.ac/problem/10166
不要想太多，直接爆搜就好

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int f[32][110],a[32],n;//f[位数][余数]
int dfs(int len,int mod,int lim)
{
	if (len==0) 
	{
		if (mod==0) return 1; else return 0;
	}
	if (f[len][mod]!=-1&&lim==0) return f[len][mod];
	int ans=0,ed;
	if (lim==1) ed=a[len]; else ed=9;
	for (int i=0;i<=ed;i++)
	{
		if (lim==1&&i==ed) ans+=dfs(len-1,(mod+i)%n,1);
		else ans+=dfs(len-1,(mod+i)%n,0);
	}
	if (lim==0) f[len][mod]=ans;
	return ans;
}
int solve(int x)
{
	memset(f,-1,sizeof(f));
	memset(a,0,sizeof(a));
	int len=0;
	while (x)
	{
		a[++len]=x%10;
		x/=10;
	}
	return dfs(len,0,1);
}
int main()
{
	int a,b;
	while (scanf("%d%d%d",&a,&b,&n)!=EOF)
		printf("%d\n",solve(b)-solve(a-1));
	return 0;
}
loj#10167. 「一本通 5.3 练习 2」不要 62

https://loj.ac/problem/10167

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int f[7][10],a[7];//f[位数i][第i位数字]
int dfs(int len,int last,int lim)
{
	if (len==0) return 1;
	if (lim==0&&f[len][last]!=-1) return f[len][last];
	int ans=0,ed;
	if (lim==1) ed=a[len]; else ed=9;
	for (int i=0;i<=ed;i++)
	{
		if (i==4) continue;
		if (last==6&&i==2) continue;
		if (lim==1&&i==ed) ans+=dfs(len-1,i,1);
		else ans+=dfs(len-1,i,0);
	}
	if (lim==0) f[len][last]=ans;
	return ans;
}
int solve(int x)
{
	memset(f,-1,sizeof(f));
	memset(a,0,sizeof(a));
	int len=0;
	while (x)
	{
		a[++len]=x%10;
		x/=10;
	}
	return dfs(len,0,1);
}
int main()
{
	int a,b;
	while (scanf("%d%d",&a,&b)!=EOF&&a&&b)
	{
		printf("%d\n",solve(b)-solve(a-1));
	}
	return 0;
}
loj#10168. 「一本通 5.3 练习 3」恨 7 不成妻

https://loj.ac/problem/10168
完全平方公式2333（没想到 装傻.jpg）
题解&分析：https://blog.csdn.net/ten_three/article/details/19698055
&https://www.cnblogs.com/kuangbin/archive/2013/05/01/3053233.html

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define LL long long
#define mod 1000000007
struct node
{
	LL s,sum,sqrsum;//个数N 各数之和 各数平方和 
	node()
	{
		s=sum=sqrsum=-1;
	}
}f[20][10][10];//f[位数][各数位之和mod7][这个数mod7] 
int a[20];
LL pi[20];
node dfs(int len,int s,int qs,int lim)//数位 各数位之和mod7 这个数mod7 是否上限 
{
	if (len==0)
	{
		node tt;
		tt.s=tt.sum=tt.sqrsum=0;
		if (s!=0&&qs!=0) tt.s=1;
		return tt;
	}
	if (lim==0&&f[len][s][qs].sqrsum!=-1) return f[len][s][qs];
	node ans; ans.s=ans.sum=ans.sqrsum=0;
	int ed=0;
	if (lim==1) ed=a[len]; else ed=9;
	for (int i=0;i<=ed;i++)
	{
		if (i==7) continue;
		/*
		第len位填i f1,f2,f2...fN记录前len位已经填好的数字
		对答案的贡献为：  (f1+ i*10^(len-1) )^2+(f2+ i*10^(len-1) )^2+...+(fN+ i*10^(len-1) )^2
		将平方和展开 整理 提公因式:
		N*(i*10^(len-1))^2+ f1^2+f2^2+...+fN^2+ 2*(i*10^(len-1))*(f1+f2+...fN) 
		其中 N即为要求的s   2*(i*10^(len-1))*(f1+f2+...fN) 即为要求的sum   整一个即为sqrsum 
		*/
		node tt;
		if (lim==1&&i==ed) tt=dfs(len-1,(s+i)%7,(qs*10+i)%7,1);
		else tt=dfs(len-1,(s+i)%7,(qs*10+i)%7,0);
		LL si=(i*pi[len-1])%mod;
		ans.s=(ans.s+tt.s)%mod; 
		ans.sum=(ans.sum+(tt.sum+ (si*tt.s)%mod ) %mod )%mod;
		ans.sqrsum=(ans.sqrsum+tt.sqrsum+( ((si*si)%mod*tt.s)%mod+ ((2*si)%mod *tt.sum) %mod )%mod )%mod;
	}
	if (lim==0) f[len][s][qs]=ans;
	return ans;
}
LL solve(LL x)
{
	memset(a,0,sizeof(a));
	int len=0;
	while (x)
	{
		a[++len]=x%10;
		x/=10;
	}
	node ans=dfs(len,0,0,1);
	return ans.sqrsum%mod;
}
int main()
{
	pi[0]=1;
	for (int i=1;i<=18;i++) pi[i]=pi[i-1]*10;
	int t;
	scanf("%d",&t);
	while (t--)
	{
		LL x,y;
		scanf("%lld%lld",&x,&y);
		if (x>y) swap(x,y);
		printf("%lld\n",(solve(y)-solve(x-1)+mod)%mod);
	}
	return 0;
}
loj#10169. 「一本通 5.3 练习 4」数字计数

https://loj.ac/problem/10169
别人都是三维f[数位][ ][前导零]+做10次，每次找一个数
只有我是f[数位][填的数字]鬼畜结构体一次找完所有数字且不带前导零……
所以注意前导零！！！具体看代码里解释的！
不加前导零那一维就要记得预处理

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
struct node
{
	bool bk;
	long long a[10];
	node()
	{
		bk=0; memset(a,-1,sizeof(a));
	}
}f[15][10];
int a[15];
long long pi[15],c[15];
node jia(node x,node y)
{
	for (int i=0;i<=9;i++) x.a[i]+=y.a[i];
	return x;
}
node dfs(int len,int x,int ze,int lim) //ze：zero是否有前导零
{
	if (len==0) {node tt;tt.bk=1;memset(tt.a,0,sizeof(tt.a));tt.a[x]=1;return tt;}
	if (ze==0&&lim==0&&f[len][x].bk==1) return f[len][x];
	node ans; ans.bk=1; memset(ans.a,0,sizeof(ans.a));
	int ed;
	if (lim==1) ed=a[len]; else ed=9;
	for (int i=0;i<=ed;i++)
	{
		node tt;
		tt=dfs(len-1,i,i==0&&ze==1,i==ed&&lim==1);
		/*if (ze==1&&i==0) ;
		else if (len>1) tt.a[i]+=tt.t;*/
		if (ze==1&&i==0) ;
		else
		{
			if (len>1) 
			{
				if (lim==1&&i==ed) tt.a[i]+=c[len-1]+1;
				else tt.a[i]+=pi[len-1];
			}
		}
		ans=jia(ans,tt);
	}
	if (ze==0&&lim==0) f[len][x]=ans;
	//关于前导零不能用/存答案 例1010 当len=3,往下找len=2的情况时 由于此时f[2][ ].a[0]=9 (10,20,30,...,90) 而答案应为18 (10,20,...,90; 01,02,...,09) 
	return ans;
}
node solve(long long x)
{
	memset(a,0,sizeof(a));
	memset(c,0,sizeof(c));
	/*for (int i=0;i<15;i++)
		for (int j=0;j<10;j++) memset(f[i][j].a,-1,sizeof(f[i][j].a)),f[i][j].bk=0;*/
	int len=0; 
	while (x)
	{
		a[++len]=x%10;
		x/=10;
		c[len]=a[len]*pi[len-1]+c[len-1]; //预处理最多有多少个数在这里面（当前枚举位出现次数）
	}
	return dfs(len,0,1,1);
}
int main()
{
	pi[0]=1;
	for (int i=1;i<15;i++) pi[i]=pi[i-1]*10;
 	long long x,y;
	scanf("%lld%lld",&x,&y);
	node r=solve(y),l=solve(x-1); 
	for (int i=0;i<=9;i++)
	{
		printf("%lld ",r.a[i]-l.a[i]);
	}
	printf("\n");
	return 0;
}

完结撒花！
（sd记忆化搜索真好用）

kkkGIGi_qtt
关注
  • 0
    点赞
  • 踩
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
数位DP小结_记忆化搜索版
pikachu_777的专栏
07-28 902
之前写数位dp是用dp方程递推形式写，接触记忆化搜索版后，发现使用dfs更加清晰明了，思路很顺，代码简洁的多W。 首先上模板，以HDU 2089 不要62为例     #define _CRT_SECURE_NO_DEPRECATE #include &lt;stdio.h&gt; #include &lt;string.h&gt; #include &lt;stdlib.h&gt; #...
数位dp(一)——Loj #10166. 「一本通 5.3 练习 1」数字游戏
hhhcbw的博客，欢迎各位来访
08-21 217
题目链接：https://loj.ac/problem/10166 题目大意 求区间[a,b]有多少数字满足数位加和%N等于0。 解题思路 看范围，很明显的数位dp。我们设定状态dp[pos][N][i]表示在%N的情况下，到第pos位，数位加和%N是i，会有多少数字满足数位加和%N等于0。 之后就是很easy的套模板环节。 AC代码 #include <iostream> #include <stdio.h> #include <cstring> #include &l
在谈数位dp_用记忆优化解决无法递推的题目
ACM之梦
12-05 645
有些数位dp题目真心难用递推，递推貌似很难想，不过搜索相对来讲是比较好设计的，现在可以考虑用搜索，如何搜索？，记忆优化肯定必须的，像dp一样设置状态。 两个例题解答如何用记忆优化解决此类问题。 hdu 3709 平衡数 故名思意，这题要求以某个点为中心点其他各个位到这个为力矩和为0的数。 不多说，举个例子 4139 以3为中心，那么左边2*4+1*1 。右边：1*9 刚好左右相等。如何设计
P1070 道路游戏 暴力和dp【1】
lattice的博客
01-24 229
P1070 道路游戏 https://www.luogu.org/problemnew/show/P1070 题目描述 小新正在玩一个简单的电脑游戏。 游戏中有一条环形马路，马路上有 n个机器人工厂，两个相邻机器人工厂之间由一小段马路连接。小新以某个机器人工厂为起点，按顺时针顺序依次将这 n个机器人工厂编号为1−n，因为马路是环形的，所以第n个机器人工厂和第1个机器人工厂是由一段马路连接...
数位dp(四)——#10169. 「一本通 5.3 练习 4」数字计数
hhhcbw的博客，欢迎各位来访
08-22 205
题目链接：https://loj.ac/problem/10169 题目大意 给定两个正整数 a 和 b，求在[a,b]中的所有整数中，每个数码 (digit) 各出现了多少次。 解题思路 数位dp，维护每个数字出现次数和数字个数。设当前状态为ans，子状态为tmp。状态转移如下。 注意有前导零时，不能把之前的零算进去，所以要特判一下。 AC代码 #include <iostream> #include <stdio.h> #include <cstring> #inc
数位dp算法详细图文解答
02-13
PPT图片结合文字讲述数位dp
数位dp——某次队内分享任务
07-21
数位动态规划(Digit Dynamic Programming, 简称数位DP)是一种在处理与数字相关的动态规划问题时，通过分析每个数字位上的状态来优化算法的技术。它通常用于解决涉及数字序列特性的数学问题，例如计数特定类型的数...
数位DP.pdf
07-12
数位DP整理的模板和一些题目的题解。QAQ没钱下载资源了只好上传点东西。
初探数位dp.pdf
09-16
lazycal的集训队报告：初探数位DP 以HDU 2089，HDU 3652, URAL 1057等题目为例，介绍了数位DP的算法
数位DP入门
08-22
数位DP入门思想，包括详细的思考过程以及一些习题的讲解
#10166. 「一本通 5.3 练习 1」数字游戏
TJ
09-09 766
【题目描述】 由于科协里最近真的很流行数字游戏，某人又命名了一种取模数，这种数字必须满足各位数字之和 modN 为 0。现在大家又要玩游戏了，指定一个整数闭区间 [a,b]，问这个区间内有多少个取模数。 【输入格式】 题目有多组测试数据。每组只含三个数字 a,b,N。 【输出格式】 对于每个测试数据输出一行，表示各位数字和 modN 为 0 的数的个数。 【样例输入】 1 19 9...
「一本通 5.3 练习 2」不要 62 --数位DP
小元勋的博客
08-20 183
loj 10167 题目分析： 数位DP就是在爆搜剪枝的基础上加了个**记忆化* Code： #include <bits/stdc++.h> using namespace std; #define LL long long #define re register int n,m,a[15]; LL f[15][2][2]; LL dfs_(int pos,int six...
#10165. 「一本通 5.3 例 3」Windy 数
zsjzliziyang的博客
11-05 169
题目描述 原题来自：SCOI 2009 Windy 定义了一种 Windy 数：不含前导零且相邻两个数字之差至少为2 的正整数被称为 Windy 数。 Windy 想知道，在A 和 B 之间，包括 A 和 B，总共有多少个 Windy 数？ 输入格式 一行两个数，分别为。 输出格式 输出一个整数，表示答案。 样例 样例输入 1 1 10 样例输出 1 9 样例输入 2 25 50 样例输出 2 20 数据范围与提示 20% 的数据，满足 ； 10...
#10163. 「一本通 5.3 例 1」Amount of Degrees
TJ
09-07 687
【题目描述】 原题来自：NEERC 2000 Central Subregional，题面详见 Ural 1057。题目 求给定区间 [X,Y] 中满足下列条件的整数个数：这个数恰好等于 K个互不相等的 B的整数次幂之和。例如，设  X=15, Y=20, K=2, B=2 则有且仅有下列三个数满足题意： 17=2^4+2^0 18=2^4+2^1 20=2^4+2^2 【输...
【题解】吉哥系列故事——恨7不成妻
ezoixx130's Blog
07-14 2889
Description 单身! 依然单身！ 吉哥依然单身！ DS级码农吉哥依然单身！ 所以，他生平最恨情人节，不管是214还是77，他都讨厌！ 吉哥观察了214和77这两个数，发现： 2+1+4=7 7+7=7*2 77=7*1     最终，他发现原来这一切归根到底都是因为和7有关！所以，他现在甚至讨厌一切和7有关的数！ 什么样的数和7有关呢？ 如果一个整数符合下面3
loj#10168恨 7 不成妻
quantum11
08-28 291
这一题的题目描述很…无法描述，但还是能马上看出是一道数位dp题。限制的条件有点多，细节处理起来确实很烦。只要记录合法的数的个数，数的和以及数的平方和。 #include&lt;bits/stdc++.h&gt; #define ll long long #define il inline using namespace std; const int mo=1e9+7; #define getch...
LOJ #10165. 「一本通 5.3 例 3」Windy 数
dingjinhui4074的博客
07-25 141
题目描述 原题来自：SCOI 2009 Windy 定义了一种 Windy 数：不含前导零且相邻两个数字之差至少为2的正整数被称为 Windy 数。 Windy 想知道，在A和B之间，包括A和B，总共有多少个 Windy 数？ 输入格式 一行两个数，分别为A,B。 输出格式 输出一个整数，表示答案。 ...
HDU4507吉哥系列故事——恨7不成妻(数位dp)
Rain
01-26 508
题目链接：点击打开链接 题意： 求区间[l,r]内所有与7无关的数的平方和（取模） 定义与7无关的数：                                       1.数字的数位上不能有7                                      2.数字的数位和不能是7的倍数                                      3
数位dp模板
最新发布
05-27
数位 DP 是一类动态规划算法，用于计算某个数的某种性质。通常情况下，数位 DP 需要将数分解成某种形式，然后逐位进行计算。下面是一个简单的数位 DP 模板： ```python def digit_dp(pos, sum, limit): if pos == ...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值