Taxicabnumbers.jl

# Ramanujanista ja G.H. Hardysta kertovan elokuvan " The Man Who Knew Infinity" innoittamana
# Julia-ohjelma Taxicabnumber.jl
# Juhani Kaukoranta 19.9.2019
# Lisäinfoa https://en.wikipedia.org/wiki/Taxicab_number

println("Taxicab luku T(n) on pienin positiivinen kokonaisluku, joka voidaan esittää kahden positiivisen kokonaisluvun kuutioiden summana n eri tavalla")
println("Ta(1)= 2")
println("Ta(2)= 1729 , Hardy-Ramanujan luku, liittyy matemaatikkojen Ramanujan ja G.H.Hardy tapaamiseen, oli Hardyn taksin numero")
println("Ta(3) = 87539319, löydettiin vuonna 1957 tietokoneen avulla")
println("Ta(4) = 6963472309248, löydettiin vuonna 1991 tietokoneen avulla")
println("Ta(5) = 48988659276962496, löydettiin vuonna 1994 tietokoneen avulla")
println("Ta(6) = 24153319581254312065344, löydettiin vuonna 2008 tietokoneen avulla")
println("Nyt voit etsiä annettua taxicab lukua vastaavan lukuparien esitykset")
println("Anna taxicab-luku; ")
inputtxt = readline()
taxicabnumber = parse(BigInt,inputtxt)

raja = 1+trunc(cbrt(taxicabnumber))

function Taxicab()
    temp = 0
    for i = 1:raja
        j = trunc(cbrt(taxicabnumber-i^3)+0.5)
        if i^3 + j^3 == taxicabnumber
              if i == temp
                  break
              end
              println(taxicabnumber," on lukujen ",convert(BigInt,i)," ja ",convert(BigInt,j), " kuutioiden summa ")
              temp = j
        end
    end
end

time0 = time()
Taxicab()
time1 = time()
println("Aikaa kului ",time1 - time0, "sekuntia")