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

题意：求区间中能被每个数位上的数整除的数的数量。

题解：数位dp
考虑直接的做法，既然能被每一位整除，即可以被lcm整除，递归每一位并记忆化，判断最终的数是否满足，但dp显然不能保存这么大的数，所以要进行取模。
根据X % a = [ X % ( b * a ) ] % a，显然模数mod = lcm(1,2,3,···,8,9) = 2520。

这样做还是会超内存，考虑离散化，用dp[][][]进行记忆的时候发现存lcm的时候浪费了很多空间，我们将mod的48个因数编号进idx[]，用索引代替lcm即可。

然后就是正常数位dp的做法。

#define _CRT_SECURE_NO_WARNINGS
#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<algorithm>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#include<fstream>
#include<set>
#include<map>
#include<sstream>
#include<iomanip>
#define ll long long
using namespace std;
const int mod = 2520;
ll dp[20][2525][50], x, y;
int t, a[20], idx[2525];
ll gcd(ll a, ll b) { return b ? gcd(b, a % b) : a; }
ll lcm(ll a, ll b) { return a / gcd(a, b) * b; }
void init() {
	int cnt = 0;
	for (int i = 1; i <= mod; i++) {
		if (mod % i == 0) idx[i] = cnt++;
	}
}
ll dfs(int pos, int sta, int lc, int limit) {
	if (pos == -1) return sta % lc == 0;
	if (!limit && dp[pos][sta][idx[lc]] != -1) return dp[pos][sta][idx[lc]];
	int up = limit ? a[pos] : 9;
	ll ans = 0;
	for (int i = 0; i <= up; i++) {
		if (i == 0) ans += dfs(pos - 1, sta * 10 % mod, lc, limit && i == a[pos]);
		else ans += dfs(pos - 1, (sta * 10 + i) % mod, lcm(lc, i), limit && i == a[pos]);
	}
	if (!limit) dp[pos][sta][idx[lc]] = ans;
	return ans;
}
ll solve(ll x) {
	int pos = 0;
	while (x) {
		a[pos++] = x % 10;
		x /= 10;
	}
	return dfs(pos - 1, 0, 1, 1);
}
int main() {
	init();
	scanf("%d", &t);
	memset(dp, -1, sizeof(dp));
	while (t--) {
		scanf("%I64d%I64d", &x, &y);
		printf("%I64d\n", solve(y) - solve(x - 1));
	}
	return 0;
}