CodeForces 55D-Beautiful numbers

题目链接

题目分析

题意

输入t组数据，每组数据包含2个整数l和r，输出从l到r范围的beautiful numbers（即能被自身各个数位上的非零数字整除的数）的数量。

思路

一个数能被它各个数位上的非零数字整除，转化为可以被它们的最小公倍数整除。1到9的最小公倍数是2520。
2 =1*2
3=1*3
6=2*3
4=2*2
5=1*5
7=7*1
8=4*2=2*2*2
9=3*3
最小公倍数为 1*2*2*2*3*3*5*7=2520
数位DP，dfs(pos,mod,lcm,f),pos为当前位，mod为前面那些位对2520的模，lcm为前面那些数位的最小公倍数，f标记前面那些位是否达到上限，这样一来dp数组就要开到19*2520*2520，明显超内存了，考虑到最小公倍数是离散的，1-2520中可能是最小公倍数的其实只有48个（还是不太懂，求指教），经过离散化处理后，dp数组的最后一维可以降到48。

代码

#include<stdio.h>
#include<string.h>
#include<iostream>
using namespace std;
#define LL long long
#define MOD 2520
#define maxn 25
int gcd(int m,int n){//求m和n的最小公倍数 
    return n==0?m:gcd(n,m%n);
}
int N,dig[maxn],Hash[MOD+10];
LL l,r;
LL dp[maxn][MOD+1][50];
LL dfs(int cur,int mod,int lcm,int e){
    if(cur<0)
        return mod%lcm==0;
    if(!e&&dp[cur][mod][Hash[lcm]]!=-1)
        return dp[cur][mod][Hash[lcm]];
    int end=(e?dig[cur]:9);
    LL ans=0;
    for(int i=0;i<=end;i++){
        int tmp=(mod*10+i)%MOD;
        int tmplcm=lcm;
        if(i)
            tmplcm=lcm/gcd(lcm,i)*i;
        ans+=dfs(cur-1,tmp,tmplcm,e&&i==end);
    }
    if(!e)
        dp[cur][mod][Hash[lcm]]=ans;
    return ans;
}
LL solve(LL n){
    int len=0;
    while(n){
        dig[len++]=n%10;
        n/=10;
    }
    return dfs(len-1,0,1,1);
}
void init(){
    int cnt=0;
    for(int i=1;i<=MOD;i++)
        if(MOD%i==0)
            Hash[i]=++cnt;
    memset(dp,-1,sizeof(dp));   
}
int main(){
    int t;
    cin>>t;//数据组数 
    init();
    while(t--){
        cin>>l>>r;//每组数据的范围 
        cout<<solve(r)-solve(l-1)<<endl;//范围内的beautiful numbers数目 
    }
    return 0;
}