Description
小 D 喜欢的数有这样的性质:
令 n 为正整数,S(n)为 n 的各位数字之和,令
D(n)={S(n)D(S(n))S(n)<10S(n)≥10
小 D 喜欢的数一定能表示成 x×D(x) 这种形式(即若一个数 A 是被喜欢的,则存在一个数 x,使得 A=x*D(x))。
小 D 想知道在区间[L, R]中,有多少个数是他喜欢的。
Input
第一行一个整数 T,表示数据组数。
以下每一行两个数 L、R(保证区间合法),代表询问[L, R]。
Output
输出 T 行,每行一个数,表示在这个区间内小 D 喜欢的数出现了多少次。
你的输出当且仅当和标准输出一样才能得该测试点满分。
Sample Input
3
1 5
3 9
8 8
Sample Output
2
2
0
Hint
【样例说明】
在[1,5]中,总共有 5 个数,但只有 、 两个小 D 喜欢的数。
同理,[3,9]中,有 4、9 两个数;
[8,8]中,没有小 D 喜欢的数。
对于30%的数据,
L,R≤1000000
;
对于100%的数据,
L,R≤1018,T≤20
。
题解
发现x被喜欢,当且仅当存在k属于[1, 9]使
xmod(k×9)=k2
。
那么可以发现x被不被喜欢只与
xmod(k×9)
有关。
设
C=9×LCM(1,2,3,4,5,6,7,8,9)=22680
。
必然x 与 x+22680一定要么同时被喜欢,要么都不被喜欢。
这样就找到了循环节,预处理出22680以内的所有数被喜欢的情况,对于每个询问就可以
O(1)
计算了。
#include<ctime>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;
inline 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;
}
const int C = 22680;
int t;
ll x, y, s[C+1];
void init(){
t = read();
int tmp;
for(int i = 1; i <= C; i++){
tmp = i * ((i-1) % 9 + 1);
if(tmp <= C) s[tmp] = 1;
}
for(int i = 1; i <= C; i++)
s[i] += s[i-1];
}
ll cal(ll x){
return x/C * s[C] + s[x%C];
}
void work(){
while(t--){
x = read(); y = read();
printf("%lld\n", cal(y) - cal(x - 1));
}
}
int main(){
freopen("num.in", "r", stdin);
freopen("num.out", "w", stdout);
init();
work();
//printf("%.6lf\n", (double)clock() / CLOCKS_PER_SEC);
return 0;
}