当前位置：网站首页>【BZOJ2393】Cirno的完美算数教室

【BZOJ2393】Cirno的完美算数教室

2022-07-19 11:47:00 【LauJiYeoung】

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

弱化版的SCOI2010幸运数字

还是考虑值的增长率很快所以暴力容斥就好了

#include<bits/stdc++.h>
using namespace std;
typedef int INT;
#define int long long
int L,R;
int Baka[3000];
int vis[3000];
int A[3000];
int cnt=0;
void DFS(int x,int sum){
	if(sum>R)return;
	if(sum)Baka[++cnt]=sum;
	DFS(x+1,sum*10+2);
	DFS(x+1,sum*10+9);
}
int ans=0;
int GCD(int x,int y){
	while(y){
		int tmp=y;
		y=x%y;
		x=tmp;
	}
	return x;
}
void Solve(int x,int sum,int flag){
	if(sum<0||sum>R)return;
	if(x==A[0]+1){
		if(!sum)return;
		ans=(ans+(R/sum-L/sum)*flag);
		return;
	}
	int tmp;
	if(!sum)tmp=A[x];
	else tmp=sum/GCD(A[x],sum)*A[x];
	Solve(x+1,tmp,flag*-1);
	Solve(x+1,sum,flag);
}
bool cmp(int A,int B){
	return A>B;
}
signed main(){
	cin>>L>>R;
	L--;
	DFS(1,0);
	sort(Baka+1,Baka+1+cnt);
	for(int i=1;i<=cnt;++i){
		if(!vis[i])A[++A[0]]=Baka[i];
		for(int j=i+1;j<=cnt;++j){
			if(Baka[j]%Baka[i]==0)vis[j]=1;
		}
	}	
	sort(A+1,A+1+A[0],cmp);
	Solve(1,0,-1);
	cout<<ans;
}

原网站

版权声明
本文为[LauJiYeoung]所创，转载请带上原文链接，感谢
https://blog.csdn.net/fcb_x/article/details/82990179

边栏推荐

猜你喜欢

随机推荐