[COJ0528]BJOI幸运数

试题描述

输入

见“试题描述

输出

见“试题描述

输入示例

见“试题描述

输出示例

见“试题描述

数据规模及约定

见“试题描述

题解

首先想到一个比较暴力的数位 dp 方法:设 f[i][j][s1][s2] 表示前 i 位中,最高位为 j,其中奇数位的数字和为 s1,偶数位的数字和为 s2,满足这个条件的数的个数。

然而,这样的话每次询问都不得不每位处理时都枚举一下 s1 和 s2,受不了 T = 1000 的数据。

先贴一下这个代码:

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <algorithm>
using namespace std;
#define LL long long LL read() {
LL x = 0, f = 1; char c = getchar();
while(!isdigit(c)){ if(c == '-') f = -1; c = getchar(); }
while(isdigit(c)){ x = x * 10 + c - '0'; c = getchar(); }
return x * f;
} #define maxn 20
#define maxs 90
#define maxk 110
LL g[maxn][10][maxs][maxs];
int Gcd[maxs][maxs]; int gcd(int x, int y) { return !y ? x : gcd(y, x % y); } int K, num[maxn];
LL sum(LL x) {
int cnt = 0; LL tx = x;
while(x) num[++cnt] = x % 10, x /= 10;
LL ans = 0;
int s1 = 0, s2 = 0;
for(int i = cnt; i; i--) {
for(int j = 0; j < num[i]; j++)
for(int t1 = max(s1, 1); t1 < maxs; t1++)
for(int t2 = max(s2, 1); t2 < maxs; t2++)
if(Gcd[t1][t2] <= K && g[i][j][t1-s1][t2-s2]) {
ans += g[i][j][t1-s1][t2-s2];
// printf("%d %d %d %d: %d\n", i, j, t1-s1, t2-s2, g[i][j][t1-s1][t2-s2]);
}
if(i & 1) s1 += num[i]; else s2 += num[i];
}
if(Gcd[s1][s2] <= K && s1 && s2) ans++;
// printf("%lld: %lld\n", tx, ans);
return ans;
} int main() {
g[0][0][0][0] = 1;
for(int i = 0; i < maxn - 1; i++)
for(int j = 0; j <= 9; j++)
for(int s1 = 0; s1 < maxs; s1++)
for(int s2 = 0; s2 < maxs; s2++)
if(g[i][j][s1][s2]) {
for(int k = 0; k <= 9; k++) {
int S1 = s1, S2 = s2;
if(i + 1 & 1) S1 += k; else S2 += k;
g[i+1][k][S1][S2] += g[i][j][s1][s2];
}
// printf("%d %d %d %d: %d\n", i, j, s1, s2, g[i][j][s1][s2]);
}
for(int i = 0; i < maxs; i++)
for(int j = 0; j < maxs; j++)
Gcd[i][j] = gcd(i, j); int T = read();
while(T--) {
K = read(); LL l = read(), r = read();
printf("%lld\n", sum(r) - sum(l - 1));
} return 0;
}

对了记得预处理 Gcd[i][j],否则会更慢。。。

现在我们主要目的是优化每次询问所消耗的代价,核心在于不能再每次枚举 s1 和 s2 了。

我们考虑一下求 [1, x] 中答案的过程,我们是一位位确定数字的,所以当前状态可以形象地表示为一下形式:

XXXX____

其中“X”表示已填数字,“_”表示未填数字,我们发现刚才之所以要枚举 s1 和 s2 的原因在于前面填了数字之后,奇偶位上已经有一定的数字和了,所以在这里我们不妨将其设为状态,于是就变成了:

f[i][k][s1][s2] 表示还有 i 位未填,且前面已经填写的数字中,奇数位上数字和为 s1,偶数位上数字和为 s2,现在我们要填满后面那 i 个空位,并且满足填完之后整个数奇数位之和与偶数位之和的最大公约数不超过 k,那么满足上面这一系列条件的填法有多少种。

显然初始状态 f[0][1~maxk][1~maxs][1~maxs] = 1,其中 maxk 表示输入的 k 的最大值,maxs 表示可能出现的最大数字和;

转移:

1. f[i][k][s1][s2] -> f[i+1][k][s1+j][s2] (0 ≤ j < 10),i+1 为奇数

2. f[i][k][s1][s2] -> f[i+1][k][s1][s2+j] (0 ≤ j < 10),i+1 为偶数

这样询问的负担就减轻不少了。

这题还得卡卡常数。

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <algorithm>
using namespace std;
#define LL long long LL read() {
LL x = 0, f = 1; char c = getchar();
while(!isdigit(c)){ if(c == '-') f = -1; c = getchar(); }
while(isdigit(c)){ x = x * 10 + c - '0'; c = getchar(); }
return x * f;
} #define maxn 19
#define maxs 90
#define maxk 101
LL f[maxn][maxk][maxs][maxs];
int Gcd[maxs][maxs]; int gcd(int x, int y) { return !y ? x : gcd(y, x % y); } int num[maxn];
LL sum(LL x, int k) {
int cnt = 0;
while(x) num[++cnt] = x % 10, x /= 10;
int s1 = 0, s2 = 0;
LL ans = 0;
for(int i = cnt; i; i--) {
for(int j = 0; j < num[i]; j++) {
ans += f[i-1][k][s1][s2];
if(i & 1) s1++; else s2++;
}
}
ans += f[0][k][s1][s2];
return ans;
} int main() {
for(int i = 0; i < maxs; i++)
for(int j = 0; j < maxs; j++) Gcd[i][j] = gcd(i, j);
for(int k = 1; k < maxk; k++)
for(int i = 1; i < maxs; i++)
for(int j = 1; j < maxs; j++)
if(Gcd[i][j] <= k) f[0][k][i][j]++;
for(int i = 0; i < maxn - 1; i++)
for(int k = 1; k < maxk; k++) {
int mxs = (maxn - i >> 1) * 9;
for(int s1 = 0; s1 <= mxs; s1++)
for(int s2 = 0; s2 <= mxs; s2++) if(f[i][k][s1][s2])
for(int j = 0; j <= 9; j++) {
int t1 = s1, t2 = s2;
if(i + 1 & 1) t1 -= j; else t2 -= j;
if(t1 < 0 || t2 < 0) break;
f[i+1][k][t1][t2] += f[i][k][s1][s2];
}
} int T = read();
while(T--) {
int k = read(); LL l = read(), r = read();
printf("%lld\n", sum(r, k) - sum(l - 1, k));
} return 0;
}

