【QED】不想被排除

想要AC的dly

已于 2023-12-25 14:51:08 修改

阅读量891

点赞数 13

分类专栏： QEDfalse Online Judge 文章标签： c++ 算法数论

于 2023-12-09 19:34:43 首次发布

本文链接：https://blog.csdn.net/m0_67724631/article/details/134868424

版权

QEDfalse Online Judge 专栏收录该内容

13 篇文章 0 订阅

订阅专栏

题目描述

给出一个数 $n$ 以及 $x$ , $y$ , $z$ ，求 $1$ 到 $n$ 中，有多少个数不是 $x$ , $y$ , $z$ 中任意一个数的倍数。

输入格式

第一行输入一个整数 $T$ ，表示测试组数（ $\leq T \leq 10^{5}$ ）.

接下来会有 $T$ 行输入，每行 $4$ 个整数 $n$ , $x$ , $y$ , $z$ （ $\leq n \leq 10^{18}$ , $\leq x, y, z \leq 10^6$ ）.

输出格式

共 $T$ 行，每行输出一个整数，表示 $1$ 到 $n$ 中，不是 $x$ , $y$ , $z$ 中任意一个数倍数的整数个数。

测试样例

2
10 3 4 5
100 2 4 6

3
50

温馨提示

递归求解int范围内两数的最大公约数的代码模板如下。

// C/C++ code
int gcd(int a, int b) {
    return b == 0 ? a : gcd(b, a % b);
}

// Java code
public int gcd(int a, int b){
    return b == 0 ? a : gcd(b, a % b);
}

# Python code
def gcd(a, b):
    if b == 0:
        return a
    return gcd(b, a % b)

Idea from HGG.

思路

直接求“不是”的个数比较难求，但是我们这里可以转化成求“是x、y、z中任意一个数的倍数”的个数，那么现在就是一个容斥问题
在这里插入图片描述
注意个数的取值范围，求个数的时候向下取整

核心代码

#include <iostream>
using namespace std;

inline long long gcd(long long a,long long b) 
{    
    if(b) while((a%=b) && (b%=a));    
    return a+b;
}

void solve()
{
	long long n,x,y,z;
	scanf("%lld%lld%lld%lld",&n,&x,&y,&z);
	long long a=x/gcd(x,y)*y;
	long long b=x/gcd(x,z)*z;
	long long c=y/gcd(y,z)*z;
	long long d=a/gcd(a,z)*z;
//	cout<<a<<"\n"<<b<<"\n"<<c<<"\n"<<d<<"\n"; 
	cout<<n-n/x-n/y-n/z+n/a+n/b+n/c-n/d<<"\n";
}

int main() {
	int t=1;
	cin>>t;
	while(t--)
	{
		solve();
	}
	return 0;
}