热门标签 | HotTags
当前位置:  开发笔记 > 编程语言 > 正文

2022杭电多校七FSumire(数位DP+实用技巧)

题目链接:杭电多校7-VirtualJudgevjudge上题目显示的有问题,我下面附上官方题目:样例输入:32201

题目链接:杭电多校7 - Virtual Judge

vjudge上题目显示的有问题,我下面附上官方题目:

样例输入: 

3
2 2 0 1 5
1 4 3 11 45
10 14 11 19 198

样例输出:

6
19
1049

题意:多组样例,每组样例给出五个数k,b,d,l,r,其中b是代表b进制,f(i,b,d)代表在b进制下i的数位中出现数字d的个数。求

 分析:定义f[i][j][1/0][1/0]表示已经统计了前i位有j位是d且当前位有/无限制,有/无前导0的方案数

一开始我按照普通的数位DP思路求解,但是总是超时,然后仔细分析了一下,发现了原来的数位DP还可以具体根据题目优化一下。

比如这道题,我一开始是分两种情况来进行动态转移:

枚举当前填的数字

一种情况是当前有前导0且当前位置填0。

另一种情况就是当前位置填无前导0或者当前位置不填0,这两个合成一种情况即可。

dfs代码是这样写的:

long long dfs(int pos,int cnt,int limit,int lead)
{if(pos&#61;&#61;0) return (cnt&#61;&#61;tt);if(!limit&&!lead&&f[pos][cnt]!&#61;-1) return f[pos][cnt];int up&#61;limit?a[pos]:(b-1);long long ans&#61;0;for(int i&#61;0;i<&#61;up;i&#43;&#43;){if(lead&&i&#61;&#61;0) ans&#43;&#61;dfs(pos-1,cnt,limit&&(i&#61;&#61;a[pos]),lead);else ans&#43;&#61;dfs(pos-1,cnt&#43;(i&#61;&#61;d),limit&&(i&#61;&#61;a[pos]),0); }if(!limit&&!lead) f[pos][cnt]&#61;ans;return ans;
}

后来对这样的代码进行了很多优化发现还是超时&#xff0c;仔细分析了一下才发现当前位置如果不填d或者up或者0&#xff0c;那么填其他任何数字都是一样的&#xff0c;因为只有d才会对最后的答案造成影响&#xff0c;而up会对下一次可以取的数字造成影响&#xff0c;而0可以对前导0造成影响&#xff0c;所以对于一种情况&#xff0c;如果我们可以在0~up中做出选择&#xff0c;去掉d、up、0之外的数选哪个都是一样的&#xff0c;我们可以直接计算一下这样的数字的个数&#xff0c;然后直接令ans&#43;可选数字个数*dfs(pos-1&#xff0c;cnt&#xff0c;0&#xff0c;0)即可&#xff0c;因为选的数既不是上限也不是0&#xff0c;所以会去掉限制和前导0&#xff0c;而且由于选的数不是d&#xff0c;那么也不会对答案造成贡献&#xff0c;所以我们就可以把很多次重复的计算直接一次算出来&#xff0c;这样的话就会极大地优化时间&#xff0c;从而就可以ac&#xff0c;有同学可能会问&#xff1a;我们都记忆化了&#xff0c;算一次和算十次有什么本质区别吗&#xff0c;算完一次之后再算就直接访问数组就可以了啊&#xff0c;为什么这样还会节省时间呢&#xff1f;可以这么想假如说我们刚才讨论的那种不是up和d和0的情况有10000种&#xff0c;那么我们就需要访问10000次数组&#xff0c;这样算下来复杂度会高很多&#xff0c;而我们之前的数位DP题目一般进制数都很少&#xff0c;或者说一般都是10进制的&#xff0c;而这道题b的上界是1e9&#xff0c;这个数还是很大的&#xff0c;所以我们这种优化对于进制数比较大的问题还是很有用的。

这个技巧近适用于贡献与某位数字有直接关系的题目中。

细节见代码&#xff1a;

#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;
long long k,b,d;
const int N&#61;65,mod&#61;1e9&#43;7;
long long f[N][N][2][2];//f[i][j][1/0][1/0]表示已经统计了前i位有j位是d且当前位有/无限制&#xff0c;有/无前导0的方案数
long long a[N],tt;
long long qpow(long long a,long long b)
{if(a&#61;&#61;0) return 0;//题目中要求0^0&#61;0long long ans&#61;1;while(b){if(b&1) ans&#61;ans*a%mod;a&#61;a*a%mod;b>>&#61;1;}return ans;
}
long long dfs(int pos,int cnt,int limit,int lead)
{if(!pos) return qpow(cnt,k);if(f[pos][cnt][limit][lead]!&#61;-1) return f[pos][cnt][limit][lead];int up&#61;limit?a[pos]:(b-1);long long ans&#61;0;if(d&#61;&#61;up){if(lead){ans&#43;&#61;dfs(pos-1,cnt,0,1);//有前导0且填0ans&#43;&#61;dfs(pos-1,cnt&#43;1,limit,0);//有前导0填d ans&#43;&#61;(up-1)*dfs(pos-1,cnt,0,0)%mod;//有前导0但不填up和0 }else{ans&#43;&#61;dfs(pos-1,cnt&#43;1,limit,0);//无前导0填dans&#43;&#61;up*dfs(pos-1,cnt,0,0)%mod;//无前导0但不填d} ans%&#61;mod;}else if(d&#61;&#61;0){if(lead){ans&#43;&#61;dfs(pos-1,cnt,0,1);//有前导0且填0ans&#43;&#61;dfs(pos-1,cnt,limit,0);//有前导0且填upans&#43;&#61;(up-1)*dfs(pos-1,cnt,0,0)%mod;//有前导0但不填up和d }else{ans&#43;&#61;dfs(pos-1,cnt&#43;1,0,0);//无前导0且填dans&#43;&#61;dfs(pos-1,cnt,limit,0);//无前导0且填upans&#43;&#61;(up-1)*dfs(pos-1,cnt,0,0)%mod;//无前导0且不填up和d } ans%&#61;mod;}else if(up>d)//d不等于up和0,但是要保证up大于d {if(lead){ans&#43;&#61;dfs(pos-1,cnt,0,1);//有前导0且填0ans&#43;&#61;dfs(pos-1,cnt&#43;1,0,0);//有前导0且填dans&#43;&#61;dfs(pos-1,cnt,limit,0);//有前导0且填upans&#43;&#61;(up-2)*dfs(pos-1,cnt,0,0)%mod;//有前导0但不填up和d和0 }else{ans&#43;&#61;dfs(pos-1,cnt&#43;1,0,0);//无前导0且填dans&#43;&#61;dfs(pos-1,cnt,limit,0);//无前导0且填upans&#43;&#61;(up-1)*dfs(pos-1,cnt,0,0)%mod;//无前导0但不填up和d} ans%&#61;mod;}else//up}
long long solve(long long x)
{memset(f,-1,sizeof f);int pos&#61;0;while(x){a[&#43;&#43;pos]&#61;x%b;x/&#61;b;}return dfs(pos,0,1,1);
}
int main()
{int T;cin>>T;long long l,r;while(T--){scanf("%lld%lld%lld%lld%lld",&k,&b,&d,&l,&r);printf("%lld\n",(solve(r)-solve(l-1)&#43;mod)%mod);}return 0;
}


推荐阅读
  • 在1995年,Simon Plouffe 发现了一种特殊的求和方法来表示某些常数。两年后,Bailey 和 Borwein 在他们的论文中发表了这一发现,这种方法被命名为 Bailey-Borwein-Plouffe (BBP) 公式。该问题要求计算圆周率 π 的第 n 个十六进制数字。 ... [详细]
  • 深入解析 C++ 中的 String 和 Vector
    本文详细介绍了 C++ 编程语言中 String 和 Vector 的使用方法及特性,旨在帮助开发者更好地理解和应用这两个重要的容器。 ... [详细]
  • 本文详细介绍了在Luat OS中如何实现C与Lua的混合编程,包括在C环境中运行Lua脚本、封装可被Lua调用的C语言库,以及C与Lua之间的数据交互方法。 ... [详细]
  • 想把一组chara[4096]的数组拷贝到shortb[6][256]中,尝试过用循环移位的方式,还用中间变量shortc[2048]的方式。得出的结论:1.移位方式效率最低2. ... [详细]
  • Hanks博士是一位著名的生物技术专家,他的儿子Hankson对数学有着浓厚的兴趣。最近,Hankson遇到了一个有趣的数学问题,涉及求解特定条件下的正整数x,而不使用传统的辗转相除法。 ... [详细]
  • 本文通过C++语言实现了一个递归算法,用于解析并计算数学表达式的值。该算法能够处理加法、减法、乘法和除法操作。 ... [详细]
  • 洛谷 P4009 汽车加油行驶问题 解析
    探讨了经典算法题目——汽车加油行驶问题,通过网络流和费用流的视角,深入解析了该问题的解决方案。本文将详细阐述如何利用最短路径算法解决这一问题,并提供详细的代码实现。 ... [详细]
  • 本文详细介绍了如何在 Ubuntu 14.04 系统上搭建仅使用 CPU 的 Caffe 深度学习框架,包括环境准备、依赖安装及编译过程。 ... [详细]
  • 本文探讨了异步编程的发展历程,从最初的AJAX异步回调到现代的Promise、Generator+Co以及Async/Await等技术。文章详细分析了Promise的工作原理及其源码实现,帮助开发者更好地理解和使用这一重要工具。 ... [详细]
  • 本文详细介绍了 Redis 中的主要数据类型,包括 String、Hash、List、Set、ZSet、Geo 和 HyperLogLog,并提供了每种类型的基本操作命令和应用场景。 ... [详细]
  • 本文详细探讨了在Java中如何将图像对象转换为文件和字节数组(Byte[])的技术。虽然网络上存在大量相关资料,但实际操作时仍需注意细节。本文通过使用JMSL 4.0库中的图表对象作为示例,提供了一种实用的方法。 ... [详细]
  • 本文详细介绍了Linux系统中信号量的相关函数,包括sem_init、sem_wait、sem_post和sem_destroy,解释了它们的功能和使用方法,并提供了示例代码。 ... [详细]
  • 问题描述现在,不管开发一个多大的系统(至少我现在的部门是这样的),都会带一个日志功能;在实际开发过程中 ... [详细]
  • 本问题涉及在给定的无向图中寻找一个至少包含三个节点的环,该环上的节点不重复,并且环上所有边的长度之和最小。目标是找到并输出这个最小环的具体方案。 ... [详细]
  • 字符串中特定模式出现次数的计算方法
    本文详细探讨了如何高效地计算字符串中特定模式(如'pat')的出现次数,通过实例分析与算法解析,帮助读者掌握解决此类问题的方法。 ... [详细]
author-avatar
捕风的默小墨
这个家伙很懒,什么也没留下!
PHP1.CN | 中国最专业的PHP中文社区 | DevBox开发工具箱 | json解析格式化 |PHP资讯 | PHP教程 | 数据库技术 | 服务器技术 | 前端开发技术 | PHP框架 | 开发工具 | 在线工具
Copyright © 1998 - 2020 PHP1.CN. All Rights Reserved | 京公网安备 11010802041100号 | 京ICP备19059560号-4 | PHP1.CN 第一PHP社区 版权所有