PROG
(MAGMAMagma) [&+[d: d in Divisors(n) | IsIntegral(SumOfDivisors(d) / d)] : n in [1..100]]
(MAGMAMagma) [&+[d: d in Divisors(n) | IsIntegral(SumOfDivisors(d) / d)] : n in [1..100]]
reviewed
approved
proposed
reviewed
editing
proposed
Antti Karttunen, <a href="/A325470/b325470.txt">Table of n, a(n) for n = 1..20000</a>
proposed
editing
editing
proposed
a[n_] := DivisorSum[n, # &, Divisible[DivisorSigma[1, #], #] &]; Array[a, 100] (* Amiram Eldar, Aug 17 2019 *)
proposed
editing
editing
proposed
(PARI) a(n)={sumdiv(n, d, if(sigma(d) % d == 0, d))} \\ Andrew Howroyd, Aug 16 2019
proposed
editing
editing
proposed