Gọi ƯCLN(a, b) là hàm ƯCLN của hai số tự nhiên a, b. Dễ thấy ta có

def UCLN(a,b):

    r = a % b

    while r != 0:

        a = b

        b = r

        r = a % b

    return b

a=int(input("Nhập số tự nhiên thứ nhất"))

b=int(input("Nhập số tự nhiên thứ hai"))

print("ƯCLN của a và b là: ",UCLN(a,b))