题目分析
cf上的D题还是很有水准,出的题不难,一眼就能知道什么意思,也很容易想到一些东西,但是如果不多加一下处理,就不能成功的解决这道题。这道题真的是一道好题,感觉特别开动人的脑经。
首先我们要求Beautiful numbers,那么我们怎么求呢,肯定要求已经加进去的各位数字的最小公倍数了,只要最后处理到头得到的那个数能够整除这个数的所有位数的最小公倍数就可以,但是你以为这样就够了吗?很明显不行。因为我们不能直接把最后得到的那个数开一维数组直接存下来,大家想想,18位数能存的下吗??肯定不行,那么我们就有一个很奇妙的东西,我们可以对1~9这9个数的最小公倍数(2520直接手算得到)取余就行了,是不是很nice?这时候的数组可以这样写dp[20][2520][2520],一算好像过亿了,数组还是开不下,自己立马傻了,但是这时候又有一个方法,我们可以把1~9这9个数能组成的最小公倍数好像没有2520这么多,实际上很少,我们用一个数组存下来即可。要找的时候直接二分查找就可以了。这样我们dfs的4个参数分别为当前处理的数的位置,前面处理的各个位数对2520的模,前面处理的各个位数上的数字的最小公倍数对应数组中的位置,是否有限制。
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long LL;
LL dp[20][2520][50];
int bit[20],temp[100],length;
int gcd(int a,int b){
return b == 0?a:gcd(b, a%b);
}
int lcm(int a,int b){
return a/gcd(a,b)*b;
}
int Binary_Search(int x){
int left = 0,right = length-1;
while(left <= right){
int mid = (left+right)/2;
if(temp[mid] == x) return mid;
else if(temp[mid] > x) right = mid-1;
else left = mid+1;
}
return -1;
}
LL dfs(int pos,int mod, int rest,int limit){
if(pos < 1)
return (mod%temp[rest]) == 0;
if(!limit && dp[pos][mod][rest] != -1)
return dp[pos][mod][rest];
LL ret = 0;
int len = limit?bit[pos]:9;
for(int i = 0; i <= len; i++){
int pre = rest;
if(i) pre = Binary_Search(lcm(temp[rest], i));
ret += dfs(pos-1, (mod*10+i)%2520, pre, limit&&i==len);
}
if(!limit) dp[pos][mod][rest] = ret;
return ret;
}
void Init(){
length = 0;
for(int i = 1; i <= 2520; i++)
if(2520 % i == 0)
temp[length++] = i;
}
LL solve(LL n){
int len = 0;
while(n){
bit[++len] = n%10;
n /= 10;
}
return dfs(len, 0, 0, 1);
}
int main(){
int t;
scanf("%d", &t);
Init();
memset(dp, -1, sizeof(dp));
while(t--){
LL a,b;
scanf("%I64d %I64d", &a, &b);
printf("%I64d\n", solve(b) - solve(a-1));
}
return 0;
}