Codeforces Beta Round #51 D Beautiful numbers(数位dp+离散化)

题目链接：https://codeforces.com/contest/55/problem/D

 

题目大意：求l到r中所有满足能被自己非0的数位整除的数字个数

 

题目思路：数位dp入门，准备从0开始学一下数位dp。这个题目感觉非常巧妙，首先被所有数位整除，就需要想到一个事情，最小的满足能被多个数字整除的数字，是什么？没错，就是他们的最小公倍数，1~9的最小公倍数是2520，也就是所有的2520的倍数，都一定能被自己的数位整除。那么也就是说，9个数字里面任意取出构成的lcm，一定能被2520整除，也就是一定是2520的因数，2520=2^3 * 3^2 * 5 * 7，那么一共只有4*3*2*2=48种不同可能。

然后我们可以意识到，2520倍数的那一部分，无论怎么弄都能被消掉，所以直接模一发，剩下来的那一部分要被当前已经用的数位的lcm整除，就能满足题目要求，所以我们只需要记录一路下来的模2520的结果和lcm即可，但是这样数组开不下，所以需要对lcm来个离散化就美滋滋了

 

以下是代码：

#include<bits/stdc++.h>
using namespace std;
#define inf 0x3f3f3f3f
#define rep(i,a,b) for(ll i=a;i<=b;i++)
#define per(i,a,b) for(ll i=a;i>=b;i--)
#define ll long long
const ll MAXN=2525;
const ll MOD =1e9+7;
ll l,r,a[MAXN],Hash[MAXN],dp[100][2525][50];
ll dfs(ll pos,ll mod,ll lcm,ll limit){
    if(pos==0)return mod%lcm==0;
    if(!limit&&dp[pos][mod][Hash[lcm]]!=-1)return dp[pos][mod][Hash[lcm]];
    ll up=limit?a[pos]:9;
    ll ans=0;
    rep(i,0,up){
        ans+=dfs(pos-1,(mod*10+i)%2520,i==0?lcm:lcm*i/__gcd(lcm,i),limit&&i==a[pos]);
    }
    if(!limit)dp[pos][mod][Hash[lcm]]=ans;
    return ans;
}
ll solve(ll x){
    ll pos=0;
    while(x){
        a[++pos]=x%10;
        x/=10;
    }
    ll ans=dfs(pos,0,1,1);
    return ans;
}
int main(){
    ll t;
    scanf("%I64d",&t);
    memset(Hash,0,sizeof(Hash));
    ll pos=0;
    rep(i,1,2520){
        if(2520%i==0){
            Hash[i]=++pos;
        }
    }
    memset(dp,-1,sizeof(dp));
    while(t--){
        scanf("%I64d%I64d",&l,&r);
        cout<<solve(r)-solve(l-1)<<endl;
    }
    return 0;
}