蓝桥杯每日一题:最大公约数(欧拉函数)

给定两个正整数 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≤1010。

输入样例:
复制代码
3
4 9
5 10
42 9999999967
输出样例:
复制代码
6
1
9999999966

解题思路:

gcd(a,m)== gcd(a+x,m) = d;

由于d为最大公约数所哟 a /= d,x/=d,m/=d后: a+x 与m互质(且x<m)可以转换为a~a+m-1中与m互质的个数。用数轴表述:

(m~a+m-1) % m == (0~a-1) % m;

问题转换为0-m-1中与m互质的个数,又因为0,m-1都与m不互质所以相当于去m的欧拉函数。

欧拉函数:

参考代码:

cpp 复制代码
#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;
typedef long long LL;
LL T,a,m;

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

LL ula(LL n)
{
    LL ans = n;
    for(int i=2;i<=n/i;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()
{
    cin>>T;
    while(T -- )
    {
        cin>>a>>m;
        
        m /= gcd(a,m);
        
        cout<<ula(m)<<endl;    
    }
    return 0;
}
相关推荐
一次旅行14 小时前
RLHF全链路深度解析:Reward Model数学推导+PPO完整实战,对比GRPO轻量化方案
人工智能·算法·机器学习
Wang's Blog14 小时前
AI Agent白手起家29: Few Shot 提示词工程实战
人工智能·算法
June`15 小时前
warp shuffle指令
c++·人工智能·算法·cuda
geovindu15 小时前
go:Bit Operation Algorithm
开发语言·后端·算法·golang·位运算法
乱七八糟的屋子15 小时前
【C++数值计算】NumCpp超详细入门教程(纯正Numpy语法、零学习成本、C++首选科学计算库、性能压测对比)
c++·数值计算·numpy·高性能计算·科学计算·ndarray
choumin15 小时前
C程序能调用extern “C“声明的带默认参数的C++函数吗?
c语言·c++·extern c·默认参数
AmyLin_200115 小时前
PDF 色彩保真工程实践【3】上手实践:工程结构、依赖集成与一键构建运行
c++·visualstudio·pdf·zlib·vcpkg·libtiff·cmyk
并不喜欢吃鱼16 小时前
从零开始 C++------ 十六.深度拆解 C++ 智能指针:unique_ptr/shared_ptr/weak_ptr 底层模拟、内存泄漏根治
开发语言·c++
Wang's Blog16 小时前
AI Agent白手起家31: 根据语义相似度动态选择 Few Shot 示例
人工智能·算法
不相心 -w-16 小时前
Cmake的基础用法
linux·开发语言·c++