[BZOJ1853][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]内“近似幸运号码”的个数
Sample Input
【样例输入1】
1 10
【样例输入2】
1234 4321
Sample Output
【样例输出1】
2
【样例输出2】
809
HINT
【数据范围】
对于30%的数据，保证1 < =a < =b < =1000000
对于100%的数据，保证1 < =a < =b < =10000000000

很明显的容斥裸题，预处理出所有‘68’数字后，利用容斥求解即可
只是这道题有两个地方需要注意
一、幸运数字中可能会存在a是b的倍数，所以这种情况要把a先剔除
二、数据范围为10^10，注意求lcm乘的时候不要爆long long，可先开个double检验一下，
哦对还有，dfs的时候先搜大的数，这样能够比较早的剔除掉大于R的情况

#include<cstdio>#include<cstring>#include<algorithm>#include<iostream>#ifdef WIN32#define AUTO "%I64d"#else#define AUTO"%lld"#endif#define LL long long using namespace std;const int maxn=10000+5;int cnt,tot;LL L,R,num[maxn],ans;void init(LL x){    if (x>R) return ;    if (x) num[++cnt]=x;    init(x*10+6);    init(x*10+8);}bool cmp(LL x,LL y){    return x>y;}LL gcd(LL a,LL b){    return b?gcd(b,a%b):a;}void dfs(int ID,LL sum,LL lcm){    if (ID>tot)    {        if (sum&1) ans+=R/lcm-(L-1)/lcm;        else if (sum) ans-=R/lcm-(L-1)/lcm;        return ;    }     dfs(ID+1,sum,lcm);    LL temp=lcm/gcd(num[ID],lcm);    if(((double)num[ID]*temp)<=R) dfs(ID+1,sum+1,num[ID]*temp);}int main(){//  freopen("scoi.in","r",stdin);    scanf(AUTO AUTO,&L,&R);    init(0);    sort(num+1,num+cnt+1,cmp);    for (int i=cnt;i>0;i--) if(num[i])      for (int j=i-1;j>0;j--) if (!(num[j]%num[i]))        num[j]=0;    for (int i=1;i<=cnt;i++) if (num[i]) num[++tot]=num[i];    dfs(1,0,1);    printf(AUTO,ans);    return 0;}