一个数x各个数位上的数之积记为f(x) <不含前导零> 求[L,R)中满足0< f(x)<= n的数的个数
(我的做法应该在这道题里面是最差的了,并且代码应该是最丑的了)
这道题的新奇的地方实际是n的范围,不然其实上是一道大水题了。但其实也只需要改动一点小地方,因为我们发现数位之积是2,3,5,7的倍数,不会有其他的质因子参与进来。所以我们定义一个五维数组,f[i][s1][s2][s3][s4],i为枚举到的位数(从后往前),s1~s4则分别表示枚举到现在的2,3,5,7各个数的指数。然后这道题就迎刃而解了。(数组一定不要因为保险而开过大,不然会MLE,不要问我怎么知道的)
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
using namespace std;
int n,a[19];
long long n1[55],n2[37],n3[19],n4[19],f[19][55][37][19][19];
bool wy(int s1,int s2,int s3,int s4)
{
long long ans=n1[s1]*n2[s2]*n3[s3]*n4[s4];
if(ans<=n)return true;
else return false;
}
long long dfs(int pos,int s1,int s2,int s3,int s4,bool limt,bool lead)
{
if(pos==0)
{
if(wy(s1,s2,s3,s4)==true)return 1;
else return 0;
}
if(wy(s1,s2,s3,s4)==false)return 0;
if(limt==false && lead==false && f[pos][s1][s2][s3][s4]!=-1)return f[pos][s1][s2][s3][s4];
int up=9;long long ans=0;
if(limt==true)up=a[pos];
for(int i=0;i<=up;i++)
{
if(lead==false && i==0)continue;
int ss1=s1,ss2=s2,ss3=s3,ss4=s4;bool bk1=false,bk2=false;
if(i==2)ss1++;
if(i==3)ss2++;
if(i==4)ss1+=2;
if(i==5)ss3++;
if(i==6)ss1++,ss2++;
if(i==7)ss4++;
if(i==8)ss1+=3;
if(i==9)ss2+=2;
if(limt==true && i==a[pos])bk1=true;
if(lead==true && i==0)bk2=true;
ans+=dfs(pos-1,ss1,ss2,ss3,ss4,bk1,bk2);
}
if(limt==false && lead==false)f[pos][s1][s2][s3][s4]=ans;
return ans;
}
long long solve(long long x)
{
int pos=0;
while(x)
{
a[++pos]=x%10;
x/=10;
}
memset(f,-1,sizeof(f));
return dfs(pos,0,0,0,0,true,true);
}
int main()
{
n1[0]=n2[0]=n3[0]=n4[0]=1;
for(int i=1;i<=54;i++)n1[i]=n1[i-1]*2;
for(int i=1;i<=36;i++)n2[i]=n2[i-1]*3;
for(int i=1;i<=18;i++)n3[i]=n3[i-1]*5;
for(int i=1;i<=18;i++)n4[i]=n4[i-1]*7;
long long l,r;
scanf("%d",&n);
scanf("%lld%lld",&l,&r);
printf("%lld\n",solve(r-1)-solve(l-1));
return 0;
}