2393: Cirno的完美算数教室
Time Limit: 10 Sec
Memory Limit: 128 MB
Submit: 406
Solved: 242
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组成的数,一共不超过2^11,然后把一些大的是小的的倍数的baka数去掉之后,还剩余大概400个。
容斥原理就可以了,容斥原理选的数>=3的时候,这三个数的最小公倍数会非常大而且超过R,我们就可以这样剪枝了,剪枝后的时间复杂度很理想
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define LL long long
const int M=10005;
LL L,R,s[15],f[M],tot=1,cnt,g[M],ans;
LL len(LL x)
{
int cnt=0;
while(x) x/=10,cnt++;
return cnt;
}
LL gcd(LL a,LL b)
{
if(b==0) return a;
return gcd(b,a%b);
}
void dfs(int deep,int flag,LL t)
{
if(deep>cnt)
{
if(t!=1) ans+=(R/t-(L-1)/t)*flag;
return;
}
dfs(deep+1,flag,t);
LL tmp=t/gcd(t,g[deep])*g[deep];
if(tmp>R) return;
dfs(deep+1,-flag,tmp);
}
int main()
{
scanf("%lld%lld",&L,&R);
s[1]=1;
for(int i=2;i<=10;i++) s[i]=s[i-1]*10;
f[1]=2,f[2]=9;
tot=2;
for(int i=2;i<=len(R);i++)
{
int now=tot;
for(int k=now-(1<<(i-1))+1;k<=now;k++)
f[++tot]=f[k]+s[i]*2,f[++tot]=f[k]+s[i]*9;
}
g[1]=2,g[2]=9;
cnt=2;
for(int i=3;i<=tot;i++)
{
bool flag=false;
for(int j=1;j<=cnt;j++)
if(f[i]%g[j]==0) {flag=true;break;}
if(!flag) g[++cnt]=f[i];
}
dfs(1,-1,1);
printf("%lld",ans);
}