动态规划——数位

已于 2023-12-24 19:23:24 修改
阅读量989 收藏 14
点赞数 26
分类专栏： ACM 文章标签： 动态规划 代理模式 算法
于 2023-12-24 19:19:55 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_71819874/article/details/135185254
版权

动态规划——数位DP

”某一区间“满足某种性质

技巧1：[x,y]=f[y]-f[x-1];

技巧2：树;

模板如下：

//框架模板数位dp
#include<bits/stdc++.h>
using namespace std;
int a,b;
const int N=35;
int f[N][10];  //f[i][j]:长度为i，以j开头的满足条件的数的个数
void init()//解决实质方法 初始化f
{

}

int dp(int x)//解决区间问题 
{
	if(x==0)return 0;
	vector<int> nums;
	while(x)nums.push_back(x%10),x=x/10;
	
	int res=0;
	int last=0;//用来记录前面有关系的一个性质 
	
	for(int i=nums.size()-1;i>=0;i--)
	{
		int x=nums[i];
        //先处理左边分支
        
        //继续讨论处理右边分支，条件判断，更新last，最后i=0时候res更新
		
	}
	
	return res;
	
}
int main()
{
	init();
	cin>>a>>b; 
	cout<<dp(b)-dp(a-1)<<endl;
	return 0;
 }

例1：度的数量

在这里插入图片描述
分析：
在这里插入图片描述

#include<bits/stdc++.h>
using namespace std;
int x,y,k,b;
const int N=35;
int f[N][N];
void init()//组合数——杨辉三角形
{
	for(int i=0;i<N;i++)
	{
		for(int j=0;j<=i;j++)
		{
			if(!j)f[i][j]=1;
			else f[i][j]=f[i-1][j-1]+f[i-1][j];
		}
	}
}
int dp(int x)
{
	if(!x)return 0;
	
	vector<int> nums;
	while(x)nums.push_back(x%b),x=x/b;
	
	int res=0;
	int last=0;//记录前面的位数占用的1
	
	for(int i=nums.size()-1;i>=0;i--)
	{
		int x=nums[i];
		if(x)//为了防止超出区间只有x不为0才有选择 
		{
			res+=f[i][k-last];//如果当前位为0 
			if(x>1)//如果当前位选择为1
			{
				res+=f[i][k-last-1];
				break; 
			} 
			else//当前本来就为1，就需要继续讨论 
			{
				last++;
				if(last>k)break; 
			}
		}
		if(i==0&&last==k)res++;//最右侧分支上的1 
	 } 
	 return res; 
}
int main()
{
	cin>>x>>y>>k>>b;
	init();
	cout<<dp(y)-dp(x-1)<<endl;
	return 0;
}

例2、数字游戏

在这里插入图片描述
分析：

在这里插入图片描述

#include<bits/stdc++.h>
using namespace std;
int a,b;
const int N=35;
int f[N][N];//f[i][j]为长度为i，且以j开头的的有多少个不降数 
void init()//动态规划预处理
{
	for(int i=0;i<=9;i++)f[1][i]=1;
	for(int i=2;i<=N;i++) 
	{
		for(int j=0;j<=9;j++)
		{
			for(int k=j;k<=9;k++)
			{
				f[i][j]+=f[i-1][k];
//				cout<<f[i][j]<<endl;
			}
		}
	}
}

int dp(int x)
{
	if(x==0)return 1;
	vector<int> nums;
	while(x)nums.push_back(x%10),x=x/10;
	
	int res=0;
	int last=0;//用来记录上一位数字 
	
	for(int i=nums.size()-1;i>=0;i--)
	{
		int y=nums[i];
		for(int j=last;j<y;j++)//左边分支的处理
		{
			res+=f[i+1][j]; 
		} 
		if(y<last)break;
		last=y;//右边分支的处理
		if(i==0)res++;//右边最后分支的处理
	}
	return res;
	
}
int main()
{e
	init();
	while(cin>>ea>>b)
	{
		cout<<dp(b)-dp(a-1)<<endl;
	}
	return 0;
 }

例3、 Windy数

在这里插入图片描述

在这里插入图片描述

如果要求dp（1567），由于第一个根节点的分叉左枝干是从1开始，那么就等于0-999都没有算进res里面

好好理解f的含义，例如f[3] [9]表示最高位为9长度为3 满足条件的windy数的个数，那么是9_ _,900–999中的个数

