BZOJ 1853 【SCOI2010】 幸运数字

Description

在中国，很多人都把6和8视为是幸运数字！lxhgww也这样认 为，于是他定义自己的“幸运号码”是十进制表示中只包含数字6和8的那些号码，比如68，666，888都是“幸运号码”！但是这种“幸运号码”总是太少 了，比如在[1,100]的区间内就只有6个（6，8，66，68，86，88），于是他又定义了一种“近似幸运号码”。lxhgww规定，凡是“幸运号 码”的倍数都是“近似幸运号码”，当然，任何的“幸运号码”也都是“近似幸运号码”，比如12，16，666都是“近似幸运号码”。 现在lxhgww想知道在一段闭区间[a, b]内，“近似幸运号码”的个数。

Input

输入数据是一行，包括2个数字a和b

Output

输出数据是一行，包括1个数字，表示在闭区间[a, b]内“近似幸运号码”的个数

HINT

【数据范围】
对于$30\%$的数据，保证$1 \leqslant a \leqslant b \leqslant1000000$
对于$100\%$的数据，保证$1 \leqslant a \leqslant b \leqslant 10000000000$

 

　　这道题一开始我还以为需要用到什么神奇的数学推导，或者一些什么奇妙的数学公式，然后看了题解之后发现是一道搜索题……

　　一个非常显然的事实就是幸运号码不会太多。把表打出来，就会发现在$10^{10}$以内的幸运数只有$2000$多一点……

　　这个时候一个非常显然的想法就是对这些数进行容斥，即加上每个数的倍数个数，减去两个数的倍数个数，加上三个的，……以此类推。

　　这样的复杂度显然是不对的，理论上可达$O(2^x)$，其中$x$为幸运数个数。但是由于多个数的倍数不能超过右边界，就可以减一大刀，实际复杂度低了不知道多少。

　　但是这样任然不够。我们还可以对幸运数进行处理，将其中是另外的幸运数的倍数的数给去掉。这样可以将需要考虑的数的个数减掉一半左右。

　　最后还有一个小优化，那就是将最后需要处理的幸运数按从大到小排好序。这样可以让乘积尽早变得更大，可以减掉许多不必要的计算。

　　加了上述优化，就差不多可以$AC$了。

　　下面贴代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define File(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)
#define maxn 10010

using namespace std;
typedef long long llg;

int la,lb,ci;
llg now,ans,mi[15],l,r;
llg a[maxn],b[maxn];

void search(int d){
	if(now>r) return;
	llg xx=now;
	now=xx+6*mi[d]; search(d+1);
	now=xx+8*mi[d]; search(d+1);
	if(xx) a[++la]=xx; now=xx;
}

llg gcd(llg a,llg b){
	llg r=a%b;
	while(r) a=b,b=r,r=a%b;
	return b;
}

void dfs(int j){
	if(j==lb+1){
		if(!ci) return;
		if(ci&1) ans+=r/now-(l-1)/now;
		else ans-=r/now-(l-1)/now;
		return;
	}
	dfs(j+1); llg xx=now; ci++;
	now=xx/gcd(xx,b[j]);
	if((double)now*b[j]<=r){
		now*=b[j];
		if(now<=r) dfs(j+1);
	}
	ci--; now=xx;
}

int main(){
	File("a");
	mi[0]=1;
	for(int i=1;i<=10;i++) mi[i]=mi[i-1]*10;
	scanf("%lld %lld",&l,&r);
	search(0); sort(a+1,a+la+1);
	for(int i=1;i<=la;i++){
		b[++lb]=a[i];
		for(int j=1;j<lb;j++)
			if(a[i]%b[j]==0){lb--; break;}
	}
	for(int i=1;i<=lb/2;i++) swap(b[i],b[lb-i+1]);
	now=1; dfs(1);
	printf("%lld",ans);
	return 0;
}