题目链接:点击查看
题目大意:给出两个正整数 a 和 m ,再给出 x 的范围为[ 0 , m ),现在要求满足 gcd ( a , m ) = gcd ( a + x , m ) 时 x 的个数
题目分析:因为辗转相除法的定义:gcd(a,b) = gcd(b,a mod b)可知,当 i ∈ [ a , m ] 和 i ∈ [ m , 2 * m ] 时,gcd( i ,m ) 相等的 i 的个数相同,因为每次都是对第二个数取模,所以具有周期性,而题目要求我们求得是 i ∈ [ a , m + a -1 ] 时的个数,因为长度为一个完整的周期,所以我们不妨将其问题转换为 i ∈ [ 0 , m ] 时,满足gcd( i , m ) = = gcd ( a , m ) 的个数
剩下的就能用欧拉函数解决了:
也就是求m / gcd ( a , m ) 的欧拉函数
代码:
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
//ll eular(ll n)
//{
// ll ans=n;
// for(int i=2;i*i<=n;i++)
// {
// if(n%i==0)
// {
// ans=ans/i*(i-1);
// while(n%i==0)
// {
// n/=i;
// }
// }
// }
// if(n>1)
// {
// ans=ans/n*(n-1);
// }
// return ans;
//}
ll eular(ll n)
{
ll ans = n;
for(ll i=2; i*i <= n; ++i)
{
if(n%i == 0)
{
ans = ans/i*(i-1);
while(n%i == 0)
n/=i;
}
}
if(n > 1)
ans = ans/n*(n-1);
return ans;
}
int main()
{
int t;
cin>>t;
while(t--)
{
ll a,m;
scanf("%lld%lld",&a,&m);
m/=__gcd(a,m);
printf("%lld\n",eular(m));
}
return 0;
}
需要把欧拉函数内for循环里的int变为long long
原文链接:https://blog.csdn.net/qq_45458915/article/details/104114132
欧拉函数:
就是对于一个正整数n,小于n且和n互质的正整数(包括1)的个数,记作φ(n) 。
欧拉函数的通式:φ(n)=n*(1-1/p1)(1-1/p2)(1-1/p3)*(1-1/p4)……(1-1/pn)
其中p1, p2……pn为n的所有质因数,n是不为0的整数。φ(1)=1(唯一和1互质的数就是1本身)。
所以,根据通式我们可以打出以下代码:
ll eular(ll n)
{
ll ans = n;
for(ll i=2; i*i <= n; ++i)
{
if(n%i == 0)
{
ans = ans/i*(i-1);
while(n%i == 0)
n/=i;
}
}
if(n > 1) ans = ans/n*(n-1);
return ans;
}