CF55D Beautiful numbers （数位dp+离散化）

题目大意：

   一个数能被他的各个数位相除，则是一个$beautifulnumber$，$T$组数据，询问每个区间有多少个$beautifulnumber$

解题思路：（来自：https://blog.csdn.net/u013081425/article/details/38908893）

如果被各个数位整除，那么一定被其最小公倍数整除

因此记忆化搜索中的参数除了pos(当前位)和up（是否达到上界），有一个prelcm表示前面的数的最小公倍数，判断这个数是否是Beautiful Numbers，还要有一个参数sum表示前面数，但是这个数太大，需要缩小它的范围。

因为1~9的最小公倍数，为2520，所以将前面的数sum对2520取余来保存，进而缩小范围

why?
可以发现所有个位数的最小公倍数是2520，假设当前的Beautiful Numbers是x，
那么 x % lcm{dig[i]} = 0,
又 2520%lcm{dig[i]} = 0，
那么x%2520%lcm{ dig[i] } = 0，x范围由9*10^18变为2520。

此时就可以设$dp[100][2520][2520]$，来表示当前状态，可是内存超了，那该咋玩？

可以发现真正需要记录的公倍数只有几十个远远到不了2520个所以我们只需要将其离散化即可进而变为dp[100][2525][50]

之后按照模板来套即可

AC代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 100;
const ll Mod = 2520;
ll T, A, B, digit[maxn], dp[maxn][Mod + 5][50], mp[Mod + 5];
ll dfs(int pos, bool limit, int div, int sum) {
	if (!pos) return (sum % div == 0);
	if (!limit && dp[pos][sum][mp[div]] != -1) return dp[pos][sum][mp[div]];
	int up = limit ? digit[pos] : 9;
	ll ans = 0;
	for (int i = 0; i <= up; i++) {
		int div_tmp = div;
		if (i) div_tmp = div_tmp * i / __gcd(div_tmp, i); //求最小公倍数
		ans += dfs(pos - 1, limit && i == digit[pos], div_tmp, (sum * 10 + i) % Mod);
	}
	if (!limit) dp[pos][sum][mp[div]] = ans;
	return ans;
}
ll solve(ll x) {
	int cnt = 0;
	while (x) digit[++cnt] = x % 10, x /= 10;
	return dfs(cnt, true, 1, 0); 
}

int main() {
	int tmp = 0;
	for (int i = 1; i <= Mod; i++) {
		if (Mod % i == 0) mp[i] = ++tmp;//离散化
	}
	memset(dp, -1, sizeof(dp));
	scanf("%lld", &T);
	while (T--) {
		scanf("%I64d%I64d", &A, &B);
		//cout << solve(B); //<< " " << solve(A - 1) << endl; 
		printf("%I64d\n", solve(B) - solve(A - 1));
	}
}