Răspuns:
Funcție sumă_divizori(n):
Suma <- 1 // Adăugăm 1 pentru că 1 este întotdeauna divizorul lui n
Pentru fiecare divizor d de la 2 până la √n:
Dacă n este divizibil cu d:
Suma <- Suma + d
Altfel dacă n / d este diferit de d: // Pentru a evita dublarea divizorilor
Suma <- Suma + n / d
Întoarce Suma
Funcție numere_prietene(a, b):
Pentru fiecare x de la a până la b:
y <- suma_divizori(x) - x // Calculăm suma divizorilor lui x și scădem x pentru a obține y
Dacă y > x și suma_divizori(y) - y == x:
Afișează x și y ca numere prietene