码迷,mamicode.com
首页 > 其他好文 > 详细

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

时间:2015-09-07 22:58:24      阅读:475      评论:0      收藏:0      [点我收藏+]

标签:codeforces   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 ) );
    }
}

版权声明:本文为博主原创文章,未经博主允许不得转载。

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

标签:codeforces   dp   数论   

原文地址:http://blog.csdn.net/qq_24451605/article/details/48274605

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!