2019-2020 ICPC Asia Hong Kong Regional Contest J. Junior Mathematician 题解(数位dp)
题目链接
题目大意
要你在[l,r]中找到有多少个数满足\(x\equiv f(x)(mod\; m)\)
\(f(x)=\sum_{i=1}^{k-1} \sum_{j=i+1}^{k}d(x,i)*d(x,j)\)
\(d(x,i)表示x的第i位数\)
题目思路
显然是数位dp,然而这个数位dp不能同时存x%m 和f(x)%m
这样会内存太大存不了,所以存差值即可
还有这个dfs的时候取模只取一次,不然会t,卡常严重
代码
#include<set>
#include<map>
#include<queue>
#include<stack>
#include<cmath>
#include<cstdio>
#include<vector>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<unordered_map>
#define fi first
#define se second
#define debug printf(" I am here\n");
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
const ll INF=0x3f3f3f3f3f3f3f3f;
const int maxn=5e3+5,inf=0x3f3f3f3f,mod=1e9+7;
const double eps=1e-10;
char l[maxn],r[maxn];
int m,base[maxn];
ll dp[maxn][100][100];
int lenl,lenr,num[maxn];
ll dfs(int pos,int pre,int dif,bool flag){
if(pos==0) return dif==0;
if(!flag&&dp[pos][pre][dif]!=-1){
return dp[pos][pre][dif];
}
int lim=flag?num[pos]:9;
ll ans=0;
for(int i=0;i<=lim;i++){
ans=ans+dfs(pos-1,(pre+i)%m,((dif+i*base[pos]-i*pre)%m+m)%m,flag&&i==lim);
}
ans%=mod;
if(!flag){
dp[pos][pre][dif]=ans;
}
return ans;
}
ll solve1(){
for(int i=1;i<=lenr;i++){
num[i]=r[lenr-i+1]-'0';
}
return dfs(lenr,0,0,1);
}
ll solve2(){
for(int i=1;i<=lenl;i++){
num[i]=l[lenl-i+1]-'0';
}
return dfs(lenl,0,0,1);
}
bool check(){
int sum1=0,sum2=0,pre=0;
for(int i=1;i<=lenl;i++){
sum1=(sum1*10+l[i]-'0')%m;//x%m
sum2=(sum2+(l[i]-'0')*pre)%m;// f(x)%m
pre=(pre+l[i]-'0')%m;
}
return sum1==sum2;
}
signed main(){
int _;scanf("%d",&_);
while(_--){
scanf("%s%s",l+1,r+1);
scanf("%d",&m);
lenl=strlen(l+1),lenr=strlen(r+1);
for(int i=1;i<=max(lenl,lenr);i++){
for(int j=0;j<m;j++){
for(int k=0;k<=m;k++){
dp[i][j][k]=-1;
}
}
}
base[1]=1;
for(int i=2;i<=max(lenl,lenr);i++){
base[i]=base[i-1]*10%m;
}
ll ans=((solve1()-solve2()+check())%mod+mod)%mod;
printf("%lld\n",ans);
}
return 0;
}
2019-2020 ICPC Asia Hong Kong Regional Contest J. Junior Mathematician 题解(数位dp)的更多相关文章
- 2019-2020 ICPC Asia Hong Kong Regional Contest
题解: https://files.cnblogs.com/files/clrs97/19HKEditorial-V1.zip Code:(Part) A. Axis of Symmetry #inc ...
- Asia Hong Kong Regional Contest 2019
A. Axis of Symmetry B. Binary Tree n 的奇偶性决定胜负. C. Constructing Ranches 路径上点权之和大于,极大值两倍,这是路径上点能拼出多边形的 ...
- Asia Hong Kong Regional Contest 2016
A. Colourful Graph 可以在$2n$步之内实现交换任意两个点的颜色,然后就可以构造出方案. #include <bits/stdc++.h> using namespace ...
- 2020.5.16-ICPC Central Europe Regional Contest 2019
A. ABB #include <bits/stdc++.h> using namespace std; #define PB push_back #define ZERO (1e-10) ...
- 2016 ACM/ICPC Asia Regional Shenyang Online 1007/HDU 5898 数位dp
odd-even number Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)T ...
- 2020.5.4-ICPC Pacific Northwest Regional Contest 2019
A. Radio Prize All boring tree-shaped lands are alike, while all exciting tree-shaped lands are exci ...
- Gym - 101981J The 2018 ICPC Asia Nanjing Regional Contest J.Prime Game 计数
题面 题意:1e6的数组(1<a[i]<1e6), mul (l,r) =l × (l+1) ×...× r, fac(l,r) 代表 mul(l,r) 中不同素因子的个数,求s ...
- 2019 ACM/ICPC Asia Regional shanxia D Miku and Generals (二分图黑白染色+01背包)
Miku is matchless in the world!” As everyone knows, Nakano Miku is interested in Japanese generals, ...
- 2020 ICPC Asia Taipei-Hsinchu Regional Problem H Optimization for UltraNet (二分,最小生成树,dsu计数)
题意:给你一张图,要你去边,使其成为一个边数为\(n-1\)的树,同时要求树的最小边权最大,如果最小边权最大的情况有多种,那么要求总边权最小.求生成树后的所有简单路径上的最小边权和. 题解:刚开始想写 ...
随机推荐
- 浅谈分布式共识算法raft
前言:在分布式的系统中,存在很多的节点,节点之间如何进行协作运行.高效流转.主节点挂了怎么办.如何选主.各节点之间如何保持一致,这都是不可不面对的问题,此时raft算法应运而生,专门 用来解决上述问题 ...
- Java动态代理——框架中的应用场景和基本原理
前言 之前已经用了5篇文章完整解释了java动态代理的原理,本文将会为这个系列补上最后一块拼图,展示java动态代理的使用方式和应用场景 主要分为以下4个部分 1.为什么要使用java动态代理 2.如 ...
- 简单Emacs配置
(global-set-key [f9] 'compile-file) (global-set-key [f10] 'gud-gdb) (global-set-key (kbd "C-s&q ...
- print( "Hello,NumPy!" )
print( "Hello,NumPy!" ) 学习痛苦啊,今天学,明天丢.这种天气,还是睡觉最舒服了. 咱说归说,闹归闹,但还是得学才行啊. 之前在学习的过程中一直都有记录笔记的 ...
- 专业之旅——GitHub 热点速览 Vol.45
作者:HelloGitHub-小鱼干 从入门到精通需要什么?AI-Expert-Roadmap 带你开启专业之旅,和 135k+ 高星项目 developer-roadmap 一样, AI-Exper ...
- 9、Django之模型层第四篇:进阶操作
一 QuerySet对象 1.1可切片 使用Python 的切片语法来限制查询集记录的数目 .它等同于SQL 的LIMIT 和OFFSET 子句. Entry.objects.all()[:5] # ...
- 【填坑往事】Android手机锁屏人脸解锁优化过程实录
背景 写这篇文章,主要是为了以后面试方便.因为我简历上写了,上一份工作的最大亮点是将人脸解锁的速度由1200ms优化到了600ms,所以这些内容已经回答无数遍了.但每次总觉得回答的不完整,或者说总感觉 ...
- 时间序列数据的定义,读取与指数平滑(Java)
应上头的要求,需要实现以下指数平滑进行资源调度负载的预测,那就是用我最喜欢的Java做一下吧. 引用<计量经济学导论>的一句话:时间序列数据区别于横截面数据的一个明显特点是,时间序列数据集 ...
- Docker学习笔记:Alpine镜像+Python3安装+http服务器
编写Dockerfile文件使用最新的Alpine镜像并安装Python3环境,如下: 因为python高于3.4则不会默认安装pip,需要手动安装. 试了很多其他办法都没安装上,唯有下载get-pi ...
- php判断用户设备类型
最近做的一个需求里面希望能判断用户访问页面的设备类型,根据不同的类型去加载不同的数据和页面样式.由于技术栈是使用的php,于是考虑在php层面去做这个判断. 假设主要判断的设备有平板和手机为主,分两个 ...