1853: [Scoi2010]幸运数字

最新推荐文章于 2024-04-24 16:06:58 发布

ndsffx501ccy

最新推荐文章于 2024-04-24 16:06:58 发布

阅读量445

点赞数

CC 4.0 BY-SA版权

分类专栏：数学

本文链接：https://blog.youkuaiyun.com/ndsffx501ccy/article/details/31423579

数学专栏收录该内容

51 篇文章

订阅专栏

本文深入探讨了利用C++编程解决数论问题的方法，具体涉及到最大公约数计算、排列组合以及递归算法的应用。通过实例展示了如何运用这些数学概念解决实际问题，为读者提供了宝贵的编程技巧和理论知识。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<algorithm>
#include<iostream>
#include<cstdio>
#include<cmath>
#define ll long long
#define N 10001
using namespace std;
ll l,r,ans,a[N],b[N];
int t,n,m;
bool vis[N];
inline ll gcd(ll x,ll y){return y==0?x:gcd(y,x%y);}
void pre(int x,ll y){
	if(y>r)return;
	if(x>0)a[++m]=y;
	pre(x+1,y*10+6);
	pre(x+1,y*10+8);
}
void dfs(int x,int y,ll z){
	if(x>n){
		if(y&1)ans+=r/z-(l-1)/z;
		else if(y)ans-=r/z-(l-1)/z;
		return;
	}
	dfs(x+1,y,z);
	ll next=z/gcd(a[x],z);
	if((double)a[x]*next<=r)dfs(x+1,y+1,a[x]*next);
}
int main(){
	//scanf("%lld%lld",&l,&r);
	cin>>l>>r;
	pre(0,0);
	sort(a+1,a+m+1);
	for(int i=1;i<=m;i++)
		if(!vis[i]){
			b[++n]=a[i];
			for(int j=i+1;j<=m;j++)
				if(!(a[j]%a[i]))vis[j]=1;
		}
	for(int i=1;i<=n;i++)
		a[n-i+1]=b[i];
	dfs(1,0,1);
	cout<<ans;
	return 0;
}