codeforces 55D D. Beautiful numbers(数位dp+数论)

题目链接：

codeforces 55D

---

题目大意：

求在[l,r]中能够整除自己每个数位上的数字的数的个数。

---

题目分析：

• 首先我们能够知道如果这个数能够整除它的每个数位上的数字，那么它一定能够整除他们的最小公倍数，是充要的。
• 那么我们定义状态dp[i][j][k]代表i位在任意组合下得到的所有数位的数字的最小公倍数为j的每个数位上的数字之积%2520为k的方案数。
• 我们可以知道所有的公倍数最大不会超过2520，而且他们都是2520的约数，所以如果他们能够整除2520的余数，那么证明他们能够整除所有数位上的数字之积，是充要的。
• 所以我们利用这样状态记录进行记忆化搜索，对于给定的数字求取比它小的所有的符合要求的数，间接的求出区间内合法的数的个数，具体求法就是对于给定数的每一位，枚举它可能的数字，如果小于当前位i的数字，那么后面方案数就是dp[n-i][j][k]，找出满足条件的j和k即可，利用深搜来实现。如果等于当前位，那么固定当前位，更新lcm，然后转移到下一位去寻找方案数。
• 清楚了状态就是一个简单的枚举和统计的操作了。具体看代码吧….写的不清楚的地方欢迎评论指正。

AC代码：

#include <iostream>#include <cstring>#include <cstdio>#include <algorithm>using namespace std;typedef long long LL;LL dp[20][50][3000];int d[20],hash[3000];const int mod = 2520;int gcd ( int a ,int b ){    return !b?a:gcd(b,a%b);}void init ( ){    memset ( dp , -1 , sizeof ( dp ) );    int cc = 0;    for ( int i = 1 ; i <= 8 ; i *= 2 )        for ( int j = 1; j <= 9 ; j*= 3 )            for ( int k = 1 ; k <= 5 ; k *= 5 )                for ( int t = 1; t <= 7 ; t *= 7 )                    hash[i*j*k*t] = ++cc;}LL dfs ( int n , int lcm , int f , int r ){    if ( !n ) return r%lcm==0;    int ll = hash[lcm];    if ( f && dp[n][ll][r] != -1 )         return dp[n][ll][r];    int x = f?9:d[n];    LL ret = 0;    for ( int i = 0 ; i <= x ;i++ )        ret += dfs ( n-1 , i==0?lcm:lcm*i/gcd(lcm,i) , i==x?f:1 , (r*10+i)%mod );    return f? dp[n][ll][r] = ret : ret;}LL solve ( LL x ){    int cc = 0;    while ( x )    {        d[++cc] = x%10;        x /= 10;    }    return dfs ( cc , 1 , 0 , 0 );}int main ( ){    int t;    LL a,b;    scanf ( "%d" , &t );    init ( );    while ( t-- )    {        scanf ( "%lld%lld" , &a , &b );        printf ( "%lld\n" , solve ( b ) - solve ( a-1 ) );    }}