Fermatabc.jl

#Fermat problem
# a^2 + b^3 = c^7
# a,b,c positive coprime integers
function Fermatabc()
    println("Anna ylärajaluku, esim 200, 500, 1000: ")
    inputtxt = readline()
    time0 = time()
	n = parse(Int128,inputtxt)
    bn = Int128(trunc(1+n^(7.0/3.0)))
    a = Int128(0)
    b = Int128(0)
    c = Int128(0)
    a2 = Int128(0)
    for c = Int128(2) : n
        for b = Int128(3) : bn
            a2 = c*c*c*c*c*c*c-b*b*b
            if a2 > 0
                a = Int128(trunc(sqrt(a2)))
                if (a*a == a2 && gcd(a,gcd(b,c))==1)
                    println(a," ",b," ",c)
                end
            end
        end
    end
    time1 = time()
    println("Aikaa kului ",time1-time0," sekuntia")
end
Fermatabc()