Beautiful numbers CodeForces - 55D 数位dp 求能被各位数字整除的个数

题目链接：https://vjudge.net/problem/CodeForces-55D
转自：https://blog.csdn.net/codeswarrior/article/details/80434276
题意：求l~r之间能被除0外各位数字整除的数的个数。
思路： 一个数能被多个数整除，那么这个数一定也能整除多个数的最小公倍数。 dp【i】【j】【k】表示到第i位已经枚举了数字j，以及j的各位数字的最小公倍数k。正常情况下dp数组需要开20 * 1e18 * 2520，因为只需要判断j是否是k的倍数，j可以对2520取模，因此j只需要开到2520。这样k数组还是太大。
设集合s为1~9中任意1个、2个……9个数的最小公倍数，则s含有1、2、3、……2 * 6=12、2 * 7=14、2 * 8=16等等，也就是说s的大小肯定是远小于2520个数的，再加上它们都能被2520整除，计算得出s的大小=48，也就可以进一步缩小k数组的大小。

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <vector>
#include <queue>
#include <stack>
using namespace std;
long long dp[50][2530][50],a[50],hsh[3000];
void init()
{
    int cnt=0;
    for(int i=1; i<=2520; i++)
    {
        if(2520%i==0)
            hsh[i]=cnt++;
    }
}
long long Lcm(long long a,long long b)
{
    return a/__gcd(a,b)*b;
}
long long dfs(long long pos,long long num,long long lcm,bool limit)
{
    if(pos==-1)
        return num%lcm==0;
    long long &dpnow=dp[pos][num][hsh[lcm]];
    if(!limit&&dpnow!=-1)
        return dpnow;
    int up=limit?a[pos]:9;
    long long ans=0;
    for(int i=0;i<=up;i++)
    {
        long long cnt;
        if(!i)
            cnt=lcm;
        else
            cnt=Lcm(i,lcm);
        ans+=dfs(pos-1,(num*10+i)%2520,cnt,limit&&i==a[pos]);
    }
    if(!limit)
        dpnow=ans;
    return ans;
}
long long solve(long long x)
{
    long long pos=0;
    while(x)
    {
        a[pos++]=x%10;
        x/=10;
    }
    return dfs(pos-1,0,1,1);
}
int main()
{
    ios::sync_with_stdio(false);
    int T;
    cin>>T;
    init();
    memset(dp,-1,sizeof(dp));
    while(T--)
    {
        long long l,r;
        cin>>l>>r;
        long long ans=solve(r);
        ans-=solve(l-1);
        cout<<ans<<"\n";
    }
    return 0;
}

错误代码
dp【i】【j】表示第i位的数字是j的情况下的答案。会导致重复，如10x和00x就会状态重复。

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <vector>
#include <queue>
#include <stack>
using namespace std;
long long dp[50][10],a[50],b[50];
long long dfs(long long pos,bool limit,long long cur)
{
    if(pos==-1)
    {
        long long p=0;
        long long cc=cur;
        while(cc)
        {
            b[p++]=cc%10;
            cc/=10;
        }
        bool flag=true;
        for(long long i=0; i<p; i++)
        {
            if(b[i]==0)
                continue;
            if(cur%b[i]!=0)
            {
                flag=false;
                break;
            }
        }
        if(flag)
            return 1;
        else
            return 0;
    }
    if(!limit&&dp[pos][cur%10]!=-1)
        return dp[pos][cur%10];
    long long up=limit?a[pos]:9;
    long long tmp=0;
    for(int i=0; i<=up; i++)
    {
        tmp+=dfs(pos-1,limit&&i==a[pos],cur*10+i);
    }
    if(!limit)
        dp[pos][cur%10]=tmp;
    return tmp;
}
long long solve(long long x)
{
    long long pos=0;
    while(x)
    {
        a[pos++]=x%10;
        x/=10;
    }
    return dfs(pos-1,true,0);
}
int main()
{
    int T;
    scanf("%d",&T);
    memset(dp,-1,sizeof(dp));
    while(T--)
    {
        long long l,r;
        scanf("%I64d%I64d",&l,&r);
        long long ans=solve(r)-solve(l-1);
        printf("%I64d\n",ans);
    }
    return 0;
}