【201703 HNSDFZ集训 R5T1 】数字

Description

小 D 喜欢的数有这样的性质：
令 n 为正整数，S(n)为 n 的各位数字之和，令

$$D(n)=\left\{ \begin{array}{rcl} S(n) & & {S(n) < 10}\\ D(S(n)) & & {S(n) \geq 10} \end{array} \right.$$
小 D 喜欢的数一定能表示成 $x\times 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 \leq 1000000$；
对于100%的数据，$L,R\leq 10^{18}，T\leq20$。

题解

发现x被喜欢，当且仅当存在k属于[1, 9]使 $x \mod (k \times 9) = k ^ 2$。
那么可以发现x被不被喜欢只与$x \mod (k \times 9)$有关。
设$C = 9 \times 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;
}