2018.10.05 bzoj2393: Cirno的完美算数教室（容斥原理+搜索）

最新推荐文章于 2021-12-21 17:05:17 发布

原创最新推荐文章于 2021-12-21 17:05:17 发布 · 315 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#数论

收录于

容斥原理同时被 2 个专栏收录

32 篇文章

订阅专栏

搜索

22 篇文章

订阅专栏

本文解析了一道经典的算法竞赛题目，通过分析题目特性，指出了解题关键在于处理数的倍数关系，并利用容斥原理与爆搜剪枝技巧求解。文章详细介绍了算法思路与实现代码，适合对算法竞赛感兴趣的学习者。

传送门
经典题目。
显然满足题意的数最多不超过1024个。
然后对于两个数 $a, b$ ，如果 $a$ 是 $b$ 的倍数，那么就不必计算 $a$ 的贡献。
处理出来之后容斥原理+爆搜剪枝就能过了。
代码：

#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll L,R,divv[1050],Div[1050],ans=0,f=1;
bool BK;
int tot=0,cnt=0;
inline void dfs(ll sum){
	if(sum>R)return;
	if(sum)divv[++tot]=sum;
	dfs(sum*10+2),dfs(sum*10+9);
}
inline ll gcd(ll a,ll b){while(b){ll t=a;a=b,b=t%a;}return a;}
inline void dfs(int pos,int dep,ll mul){
	if(mul>R)return;
	if(!dep){if(mul)ans+=f*(R/mul-(L-1)/mul),BK=0;return;}
	if(pos==cnt+1)return;
	dfs(pos+1,dep,mul),dfs(pos+1,dep-1,mul/gcd(mul,Div[pos])*Div[pos]);
}
int main(){
	cin>>L>>R,dfs(0),sort(divv+1,divv+tot+1);
	for(int i=1;i<=tot;++i){
		bool f=true;
		for(int j=i-1;j;--j)if(divv[i]%divv[j]==0){f=false;break;}
		if(f)Div[++cnt]=divv[i];
	}
	for(int i=1;i<=cnt;++i,f*=-1){
		BK=1,dfs(1,i,1);
		if(BK)break;
	}
	cout<<ans;
	return 0;
}