Exponentiation modulaire tres rapide

Bonjour,
Savez vous comment cette algorithme fonctionne réellement?
Je cherche a savoir comment ça se fait que après que la fonction soit appelé jusqu'à k=1 alors elle arrive a revenir en arrière et retrouver son chemin pour donner la bonne solution (je cherche donc des caractéristiques propres a python sur les fonctions récursives)
Si quelqu'un peut m'expliquer le principe de cette pile d'éxécution (si c'est bien cela), j'en serait ravi.

Voici le code:

import time
a=int(input())
k=int(input())
n=int(input())
hey=time.time() #rajouter pour obtenir le temps
def puis(a,k,n):
    if k==0:
       return 1
    b=puis(a,(k//2),n)
    if k%2==0:
        return ((b*b)%n)
    return ((b*b*a)%n)
print(puis(a,k,n))
ho=time.time()
print(ho-hey)

Merci d'avance