【bzoj1853】: [Scoi2010]幸运数字 数论-容斥原理

【bzoj1853】: [Scoi2010]幸运数字

预处理出所有幸运数字然后容斥原理

但是幸运数字是2^logn个数的 直接搞会炸

所以把成倍数的处理掉

然后发现还是会T 所以数字要从大到小处理会快很多 我也不知道为什么

还有就是剪枝的时候要开double 不然爆longlong

/* http://www.cnblogs.com/karl07/ */
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <map>
#include <algorithm>
using namespace std;
#define ll long long

const ll N=10000000005;
ll num[5000];
ll a,b,cnt,ans,cc;

ll gcd(ll a,ll b){ return b==0 ? a : gcd(b,a%b);}
ll lcm(ll a,ll b){ return a/gcd(a,b)*b;}

void dfs(ll x,ll n){
    if (x>n) return;
    num[++cnt]=x;
    dfs(x*10+6,n);
    dfs(x*10+8,n);
}

void rc(ll c,ll x,ll n,ll nw){
    if (x==0){
        if (c==0) return;
        if (c&1) ans+=n/nw; else ans-=n/nw;
        return;
    }
    if ((double)nw/gcd(num[x],nw)*num[x] <= n) rc(c+1,x-1,n,lcm(nw,num[x]));
    rc(c,x-1,n,nw);
}

ll work(ll n){
    if (n<6) return 0;
    cnt=0;ans=0;
    dfs(6,n);dfs(8,n);cc=cnt;
    num[0]=1;
    sort(num+1,num+cnt+1);
    for (int i=1;i<=cnt;i++){ 
        for (int j=i+1;j<=cnt && num[i]!=N;j++){
            if (num[j]%num[i]==0) num[j]=N,cc--;
        }
    }
    sort(num+1,num+cnt+1);
    rc(1,cc-1,n,num[cc]);
    rc(0,cc-1,n,1);
    return ans;
}

int main(){
    scanf("%lld%lld",&a,&b);
    printf("%lld
",work(b)-work(a-1));
    return 0;
}