BZOJ2393: Cirno的完美算数教室 容斥原理

2393: Cirno的完美算数教室

Time Limit: 10 Sec   Memory Limit: 128 MB

Submit: 406   Solved: 242

[Submit][Status][Discuss]

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