[bzoj 3679]数字之积

124 篇文章 2 订阅
9 篇文章 0 订阅

一个数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;
}
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值