Terras_modified_Collatz.py

# Collatzin konjektuuri
# Terras-modified algoritm : not n = 3*n+1, but n = (3*n+1)//2.
# Terras-modified algoritm shortens collatz number steps
# Juhani Kaukoranta 17.3.2019
# Jonon luvut tallennetaan numpy-taulukkoon

import numpy as np

def Collatz(n):
    x = np.array([n])
    while n > 1:
        if n % 2 == 0:
            n = n//2
        else:
            n = (3*n + 1) // 2
 # lisätään saatu luku n numpy-taulukkoon x
        x = np.append(x,n)
#  Collatz-jono on  valmiina palautettavaksi
    return x

luku = int(input("Anna joku positiivinen kokonaisluku "))
# numpy-array collatzjono sisältää Collatzin konjektuurin luvut
terras_collatzjono = Collatz(luku)
print("Terras-modified Collatz-jonossa alkioita ",terras_collatzjono.size,", jonon alkiot ovat:")
print(terras_collatzjono)