Codeforces55D_Beautiful numbers_记忆化搜索版数位DP

题意

求 [a, b] 中满足以下条件的数字的个数：
这个数字能被它各个位上的非零数字整除。

思路

数位dp。

dp第一维参数自然是位数。
这个数字能整除各个位上的数字，也就是能整除它们的最小公倍数，第二维参数便是前 i 维数字的最小公倍数。
因为要考察这个数字能否整除这个最小公倍数，因此第三位参数是这个数的值。

这样就写出了状态 dp[i][j][k]，但是后两维的规模是在太大，MLE 成狗。

考虑到 1 ~ 9 的最小公倍数是 2520 ，因此第三位只需要维护 %2520 下的值即可。
于是此时的规模变成 dp[20][2520][2520], 还是太大。

又考虑到，1 ~ 9 的最小公倍数并不是 2520 个，实际上这个数字只有 48。因此可以利用哈希将第二维离散化，于是规模变成了 dp[20][50][2520]。至此这个题解决。

题目链接

Vjudge链接
https://vjudge.net/contest/169049#problem/A

AC代码

#include<cstdio>
#include<iostream>
#include<cstring>

using namespace std;

typedef long long LL;

const int maxlcm = 2520;

int t;
LL a, b;
int A[25];
int Hash[2525];
LL dp[25][55][2525];

int gcd(int a, int b)
{
    return b == 0 ? a : gcd(b, a % b);
}

int Lcm(int a, int b)
{
    return a / gcd(a, b) * b;
}

LL dfs(int p, int lcm, int num, bool ceil)
{
    if(p == 0) return num % lcm == 0;
    if(!ceil && dp[p][Hash[lcm]][num] != -1) return dp[p][Hash[lcm]][num];

    LL res = 0;
    int ub = ceil ? A[p] : 9;

    for(int i= 0; i<= ub; i++)
    {
        int new_lcm = i ? Lcm(lcm, i) : lcm;
        int new_num = (num * 10 + i) % maxlcm;
        res += dfs(p - 1, new_lcm, new_num, ceil && i == ub);
    }

    if(!ceil) dp[p][Hash[lcm]][num] = res;

    return res;
}

LL solve(LL x)
{
    int len = 0;
    while(x)
    {
        A[++ len] = x % 10;
        x /= 10;
    }

    return dfs(len, 1, 0, true);
}

int main()
{
    int len = 0;
    for(int i= 1; i<= maxlcm; i++)
        if(maxlcm % i == 0) Hash[i] = len ++;

    memset(dp, -1, sizeof dp);

    scanf("%d", &t);
    while(t --)
    {
        scanf("%lld %lld", &a, &b);
        cout << solve(b) - solve(a - 1) << endl;
    }

    return 0;
}