Gausstotientit.py

# Totientit(1...n) käyttäen Gaussin
import numpy as np # käytetään numpy matriisia(vektoria)
from numba import njit #numba-kääntäjän kiihdytys
@njit
def Gausstotientit(n):
    phi =np.zeros(n+1,dtype=np.int32)
    phi[0]=0
    phi[1]=1
    for i in range(2,n+1):
        phi[i] = i-1
    for i in range(2,n+1):
        for j in range(2*i,n+1,i):
            phi[j] -= phi[i]
    return phi[1:n+1]