CodeForces - 55D Beautiful numbers（数位dp 离散化）

题目地址：点击打开链接

题意：统计某段区间内满足它所有非零数位能整除这个数本身的数的个数。

思路：要求这个数能被每一位整除，也就是这个数能被所有位数的最小公倍数整除。2-9的最小公倍数是2520，可知

dp[pos][preMod][lcm]需要的开的大小dp[20][2520][2520], 显然是开不下的而且会超时。所以需要想办法优化。 我们可

以发现第三维，2-9的最小公倍数最大为2520，但是组合形成的最小公倍数却不多，只有不到50个，所以我们可以对

这一维进行离散化，能形成的每个最小公倍数对应一个哈希值。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn = 2530;
typedef long long ll;
ll dp[25][50][maxn];
int a[25], has[maxn];

void init()
{
    memset(dp, -1, sizeof(dp));
    int k = 1;
    for(int i = 1; i <= 2520; i++)
        if(2520%i == 0)
            has[i] = k++;
}

int __lcm(int x, int y)
{
    return x/__gcd(x, y)*y;
}

ll dfs(int pos, int lcm, int mod, int limit)
{
    if(pos < 0) return (mod%lcm==0);
    if(!limit && dp[pos][has[lcm]][mod] != -1) return dp[pos][has[lcm]][mod];
    int up = limit ? a[pos] : 9;
    ll tmp = 0;
    for(int i = 0; i <= up; i++)
    {
        if(!i) tmp += dfs(pos-1, lcm, mod*10%2520, limit && i==a[pos]);
        else tmp+= dfs(pos-1, __lcm(lcm, i), (mod*10+i)%2520, limit && i==a[pos]);
    }
    if(!limit) dp[pos][has[lcm]][mod] = tmp;
    return tmp;
}

ll solve(ll x)
{
    int pos = 0;
    while(x)
    {
        a[pos++] = x%10;
        x /= 10;
    }
    return dfs(pos-1, has[1], 0, 1);
}

int main(void)
{
    int t;
    init();
    cin >> t;
    while(t--)
    {
        ll l, r;
        scanf("%I64d%I64d", &l, &r);
        printf("%I64d\n", solve(r)-solve(l-1));
    }
    return 0;
}