Calcul de suites python

Bonjour,
J'essaie de résoudre un calcul en python qui est le suivant:

U1 = (X**2 - 2)%N
U2 = (X*U1 -X)%N
Un = (X* Un-1 - Un-2)%N

Je connais la valeur de Un pour n=45000 par exemple et je dois trouver la valeur de X correspondante.
La valeur N est connue.
je souhaite simplement savoir s'il est possible de résoudre cela en python, car je n'ai trouvé aucune manière mathématique de la faire à part l'utilisation d'une boucle (mais X étant trop grand, ça ne donne rien)?

merci d'avance