询问一个区间内能被它的所有非零位整除的数的个数。
数位dp,1-9的最小公倍数是2520,dp[pos][mod][lcm][f],pos是当前位置,mod取的是到当前位置%2520,实际上1-2520只有48个可取到的公倍数,因此lcm=48,表示到当前位置数位的最小公倍数,若到0位,mod%lcm[lcm]==0,则返回1,否则返回0,f表示的是之后取值的范围,当可以取任意值时这个结果对所有数都是成立的,也就是不受当前询问的影响,因此每次询问不必重新初始化dp数组。
#include<iostream>
#include<string>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<iomanip>
#include<map>
#include<algorithm>
#include<queue>
#include<set>
#define inf 1000000000
#define pi acos(-1.0)
#define eps 1e-8
#define seed 131
using namespace std;
typedef pair<int,int> pii;
typedef unsigned long long ull;
typedef long long ll;
const int maxn=100005;
int MOD=2520;
ll dp[20][2520][48];
int lcm[48];
int digit[20];
void init()
{
int e=0;
for(int i=1;i<=2520;i++)
{
if(MOD%i==0)
lcm[e++]=i;
}
}
int gcd(int a,int b)
{
return b==0?a:gcd(b,a%b);
}
int so(int lc,int i)
{
int p=lcm[lc]*i/gcd(lcm[lc],i);
return lower_bound(lcm,lcm+48,p)-lcm;
}
ll dfs(int pos,int mod,int lc,int f)
{
if(pos==0)
{
if(mod%lcm[lc]==0)
return 1;
return 0;
}
if(f==1&&dp[pos][mod][lc]!=-1)
return dp[pos][mod][lc];
//dp[pos][mod][lc][f]=0;
ll ans=0;
int mx=(f?9:digit[pos-1]);
//cout<<mx<<endl;
for(int i=0;i<=mx;i++)
{
int m=(mod*10+i)%MOD;
int lcc;
if(i!=0)
lcc=so(lc,i);
else
lcc=lc;
int ff;
if(f==1)
ff=1;
else
{
if(i==mx)
ff=0;
else
ff=1;
}
ans+=dfs(pos-1,m,lcc,ff);
}
if(f==1)
dp[pos][mod][lc]=ans;
return ans;
}
ll cal(ll u)
{
int e=0;
while(u)
{
digit[e++]=u%10;
u/=10;
}
//memset(dp,-1,sizeof(dp));
return dfs(e,0,0,0);
}
int main()
{
memset(dp,-1,sizeof(dp));
init();
ll l,r;
int t;
scanf("%d",&t);
while(t--)
{
cin>>l>>r;
cout<<cal(r)-cal(l-1)<<endl;
}
return 0;
}