【数位dp && 优化数组范围】CodeForces - 55D Beautiful numbers

Problem Description

输入T组测试数据，每组测试数据输入一组区间范围。让你求区间范围内有多少个数是美丽的数。
美丽的数：数 能被 自身各位数整除（不包括0）那么这个数就是美丽的数。
例如：
120 能被1整除 同时 也能被2整除，所以它是一个美丽的数。

思路：

这题得用到的前置技能，大整数取模
123%5 = ( (1*10+2)%5 )*10+3%5
这题用到的一些小技巧都是为了优化记忆化数组的大小。数位dp的时间复杂度，跟数组大小有着密切关系。
1-9的最小公倍数2520。如果一个数能够被2520整除，那么它一定是美丽的数。
我们搜索到一个数的时候，我们都判断该数是不是美丽的数，我们肯定得求出这个数所有位的lcm。
lcm最大是2520，但是1-9能组成的lcm个数 不重复的只有50个不到（可以打表验证）。
有了上面这些 你就可以确定dp数组的范围了。然后就是实现它了

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int MOD = 2520;
ll dp[50][2550][50], a[50];//核心点就是如何去压缩dp数组的范围
ll Hash[2550];
ll _gcd(ll x, ll y)
{
    return y ? _gcd(y, x%y) : x;
}
ll dfs(int pos, int num, int lcm, int limit)
{
    //printf("%d\n", num);
    if(pos < 0) return num%lcm?0:1;
    if(!limit && dp[pos][num][Hash[lcm]] != -1) return dp[pos][num][Hash[lcm]];
    int up = limit ? a[pos] : 9;
    ll res = 0;
    for(int i = 0; i <= up; i++)
    {
        if(!i) res += dfs(pos-1, (num*10+i)%MOD, lcm, limit && i == a[pos]);
        else res += dfs(pos-1, (num*10+i)%MOD, lcm*i/_gcd(lcm, i), limit && i == a[pos]);
    }
    if(!limit) return dp[pos][num][Hash[lcm]] = res;
    return res;
}
ll solve(ll x)
{
    int pos = 0;
    while(x)
    {
        a[pos++] = x%10;
        x /= 10;
    }
    return dfs(pos-1, 0, 1, 1);
}
void init()
{
    memset(dp, -1, sizeof(dp));
    int cnt = 0;
    for(int i = 1; i <= MOD; i++)//如果MOD能被i整除，那么i肯定是其中一个能组成的lcm
        if(MOD % i == 0) Hash[i] = cnt++;
}
int main()
{
    int T;
    cin >> T;
    init();
    ll lt, rt;
    while(T--)
    {
        cin >> lt >> rt;
        cout << solve(rt) - solve(lt-1) << endl;
    }
    return 0;
}