题意:求区间中能被每个数位上的数整除的数的数量。
题解:数位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;
}