#include<bits/stdc++.h>
using namespace std;
int a,b;
const int N=35;
int f[N][10];   
void init()//解决实质方法 
{
	for(int i=0;i<=9;i++)f[1][i]=1;
	for(int i=2;i<=N;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 dp(int x)//解决区间问题 
{
	if(x==0)return 0;
	vector<int> nums;
	while(x)nums.push_back(x%10),x=x/10;
	
	int res=0;
	int last=-2;//用来记录前面有关系的一个性质 
	
	for(int i=nums.size()-1;i>=0;i--)
	{
		int y=nums[i];
		for(int j=i==nums.size()-1;j<y;j++)
		{
			if(abs(j-last)>=2)
			res+=f[i+1][j];
		}
		if(abs(y-last)>=2)last=y;
		else break;
		if(!i)res++;
		
	}
	//特殊处理000——的情况 ，也就是一开始少算了前导0的部分区间
	for(int i=1;i<nums.size();i++)
	{
		for(int j=1;j<=9;j++)
		res+=f[i][j];
	}

	return res;
	
}
int main()
{
	init();
	cin>>a>>b; 
	cout<<dp(b)-dp(a-1)<<endl;
	return 0;
 }

4、数字游戏

在这里插入图片描述

#include<bits/stdc++.h>
using namespace std;
int a,b;
const int N=35;
int M;
int f[N][10][100];  //数长为i，以j为开头，每个数字之和模为t 
void init()//解决实质方法 
{
	memset(f,0,sizeof(f));
	for(int i=0;i<=9;i++)f[1][i][i%M]++;
	
	for(int i=2;i<N;i++)
	{
		for(int j=0;j<=9;j++)
		{
			for(int t=0;t<=M-1;t++)
			{
				for(int k=0;k<=9;k++)
				{
					f[i][j][t]+=f[i-1][k][((t-j)%M+M)%M];
				//	cout<<f[i][j][t]<<endl;
				}
			}
		}
	}
}

int dp(int x)//解决区间问题 
{
	if(x==0)return 1;
	vector<int> nums;
	while(x)nums.push_back(x%10),x=x/10;
	
	int res=0;
	int last=0;//用来记录前面有关系的一个性质 
	
	for(int i=nums.size()-1;i>=0;i--)
	{
		int y=nums[i];
		for(int j=0;j<y;j++)//处理左边的分支 
		{
			res+=f[i+1][j][(M+-last)%M];
		}
		last=(last+y)%M;
		if(i==0&&last==0)res++;
		
	}

	return res;
	
}
int main()
{
	while(cin>>a>>b>>M)
	{
	
	init();
	cout<<dp(b)-dp(a-1)<<endl;
}
	return 0;
 }

5、 不要62

在这里插入图片描述

有些该有的条件不可省略

#include<bits/stdc++.h>
using namespace std;
int a,b;
const int N=10;
int f[N][10];  //数长为i，以j为开头
//f[i][j]+=f[i-1][k],如果j为2,k就不会为6 
void init()//解决实质方法 
{
	for(int i=0;i<=9;i++)f[1][i]=1;
	f[1][4]=0;
	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;
				f[i][j]+=f[i-1][k];
			}
		}
	 } 
}

int dp(int x)//解决区间问题 
{
	if(x==0)return 1;
	vector<int> nums;
	while(x)nums.push_back(x%10),x=x/10;
	
	int res=0;
	int last=0;//用来记录前面有关系的一个性质 
	
	for(int i=nums.size()-1;i>=0;i--)
	{
		int y=nums[i];
		for(int j=0;j<y;j++)
		{
			if(last==6&&j==2)continue;
			res+=f[i+1][j];
		}
		if(last==6&&y==2)break;
		if(y==4)break;
		last=y;
		if(!i)res++;
	} 
	return res;
	
}
int main()
{
	init();
	while(~scanf("%d %d",&a,&b),a&&b)
	{ 
		cout<<dp(b)-dp(a-1)<<endl;
	}
	return 0;
 }

6、

[USACO06NOV] Round Numbers S

如果一个正整数的二进制表示中， 0 0 0 的数目不小于 1 1 1 的数目，那么它就被称为「圆数」。

例如， 9 9 9 的二进制表示为 1001 1001 1001，其中有 2 2 2 个 0 0 0 与 2 2 2 个 1 1 1。因此， 9 9 9 是一个「圆数」。

请你计算，区间 [ l , r ] [l,r] [l,r] 中有多少个「圆数」。

一行，两个整数 l , r l,r l,r。

一行，一个整数，表示区间 [ l , r ] [l,r] [l,r] 中「圆数」的个数。

2 12

6

【数据范围】

