Prolog - suma de divisores

hola amigos :

como puedo realizar una aplicacion para calcular la suma de divisores de un numero dado

sumdiv(10,X)

X=18

gracias

cual es la formula de la suma de divisores de un numero

suma(a,b,c) :- b>a, c = 0,!
suma(a,b,c) :- a mod b = 0, b1 = b + 1, suma (a,b1,c1), c = c1 + b, !
suma(a,b,c) :- b1 = b + 1, suma (a,b1,c1)

Supongamos el nº N factorizado:
N = (p1)^e1·(p2)^e2·(p3)^e3·...·(pn)^en
La suma de sus divisores viene dada por la función Sigma(N) =
= [[(p1)^(e1+1)]-1]/(p1-1)]]·[[(p2)^(e2+1)]-1]/(p2-1)]]·...·[[(pn)^(en+1)]-1]/(pn-1)]]
Ejemplo:
18=2·(3^2)
Sigma(18)=[[(2^(1+1))-1]/(2-1)]·[[(3^(2+1))-1]/(3-1)] = (3/1)·(26/2) = 39
Puedes comprobarlo:
Los divisores de 18 son: 1,2,3,6,9 y 18. Si los sumas: 1+2+3+6+9+18 = 39
Ejemplo:
12=(2^2)·3
Sigma(12)=[[(2^(2+1))-1]/(2-1)]·[[(3^(1+1))-1]/(3-1)] = (7/1)·(8/2) = 28
Puedes comprobarlo:
Los divisores de 12 son: 1,2,3,4,6 y 12. Si los sumas: 1+2+3+4+6+12 = 28

necesito hacer un programa que lea 2 numeros y me indique que son amigos por la sumaa de sus divisores

okeii laa formulaa utilizada generalmentee es: pichi_cacaa