数位dp训练笔记

已于 2023-06-25 18:39:01 修改
阅读量213 收藏
点赞数
分类专栏： 学习笔记 搜索 文章标签： 笔记
于 2023-06-04 10:03:13 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/sophilex/article/details/131026042
版权

数位dp基本上是处理位数相关的问题，不过也不一定

大部分题目存在一定套路，但是有一些也不好想

其dp状态的设置大致是由枚举的位数，以及每一位的状态（是否前导0，是否顶到上界），以及一些依题目而定的性质

时间复杂度基本上就是dp枚举的状态数

没了

[AHOI2009]同类分布

大意：

给出两个数a,b，求出[a,b]中各位数字之和能整除原数的数的个数。

a,b<=1e18

思路：

考虑dp状态dp[i][j][k]表示枚举到前i位，前i位数字的和是j，前i位组成的数字是k

这是一个比较niave的想法，但是数据范围不支持我们这样处理第三维。

考虑到最后只要求整除，所以可以考虑用k%j来代替k

但是这样还涉及到一个问题，就是如果j是不断变化的，就很难实现状态转移

所以我们可以在外层枚举j，然后dfs的时候就保持j不变

就好了

由于这里是考虑每一位数字的和，所以我们不用考虑前导0的问题

code

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define endl '\n'
const ll N=1e5+10;
ll mod;
ll n,m;
ll a[20];
ll cnt=0;
ll dp[20][200][200];
ll dfs(ll x,ll sum,ll rel,ll op)
{
	if(x==0) return sum==0&&rel==0;
	if(!op&&dp[x][rel][sum]!=-1) return dp[x][rel][sum]; 
	ll lim=op?a[x]:9;
	ll tot=0;
	for(int i=0;i<=lim&&i<=sum;++i)
	{	
		tot+=dfs(x-1,sum-i,(rel*10%mod+i)%mod,op&&i==lim);	
	}
	if(!op) dp[x][rel][sum]=tot;
	return tot;
}
ll f(ll x)
{
	cnt=0;
	while(x)
	{
		a[++cnt]=x%10;
		x/=10;
	}
	ll det=0;
	for(int i=1;i<=9*cnt;++i)
	{
		mod=i;
		memset(dp,-1,sizeof dp);
		det+=dfs(cnt,i,0,1);
	}
	return det;
}
void solve()
{
	
	cin>>n>>m;
	cout<<f(m)-f(n-1)<<endl;
}
int main()
{
	ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
//	ll t;cin>>t;while(t--)
	solve();
	return 0;
}

虽然是紫题，但是理解了套路就很水

Classy Numbers

大意：

给出两个数a,b，求出[a,b]中各位数字中非0数不大于3的数字个数。

a,b<=1e18

思路：
板子题

这里显然需要考虑前导0

code

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define endl '\n'
const ll N=1e5+10;
ll b,l,r;
int a[20];
ll cnt=0;
ll dp[20][200];
ll dfs(ll x,int sum,ll op)
{
	if(x==0) return sum<=3;
	if(!op&&dp[x][sum]!=-1) return dp[x][sum];
	int lim=op?a[x]:9;
	ll tot=0;
	for(int i=0;i<=lim;++i)
	{
		tot+=dfs(x-1,sum+(i!=0),op&&i==lim);
	}
	if(!op) dp[x][sum]=tot;
	return tot;
}
ll f(ll x)
{
	cnt=0;
	while(x)
	{
		a[++cnt]=x%10;
		x/=10;
	}
	return dfs(cnt,0,1);
}
void solve()
{
	cin>>l>>r;
	cout<<f(r)-f(l-1)<<endl;
}
int main()
{
	ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
	memset(dp,-1,sizeof dp);
	ll t;cin>>t;while(t--)
	solve();
	return 0;
}

Segment Sum

大意：

给定K,L,R，求L~R之间最多不包含超过K种数码的数的和。

K<=10，L,R<=1e18

思路：
如果只是求满足条件的数字的个数的话就是上面提到的板子题了

这里要求和，我们同样可以考虑对相同状态进行合并求和

f[i][j]表示当前枚举到第i位，出现过的数码种类的状态为j，也就是我们要求的答案数组

