bzoj1853: [Scoi2010]幸运数字（容斥原理）

题目传送门
这道题。。
数据有点强啊。。。。

解法：
我好垃圾我一开始的想法如下：
先把幸运数字都求出来。。
然后用容斥原理去搜。。
奇数个的时候就加上答案，偶数个的时候就减去答案。

结果。。
数据规模十的十次方。
那么有两千个幸运数字诶。
容斥原理2的两千次方。
T到死。。。。。

想了很久，从大到小会不会常数稍微小一点呢。。
（以前做过那种爆搜填木桶的题，从大到小可以AC从小到大要T）
然后打了个代码。
T了。
T了很多次。。
找不出原因啊上网%题解呀。。。。
发现还有这种操作？！
dfs里面求最小公倍数会炸。。。
要把它转成double类型然后再判断？！
不懂，求神犇指教。

代码实现：

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
typedef long long ll;
ll st,ed;int len1,len2;
ll a[11000],b[11000];
ll s[11000];int len;
void dfs2(ll x) {if(x>ed)return ;if(x!=0)a[++len1]=x;dfs2(x*10+6);dfs2(x*10+8);}
bool cmp(ll n1,ll n2) {return n1<n2;}
bool cmp1(ll n1,ll n2) {return n1>n2;}
ll get_sum(ll x) {return ed/x-(st-1)/x;}
ll gcd(ll x,ll y) {if(x==0)return y;return gcd(y%x,x);}
ll lcu(ll x,ll y) {if(x>y)swap(x,y);return x/gcd(x,y)*y;}
ll ans=0;bool v[11000];
void dfs(ll Gcd,int x,int k) {
    if(x==len+1) {
        if(k==0)return ;
        if(k%2==1) ans+=get_sum(Gcd);
        else ans-=get_sum(Gcd);
        return ;
    }
    dfs(Gcd,x+1,k);
    ll t=Gcd/gcd(Gcd,s[x]);
    if((double)t*s[x]<=ed)dfs(lcu(s[x],Gcd),x+1,k+1);
}
int main() {
    scanf("%lld%lld",&st,&ed);
    len1=0;dfs2(0);
    sort(a+1,a+1+len1,cmp);
    memset(v,true,sizeof(v));
    for(int i=1;i<=len1;i++)if(v[i]==true)for(int j=i+1;j<=len1;j++)if(a[j]%a[i]==0)v[j]=false;   //如果大的包含小的那我们肯定保留小的呀，因为小的迟早会覆盖到大的。
    len=0;for(int i=1;i<=len1;i++)if(v[i]==true)s[++len]=a[i];
    sort(s+1,s+1+len,cmp1);
    ans=0;dfs(1,1,0);printf("%lld\n",ans);
    return 0;
}