Problem Description
输入T组测试数据,每组测试数据输入一组区间范围。让你求区间范围内有多少个数是美丽的数。
美丽的数:数 能被 自身各位数整除(不包括0)那么这个数就是美丽的数。
例如:
120 能被1整除 同时 也能被2整除,所以它是一个美丽的数。
思路:
这题得用到的前置技能,大整数取模
123%5 = ( (1*10+2)%5 )*10+3%5
这题用到的一些小技巧都是为了优化记忆化数组的大小。数位dp的时间复杂度,跟数组大小有着密切关系。
1-9的最小公倍数2520。如果一个数能够被2520整除,那么它一定是美丽的数。
我们搜索到一个数的时候,我们都判断该数是不是美丽的数,我们肯定得求出这个数所有位的lcm。
lcm最大是2520,但是1-9能组成的lcm个数 不重复的只有50个不到(可以打表验证)。
有了上面这些 你就可以确定dp数组的范围了。然后就是实现它了
#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int MOD = 2520;
ll dp[50][2550][50], a[50];//核心点就是如何去压缩dp数组的范围
ll Hash[2550];
ll _gcd(ll x, ll y)
{
return y ? _gcd(y, x%y) : x;
}
ll dfs(int pos, int num, int lcm, int limit)
{
//printf("%d\n", num);
if(pos < 0) return num%lcm?0:1;
if(!limit && dp[pos][num][Hash[lcm]] != -1) return dp[pos][num][Hash[lcm]];
int up = limit ? a[pos] : 9;
ll res = 0;
for(int i = 0; i <= up; i++)
{
if(!i) res += dfs(pos-1, (num*10+i)%MOD, lcm, limit && i == a[pos]);
else res += dfs(pos-1, (num*10+i)%MOD, lcm*i/_gcd(lcm, i), limit && i == a[pos]);
}
if(!limit) return dp[pos][num][Hash[lcm]] = res;
return res;
}
ll solve(ll x)
{
int pos = 0;
while(x)
{
a[pos++] = x%10;
x /= 10;
}
return dfs(pos-1, 0, 1, 1);
}
void init()
{
memset(dp, -1, sizeof(dp));
int cnt = 0;
for(int i = 1; i <= MOD; i++)//如果MOD能被i整除,那么i肯定是其中一个能组成的lcm
if(MOD % i == 0) Hash[i] = cnt++;
}
int main()
{
int T;
cin >> T;
init();
ll lt, rt;
while(T--)
{
cin >> lt >> rt;
cout << solve(rt) - solve(lt-1) << endl;
}
return 0;
}