对于 100 % 100\% 100% 的数据， 1 ≤ l , r ≤ 2 × 1 0 9 1\le l,r\le 2\times 10^9 1≤l,r≤2×109。

【样例说明】

区间 [ 2 , 12 ] [2,12] [2,12] 中共有 6 6 6 个「圆数」，分别为 2 , 4 , 8 , 9 , 10 , 12 2,4,8,9,10,12 2,4,8,9,10,12。

#include<bits/stdc++.h>
using namespace std;
int a,b;
const int N=50;
int f[N][2][N];//数长为i,圆数的个数
//f[i][0][j]+=f[i][0][j-1] 
//f[i][j][p][t]=f[i-1][k1][p-j][t-t*i0]
//第三维为数字之和的模，第四维为证书之和的模 

void init()//解决实质方法 
{
	f[1][0][1]=1;
	f[1][1][0]=1; 
	for(int i=2;i<N;i++)
	{
		for(int j=0;j<=1;j++)
		{
			for(int k=i;k>=0;k--)//表示的是0的个数 
			{
				if(j)f[i][j][k]=f[i-1][0][k]+f[i-1][1][k];
				if(k&&!j) f[i][j][k]=f[i-1][0][k-1]+f[i-1][1][k-1];
			
			}
		}
	}
}

int dp(int x)//解决区间问题 
{
//	if(x==0)return 1;
	vector<int> nums;
	while(x)nums.push_back(x%2),x=x/2;
	
	int res=0;
	int last=0;//用来记录前面有关系的一个性质 
	int l=nums.size();
	
	
	for(int i=nums.size()-2;i>=0;i--)
	{
		int y=nums[i];
		if(y)//y为1———
		{
			for(int j=l-last;j>=(l+1)/2-last;j--)
			{
				if(j<=i+1)res+=f[i+1][0][j];
			}
		}
		if(!y)last++;
		if(!i&&last>=(l+1)/2)res++;
		
	} 
	for(int i=1;i<l;i++){
		for(int j=i;j>=(i+1)/2;j--)
			res+=f[i][1][j];
	}
	return res;
	
}
int main()
{
	init();
	cin>>a>>b;
	cout<<dp(b)-dp(a-1)<<endl;

	return 0;
 }

7、花神花神的数论题

众所周知，花神多年来凭借无边的神力狂虐各大 OJ、OI、CF、TC …… 当然也包括 CH 啦。

话说花神这天又来讲课了。课后照例有超级难的神题啦…… 我等蒟蒻又遭殃了。 花神的题目是这样的：设 sum ( i ) \text{sum}(i) sum(i) 表示 i i i 的二进制表示中 1 1 1 的个数。给出一个正整数 N N N ，花神要问你 ∏ i = 1 N sum ( i ) \prod_{i=1}^{N}\text{sum}(i) ∏i=1N​sum(i) ，也就是 sum ( 1 ) ∼ sum ( N ) \text{sum}(1)\sim\text{sum}(N) sum(1)∼sum(N) 的乘积。

一个正整数 N N N。

一个数，答案模 10000007 10000007 10000007 的值。

3

2

对于 100 % 100\% 100% 的数据， 1 ≤ N ≤ 1 0 15 1\le N\le 10^{15} 1≤N≤1015。

运用快速幂的方法

#include<bits/stdc++.h>
using namespace std;
typedef long long int ll;
ll x;
const ll N=55;
ll f[N][2][N];//长度为i，以j开头的，1的个数为k的数有多少个 
ll mod=10000007;
void init()//解决实质方法 
{
	f[1][0][0]=1;
	f[1][1][1]=1; 
	for(int i=2;i<N;i++)
	{
		for(int j=0;j<=1;j++)
		{
			for(int k=i;k>=0;k--)//表示的是1的个数 
			{
				if(j&&k)f[i][j][k]=f[i-1][0][k-1]+f[i-1][1][k-1];
				if(!j) f[i][j][k]=f[i-1][0][k]+f[i-1][1][k];
			}
		
		}                            
	}
}
ll quick_mi(ll a,ll b)
{
	ll res=1;
	while(b)
	{
		if(b&1) res=res*a%mod;
		b=b/2;
		a=a*a%mod;
	}
	return res;
}
ll dp(ll x)//解决区间问题 
{
//	if(x==0)return 1;
	vector<int> nums;
	while(x)nums.push_back(x%2),x=x/2;
	
	ll res=1;
	int last=1;//用来记录前面有关系的一个性质 
	int l=nums.size(); 
	for(int i=nums.size()-2;i>=0;i--)//从第二位分支开始  
	{
		int y=nums[i];
		if(y)//y为1———
		{
			for(int j=0;j<=i;j++)
			{
				res=(res*quick_mi(j+last,f[i+1][0][j])%mod)%mod;
			}
		}
		if(y)last++;//增加1的个数 
		
		if(!i)res=(res*last)%mod;
		
	} 
	for(int i=1;i<l;i++){
		for(int j=1;j<=i;j++)
			res=(res*quick_mi(j,f[i][1][j])%mod)%mod;
	}
	return res;
	
}
int main()
{
	init();
	cin>>x;
	cout<<dp(x)<<endl;

	return 0;
 }
