现在的位置: 首页 > 综合 > 正文

2393: Cirno的完美算数教室

2018年01月13日 ⁄ 综合 ⁄ 共 694字 ⁄ 字号 评论关闭
#include<algorithm>
#include<iostream>
#include<cstdio>
#include<cmath>
#define ll long long
#define N 2001
using namespace std;
int l,r,t,m,n,ans;
ll a[N],b[N];
bool vis[N];
inline ll gcd(ll x,ll y){return y==0?x:gcd(y,x%y);}
void pre(int x,int y){
	if(x>t||y>r)return;
	if(x>0)a[++m]=y;
	pre(x+1,y*10+2);
	pre(x+1,y*10+9);
}
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=a[x]*z/gcd(a[x],z);
	if(next<=r)dfs(x+1,y+1,next);
}
int main(){
	scanf("%d%d",&l,&r);
	t=(int)(log(r)/log(10))+1;
	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);
	printf("%d",ans);
	return 0;
}

容斥原理

抱歉!评论已关闭.