[COJ0528]BJOI幸运数的更多相关文章

  1. 京东2017校园招聘笔试题 【第K个幸运数】

    题目描述 4和7是两个幸运数字,我们定义,十进制表示中,每一位只有4和7两个数的正整数都是幸运数字. 前几个幸运数字为:4,7,44,47,74,77,444,447... 现在输入一个数字K,输出第 ...

  2. [51NOD1230]幸运数(数位DP)

    题目链接:http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1230 dp(l,s,ss)表示长度为l的数各位和为s,各位平方 ...

  3. Problem 1007 幸运数 线段树成段更新

    题目链接: 题目 Problem 1007 幸运数 Time Limit: 2000 mSec Memory Limit : 131072 KB 问题描述 皮特的幸运数是2和5.只由幸运数字2和5组成 ...

  4. pongo英雄会-幸运数题解

    显然我们只要知道1~x范围有多少幸运数(用f(x)表示),lucky(x,y)=f(y)-f(x-1). 解法1. 计算排列数 由于y<=1000000000这个规模,我们不能暴力验证每个数是否 ...

  5. Python 计算当真因子个数为偶数个时为幸运数,计算区间内幸运数之和

    晚饭后朋友发来个问题,正好无事做,动手写了一下 若一个正整数有偶数个不同的真因子,则称该数为幸运数.如4含有2个真因子为 1 和 2 .故4是幸运数.求[2,100]之间的全部幸运数之和. 常规思路 ...

  6. 蓝桥杯  历届试题 幸运数  dfs

    历届试题 幸运数 时间限制:1.0s   内存限制:256.0MB 问题描述 幸运数是波兰数学家乌拉姆命名的.它采用与生成素数类似的"筛法"生成 . 首先从1开始写出自然数1,2, ...

  7. C. cltt的幸运数LCAtarjan

    /*C: cltt的幸运数 Time Limit: 1 s      Memory Limit: 128 MB Submit Problem Description 一棵树有n个节点,共m次查询,查询 ...

  8. C. cltt的幸运数LCAdfs

    /*C: cltt的幸运数 Time Limit: 1 s      Memory Limit: 128 MB Submit Problem Description 一棵树有n个节点,共m次查询,查询 ...

  9. [蓝桥杯]PREV-10.历届试题_幸运数

    问题描述 幸运数是波兰数学家乌拉姆命名的.它采用与生成素数类似的“筛法”生成 . 首先从1开始写出自然数1,,,,,,.... 就是第一个幸运数. 我们从2这个数开始.把所有序号能被2整除的项删除,变 ...

随机推荐

  1. for 循环的关键字 break和continue

    for(int i=0;i<100;i++){ if(i==50){ continue:// 跳出当前循环,执行下面的循环,就是说,当i=50的时候,跳出循环,从i=51开始继续循环 //业务逻 ...

  2. SQL Server编程(03)自定义存储过程

    存储过程是一组预编译的SQL语句,它可以包含数据操纵语句.变量.逻辑控制语句等. 存储过程允许带参数: 输入参数:可以在调用时向存储过程传递参数,此类参数可用来向存储过程中传入值(可以有默认值) 输出 ...

  3. UIProgressView(进度条控件)

    UIProgressView *pr=[[UIProgressView alloc]init]; pr.frame=CGRectMake(150.0, 190.0, 130.0, 30.0);//进度 ...

  4. 根据wsdl文件生成webservice 的.cs文件 及 生成dll C#调用

               Visual Studio 2013->Visual Studio Tools->VS2013 开发人员命令提示 命令行输入 wsdl E:\WS.wsdl /out ...

  5. Hadoop-2.7.3 问题收集

    问题1:在DFS Lcation 上不能多文件进行操作: 在Hadoop上的每个节点上修改该文件       conf/mapred-site.xml 增加: <property> < ...

  6. C#中导入Win32 API函数

    C#中导入Win32 API的方法: 1.引用命名空间 using System.Net.Security; using System.Runtime.InteropServices; 2. [Dll ...

  7. JavaScript中变量和函数声明的提升

    现象: 1.在JavaScript中变量和函数的声明会提升到最顶部执行. 2.函数的提升高于变量的提升. 3.函数内部如果用var声明了相同名称的外部变量,函数将不再向上寻找. 4.匿名函数不会提升. ...

  8. 用IIS配置反向代理

    http://my.oschina.net/tanyixiu/blog/123832 目标服务器:targetServer 配置反向代理的服务器:reveseProxServer 1.确定最终访问的网 ...

  9. Exception Handling in ASP.NET Web API webapi异常处理

    原文:http://www.asp.net/web-api/overview/error-handling/exception-handling This article describes erro ...

  10. maven可选依赖(Optional Dependencies)和依赖排除(Dependency Exclusions)

    我们知道,maven的依赖关系是有传递性的.如:A-->B,B-->C.但有时候,项目A可能不是必需依赖C,因此需要在项目A中排除对A的依赖.在maven的依赖管理中,有两种方式可以对依赖 ...