第十一届蓝桥杯大赛第二次模拟 5.反倍数

反倍数

【问题描述】

给定三个整数 a, b, c，如果一个整数既不是 a 的整数倍也不是 b 的整数倍还不是 c 的整数倍，则这个数称为反倍数。
请问在 1 至 n 中有多少个反倍数。

【输入格式】

输入的第一行包含一个整数 n。
第二行包含三个整数 a, b, c，相邻两个数之间用一个空格分隔。

【输出格式】

输出一行包含一个整数，表示答案。

【样例输入】

30

2 3 6

【样例输出】

10

【样例说明】

以下这些数满足要求：1, 5, 7, 11, 13, 17, 19, 23, 25, 29。

【评测用例规模与约定】

对于 40% 的评测用例，1 <= n <= 10000。

对于 80% 的评测用例，1 <= n <= 100000。

对于所有评测用例，1 <= n <= 1000000，1 <= a <= n，1 <= b <= n，1 <= c <= n。

题解：

#include<iostream>
#include<algorithm>
using namespace std;

int gys(int m,int n){
	if(n>m)                  //函数求最大公约数
		swap(m,n);
	while(int r=m%n){
		m=n;
		n=r;
	}
	return n;
}

int main(){
	int n,a[3],res=0;
	cin>>n;
	for(int i=0;i<3;i++)
		cin>>a[i];
	sort(a,a+3);             //将a,b,c按从小到大排序
	if(a[1]%a[0]==0){        //分不同情况进行讨论
		if(a[2]%a[0]==0)     //b,c都是a的倍数
			res=n-n/a[0];
		else                 //b是a的倍数但c不是a的倍数
			res=n-n/a[0]-n/a[2]+n/(a[0]*a[2]/gys(a[2],a[0]));
	}
	else{
		if(a[2]%a[0]==0)     //b不是a的倍数但c是a的倍数
			res=n-n/a[0]-n/a[1]+n/(a[0]*a[1]/gys(a[1],a[0]));
		else{
			if(a[2]%a[1]==0) //b,c都不是a的倍数但b是c的倍数
				res=n-n/a[0]-n/a[1]+n/(a[0]*a[1]/gys(a[1],a[0]));
			else             //b,c都不是a的倍数而且b也不是c的倍数
				res=n-n/a[0]-n/a[1]-n/a[2]+n/(a[0]*a[1]/gys(a[1],a[0]))+n/(a[0]*a[2]/gys(a[2],a[0]))+n/(a[1]*a[2]/gys(a[2],a[1]))-n/(a[0]*a[1]*a[2]/gys(gys(a[2],a[1]),a[0]));
		}                    //这种情况最为复杂
	}
	cout<<res<<endl;
	return 0;
}