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;
}