Taxicabnumbers.py

# Ramanujanista ja G.H. Hardysta kertovan elokuvan " The Man Who Knew Infinity" innoittamana
# Juhani Kaukoranta 19.9.2019
# Lisäinfoa https://en.wikipedia.org/wiki/Taxicab_number
import math
import time
print("Taxicab luku T(n) on pienin positiivinen kokonaisluku, joka voidaan esittää kahden positiivisen kokonaisluvun kuutioiden summana n eri tavalla")
print("Ta(1)= 2")
print("Ta(2)= 1729 , Hardy-Ramanujan luku, liittyy matemaatikkojen Ramanujan ja G.H.Hardy tapaamiseen, oli Hardyn taksin numero")
print("Ta(3) = 87539319, löydettiin vuonna 1957 tietokoneen avulla")
print("Ta(4) = 6963472309248, löydettiin vuonna 1991 tietokoneen avulla")
print("Ta(5) = 48988659276962496, löydettiin vuonna 1994 tietokoneen avulla")
print("Ta(6) = 24153319581254312065344, löydettiin vuonna 2008 tietokoneen avulla")
print("Nyt voit etsiä annettua taxicab lukua vastaavan lukuparien esitykset")
taxicabnumber = int(input("Anna taxicab-luku :"))
print("taxi=  ",taxicabnumber)
raja = 1+int(taxicabnumber**(1.0/3))

def Taxicab():
    temp = 0
    for i in range(1,raja):
        j = int((taxicabnumber - i ** 3) ** (1.0 / 3) + 0.5)
        if (i**3 + j**3 == taxicabnumber):
            if i == temp:
                break
            print(taxicabnumber," on lukujen ",i," ja ",j, "kuutioiden summa")
            temp = j

time0 = time.perf_counter()
Taxicab()
time1 = time.perf_counter()
print("Aikaa kului ",time1 - time0, "sekuntia")