【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

题解：同2393只是数据比较强,几个数乘的时候会爆longlong,所以需要用long double

代码：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#define N 5010
using namespace std;
long long a[N],b[N],vis[N],ans,l,r;
long long gcd(long long a,long long b){return b?gcd(b,a%b):a;}
void pro(long long x){if (x>r) return;a[++a[0]]=x;pro(x*10+6);pro(x*10+8);} 
void dfs(int p,int f,long long v){
    if(!p){if(v!=1)ans+=(r/v-(l-1)/v)*f;return;}
    dfs(p-1,f,v);
    long long temp=b[p]/gcd(v,b[p]);if((long double)temp*v>r) return;
    dfs(p-1,-f,temp*v);
}
int main(){
  scanf("%lld%lld",&l,&r);pro(6);pro(8);sort(a+1,a+a[0]+1); 
  for (int i=1;i<a[0];i++)
   for (int j=i+1;j<=a[0];j++)
     if (a[j]%a[i]==0) vis[j]=1;
  for (int i=1;i<=a[0];i++) if (!vis[i]) b[++b[0]]=a[i]; 
  dfs(b[0],-1,1);
  cout<<ans<<endl;
}