codeforces 55D Beautiful numbers

最新推荐文章于 2024-06-17 14:54:52 发布

原创最新推荐文章于 2024-06-17 14:54:52 发布 · 316 阅读

1 ·

CC 4.0 BY-SA版权

动态规划-数位DP 专栏收录该内容

10 篇文章

订阅专栏

本文介绍了一种使用数位DP算法解决特定数学问题的方法：在一个给定区间内找出能被其所有非零位整除的数的数量。通过预处理1到2520之间的最小公倍数，利用动态规划技术来高效地解决问题。

询问一个区间内能被它的所有非零位整除的数的个数。

数位dp，1-9的最小公倍数是2520，dp[pos][mod][lcm][f]，pos是当前位置，mod取的是到当前位置%2520，实际上1-2520只有48个可取到的公倍数，因此lcm=48，表示到当前位置数位的最小公倍数，若到0位，mod%lcm[lcm]==0，则返回1，否则返回0，f表示的是之后取值的范围，当可以取任意值时这个结果对所有数都是成立的，也就是不受当前询问的影响，因此每次询问不必重新初始化dp数组。

#include<iostream>
#include<string>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<iomanip>
#include<map>
#include<algorithm>
#include<queue>
#include<set>
#define inf 1000000000
#define pi acos(-1.0)
#define eps 1e-8
#define seed 131
using namespace std;
typedef pair<int,int> pii;
typedef unsigned long long ull;
typedef long long ll;
const int maxn=100005;
int MOD=2520;
ll dp[20][2520][48];
int lcm[48];
int digit[20];
void init()
{
    int e=0;
    for(int i=1;i<=2520;i++)
    {
        if(MOD%i==0)
            lcm[e++]=i;
    }
}
int gcd(int a,int b)
{
    return b==0?a:gcd(b,a%b);
}
int so(int lc,int i)
{
    int p=lcm[lc]*i/gcd(lcm[lc],i);
    return lower_bound(lcm,lcm+48,p)-lcm;
}
ll dfs(int pos,int mod,int lc,int f)
{
    if(pos==0)
    {
        if(mod%lcm[lc]==0)
            return 1;
        return 0;
    }
    if(f==1&&dp[pos][mod][lc]!=-1)
        return dp[pos][mod][lc];
    //dp[pos][mod][lc][f]=0;
    ll ans=0;
    int mx=(f?9:digit[pos-1]);
    //cout<<mx<<endl;
    for(int i=0;i<=mx;i++)
    {
        int m=(mod*10+i)%MOD;
        int lcc;
        if(i!=0)
            lcc=so(lc,i);
        else
            lcc=lc;
        int ff;
        if(f==1)
            ff=1;
        else
        {
            if(i==mx)
                ff=0;
            else
                ff=1;
        }
        ans+=dfs(pos-1,m,lcc,ff);
    }
    if(f==1)
        dp[pos][mod][lc]=ans;
    return ans;
}
ll cal(ll u)
{
    int e=0;
    while(u)
    {
        digit[e++]=u%10;
        u/=10;
    }
    //memset(dp,-1,sizeof(dp));
    return dfs(e,0,0,0);
}
int main()
{
    memset(dp,-1,sizeof(dp));
    init();
    ll l,r;
    int t;
    scanf("%d",&t);
    while(t--)
    {
        cin>>l>>r;
        cout<<cal(r)-cal(l-1)<<endl;
    }
    return 0;
}