cf55D 数位dp记忆化搜索+状态离散

时间:2024-01-12 18:03:20
/*
漂亮数定义:可以整除任意数位上的数
求出区间[l,r]之间的漂亮数个数
因为
dp[i][j][k]:i位前模lcm的值是j,i位前lcm是k的漂亮数个数
*/
#include<bits/stdc++.h>
using namespace std;
#define ll long long ll dp[][][],has[],a[],len,tot;
void init(){//2520最多也就51个约数
for(int i=;i<=;i++)
if(%i==)has[i]=++tot;
}
ll LCM(ll a,ll b){return a*b/__gcd(a,b);}
ll dfs(ll pos,ll mod,ll lcm,ll lim){
if(pos==)return mod%lcm==;
if(!lim && dp[pos][mod][has[lcm]]!=-)return dp[pos][mod][has[lcm]];
ll res=,num=lim?a[pos]:;
for(int i=;i<=num;i++){
ll c_mod=(mod*+i)%;
ll c_lcm=lcm;
if(i)c_lcm=LCM(lcm,i);
res+=dfs(pos-,c_mod,c_lcm,lim&&i==num);
}
if(!lim)dp[pos][mod][has[lcm]]=res;
return res; }
ll calc(ll x){
len=;
memset(a,,sizeof a);
while(x){
a[++len]=x%;
x/=;
}
return dfs(len,,,);
}
int main(){
init();
int t;
cin>>t;
memset(dp,-,sizeof dp);
while(t--){
ll A,B;
cin>>A>>B;
cout<<calc(B)-calc(A-)<<endl;
}
}