Description

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


Solution

设 \(f[i][j][k][q]\) 表示 枚举到第 i 位,当前数字为 j ,各个数位上的数字和为 k,原数模 k 的余数为 q

可以发现,从 i - 1 位推到第 i 位并不好推,所以可以转变一下

用第 i 位推第 i + 1 位

显然有

\[f[i+1][s][j+s][(k+s\times 10^i)\%x]=\sum_{t=0}^9f[i][t][j][k]
\]

约束条件为 \(i\in [1,len]\),\(j\in [0,\min\{i\times 9,x\}]\),\(s\in[0,9]\),\(k\in[0,x)\)

其中 \(i\) 是枚举位数,\(s\) 是下一位上的数,\(t\) 是当前位上的数,\(x\) 是总数位和,\(k\) 是余数,\(j\) 是当前累计的数位和

然后对于答案的统计:

  • 对于位数小于当前位的,\(ans+=f_{i,10,x,0}\) \(−\) \(f_{i,0,x,0}\)
  • 对于位数相同的,计算小于最高位的 \(ans+=f_{num_i,s,x,0}\)
  • 然后对于当前不同的一位进行讨论即可,计算出后面所需要的余数,然后加和

\(num_i\) 是 \(i\) 的最高位,\(x\) 是要求的位数和

答案就是所有位数和的 \(ans\) 之和, 同时满足满足 \([l,r]\)


Code

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<queue>
#include<cmath>
#include<cstring>
#define int long long using namespace std; int A,B;
int f[20][11][172][172];
int highest[20],a[20][2];
int Mod,Pow[20]; inline int read(){
int s=0,w=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-') w=-1;ch=getchar();}
while(ch>='0' &&ch<='9') s=(s<<1)+(s<<3)+ch-'0',ch=getchar();
return w*s;
} inline int Min(int x,int y){return x<y?x:y;}
inline void F(int x,int y){highest[y]=0;while(x){a[++highest[y]][y]=x%10;x/=10;}} inline void init(int x,int len){
memset(f,0,sizeof f);
for(register int i=0;i<=9;i++) f[1][i][i][i%x]++,f[1][10][i][i%x]++;
for(register int i=1;i<len;i++)
for(register int j=0;j<=Min(i*9,x);j++)
for(register int k=0;k<x;k++){
if(f[i][10][j][k])for(register int s=0;s<=9;s++){
f[i+1][s][j+s][(k+s*Pow[i])%x]+=f[i][10][j][k];
f[i+1][10][j+s][(k+s*Pow[i])%x]+=f[i][10][j][k];
}
}
Mod=x;
} inline int solve(int s,int x){
int ans=0,num=a[highest[s]][s],now=num;
for(register int i=1;i<highest[s];i++) ans+=f[i][10][x][0]-f[i][0][x][0];//处理最高位之前
for(register int i=1;i<a[highest[s]][s];i++) ans+=f[highest[s]][i][x][0];//处理最高位
if(x<now) return ans;
for(register int i=highest[s]-1;i;i--){
int tot=(x-num*Pow[i]%x)%Mod;
for(register int j=0;j<a[i][s];j++)ans+=f[i][j][x-now][tot];//操作最高位
now+=a[i][s];num=(num*10+a[i][s])%Mod;
if(x<now) break;
}
return ans;
} signed main(){
int ans=0;A=read();B=read();F(A,0),F(B+1,1);Pow[0]=1;
for(register int i=1;i<=18;i++) Pow[i]=Pow[i-1]*10;
for(register int i=1;i<=highest[1]*9;i++){init(i,highest[1]);ans+=solve(1,i)-solve(0,i);}
printf("%lld",ans);
return 0;
}

