python蓝桥杯 既约分数

学习点：辗转相除法求最大公约数 

如果一个分数的分子和分母的最大公约数是 11，这个分数称为既约分数。

例如 3/4 ,1/8, 7/1​ 都是既约分数。

请问，有多少个既约分数，分子和分母都是 1到 2020 之间的整数（包括 11 和 20202020）？

运行限制

• 最大运行时间：2s
• 最大运行内存: 128M

自己写代码的时候，求两个数的最小公约数是不是1总是报错还因为循环会超时，后来就干脆改了写求两个数的最小公约数，用辗转相除法求两个数的最大公约数。

参考python 最大公约数与最小公倍数_l_vip的博客-CSDN博客

def gcd(a,b):
    # a作为除数 必须大于b
    a, b = max(a,b),min(a,b) #或者写为 a, b = (a, b) if a >=b else (b, a)
    while b:
        a,b = b,a%b
    return a

 完整代码：

def gcd(a,b):
    # a作为除数 必须大于b
    a, b = (a, b) if a >=b else (b, a)
    while b:
        a,b = b,a%b
    return a

def num(x):
    m=0
    for i in range(x+1,2021): #从比x打的数开始算，防止重复计算
        if(gcd(i,x)==1):  #如果最大公约数是1，就加1（既约分数对+1）
            m+=1
    return m

M=[2019]   #1和后面所有的数都可以组成既约分数对，所以就放一个2019在列表里
for i in range(2,2021):  #从2开始循环求每个数和比他大的数可以组成的对
    m=num(i)
    M.append(m)
print(sum(M)*2)

但是，这里输出的结果报错，百思不得其解，上网搜了一下正确答案，发现就差1，知道了，是因为没有算 1/1 这个分数，最后再把1/1加进去。（2/2 3/3并不是既约分数）

def gcd(a,b):
    # a作为除数 必须大于b
    a, b = (a, b) if a >=b else (b, a)
    while b:
        a,b = b,a%b
    return a

def num(x):
    m=0
    for i in range(x+1,2021):
        if(gcd(i,x)==1):
            m+=1
    return m

M=[2019]
for i in range(2,2021):
    m=num(i)
    M.append(m)
print(sum(M)*2+1)