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

我们发现纯幸运数字只有2046个，而去掉倍数之后只剩了943个。
我们直接容斥，暴力枚举就好了qaq
注意剪枝，从大到小枚举，求lcm避免爆ll等玄学技巧就好了qaq

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 2100
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(T==S){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline ll read(){
    ll x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
ll l,r,a[N],ans=0;
int n=-1;
bool ban[N];
inline ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}
inline void dfs1(ll x){
    if(x>r) return;a[++n]=x;
    dfs1(x*10+6);dfs1(x*10+8);
}
inline void dfs(int x,ll lcm,int cnt){
    if(lcm>r) return;
    if(cnt) ans+=(cnt&1)?r/lcm-(l-1)/lcm:-(r/lcm-(l-1)/lcm);
    for(int i=x;i<=n;++i){
        ll g=gcd(a[i],lcm);
        if(lcm/g>r/a[i]) continue;
        dfs(i+1,lcm/g*a[i],cnt+1);
    }
} 
int main(){
//  freopen("a.in","r",stdin);
    l=read();r=read();dfs1(0);
    sort(a+1,a+n+1);
    for(int i=1;i<=n;++i){
        if(ban[i]) continue;
        for(int j=i+1;j<=n;++j)
            if(!ban[j]&&a[j]%a[i]==0) ban[j]=1;
    }int m=n;n=0;for(int i=1;i<=m;++i) if(!ban[i]) a[++n]=a[i];
    reverse(a+1,a+n+1);dfs(1,1,0);
    printf("%lld\n",ans);
    return 0;
}