【AcWing】蓝桥杯集训每日一题Day30|欧拉函数|3999.最大公约数(C++)

3999.最大公约数

3999. 最大公约数 - AcWing题库
难度：困难
时/空限制：1s / 256MB
总通过数：1031
总尝试数：1671
来源： AcWing第21场周赛
算法标签 数论欧拉函数

---

题目内容

给定两个正整数 a,m，其中 a<m。
请你计算，有多少个小于 m 的非负整数 x 满足：
$$gcd(a,m)=gcd(a+x,m)$$

输入格式

第一行包含整数 T，表示共有 T 组测试数据。
每组数据占一行，包含两个整数 a,m。

输出格式

每组数据输出一行结果，一个整数，表示满足条件的非负整数 x 的个数。

数据范围

前三个测试点满足，1≤T≤10。
所有测试点满足，1≤T≤50，1≤a<m≤10^10。

输入样例：

3
4 9
5 10
42 9999999967

输出样例：

6
1
9999999966

题目解析

数论问题
设最大公约数是d
$$d=gcd(a,m)=gcd(a+x,m)$$
由于d是它们的最大公约数
d能够整除a，也能够整除m，也能够整除a+x
因为d可以整除a，也可以整除a+x
所以d可以整除x

技巧

可以把每一个数除一个最大公约数，除完之后，剩余的部分就互质了
$$a^{\prime }=\frac{a}{d}$$
在a当中把最大公约数除干净
$$x^{\prime }=\frac{x}{d}$$
$$m^{\prime }=\frac{m}{d}$$
等价于a’和m‘的最大公约数是1
$$(a^{\prime },m^{\prime })=(a^{\prime }+x^{\prime },m^{\prime })=1$$
从0到m之间有多少个x是满足这个式子的
由于x和x’是一一对应的

问题就变成了

有多少个x‘满足$0\le x^{\prime }<m^{\prime }$且$(a^{\prime }+x^{\prime },m^{\prime })=1$
如果a+x的最大公约数等于a和m的最大公约数，等价于a’+x‘和m’是互质的

如果x‘在这个范围内，并且x满足这两组数是互质的
x就会在0~m-1之间，并且a+x和m的最大公约数就是d

如何求有多少个x’满足它们是互质的

求最大公约数的时候用的是辗转相除法
$$(a,b)=(aMODb,b)$$

任取一个a’+x’标记位t，t和m‘互质
等价于t%m’和m‘互质

要求有多少个x’满足互质的话，等价于有多少个x‘满足$((a^{\prime }+x^{\prime })MOD{m}^{\prime },m^{\prime })=1$
模完m’之后，对于所有的x‘，形式比较简单

x’是在0到m-1当中
a’+x’就应该在a’到a’+m-1这个范围内
a’+x’每个数模上m‘之后，余数恰好就是0到m’-1

任意连续m‘个数，在模m的意义下，都和0~m’-1是同余的

问题就变成了

在0~m‘-1当中，有多少个数和m’互质
等价于1到m‘之间有多少个数和m’互质
就是m’的欧拉函数

时间复杂度$O(\sqrt{n})$，500万的计算量

代码

#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;

typedef long long LL;

LL gcd (LL a, LL b)
{
	return b ? gcd(b, a % b) : a;
}

LL phi(LL m)
{
	LL res = m;
	//先把m分解质因数
	for (LL i = 2; i <= m / i; i ++)
		//能整除，说明i是质因子
		if (m % i == 0)
		{
			//把质因子整除干净
			while (m % i == 0) m /= i;
			res = res / i * (i - 1);
		}
	//如果m大于1的话，表示有一个大于根号m的质因子
	if (m > 1)
		res = res / m * (m - 1);
	return res;
}

int main()
{
	int T;
	cin >> T;

	while (T --)
	{
		LL a, m;
		cin >> a >> m;
		//求一下最大公约数
		LL d = gcd(a, m);
		//m/d就是m'
		cout << phi(m / d) << endl;
	}

	return 0;
}