CodeForces - 55D A - Beautiful numbers 离散化+数位dp+取余优化

数位dp

怎么说呢，不愧是CF的题目。很有思维的感觉，这题主要算法就是数位dp+离散化+mod优化。
1，题目要求被所有位上的数字都能整除，我第一个冒出来的想法就是9维dp，每一位一个数，然而是不可能的，想了很久都没有结果，最后去看了dalao的想法，一下就悟了，其实就是优化了mod，把除以9个数化成整除9个数的最小公倍数2520，然后在判断的时候return left%p==0，p是经过的数位上所有数的最小公倍数。
2，但是这样还是不够，因为即使我们解决了dp是何种状态的问题，也无法定义一个dp[20][2520][2520],分别是位数，left（余数），p，这个时候应该想到只能对p优化，以为p是最小公倍数，所以有很多空间是浪费的，我们通过离散化得到p的对应hash数组就能构造一个dp[20][50][2520]了。

//
// Created by SANZONG on 2020/10/23.
//

#include "cstdio"
#include "cstring"
#include "iostream"
#include "algorithm"
#include "map"
#include "cmath"

typedef long long ll;
using namespace std;
int dig[50];
ll dp[25][2530][50];
int hashnum[2530];

/*
ll gcd(ll a,ll b)
{
    if(b == 0)
        return a;
    else
        return gcd(b,a%b);
}
*/

int Hash()
{
    int pos = 0;
    for(int i = 1;i <= 2520;i++) {
        if (2520 % i == 0)
            hashnum[i] = ++pos;
    }
}


ll dfs(ll pos, ll left, bool limit,ll p) {
    if (pos < 0 )
        return left%p == 0;
    if (!limit && ~dp[pos][left][hashnum[p]] ) return dp[pos][left][hashnum[p]];
    int ub = limit ? dig[pos] : 9;
    ll tmp = 0;
    for (int i = 0 ; i <= ub; ++i) {
        tmp += dfs(pos-1,(left * 10 + i) %2520 ,limit && i==ub,i?(p * i/__gcd((int)p,(int)i)):p);
    }
    if (!limit) dp[pos][left][hashnum[p]] = tmp;
    return tmp;
}

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


int main() {
//    freopen("in.txt","r",stdin);
//    freopen("out.txt","w",stdout);
    int T;
    Hash();
    cin >> T;
    memset(dp,-1,sizeof(dp));
    while (T--) {
        ll a, b;
        cin >> a >> b;
        cout << solve(b) - solve(a-1) << endl;
    }
}