bzoj2396 Cirno的完美算数教室 容斥

Description

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

一行正整数L R
( 1 < L < R < 10^10)
Output

一个正整数，代表所求的答案
Sample Input

1 100
Sample Output

58

一开始我以为只是2和9，而且一算样例也是对的，内心OS：我去这不大水题？简直堪比A+B啊。
本着对bzoj的题目质量的信任，重新看了一遍题，才发现自己看错了。。原来29,92之类的都行。这就有点坑了，感觉也不能用数论知识。
然后才发现能用容斥。。因为这样的数在10^10以内也就几百个，直接用容斥暴力算，集合内有奇数个就加，偶数就减。

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
int l,r,t,n,m;
typedef long long ll;
const int N=1e5;
int a[N],tot=0,vis[N],b[N],cnt=0;
ll ans=0;
inline void pre(int x,int y)
{
    if(x>t||y>r)return ;
    if (x>0)a[++tot]=y;
    pre(x+1,y*10+2);
    pre(x+1,y*10+9);
}
inline int gcd(int x,int y)
{
    if (x%y==0)return y;
    else return gcd(y,x%y);
}
inline void dfs(int x,int y,ll z)
{
    if (x>n)
    {
        if (y%2==1)ans+=r/z-(l-1)/z;
        else if (y)ans-=r/z-(l-1)/z;
        return;
    }
    dfs(x+1,y,z);
    ll lcm=a[x]*z/gcd(a[x],z);
    if (lcm<=r)
    dfs(x+1,y+1,lcm);
}
int main()
{
    scanf("%d%d",&l,&r);
    t=log(r)/log(10)+1;
    pre(0,0);
    sort(a+1,a+tot+1);
    fo(i,1,tot)
    {
        if (!vis[i])
        b[++n]=a[i];
        fo(j,i+1,tot)
        if (a[j]%a[i]==0)vis[j]=1;
    }
    fo(i,1,n)
    a[n-i+1]=b[i];
    dfs(1,0,1);
    printf("%lld\n",ans);
}