元气满满的小龙
关注
专栏目录
统计学——中位数、众数
weixin_30662011的博客
03-26 2119
（2）中位数 中位数是指将数据按大小顺序排列起来，形成一个数列，居于数列中间位置的那个数据。中位数用Me表示。 　 从中位数的定义可知，所研究的数据中有一半小于中位数，一半大于中位数。中位数的作用与算术平均数相近，也是作为所研究数据的代表值。在一个等差数列或一个正态分布数列中，中位数就等于算术平均数。 在数列中出现了极端变量值的情况下，用中位数作为代表值要比用算术平均数更好，...
学习记忆——记忆宫殿——编码——数字编码——三位数
stqer的博客
09-06 1195
https://tieba.baidu.com/p/5074107399
乘积最大问题——动态规划
Snitro的博客
03-19 1万+
乘积最大问题——动态规划
数位dp
zzqwtc的博客
01-25 2194
数位统计DP 给定两个整数 a 和 b，求 a 和 b 之间的所有数字中0 ~ 9的出现次数。 假设n=abcdefgn = abcdefgn=abcdefg 计算ddd位上xxx出现的次数 resresres记录答案 1、(001——abc−1)(001——abc-1)(001——abc−1) xxx (000——999)(000——999)(000——999) res+=abc∗1000res += abc * 1000res+=abc∗1000 2、abcabcabc xxx ​ 2.1、d&l
PTA——逆序的三位数
Just__Joker的博客
11-24 1万+
PTA——逆序的三位数 题目: 程序每次读入一个正3位数，然后输出按位逆序的数字。注意：当输入的数字含有结尾的0时，输出不应带有前导的0。比如输入700，输出应该是7。 输入格式 每个测试是一个3位的正整数。 输出格式 输出按位逆序的数 输入样例 123 输出样例 321 解题思路如下： 输入一个三位数，接下来运用int类型整数的特性，将三位数的百位，十位，个位都存在三个变量中，最后在输出处，将存个位的变量乘100，十位乘10，百位乘1。即可得到逆序三位数。 整体代码如下： #include<st
java —— 数字反转
sokoooo的博客
01-07 5087
java —— 数字反转 给定一个整数，请将该数各个位上数字反转得到一个新数。新数也应满足整数的常见形式，即除非给定的原数为零，否则反转后得到的新数的最高位数字不应为零（参见样例2）。 输入格式：输入共1 行，一个整数N。 输出格式：输出共1 行，一个整数，表示反转后的新数。 限制： -1,000,000,000 ≤ N ≤ 1,000,000,000 。（每个测试点1s） 样例 1 ：
递归、迭代与动态规划
Wopelo的博客
07-09 3568
在这篇博客里，我将以计算斐波那契数列指定位置的数为例介绍递归、迭代与动态规划。首先我们得弄清楚这三者的定义。 递归——程序调用自身,也就是函数自己调用自己。递归通常从顶部将问题分解，通过解决掉所有分解出来的小问题，来解决整个问题； 迭代——利用变量的原值推算出变量的一个新值。递归中一定有迭代,但是迭代中不一定有递归,大部分可以相互转换； 动态规划——通常与递归相反，其从底部开始解决问题。将所有
蓝桥杯—— 算法训练 数位分离
Coder_PWH
04-13 2536
蓝桥杯—— 算法训练 数位分离        问题描述 　　编写一个程序，输入一个1000 以内的正整数，然后把这个整数的每一位数字都分离出来，并逐一地显示。 　　输入格式：输入只有一行，即一个1000以内的正整数。 　　输出格式：输出只有一行，即该整数的每一位数字，之间用空格隔开。 　　输入输出样例 样例输入 769 样例输出 7 6 9 mport java.u
位掩码(BitMask)——介绍与使用
热门推荐
数字世界 万物互联
06-30 1万+
今天在群里看到有人讨论小白鼠试毒药的问题，感觉很有意思，生活中处处都有技术解决问题的思想，思想的提高也会带动技术水平的提高吧，转自一位网友的帖子，仅记录学习应用，多谢。 一、前言 位运算在我们实际开发中用得很少，主要原因还是它对于我们而言不好读、不好懂、也不好计算，如果不经常实践，很容易就生疏了。但实际上，位......
Python蓝桥杯练习03——洁净数
m0_48936146的博客
02-27 5050
1.题目 洁净数 2.题目描述 小明不喜欢数字2，包括那些数位上包含数字2的数，如果一个数的数位不包含数字2，小明将它成为洁净数。 请问在整数1至n中，洁净数有多少个？ 3.输入格式 输入的第一行包含一个整数 4.输出格式 输出一行包含一个整数，表示答案 5.样例输入 30 6.样例输出 18 7.字符串方法 n=int(input()) s=0 for i in range(1,n+1): if '2' not in str(i): #将遍历i转换成字符串
大一新生工作学习生活必备——数位名校校长总结.pdf
02-26
大一新生工作学习生活必备——数位名校校长总结.pdf
信息化教学设计——三位数乘以两位数.doc
11-06
信息化教学设计——三位数乘以两位数是针对小学四年级数学的一次教学活动，旨在教授学生如何进行三位数乘以两位数的笔算。教学设计基于学生已经掌握的两位数乘以两位数的基础，强调算理和算法的理解，并通过情境教学...
C语言练习——3位数各位数字的立方和等于本身.c
06-02
C语言基础的练习
动态规划算法python_动态规划——DP算法（Dynamic Programing）
weixin_39527372的博客
11-29 1722
一、斐波那契数列（递归VS动态规划）1、斐波那契数列——递归实现（python语言）——自顶向下递归调用是非常耗费内存的，程序虽然简洁可是算法复杂度为O(2^n)，当n很大时，程序运行很慢，甚至内存爆满。1 deffib(n):2 #终止条件，也就是递归出口3 if n == 0 or n == 1:4 return 15 else:6 #递...
【AGC005D】~K Perm Counting（计数抽象成图）
2301_77025310的博客
10-01 2813
注意到位置为id，权值为v ,不合法的情况，当且仅当 v = id+k或 v= id-k。dp(i,j,pd)表示考虑到第i号点， 连了j条边，是否有连接i 到 i-1号点。因此，我们把每一个位置和权值抽象成点 ，不合法的情况之间连一条边，可以构成二分图。由此可知，当选了n条边，就恰好n个位置不合法，限制条件是：连的边不能相邻，求出f(m) ，f(m)指代至少有m个位置不合法的方案数。由此总共有2n 个点 k 条链，链与链之间无边 互不干涉。把二分图展开成k条链，进行dp。简单的乘法原理罢了。
【算法】（Python）动态规划
yannan20190313的博客
11-07 894
【算法】（Python）动态规划（最小花费爬楼梯，买卖股票最佳时机，最低加油次数）
Day41 | 动态规划 ：完全背包应用 完全平方数&&单词拆分（类比爬楼梯）
最新发布
m0_74795952的博客
11-09 798
代码都不难写，
Java 代理模式详解
aaaa_1111111的博客
11-05 900
姨妈在这里就可以看作是代理你的代理对象，代理的行为（方法）是接收和回复新郎的提问。涵盖了静态代理和动态代理实战、静态代理和动态代理的区别、JDK 动态代理和 Cglib 动态代理区别等内容。， 例如 Spring 中的 AOP 模块中：如果目标对象实现了接口，则默认采用 JDK 动态代理，否则采用 CGLIB 动态代理。说到动态代理，Spring AOP、RPC 框架应该是两个不得不提的，它们的实现都依赖了动态代理。代理模式有静态代理和动态代理两种实现方式，我们 先来看一下静态代理模式的实现。
Spring框架之代理模式 (Proxy Pattern)
m0_61840987的博客
11-08 1115
代理模式是一种非常有用的设计模式，能够在不修改目标对象的情况下增强其功能。通过静态代理和动态代理的不同实现方式，可以根据具体场景选择最合适的代理模式。代理模式在实际开发中应用广泛，尤其是在权限控制、日志记录、缓存管理和事务处理等方面发挥了重要作用。理解和灵活应用代理模式，将有助于开发出更为灵活、可维护的系统。
数位动态规划：区间统计与深度解析
第3章 "数位动态规划" 是2020年4月22日的一份资料，主要探讨了在信息奥林匹克（信奥）竞赛中常见的一种问题类型——与数位相关的区间统计问题。这类问题通常具有深厚的数学背景，常规的组合数学方法可能不足以直接...
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值