bzoj2393 Cirno的完美算数教室 容斥

Description

---

~Cirno发现了一种baka数，这种数呢~只含有2和⑨两种数字~~
现在Cirno想知道~一个区间中~~有多少个数能被baka数整除~
但是Cirno这么天才的妖精才不屑去数啦
只能依靠聪明的你咯。

( 1 < L < R < 10^10)

Solution

---

第一眼以为是奇怪的数位dp
第二眼发现可以用容斥做。先跑一遍找出所有的baka数，然后搜索减去算重的lcm

Code

---

#include <stdio.h>
#include <algorithm>
#include <vector>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
typedef long long ll;
ll pre[1025],ans[1025],prt=0,l,r;
ll gcd(ll x,ll y) {return !y?x:gcd(y,x%y);}
bool vis[1025];
void dfs1(ll x,ll y) {
    if (y>r) return ;
    if (y) pre[++pre[0]]=y;
    dfs1(x+1,y*10+2);
    dfs1(x+1,y*10+9);
}
void dfs2(int dep,ll y,int v) {
    if (y>r) return ;
    if (dep==ans[0]+1) {
        if (y!=1) {
            prt+=(ll)(v*(r/y-(l-1)/y));
            // printf("%lld\n", v*(r/y-(l-1)/y));
        }
        return ;
    }
    dfs2(dep+1,y*ans[dep]/gcd(y,ans[dep]),v*-1);
    dfs2(dep+1,y,v);
}
int main(void) {
    scanf("%lld%lld",&l,&r);
    dfs1(0,0); std:: sort(pre+1,pre+pre[0]+1);
    rep(i,1,pre[0]) {
        if (!vis[i]) ans[++ans[0]]=pre[i];
        rep(j,i+1,pre[0]) if (pre[j]%pre[i]==0) vis[j]=1;
    }
    dfs2(1,1,-1);
    // puts("");
    printf("%lld\n", prt);
    return 0;
}