（CodeForces - 55D）Beautiful numbers（数位DP）

题目链接：Problem - 55D - Codeforces

简化题意：求区间[l,r]中满足能够被自身的各非零数字整除的数的个数。

这显然是一道数位DP题，我们先来看看如何判断一个数是否能够被自身的各非零数字整除，这个不是特别难判断，能够被自身的各非零数字整除那就肯定能够被自身的各非零数字的最小公倍数整除，所以在进行数位DP中我们要记录能够遍历到当前位的最小公倍数。除了这个呢？我们还需要记录遍历到当前位所形成的数的大小，比如对于数32120，假如我们遍历到了第三位，那么当前数就是321，1~10的最小公倍数是2520，而题目中给的数是9e18，这么大肯定不可能放到数组里面啊，这个时候就需要用一个性质了，因为我们最后每个数自身各非零数字的最小公倍数肯定能够整除2520，所以我们就可以在过程中保留当前数对2520的余数而不是真正的保留当前数，过程中对2520取余不会影响我们的最终对最小公倍数的取余结果的。（因为其本质是对商操作而没有对余数操作）这样我们就可以把数组开到f[20][2520][2520]，但是大家发现这样还是开不下，但是余数方面已经没有办法再进行优化了只能从当前位的各位非零数字的最小公倍数方面考虑了，其实我们容易发现，1~2520中并不是每个数都能够被1~10中随机几个数的乘积表示的，我们只需要把能够被乘积表示的数单独求出来放在一个数组中即可，这也就是本道题的核心优化方法了，下面是代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
using namespace std;
typedef long long ll;
const int N=23,mod=2520;
ll f[N][2600][50],a[N];
int book[2600];//book[i]记录i是第几个可以整除2520的数 
int gcd(int x,int y)
{
	if(x==0) return y;
	return gcd(y%x,x);
}
int lcm(int x,int y)
{
	return x*y/gcd(x,y);
}
ll dp(int pos,int re,int lcmm,int limit)//re记录当前余数 
{
	if(!pos) return (re%lcmm==0);
	if(!limit&&f[pos][re][book[lcmm]]!=-1) return f[pos][re][book[lcmm]];
	int up=limit?a[pos]:9;
	ll ans=0;
	for(int i=0;i<=up;i++)
	{
		if(i==0)
			ans+=dp(pos-1,re*10%mod,lcmm,limit&&(i==up));
		else
			ans+=dp(pos-1,(re*10+i)%mod,lcm(lcmm,i),limit&&(i==up));
	}
	if(!limit) f[pos][re][book[lcmm]]=ans;
	return ans;
}
ll solve(ll x)
{
	int pos=0;
	while(x)
	{
		a[++pos]=x%10;
		x/=10;
	}
	return dp(pos,0,1,1);
}
void init()
{
	int cntt=0;
	memset(f,-1,sizeof f);//数位dp中，记忆化数组只需要初始化一次 
	for(int i=1;i<=2520;i++)
		if(mod%i==0) book[i]=++cntt;		
}
int main()
{
	int T;
	cin>>T;
	init();
	while(T--)
	{
		ll l,r;
		cin>>l>>r;
		cout<<solve(r)-solve(l-1)<<endl;
	}
	return 0;
}