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