Divisores

s = 0
For j = 1 To k
If k / j = k \ j Then s = s + j ‘Este FOR-NEXT calcula la función sigma
de k
Next j
If s = n Then a = k: vale = False ‘comprueba si SIGMA coincide con el
argumento n
k = k + 1
Wend
mfsigma = a
End Function
Con esta función puedes determinar si un número coincide con la
función SIGMA de otro. Por ejemplo MF_SIGMA(2014)=0, luego no
existe ningún otro número cuya suma de divisores sea 2014. Si
embargo, MF_SIGMA(2012)=2011, porque este último es primo, y
MF_SIGMA(2016)=660, porque
2016=
660+330+220+165+132+110+66+60+55+44+33+30+22+20+15+12+11
+10+6+5+4+3+2+ 1
Puedes usar también la función definida en PARI
mfsigma1(n)={k=0;while(k=n,k=0); return(k)}
{print(mfsigma1(20))}
Con él, cambiando el valor de 20 por otro cualquiera, puedes encontrar
su MF_SIGMA
Las otras sigmas
Si sumamos los cuadrados de los divisores de un número nos resulta
la función SIGMA_2, con los cubos SIGMA_3 y, en general, podemos
definir toda la familia para exponentes mayores.
60