洛谷P4127的更多相关文章

  1. 洛谷 P4127 [AHOI2009]同类分布 解题报告

    P4127 [AHOI2009]同类分布 题目描述 给出两个数\(a,b\),求出\([a,b]\)中各位数字之和能整除原数的数的个数. 说明 对于所有的数据,\(1 ≤ a ≤ b ≤ 10^{18 ...

  2. 洛谷 P4127 [AHOI2009]同类分布

    题意简述 求l~r之间各位数字之和能整除原数的数的个数. 题解思路 数位DP 代码 #include <cstdio> #include <cstring> typedef l ...

  3. 洛谷P4127同类分布

    传送 我们要在dfs的板子里记录哪些量呢?当前填的所有数的和sum?当前填的数构成的数值all? sum可以留下,数值就扔掉叭.数值最大是1e18,要是留下,在g数组里有一维的大小是1e18.也许可以 ...

  4. bzoj1799(洛谷4127)同类分布(月之谜)

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1799 https://www.luogu.org/problemnew/show/P4127 ...

  5. 洛谷1640 bzoj1854游戏 匈牙利就是又短又快

    bzoj炸了,靠离线版题目做了两道(过过样例什么的还是轻松的)但是交不了,正巧洛谷有个"大牛分站",就转回洛谷做题了 水题先行,一道傻逼匈牙利 其实本来的思路是搜索然后发现写出来类 ...

  6. 洛谷P1352 codevs1380 没有上司的舞会——S.B.S.

    没有上司的舞会  时间限制: 1 s  空间限制: 128000 KB  题目等级 : 钻石 Diamond       题目描述 Description Ural大学有N个职员,编号为1~N.他们有 ...

  7. 洛谷P1108 低价购买[DP | LIS方案数]

    题目描述 “低价购买”这条建议是在奶牛股票市场取得成功的一半规则.要想被认为是伟大的投资者,你必须遵循以下的问题建议:“低价购买:再低价购买”.每次你购买一支股票,你必须用低于你上次购买它的价格购买它 ...

  8. 洛谷 P2701 [USACO5.3]巨大的牛棚Big Barn Label:二维数组前缀和 你够了 这次我用DP

    题目背景 (USACO 5.3.4) 题目描述 农夫约翰想要在他的正方形农场上建造一座正方形大牛棚.他讨厌在他的农场中砍树,想找一个能够让他在空旷无树的地方修建牛棚的地方.我们假定,他的农场划分成 N ...

  9. 洛谷P1710 地铁涨价

    P1710 地铁涨价 51通过 339提交 题目提供者洛谷OnlineJudge 标签O2优化云端评测2 难度提高+/省选- 提交  讨论  题解 最新讨论 求教:为什么只有40分 数组大小一定要开够 ...

随机推荐

  1. 记 CentOS 服务器上安装 neo4j 图数据库及本地访问

    下载 去官网下载压缩包放到服务器上.地址为neo4j 下载中心,我这里选择的是 Neo4j 3.5.25 (tar).具体如何做呢?我这里使用的是土方法,即先压缩包下载到本地电脑(win 10系统), ...

  2. env: python3: No such file or directory exit status 127 FER ESPectro Core 编译出错

    苹果电脑安装了Arduino,布置ESP8266开发环境,编译程序过程中出现错误: env: python3:No such file or directoryexit status 127为开发板 ...

  3. 在 CAP 中使用 AOP ( Castle.DynamicProxy )

    简介 本篇文章主要介绍如何在 CAP 中集成使用 Castle.DynamicProxy,Castle DynamicProxy 是一个用于在运行时动态生成轻量级.NET代理的库.代理对象允许在不修改 ...

  4. java:原子类的CAS

    当一个处理器想要更新某个变量的值时,向总线发出LOCK#信号,此时其他处理器的对该变量的操作请求将被阻塞,发出锁定信号的处理器将独占共享内存,于是更新就是原子性的了. 1.compareAndSet- ...

  5. AI算法测评事项

    前言 注:大概2017年-2018年国内人工智能热度达到顶峰,随后热度开始逐渐减少.2018年前人工智能被投资界.学术界.工业界和媒体炒的特别热,各大企业都想尝试一下深度学习技术在业务场景的应用.试水 ...

  6. Nginx解决前端访问资源跨域问题

    被前端跨域问题折磨快2天后,终于用ngnx的方式解决了,所以在此总结下. 该篇只探讨如何用Ngnx解决跨域问题,对于原理不作讨论. 1.首先介绍Windows环境下Nignx的相关命令操作 nginx ...

  7. 【项目实践】一文带你搞定Spring Security + JWT

    以项目驱动学习,以实践检验真知 前言 关于认证和授权,R之前已经写了两篇文章: [项目实践]在用安全框架前,我想先让你手撸一个登陆认证 [项目实践]一文带你搞定页面权限.按钮权限以及数据权限 在这两篇 ...

  8. WixVersionControl Wix项目版本控制

    原文链接:https://www.swack.cn/wiki/001565675133949eff0d3d5a51f48288cf6d8248905e28f000/001569821278313e6b ...

  9. mmal商城商品模块总结

    学习目标 FTP服务器的对接 SpringMVC文件上传 流读取properties配置文件 抽象POJO.BO.VO对象之间的转换关系及解决思路 joda-time快速入门 静态代码块 mybati ...

  10. Netty源码解析 -- FastThreadLocal与HashedWheelTimer

    Netty源码分析系列文章已接近尾声,本文再来分析Netty中两个常见组件:FastThreadLoca与HashedWheelTimer. 源码分析基于Netty 4.1.52 FastThread ...