最大公約数を求める

苦労した結果こうなりましたが、あまり綺麗にはなりませんでした。

a, b = input("自然数aは？:"), input("自然数bは？:")
for x in range(2, int(a) + 1):
    if int(a) % x == 0:
        if int(b) % x == 0:
            print(x)
        else:
            print("互いに素です")
            break
    else:
        print("aは素数です")
        break