【数位DP】BZOJ3629数字之积

传送门

数位DP大法好…
乘积太大保存不下来肿么办？
这乘积都是$1$到$9$的数字乘起来滴，于是乎用质因数表示就好了，特别注意对每一种质数的最大个数最好卡着开数组，要不然很容易MLE滴…

原谅蒟蒻太弱，代码为了放错就写的很丑= =

#include <iostream>
#include <cstdio>
#include <cstring>
#define LL long long int
using namespace std;

int n, w[20], len;
LL p2[31], p3[20], p5[14], p7[12];
LL l, r, f[20][31][20][14][12][2];

inline void add(int num,int l,int a,int b,int c,int d,int t,LL p)
{
    LL tmp=num*p2[a]*p3[b]*p5[c]*p7[d];
    if(tmp>n)return;
    if(num==2)++a;
    else if(num==3)++b;
    else if(num==4)a+=2;
    else if(num==5)++c;
    else if(num==6)++a, ++b;
    else if(num==7)++d;
    else if(num==8)a+=3;
    else if(num==9)b+=2;
    else if(num==10)++a, ++c;
    f[l][a][b][c][d][t]+=p;
}

LL cal(LL m)
{
    if(m==0)return 1;
    memset(f,0,sizeof f);

    len=0;
    while(m){w[++len]=m%10;m/=10;}
    add(w[len],len,0,0,0,0,1,1);
    for(int i=w[len]-1;i;--i)add(i,len,0,0,0,0,0,1);
    for(int i=len;i;--i)add(1,i,0,0,0,0,0,1);

    for(int i=len-1;i;--i)
    {
        for(int a=0;p2[a]<=n;++a)
            for(int b=0, e=n/p2[a];p3[b]<=e;++b)
                for(int c=0, e2=e/p3[b];p5[c]<=e2;++c)
                    for(int d=0, e3=e2/p5[c];p7[d]<=e3;++d)
                        for(int k=1;k<=9;++k)
                        {
                            if(k==w[i])add(k,i,a,b,c,d,1,f[i+1][a][b][c][d][1]);
                            else if(k<w[i])add(k,i,a,b,c,d,0,f[i+1][a][b][c][d][1]);
                            add(k,i,a,b,c,d,0,f[i+1][a][b][c][d][0]);
                        }
    }
    LL ans=0;
    for(int a=0;p2[a]<=n;++a)
            for(int b=0, e=n/p2[a];p3[b]<=e;++b)
                for(int c=0, e2=e/p3[b];p5[c]<=e2;++c)
                    for(int d=0, e3=e2/p5[c];p7[d]<=e3;++d)
                        if(p2[a]*p3[b]*p5[c]*p7[d]<=n)
                            ans+=f[1][a][b][c][d][0]+f[1][a][b][c][d][1];
    return ans;
}

int main()
{
    p2[0]=p3[0]=p5[0]=p7[0]=1;
    for(int i=1;i<31;++i)p2[i]=p2[i-1]<<1;
    for(int i=1;i<20;++i)p3[i]=p3[i-1]*3;
    for(int i=1;i<14;++i)p5[i]=p5[i-1]*5;
    for(int i=1;i<12;++i)p7[i]=p7[i-1]*7;
    scanf("%d",&n);
    cin>>l>>r;
    cout<<cal(r-1)-cal(l-1);
    return 0;
}