CF55D Beautiful numbers（dp+离散化+数论）

要求整除每一位上的数，就是要mod lcm(a1,a2,…an)余0
M=lcm(1,2,…,9)=2520
因为$a \%(nx)\equiv a \mod x$
所以我们只关心x%2520的余数即可。
我们有dp[i][j][x]表示前面的数的各位的lcm为j,前面的数%2520余x,还剩的i位有多少种合法填法。
然后因为这样还是太慢，我们把j离散化掉，最多只有48种（4*3*2*2，考虑每一个质数的指数）。
然后记忆化搜索来求0~x有多少个合法的数，注意处理不大于x这个限制。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define mod 2520
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(S==T) return EOF;}
    return *S++;
}
inline ll read(){
    ll x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int dig[30],hs[2522],tot=0;ll dp[20][50][2520];
//dp[i][j][x],前面的数的各位的lcm为j,前面的数%2520余x,还剩的i位有多少种合法填法。
inline int gcd(int x,int y){return y?gcd(y,x%y):x;}
inline int lcm(int x,int y){return x&&y?x*y/gcd(x,y):x+y;}
inline ll dfs(int cnt,int j,int x,bool flag){//填第cnt位,前面已经填过的数的lcm为j,余数为x,flag为1表示已经小于lim，否则还未小于
    if(!cnt) return x%j==0;
    if(flag){
        if(~dp[cnt][hs[j]][x]) return dp[cnt][hs[j]][x];dp[cnt][hs[j]][x]=0;
        for(int i=0;i<=9;++i) dp[cnt][hs[j]][x]+=dfs(cnt-1,lcm(i,j),(x*10+i)%mod,1);
        return dp[cnt][hs[j]][x];
    }ll res=0;
    for(int i=0;i<=dig[cnt];++i)
    res+=dfs(cnt-1,lcm(i,j),(x*10+i)%mod,i!=dig[cnt]);
    return res;
}
inline ll cal(ll x){
    int len=0;
    while(x){
        dig[++len]=x%10;x/=10;
    }return dfs(len,1,0,0);
}
int main(){
//  freopen("a.in","r",stdin);
    int tst=read();memset(dp,-1,sizeof(dp));
    for(int i=1;i<=mod;++i)
        if(mod%i==0) hs[i]=++tot;
    while(tst--){
        ll l=read(),r=read();
        printf("%I64d\n",cal(r)-cal(l-1));
    }return 0;
}