Codeforces Beta Round #51 D.Beautiful numbers (数位dp+记忆化搜索)

题目

一个美丽的数被定义为可以被它的每一位的数整除的数，

给定区间[l,r](1<=l<=r<=9e18)，求该区间内美丽数的个数

思路来源

https://www.cnblogs.com/kuangbin/archive/2013/05/01/3052670.html

题解

还是kuangbin爷的代码风格踏实清楚明白……

注意到被每个数位都能整除，就相当于被这些数的lcm整除，

①数论相关部分，这题维护当前数的值是不可能的，

但判整除只会判2520的因子（2520是1到9的最小公倍数），

所以在模2520意义下判整除，就将第二维限制到了2520以内

②但第三维lcm的规模也是2520的，事实上只有其因子才会被取到

那就对其离散化一下，存到一个to[]的hash数组里，就把空间降到20*2520*50以内了

③数位dp的经典记忆化搜索的写法，利用无后效性

当当前数%2520的值相同，数位的最小公倍数相同，且均不受上界前缀限制的时候，

两个dp[pos][sum][lcm]的值就是一样的，

所以状态大致只有20*2520*50加上上界前缀限制的部分

上界前缀每个位都要跑一次就大致相当于原规模*2叭

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn=20; 
const int mod=2520;//1到9的lcm 
int to[mod+5];
//dp[i][j][k]表示i位 的这些数位 和为j(模mod) 且最小公倍数为k的方案数 
int bit[maxn];//把原数拆成每一位 
ll dp[maxn][mod][50];//hash之后最后一维只有被mod整除的 
void init()
{
	int cnt=0;
	for(int i=1;i<=mod;++i)
	if(mod%i==0)to[i]=cnt++;
}
int GCD(int a,int b)
{
	return b?GCD(b,a%b):a; 
}
int LCM(int a,int b)
{
	return a/GCD(a,b)*b;
}
//dfs(当前第pos位,当前数位和,当前数位的最小公倍数,是否自由取值,即不为上界的前缀) 
ll dfs(int pos,int sum,int lcm,bool ok)
{
	if(!pos)return sum%lcm==0;
	if(ok&&~dp[pos][sum][to[lcm]])return dp[pos][sum][to[lcm]];
	ll ans=0;
	int mx=ok?9:bit[pos];
	for(int i=0;i<=mx;++i)
	{
		int nowsum=(sum*10+i)%mod;
		int nowlcm=lcm;
		if(i)nowlcm=LCM(nowlcm,i);
		//已经不为上界前缀 或者之前为上界前缀但当前位不为上界 都会导致现在不为上界前缀 
		ans+=dfs(pos-1,nowsum,nowlcm,ok||(i!=mx));
	}
	if(ok)dp[pos][sum][to[lcm]]=ans;
	return ans;
} 
ll solve(ll x)
{
	bit[0]=0;
	for(;x;x/=10)
	bit[++bit[0]]=x%10;
	return dfs(bit[0],0,1,0);
}
int t;
ll l,r;
int main()
{
	init();
	memset(dp,-1,sizeof dp);
	scanf("%d",&t);
	while(t--)
	{
		scanf("%I64d%I64d",&l,&r);
		printf("%I64d\n",solve(r)-solve(l-1));
	}
	return 0;
}