g[i][j]表示当前枚举到第i位，出出现过的数码种类的状态为j的合法数字个数，也就是板子

考虑如何用g来推f

这里j可以10位二进制状压

我们每往下走一位，如果我们枚举第i位上填的数字是t，用j'表示下一位的数码种类数

g[i][j]=\sum g[i-1][j']

f[i][j]=\sum 10^i*t*g[i-1][j']+f[i-1][j']

稍微意会一下应该就能懂了

f是数字的和，g是合法数字的个数

有了这个之后，我们就直接推就可以了。然后因为需要下一个状态的f和g，所以我们dfs要返回两个值

code

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define pii pair<ll,ll>
#define mk make_pair
#define endl '\n'
const ll N=1e5+10;
const ll mod=998244353;
ll n,m,k;
ll a[20];
ll cnt=0;
pii dp[40][1030];
ll p[40];
void init()
{
	p[0]=1;
	for(int i=1;i<=20;++i) p[i]=p[i-1]*10ll%mod;
}
bool check(ll x)
{
	ll cn=0;
	while(x)
	{
		cn+=(x%2);
		x/=2;
	}
	return cn<=k;
}
pii dfs(ll x,ll sum,ll head,ll op)
{
	if(x==0) return mk(0,1);
	if(!op&&!head&&dp[x][sum]!=mk(-1ll,-1ll)) return dp[x][sum];
	ll lim=op?a[x]:9;
	ll s1=0,s2=0;
	for(ll i=0;i<=lim;++i)
	{
		//f是数字的和，g是合法数字的个数
		pii gt=mk(0,0);
		if(head&&i==0) gt=dfs(x-1,0,1,op&&i==lim);
		else if(check(sum|(1<<i))) gt=dfs(x-1,sum|(1<<i),0,op&&i==lim);
		s1=(((s1+i*p[x-1]%mod*gt.second%mod)%mod)+gt.first)%mod;	
		s2=(s2+gt.second)%mod;
	} 
	if(!op&&!head) dp[x][sum]=mk(s1,s2);
	return mk(s1,s2);
}
ll f(ll x)
{
	cnt=0;
	while(x)
	{
		a[++cnt]=x%10;
		x/=10;
	}
	return dfs(cnt,0,1,1).first;
}
void solve()
{
	init();
	for(int i=0;i<=20;++i)
	{
		for(int j=0;j<=1025;++j)
		{
			dp[i][j]=mk(-1ll,-1ll);
		}
	}
	cin>>n>>m>>k;
	cout<<((f(m)-f(n-1))%mod+mod)%mod<<endl;
}
int main()
{
	ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
	solve();
	return 0;
}

 E. Salazar Slytherin's Locket

大意：
求l...r之间转成b进制后，每一位都是偶数的数的个数

思路：
不再是10进制，只要在预处理每一位的时候换个底数就好了

注意一下前导0

code

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define endl '\n'
const ll N=1e5+10;
ll b,l,r;
int a[70];
ll cnt=0;
ll dp[70][1026][12];
ll dfs(ll x,int sta,ll op,ll head)
{
	if(x==0) return sta==0;
	if(!op&&!head&&dp[x][sta][b]!=-1) return dp[x][sta][b];
	int lim=op?a[x]:b-1;
	ll tot=0;
	for(int i=0;i<=lim;++i)
	{
		if(head&&i==0) tot+=dfs(x-1,0,op&&i==lim,1);
		else tot+=dfs(x-1,sta^(1<<i),op&&i==lim,0);
	}
	if(!op&&!head) dp[x][sta][b]=tot;
	return tot;
}
ll f(ll x)
{
	cnt=0;
	while(x)
	{
		a[++cnt]=x%b;
		x/=b;
	}
	return dfs(cnt,0,1,1);
}
void solve()
{
	//now++;
	cin>>b>>l>>r;
	cout<<f(r)-f(l-1)<<endl;
}
int main()
{
	ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
	memset(dp,-1,sizeof dp);
	ll t;cin>>t;while(t--)
	solve();
	return 0;
}

花神的数论题

大意：
设sum(i)表示i的二进制表示中1的个数，求\prod_{i=1}^{n} sum(i)

思路：
考虑枚举贡献，枚举二进制1的个数为j，假设有k个这样的数字，答案就是\prod j^{k}

code

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const ll mod=1e7+7;
const ll N=70;
ll cnt[N];
ll n;
ll p,ans=1;
ll dp[N][N];//dp[i][j],当前枚举到第i位，前面的数中
//有j个1的情况下最终满足sum值等于p的数字的个数 
ll ksm(ll x,ll y)
{
	ll ans=1;
	while(y)
	{
		if(y&1) ans=ans*x%mod;
		x=x*x%mod;
		y>>=1;
	}
	return ans;
}

ll dfs(ll pos,bool limit,ll sum)
{
	if(pos==0) return (sum==p);//边界条件
	if(!limit&&~dp[pos][sum]) return dp[pos][sum];
	ll up=limit?cnt[pos]:1;
	ll cn=0;
	for(int i=0;i<=up;++i)
	{ 
		if(sum+(i==1)>p) continue;
		cn+=dfs(pos-1,limit&&(i==up),sum+(i==1));
	} 
	if(!limit) dp[pos][sum]=cn;
	return cn;
	 
}
void solve()
{
	ll d=0;
	while(n)
	{
		cnt[++d]=n%2;
		n/=2;
	}
	for(int i=1;i<=d;++i)
	{
		memset(dp,-1,sizeof dp);
		p=i;
		ans=ans*ksm(i,dfs(d,1,0))%mod;
	}
	cout<<ans<<endl;
}
int main()
{
	cin>>n;
	solve();
	return 0;
}

萌数

[HAOI2010]计数

[SCOI2014]方伯伯的商场之旅 题解

Magic Numbers

大意：

思路：

这题其实就跟上面的没有什么区别，只要按要求做就好了。前两个限制都很好处理，对于最后一个，我们只要在dp状态里加上一个前缀模m的值即可

code

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define endl '\n'
const ll N=1e5+10;
const ll mod=1000000007;
ll m,d;
string l,r;
ll a[2010];
ll cnt=0;
ll dp[2010][2005];
ll dfs(ll x,ll mo,ll op)
{
	if(x==0) return mo==0;
	if(!op&&dp[x][mo]!=-1) return dp[x][mo];
	ll tot=0;
	ll lim=op?a[x]:9;
	ll det=cnt-x+1;
	if(det%2==0)
	{
		if(d>lim) return 0;
		tot=(tot+dfs(x-1,(mo*10%m+d)%m,op&&d==lim))%mod;
	}
	else
	{
		for(int i=0;i<=lim;++i)
		{
			if(i==d) continue;
			tot=(tot+dfs(x-1,(mo*10%m+i)%m,op&&i==lim))%mod;
		}
	}
	if(!op) dp[x][mo]=tot;
	return tot;
}
ll f(string s,ll fl)
{
	memset(dp,-1,sizeof dp);
	cnt=0;
	ll len=s.size();
	for(int i=len-1;i>=0;--i)
	{
		a[++cnt]=s[i]-'0';
	}
//	for(int i=1;i<=cnt;++i) cout<<a[i];
//		cout<<endl;
	if(fl)
	{
		a[1]--;
		for(int i=1;i<=cnt;++i)
		{
			if(a[i]==-1)
			{
				a[i]=9;
				a[i+1]--;
			}
		}
//		for(int i=1;i<=cnt;++i) cout<<a[i];
//		cout<<endl;
	}
	return dfs(cnt,0,1);
}
void solve()
{
	cin>>m>>d>>l>>r;
	cout<<((f(r,0)-f(l,1))%mod+mod)%mod<<endl;
}
int main()
{
	ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
//	ll t;cin>>t;while(t--)
	solve();
	return 0;
}

 Beautiful numbers

大意：

思路：

 我们很难在枚举的过程中将dp状态设置为与x模2-9的值都有关，因为这样不好合并相同状态。

这里有两个小结论

1.两个整数a,b满足a|x,b|x的充要条件是lcm(a,b)|x

Proof：

->:

由a|x,b|x知，x=m*a=n*b

则我们有a|(n*b)

若(a,b)=1,则由上述结论我们有a|n,故x=n*b=(k*a)*b=k*lcm(a,b),即lcm(a,b)|x     ----------1

若(a,b)=d,则d*(a/d)=n*d*(b/d),即(a/d)=n*(b/d)，其中(a/d,b/d)=1，转为推导1

故该方向得证

<-:

显然

则命题得证

推广一下就是：若干个数ai都整除x的充要条件是lcm{ai}|x

2.考虑n个数a1,a2...an,记其LCM为L，则x%ai=(x%L)%ai

证明显然 

有了上述性质，我们不难发现，可以用x%2520(1-9的最小公倍数)来代替x的值，那么首先值域就已经大大简化了。然后由性质1，我们只需要维护出现过数字的lcm即可，最后的判断条件就是前缀%2520的值%lcm=0

那么我们记录dpi,j,k表示前i位，前缀模2520的值为j，前i位中非0位的lcm为k，就可以dp了。但是空间有点大，考虑优化。

不难发现，第三维其实并没有很多数字，2-9中任意若干个数字的lcm一定是2520的因数，而2520的因数总共只有48个，就大大简化了空间，并且我们可以预处理

这样就ok了

code

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define endl '\n'
const ll N=1e5+10;
const ll mod=2520;
ll n,m;
ll a[22];
ll cnt=0;
int mp[2522];
ll dp[22][2522][52];
int Gt[2522][2522];
ll cm(ll x,ll y)
{
	if(x==0||y==0) return x|y;
	ll Gcd=__gcd(x,y);
	if(x>=y) swap(x,y);
	if(Gt[x][y]) return Gt[x][y]; 
	return Gt[x][y]=x*y/Gcd;
}
ll dfs(ll x,ll op,ll mo,ll Lcm)
{
	if(x==0) return mo%Lcm==0;
	if(!op&&dp[x][mo][mp[Lcm]]!=-1) return dp[x][mo][mp[Lcm]];
	ll tot=0;
	ll lim=op?a[x]:9;
	for(int i=0;i<=lim;++i)
	{
		if(i==0) tot+=dfs(x-1,op&&i==lim,mo*10%mod,Lcm);
		else tot+=dfs(x-1,op&&i==lim,(mo*10%mod+i)%mod,cm(Lcm,i));
	}
	if(!op) dp[x][mo][mp[Lcm]]=tot;
	return tot;
	
}
void init()
{
	ll cn=0;
	for(int i=1;i<=2520;++i)
	{
		if(2520%i) continue;
		mp[i]=++cn;
	}
//	cout<<cn<<endl;
}
ll f(ll x)
{
	cnt=0;
	while(x)
	{
		a[++cnt]=x%10;
		x/=10;
	}
	return dfs(cnt,1,0,1);
}
void solve()
{
	
	cin>>n>>m;
	cout<<f(m)-f(n-1)<<endl;	
}
int main()
{
	ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
	memset(dp,-1,sizeof dp);
	init();
	ll t;cin>>t;while(t--)
	solve();
	return 0;
}

最后再提一嘴，有些时候会遇到题目说每一个数字是按位排序过的，也就是同一个数字内数位升序排列，然后要求处理相关问题。比如CF908G New Year and Original Order，[SDOI2010]代码拍卖会，这种题有一个不常见但是很关键的套路，就是将每一个数字拆分成若干个由1组成的数字的和，这个性质是由其数位升序带来的，做题的时候要留一个心眼。但是本人太菜，这两题还不会...留个坑。

sophilex
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
数位dp算法详细图文解答
02-13
PPT图片结合文字讲述数位dp
数位DP.pdf
07-12
数位DP整理的模板和一些题目的题解。QAQ没钱下载资源了只好上传点东西。
数位DP入门
08-22
数位DP入门思想，包括详细的思考过程以及一些习题的讲解
初探数位dp.pdf
09-16
lazycal的集训队报告：初探数位DP 以HDU 2089，HDU 3652, URAL 1057等题目为例，介绍了数位DP的算法
数位dp——某次队内分享任务
07-21
好不容易做的，放在电脑中不小心删了咋整
随心笔记，第六更
leen的博客
06-13 473
今天给大家分享的是使用idea 将xml转为JavaBean，并且取其中的数据和其他的数据进行交互。
Python学习笔记11 -- 细碎内容
m0_71291382的博客
06-12 223
记录一些琐碎点，包括如何注释、变量常量说明、字符串拼接的三种方法及转义字符的简单说明
论文阅读笔记：DepGraph: Towards Any Structural Pruning
HollowKnightz的博客
06-15 895
论文阅读笔记：DepGraph: Towards Any Structural Pruning
《雅思口语真经总纲1.0》笔记——第二章：官方评分标准真经——Lexical Resource 词汇资源——三十组同义形容词列表
Dontla的博客
06-09 619
直接提取词汇重点，不看例句了。
[ue5]建模场景学习笔记（5）——必修内容可交互的地形，交互沙（3）
qq_40128113的博客
06-09 316
5.第一想法直接clear掉，但直接clear会把想要保留的痕迹一并去除，跑起来就没有痕迹！所以不行，要做一个转载，先把痕迹保留在一个载体上，清楚原来痕迹，用sandsave表示，再把这个载体的痕迹保存回原来rt，这样之前的痕迹被清空了，只留下了保存在载体上的部分，正好符合我们想要的部分。我们现在已经能够让这片地形出现在任意地方，只要角色走在这片地形上，就能够产生痕迹，但这片区域总是需要人工指定，又无法把这片区域无限扩大（显存爆炸），因此尝试使角色无论走到哪里都能产生痕迹。否则绘制不到画板上。
Nginx学习笔记（十）如何配置HTTPS协议？（公网）
ACGkaka的博客
06-10 706
Nginx学习笔记（十）如何配置HTTPS协议？（公网）
笔记97：C++ 中 string / char 和 int 之间相互转化
weixin_44063383的博客
06-12 311
方法：使用强制类型转化进行 char 和 int 之间的转换；
802.11漫游流程简单解析与笔记_Part3
Mr_liu_666的博客
06-10 236
看过Part1的应该都知道，标准关联流程是auth*2 associate*2 key*4，但ns3里面没有与auth有关的部分，third脚本里面有cap字段，但无auth也无key，所以part3是空的
程序固化——FPGA学习笔记6
m0_69082048的博客
06-12 459
BIN:一般是由Vivado软件编译产生的，存储在特定目录下的二进制文件MCS:一般通过VivadoGUl界面操作或者TCL命令生成，MCS文件里包含了BIN文件的内容，除此之外，每行的开始有地址信息，最后一个Byte是CRC校验信息。
【perfetto分析性能学习笔记】
q_2474932226的博客
06-12 280
perfrtto性能分析
1）Java项目笔记搭建系统梳理相关知识
最新发布
WLKUYM的博客
06-15 354
工作有年头了，学到了很多技术，收获了很多。但是对与工作相关的专业技能知识的掌握杂而乱，不够全面系统。因此做了个项目笔记来辅助梳理所掌握内容，方便需要时拿出来看看，前面做的时候，反反复复重建过好几次，目前的笔记结构较为满意（看过很多资料，没有令我满意的，别人总结的再好，终究不是自己的）事先说明，咱不是卖书的，我搞得也是电子书。毕竟实体书太沉了，贵（或许买了就是块板砖，压箱底吃灰了，浪费钱），拿着也不方便，各位可以从脚本之家（宝藏哦）搜索对应书籍下载。
python实践笔记(一): 模块和包
Miracle8070
06-10 574
这次是从实践的再去补充python的相关内容，目的是要写出漂亮的python代码，增强代码的可读，可维护，灵活和高效，方便调试和监控
【电路笔记】-共集极放大器
视觉与物联智能
06-15 613
本文介绍另一种用于放大信号的双极晶体管架构，通常称为共集电极放大器 (CCA)。 CCA 有时也可以称为射极跟随器放大器，我们将在本文后面了解其原因。
Matlab笔记
BiGmowang的博客
06-04 609
xor(3<5,1<3)--logical=0---xor中当且仅当一个是正确的返回1。xor(3<5,1<3)--logical=0---xor中当且仅当一个是正确的返回1。intmax('int8')--int8能存的最大值--a=int8(17)intmax('int8')--int8能存的最大值--a=int8(17)plot(x,y);plot(x,y);mod(13,3)--13对3取余--rem(13,3)差不多。mod(13,3)--13对3取余--rem(13,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、付费专栏及课程。

余额充值