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;}