OFFSET
1,8
FORMULA
a(p^k)) = k - 1, for p prime and k >= 1.
EXAMPLE
MATHEMATICA
a[n_] := DivisorSum[n, 1 &, CompositeQ[#] && CoprimeQ[#, DivisorSigma[1, #]] &]; Array[a, 100] (* Amiram Eldar, Jul 19 2024 *)
PROG
(Magma) f:=func<n|n ne 1 and not IsPrime(n) and Gcd(n, DivisorSigma(1, n)) eq 1>; [#[d:d in Divisors(k)|f(d)]:k in [1..100]];
CROSSREFS
KEYWORD
nonn
AUTHOR
Marius A. Burtea, Jul 12 2024
STATUS
approved