题目链接:点击这里
题目大意:
求
[
l
,
r
]
[l,r]
[l,r] 有多少种漂亮数字,其定义为能被他所有非零数字整除的数字
题目分析:
数位dp
如果一个数能整除的他的所有非零数位,相当于他能整除所有非零数位的最小公倍数,所以在记忆化是我们可以维护一下数位前缀的
l
c
m
lcm
lcm 。
我们发现
l
c
m
(
1
,
2
,
.
.
.
,
9
)
=
2520
lcm(1,2,...,9)=2520
lcm(1,2,...,9)=2520 ,所以可以利用取余的性质:若
a
m
o
d
b
=
0
a\mod b=0
amodb=0 且
c
m
o
d
b
=
0
c\mod b=0
cmodb=0 那么
(
a
m
o
d
c
)
m
o
d
b
=
0
(a\mod c) \mod b=0
(amodc)modb=0 简化状态,所以我们可以设计状态:
d
p
[
i
]
[
j
]
[
k
]
dp[i][j][k]
dp[i][j][k] 表示处理了前
i
i
i 位,前面数的
l
c
m
lcm
lcm 是
j
j
j ,和是
k
k
k ,由于
1
1
1~
9
9
9 组成的
l
c
m
lcm
lcm 只有
48
48
48 个,所以可以先离散化一下,这样数组的第二维只需要开
48
48
48 左右即可;数组第三维利用前面所述的性质开
2520
2520
2520 左右即可
具体细节见代码:
//#pragma GCC optimize(2)
//#pragma GCC optimize("Ofast","inline","-ffast-math")
//#pragma GCC target("avx,sse2,sse3,sse4,mmx")
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<set>
#include<map>
#include<stack>
#include<queue>
#define ll long long
#define inf 0x3f3f3f3f
#define int ll
#define endl '\n'
#define IOS ios::sync_with_stdio(0); cin.tie(0); cout.tie(0)
using namespace std;
int read()
{
int res = 0,flag = 1;
char ch = getchar();
while(ch<'0' || ch>'9')
{
if(ch == '-') flag = -1;
ch = getchar();
}
while(ch>='0' && ch<='9')
{
res = (res<<3)+(res<<1)+(ch^48);//res*10+ch-'0';
ch = getchar();
}
return res*flag;
}
const int maxn = 1e6+5;
const int mod = 1e9+7;
const double pi = acos(-1);
const double eps = 1e-8;
int n,k,a[maxn],dp[25][50][2555],mp[maxn],cnt;
int _lcm(int a,int b)
{
return a/__gcd(a,b)*b;
}
const int maxlcm = 2520;
int dfs(int pos,int lcm,int sum,bool f1)
{
if(pos == -1) return sum%lcm==0;
if(!f1 && dp[pos][mp[lcm]][sum] != -1) return dp[pos][mp[lcm]][sum];
int up = f1 ? a[pos] : 9,res = 0;
for(int i = 0;i <= up;i++)
{
int cursum = (sum*10+i)%maxlcm;
int curlcm = i ? _lcm(lcm,i) : lcm;
res += dfs(pos-1,curlcm,cursum,f1&&i==up);
}
if(!f1) dp[pos][mp[lcm]][sum] = res;
return res;
}
int calc(int x)
{
int cnt = 0;
while(x)
{
a[cnt++] = x%10;
x /= 10;
}
return dfs(cnt-1,1,0,1);
}
signed main() {
for(int i = 1;i <= maxlcm;i++)
if(maxlcm%i == 0) mp[i] = ++cnt;
memset(dp,-1,sizeof(dp));
int t = read();
while(t--)
{
int l = read(),r = read();
printf("%lld\n",calc(r)-calc(l-1));
}
return 0;
}