栏目分类:
子分类:
返回
名师互学网用户登录
快速导航关闭
当前搜索
当前分类
子分类
实用工具
热门搜索
名师互学网 > IT > 软件开发 > 后端开发 > C/C++/C#

数论-欧几里得最大公约数- same gcd

C/C++/C# 更新时间: 发布时间: IT归档 最新发布 模块sitemap 名妆网 法律咨询 聚返吧 英语巴士网 伯小乐 网商动力

数论-欧几里得最大公约数- same gcd

same gcd 题解 题意

​ 给定a,b,问有多少个x(0<=x 思路

d = g c d ( a , b ) = g c d ( a + x , b ) → 1 = g c d ( a + x d , b d ) d=gcd(a,b)=gcd(a+x,b)to 1=gcd(frac{a+x}{d},frac{b}{d}) d=gcd(a,b)=gcd(a+x,b)→1=gcd(da+x​,db​)

a + x d ∈ [ a d , a d + b d ) frac{a+x}{d}in[frac{a}{d},frac{a}{d}+frac{b}{d}) da+x​∈[da​,da​+db​)

题目转化成,求在上述范围内与b/d互质的个数,进而转化成求b/d的欧拉函数
e u l e r s [ 1 , b d ) ≡ e u l e r s [ a d , a d + b d ) eulers[1,frac{b}{d})equiv eulers[frac{a}{d},frac{a}{d}+frac{b}{d}) eulers[1,db​)≡eulers[da​,da​+db​)

代码
#include 

using namespace std;
typedef long long LL;

LL gcd(LL a, LL b)  // 欧几里得算法
{
    return b ? gcd(b, a % b) : a;
}

LL phi(LL x){
    LL res=x;
    for(int i=2;i<=x/i;i++){
        if(x%i==0)res=res/i*(i-1);
        while(x%i==0)x/=i;
    }
    if(x!=1)res=res/x*(x-1);
    
    return res;
}

int main()
{
    int t;
    cin>>t;
    while(t--){
        LL a,b;//注意数据范围
        cin>>a>>b;
        cout<
转载请注明:文章转载自 www.mshxw.com
本文地址:https://www.mshxw.com/it/887228.html
我们一直用心在做
关于我们 文章归档 网站地图 联系我们

版权所有 (c)2021-2022 MSHXW.COM

ICP备案号:晋ICP备2021003244-6号