bzoj 3679: 数字之积
Description
一个数x各个数位上的数之积记为\(f(x)\) <不含前导零>
求[L,R)中满足\(0<f(x)<=n\)的数的个数
solution
最后\(f(x)\)可以拆分成2,3,5,7的乘积,我们就将 \(2,3,5,7\) 压进状态,然后就是基础的数位DP,分是否严格小于两种状态转移即可
具体实现需要一些技巧:
预处理出每一个数含有 \(2,3,5,7\)的个数
预处理出 \(2,3,5,7\) 的幂,方便剪枝
注意数字不能含有 \(0\),我们每DP一位,要新加入 \([1,9]\) 的状态,即前导零的情况
还有一种解法是用 \(map\) 压乘积,网上大部分都是这么做的,也能通过,且简洁很多
tips:代码实现比较简单,但我已不想再多看一眼我的代码....
#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#define RG register
#define il inline
#define iter iterator
#define Max(a,b) ((a)>(b)?(a):(b))
#define Min(a,b) ((a)<(b)?(a):(b))
using namespace std;
const int N=36;
typedef long long ll;
int m,lim[6],pri[6]={0,2,3,5,7};ll f[21][2][33][21][15][13];
char S[22];
ll m5[N],m2[N],m7[N],m3[N];int v[15][6];
ll solve(int *a,int n){
if(n==1){
int ret=0;
for(int i=1;i<a[1];i++)ret+=(i<=m);
return ret;
}
ll x,mu;memset(f,0,sizeof(f));
for(int i=1;i<n;i++){
if(i==1)for(int j=1;j<=a[i];j++)
f[i][j<a[i]][v[j][1]][v[j][2]][v[j][3]][v[j][4]]++;
else for(int j=1;j<=9;j++)
f[i][1][v[j][1]][v[j][2]][v[j][3]][v[j][4]]++;
for(int j=0;j<=lim[1];j++){
if(m2[j]<=m)
for(int k=0;k<=lim[2];k++){
if(m2[j]*m3[k]<=m)
for(int g=0;g<=lim[3];g++){
if(m2[j]*m3[k]*m5[g]<=m)
for(int b=0;b<=lim[4];b++){
mu=m2[j]*m3[k]*m5[g]*m7[b];
if(mu>m)break;
x=f[i][1][j][k][g][b];
if(x)
for(int d=1;d<=9 && mu*d<=m;d++)
f[i+1][1][j+v[d][1]][k+v[d][2]][g+v[d][3]][b+v[d][4]]+=x;
x=f[i][0][j][k][g][b];if(!x)continue;
for(int d=1;d<=a[i+1] && mu*d<=m;d++)
f[i+1][d<a[i+1]][j+v[d][1]][k+v[d][2]][g+v[d][3]][b+v[d][4]]+=x;
}
}
}
}
}
for(int i=n;i<=n;i++)
for(int j=1;j<=9;j++)
f[i][1][v[j][1]][v[j][2]][v[j][3]][v[j][4]]++;
ll ret=0;
for(int j=0;j<=lim[1];j++)
if(m2[j]<=m)
for(int k=0;k<=lim[2];k++)
if(m3[k]*m2[j]<=m)
for(int g=0;g<=lim[3];g++)
if(m2[j]*m3[k]*m5[g]<=m)
for(int b=0;b<=lim[4];b++){
if(m2[j]*m3[k]*m5[g]*m7[b]>m)break;
ret+=f[n][1][j][k][g][b];
}
return ret;
}
int s1[22],s2[22],l1,l2;
void work()
{
scanf("%d",&m);
for(int i=1;i<=4;i++)lim[i]=log(m+1)/log(pri[i]);
m2[0]=1;for(int i=1;i<=lim[1];i++)m2[i]=m2[i-1]<<1;
m3[0]=1;for(int i=1;i<=lim[2];i++)m3[i]=m3[i-1]*3;
m5[0]=1;for(int i=1;i<=lim[3];i++)m5[i]=m5[i-1]*5;
m7[0]=1;for(int i=1;i<=lim[4];i++)m7[i]=m7[i-1]*7;
v[2][1]=1;v[3][2]=1;v[4][1]=2;v[6][1]=1;v[6][2]=1;
v[5][3]=1;v[7][4]=1;v[8][1]=3;v[9][2]=2;
scanf("%s",S+1);l1=strlen(S+1);
for(int i=1;i<=l1;i++)s1[i]=S[i]-'0';
scanf("%s",S+1);l2=strlen(S+1);
for(int i=1;i<=l2;i++)s2[i]=S[i]-'0';
ll ans=solve(s2,l2)-solve(s1,l1);
printf("%lld\n",ans);
}
int main()
{
freopen("pp.in","r",stdin);
freopen("pp.out","w",stdout);
work();
return 0;
}
bzoj 3679: 数字之积的更多相关文章
- BZOJ 3679 数字之积 数位DP
思路:数位DP 提交:\(2\)次 错因:进行下一层\(dfs\)时的状态转移出错 题解: 还是记忆化搜索就行,但是要用\(map\)记忆化. 见代码 #include<cstdio> # ...
- 【BZOJ3679】数字之积 DFS+DP
[BZOJ3679]数字之积 Description 一个数x各个数位上的数之积记为f(x) <不含前导零>求[L,R)中满足0<f(x)<=n的数的个数 Input 第一行一 ...
- Bzoj 2154: Crash的数字表格(积性函数)
2154: Crash的数字表格 Time Limit: 20 Sec Memory Limit: 259 MB Description 今天的数学课上,Crash小朋友学习了最小公倍数(Least ...
- BZOJ 4816 数字表格
首先是惯例的吐槽.SDOI题目名称是一个循环,题目内容也是一个循环,基本上过几年就把之前的题目换成另一个名字出出来,喜大普奔亦可赛艇.学长说考SDOI可以考出联赛分数,%%%. 下面放解题报告.并不喜 ...
- [BZOJ 1833] 数字计数
Link: BZOJ 1833 传送门 Solution: 比较明显的数位DP 先预处理出1~9和包括前导0的0的个数:$pre[i]=pre[i-1]*10+10^{digit-1}$ (可以分为首 ...
- BZOJ3679 : 数字之积
设f[i][p2][p3][p5][p7][j][k]表示前i位,2,3,5,7的次数,前i位是否等于x,是否有数字的方案数 然后数位DP即可,ans=cal(r)-cal(l) #include&l ...
- BZOJ 1049 数字序列
Description 现在我们有一个长度为n的整数序列A.但是它太不好看了,于是我们希望把它变成一个单调严格上升的序列.但是不希望改变过多的数,也不希望改变的幅度太大. Input 第一行包含一个数 ...
- BZOJ 1833 数字计数 数位DP
题目链接 做的第一道数位DP题,听说是最基础的模板题,但还是花了好长时间才写出来..... 想深入了解下数位DP的请点这里 先设dp数组dp[i][j][k]表示数位是i,以j开头的数k出现的次数 有 ...
- BZOJ 1049 数字序列(LIS)
题目链接:http://www.lydsy.com:808/JudgeOnline/problem.php?id=1049 题意:给出一个数列A,要求:(1)修改最少的数字使得数列严格递增:(2)在( ...
随机推荐
- 使用SecureCRTP 连接生产环境的web服务器和数据库服务器
一.使用SecureCRTP 连接生产环境的web服务器 首先,需要知道以下参数信息: 1.web服务器的ip地址 2.服务器的端口号 3.会话连接的用户名和密码 4.服务器的用户名 ...
- 一个CSS简单入门网站
讲的知识简单明了,很实用: http://zh.learnlayout.com/
- slf4j 与 log4j2 实战讲解与日志分割
这两天搭建项目的时候用到log4j2在这里把自己的问题与了解拿出来与大家分享一下. 1.为什我要用 因为,使用slf4j可以很好的保证我们的日志系统具有良好的兼容性,兼容当前常见几种日志系统,而使用l ...
- 07-TypeScript的For循环
在传统的JavaScript中,关于循环,可以有两种方式,一种是forEach,一种是for. forEach的用法如下: var sarr=[1,2,3,4]; sarr.desc="he ...
- 从同步阻塞聊到Java三种IO方式
本文总结自 https://zhuanlan.zhihu.com/p/34408883, https://www.zhihu.com/question/19732473中愚抄的回答, http://b ...
- Django 使用celery任务队列的配置
celery 情景:用户发起request,并等待response返回.在本些views中,可能需要执行一段耗时的程序,那么用户就会等待很长时间,造成不好的用户体验,比如发送邮件.手机验证码等. 使用 ...
- spring mvc跨域(ajax post json)--filter方案
@RequestMapping(value = "/login.do",method = RequestMethod.POST) public Message login(Http ...
- Python 爬虫性能相关
性能相关 在编写爬虫时,性能的消耗主要在IO请求中,当单进程单线程模式下请求URL时必然会引起等待,从而使得请求整体变慢. import requests def fetch_async(url): ...
- Spring之IOC(一)
Spring之IOC(一) 对于Spring的基本概念在这里就不在赘述了. Spring的最核心的两部分则是AOP(面向切面编程)和IOC(控制反转).本篇文章主要讲我对IOC的理解. IOC:即In ...
- c语言清屏、等待、随机函数
清屏函数 #include<conio.h> system("CLS");或system(cls); 等待函数 #include<windows.h> S ...