OFFSET
1,2
COMMENTS
If p and q are different prime numbers and p + q is in A007850 (Giuga numbers) then m = p*q is a term because sigma(m) = sigma(p*q) = p*q + p + q + 1 and m + m'' = p*q + (p + q)' = p*q + p + q + 1 and sigma(m) = m + m''.
EXAMPLE
sigma(1) = 1 and 1 + 1'' = 1 so 1 is a term.
sigma(161) = 1 + 7 + 23 + 161 = 192 and 161 + 161'' = 161 + 30' = 161 + 31 = 192 so 161 is a term.
sigma(12690) = sigma(2*3^3*5*47) = 34560 and 12690 + 12690'' = 12690 + A068346(12690) = 12690 + 21870 = 34560 so 12690 is a term.
MATHEMATICA
d[0] = d[1] = 0; d[n_] := n * Plus @@ ((Last[#]/First[#]) & /@ FactorInteger[n]); Select[Range[160000], DivisorSigma[1, #] == # + d[d[#]] &] (* Amiram Eldar, Oct 18 2021 *)
PROG
(Magma) f:=func<n|n le 1 select 0 else n*(&+[Factorisation(n)[i][2] / Factorisation(n)[i][1]: i in [1..#Factorisation(n)]])>; [n:n in [1..160000]| DivisorSigma(1, n) eq n+Floor(f(Floor(f(n))))];
(PARI) ad(n) = vecsum([n/f[1]*f[2]|f<-factor(n+!n)~]); \\ A003415
isok(k) = sigma(k) == k+ad(ad(k)); \\ Michel Marcus, Oct 18 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Marius A. Burtea, Oct 18 2021
STATUS
approved