BZOJ 2393 容斥原理

Description

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

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

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

1 100
Sample Output

58
HINT

分析：范围为10^10，实际上我们就要处理10^10中的只含有2和9的数，然后在一个区间内枚举它们的倍数，看有多少个。然后不难看出如果直接算的话会有重复，比如lcm(2,9)会被算很多次，于是考虑容斥。ans=ans(2)+ans(9）-ans(lcm(2,9))+……-/+ans；但我们发现如果后面出现的2，9数如果是前面2,9数的倍数，那它肯定已经被算进去了，不需要再算一遍，于是把它们踢掉。考虑一个小优化，若lcm(a,b)已经大于了区间右端，则返回。

# include <iostream>
# include <cstdio>
# include <cstring>
# include <cmath>
# include <list>
# include <map>
# include <queue>
# include <algorithm>
using namespace std;
typedef long long ll;
ll read()
{
    ll f=1,i=0;char ch=getchar();
    while(ch<'0'||ch>'9') 
    {
        if(ch=='-') f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        i=(i<<3)+(i<<1)+ch-'0';
        ch=getchar();
    }
    return f*i;
}
vector<ll>q,tmp;
int l,r,m,ans;
bool vis[100010];
inline void initial(int x)
{
    if(x>r) return ;
    if(x) tmp.push_back(x);
    initial(x*10+2);initial(x*10+9);
}
int gcd(int a,int b)
{
    return b==0 ? a : gcd(b,a%b);
}
int lcm(int a,int b)
{
    return a/gcd(a,b)*b;
}
inline void dfs(int nowp,int pos,int f)
{
    if(nowp==q.size())
    {
        if(pos!=1) ans+=f*(r/pos-(l-1)/pos);
        return ;
    }
    dfs(nowp+1,pos,-f);
    if(lcm(pos,q[nowp])<r) dfs(nowp+1,lcm(pos,q[nowp]),f);
}
int main()
{
    l=read(),r=read();
    initial(0);
    sort(tmp.begin(),tmp.end());
    for(int i=0;i<tmp.size();++i)
    {
        if(vis[i]) continue;
        q.push_back(tmp[i]);
        for(int j=i+1;j<tmp.size();++j)
            if(tmp[j]%tmp[i]==0) vis[j]=true;
    }
    dfs(0,1,1);
    cout<<ans<<endl;    
}