D. Same GCDs(欧拉函数)

Problem - 1295D - Codeforces

给你两个整数a和m. 计算整数x的数量，使0≤x<m且gcd(a,m)=gcd(a+x,m).

注意：gcd(a,b)是a和b的最大公除数。

输入
第一行包含单个整数T（1≤T≤50）--测试用例的数量。

接下来的T行包含测试案例--每行一个。每行包含两个整数a和m（1≤a<m≤1010）。

输出
打印T个整数 - 每个测试案例一个。对于每个测试案例，打印适当的x-s的数量。

例子
输入复制
3
4 9
5 10
42 9999999967
输出拷贝
6
1
9999999966
注意
在第一个测试案例中，适当的x-s是[0,1,3,4,6,7]。

在第二个测试案例中，唯一合适的x是0。

题意:给a,m 求gcd(a,m) = gcd(a+x,m)  0 <=x <= m

题解:根据gcd(a,b) == gcd(a-b,b)

所以a+x>=m

gcd(a+x-m,m) = gcd(a+x,m)

k = (a+x)mod m

gcd(k,m) = d

可得gcd(k/d,m/d) = 1

题意转换为了求m/d的欧拉函数(就是对于一个正整数n，小于n且和n互质的正整数（包括1）的个数，记作φ(n) 。)

//#include<iostream>
//#include<cmath>
//#include<math.h>
//#include<algorithm>
//#include<cstdio>
#include<bits/stdc++.h>
using namespace std;
int main()
{
	int t;
	scanf("%d",&t);
	while(t--)
	{
		long long a,m;
		cin >> a>> m;
		long long k = __gcd(a,m);
		long long x = m/k;
		long long ans = x;
		for(long long i = 2;i * i <= x;i++)
		{
			if(x%i == 0)
			{
				ans = ans/i*(i - 1);
				while(x%i==0)
				{
					x = x/i;
				}
			}
		}
		if(x > 1)
		{
			ans = ans/x*(x-1);
		}
		cout <<ans <<"\n";
	}
}