2019.01.17 bzoj1853: [Scoi2010]幸运数字（容斥+dfs）

版权声明：随意转载哦......但还是请注明出处吧： https://blog.csdn.net/dreaming__ldx/article/details/86520491

传送门
搜索菜题，然而第一次没有注意然后爆 $long long$了。
题意：称所有数位由 $6,8$组成的数为幸运数字，问一个一个区间 $[l,r]$中所有幸运数字及其倍数的个数。

---

思路：
先把所有的幸运数字找出来并筛去那些会算重的，剩下一共不超过 $1000$个数。
即如果满足 $num_i|num_j$就删去 $num_j$。
然后考虑直接容斥+搜索算出答案。
现在就只用想怎么剪枝了。

1. 当前的 $lcm$大于b
2. 把所有数从大到小排序。
3. 由于每一次取 $lcm$至少乘三，因此只要一个数是 $\le b/3$的就不用加入最终的 $num$集合。

然后直接搜就行了。
代码：

#include<bits/stdc++.h>
#define ri register int
using namespace std;
typedef long long ll;
const ll divv=1e9;
ll a,b,val[2500],lim,tot=0,ans=0,sig=0,coef[2500];
inline void calc(ll mul){
	if(mul>b)return;
	if(mul)val[++tot]=mul;
	calc(mul*10+6),calc(mul*10+8);
}
inline ll gcd(ll a,ll b){while(b){ll t=a;a=b,b=t%a;}return a;}
inline ll lcm(ll a,ll b){return a/gcd(a,b)*b;}
inline bool check(ll x,ll y){
	if((x/divv)*(y/divv))return 0;
	return x*y<=b;
}
inline void dfs(int pos,int tim,ll mult){
	if(mult>b)return;
	if(pos==sig+1){
		if(mult==1)return;
		ans+=(b/mult-a/mult)*(tim&1?1:-1);
		return;
	}
	dfs(pos+1,tim,mult);
	ll x=mult/gcd(mult,coef[pos]),y=coef[pos];
	if(check(x,y))dfs(pos+1,tim+1,x*y);
}
int main(){
	cin>>a>>b,--a;
	calc(0),ans=0;
	sort(val+1,val+tot+1);
	for(ri i=1;i<=tot;++i){
		bool f=1;
		for(ri j=1;j<i;++j)if(val[i]==val[i]/val[j]*val[j]){f=0;break;}
		if(f){
			if(val[i]<=b/3)coef[++sig]=val[i];
			else ans+=b/val[i]-a/val[i];
		}
	}
	reverse(coef+1,coef+sig+1);
	dfs(1,0,1);
	cout<<ans;
	return 0;
}