51nod-1232: 完美数

Posted never-mind

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了51nod-1232: 完美数相关的知识,希望对你有一定的参考价值。

【传送门:51nod-1232


简要题意:

  完美数定义:一个数能够被它每一位非零的数整除,例如:11,12,15,但13就不是完美数,因为13不能被3整除

  给出T组询问,每组询问输入l,r,输出l到r中有多少个完美数


题解:

  首先很容易想到只要当前数被每一位所有非零的数的lcm整除就好了

  实际上最大的lcm只有2520,而且只能够构成48个不同lcm的值

  那么就可以利用这个来节省空间了,将lcm离散编号

  设f[i][j][k]为第i位已经填好数字后,当前构成的数%2520=j,且当前数位的数的lcm的离散编号为k

  然后直接记忆化搜索就行了


参考代码:

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
using namespace std;
typedef long long LL;
LL f[21][3100][51];
int p[3100],Mod=2520;
void pre()
{
    int k=0;p[0]=0;
    for(int i=1;i<=Mod;i++) if(Mod%i==0) p[i]=++k;
}
int cnt[21],tot;
int gcd(int a,int b)
{
    if(a==0) return b;
    else return gcd(b%a,a);
}
int getd(int a,int b)
{
    if(a==0) return b;
    if(b==0) return a;
    return a*b/gcd(a,b);
}
LL dfs(int tp,int x,int d1,int d2)
{
    if(x==1)
    {
        if(d2==0) return 0;
        return d1%d2==0;
    }
    if(tp==0&&f[x][d1][p[d2]]!=-1) return f[x][d1][p[d2]];
    if(tp==0)
    {
        f[x][d1][p[d2]]=0;
        for(int i=0;i<=9;i++)
        {
            if(i==0) f[x][d1][p[d2]]+=dfs(tp,x-1,(d1*10+i)%Mod,d2);
            else f[x][d1][p[d2]]+=dfs(tp,x-1,(d1*10+i)%Mod,getd(d2,i));
        }
        return f[x][d1][p[d2]];
    }
    else
    {
        LL ans=0;
        for(int i=0;i<=cnt[x-1];i++)
        {
            if(cnt[x-1]==0) ans+=dfs(tp,x-1,(d1*10+i)%Mod,d2);
            else if(i==0) ans+=dfs(tp^1,x-1,(d1*10+i)%Mod,d2);
            else if(i==cnt[x-1]) ans+=dfs(tp,x-1,(d1*10+i)%Mod,getd(d2,i));
            else ans+=dfs(tp^1,x-1,(d1*10+i)%Mod,getd(d2,i));
        }
        return ans;
    }
}
LL solve(LL x)
{
    if(x==0) return 0;
    LL d=x;tot=0;
    while(d!=0){cnt[++tot]=d%10;d/=10;}
    LL ans=0;
    for(int i=0;i<=cnt[tot];i++)
    {
        if(i==0) ans+=dfs(0,tot,0,0);
        else if(i==cnt[tot]) ans+=dfs(1,tot,i,i);
        else ans+=dfs(0,tot,i,i);
    }
    return ans;
}
int main()
{
    pre();
    int T;
    scanf("%d",&T);
    memset(f,-1,sizeof(f));
    while(T--)
    {
        LL l,r;
        scanf("%lld%lld",&l,&r);
        printf("%lld
",solve(r)-solve(l-1));
    }
    return 0;
}

 

以上是关于51nod-1232: 完美数的主要内容,如果未能解决你的问题,请参考以下文章

51nod 1232 完美数 数位dp

51nod 1232 完美数

51nod1780 完美序列

51nod 1623 完美消除(数位DP)

动态规划51nod1780 完美序列

使用片段着色器在特定位置绘制完美的水平线