Codeforces - 55D Beautiful numbers

Posted Pneuis

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了Codeforces - 55D Beautiful numbers相关的知识,希望对你有一定的参考价值。

题意: 

  在给定范围内(9*10^18)求出有多少个数满足被所有非零位的数整除。

思路:

  可以先想到如果被所有非零位的数整除,那么就会被他们的lcm整除。

  而lcm(1,2...9) = 2520.

  即2520可以整除1 - 9中任意组合的lcm。

  所以可以用%2520的值表示状态,最后用最终状态值%lcm(x1...xn),判断是否为0。

  1 - 9lcm组合的情况数由 2^3*3^2*5*7可以得出为:4*3*2*2 = 48 ,预处理提取出来。

技术分享图片
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MOD = 2520;
int t;
int bit[20], index[MOD+5];
ll l, r;
ll dp[20][MOD][50];
int gcd(int x, int y) {
    return y==0?x:gcd(y, x%y);
}
int lcm(int x, int y) {
    return x/gcd(x, y)*y;
}
void init() {
    int cnt = 0;
    memset(dp, -1, sizeof(dp));
    for(int i = 1; i <= MOD; i++)
    if(!(MOD % i)) index[i] = cnt++; 
}
ll dfs(int pos, int presum, int prelcm, bool limit) {
    if(pos == -1) return presum%prelcm==0;
    if((!limit) && dp[pos][presum][index[prelcm]] != -1) 
        return dp[pos][presum][index[prelcm]];
    int up = limit?bit[pos]:9;
    ll res = 0;
    for(int i = 0; i <= up; i++) {
        int nowsum = (presum*10+i)%MOD;
        int nowlcm = i==0?prelcm:lcm(prelcm, i);
        res += dfs(pos-1, nowsum, nowlcm, limit && i == bit[pos]);
    } 
    if(!limit) dp[pos][presum][index[prelcm]] = res;
    return res;
}
ll solve(ll x) {
    int pos = 0;
    while(x) {
        bit[pos++] = x%10;
        x /= 10; 
    }
    return dfs(pos-1, 0, 1, true);
}
int main() {
    init();
    scanf("%d", &t);
    while(t--) {
        scanf("%I64d%I64d", &l, &r);
        printf("%I64d\n", solve(r) - solve(l-1));
    }
}
View Code

 

以上是关于Codeforces - 55D Beautiful numbers的主要内容,如果未能解决你的问题,请参考以下文章

CodeForces - 55D Beautiful numbers

codeforces55D Beautiful numbers

CodeForces 55D Beautiful numbers

CodeForces - 55D Beautiful numbers (数位DP)

CodeForces - 55D - Beautiful numbers(数位DP,离散化)

CodeForces 55D Beautiful numbers