$CF55D [数位DP]$

题面

数位DP+状压。

首先，按照数位DP的基本套路，每个个位数的最小公倍数为2520，所以只用考虑模2520的情况。考虑一个DP。dp[i][j][k]表示当前是第i位，2~9的数的集合为j，模2520为k的方案数。然后，就是数位DP的基本套路解决这道题。

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
inline LL read () {
    LL res =  ;
    int f () ;
    char ch = getchar ();
    while (!isdigit(ch)) {
        if (ch == '-') f = - ;
        ch = getchar();
    }
    while (isdigit(ch)) res = (res << ) + (res << ) + (ch ^ ),ch = getchar();
    return res * f ;
}
const int N=;
const int mod=;
int p[]= {,,,,,,,,,,,,,,,,,,,
            ,,,,,,,,,,,,,,,,,,
            ,,,,,,,,,,,
           };
int _lcm(int m,int n) {
    return (m*n)/__gcd(m, n);
}
LL f[N][][mod+], bit[N], has[mod+];
inline LL dfs(int i,int lcm,int left, bool e) {
    if(i==) return lcm&&left%p[lcm]== ;
    if(!e && ~f[i][lcm][left])    return f[i][lcm][left];
    LL ans=;
    int u= e? bit[i]: ;
    for(register int d=; d<=u; d++) {
        int t= lcm? has[_lcm(p[lcm], max(d,))]: max(d,);
        ans+=dfs(i-, t, (left*+d)%mod, e&&d==u);
    }
    return e==true? ans: f[i][lcm][left]=ans;
}
LL cal(LL n) {
    int len=;
    while(n) {
        bit[++len]=n%;
        n/=;
    }
    return dfs(len,,,true);
}
signed main() {
    ios::sync_with_stdio(false);
    memset(f,-,sizeof(f));
    for(register int i=; i<; i++) has[p[i]]=i;
    LL t=read();
    while( t-- ) {
        LL L=read(),R=read();
        cout<<cal(R)-cal(L-)<<endl;
    }
    return ;
}

跑的好慢啊QwQ