OFFSET
1,2
EXAMPLE
sigma(1) = 1 = 1*sigma(1).
sigma(10) = 18 = 18*(1*sigma(1)) = 3*(2*sigma(2)).
sigma(6) = 12 = 12*(1*sigma(1)) = 2*(2*sigma(2)) = 1*(3*sigma(3)).
sigma(30) = 72 = 72*(1*sigma(1)) = 11*(2*sigma(2)) = 6*(3*sigma(3)) = 1*(6*sigma(6)) .
MATHEMATICA
a[n_] := Module[{k = 1}, While[Count[Divisors[k], _?(Divisible[DivisorSigma[1, k], # * DivisorSigma[1, #]] &)] != n, k++]; k]; Array[a, 25] (* Amiram Eldar, May 12 2021 *)
PROG
(Magma) sd:=func<n, d| DivisorSigma(1, n) mod (d*DivisorSigma(1, d)) eq 0>; a:=[]; for n in [1..32] do k:=1; while #[d:d in Divisors(k)|sd(k, d)] ne n do k:=k+1; end while; Append(~a, k); end for; a;
(PARI) isok(k, n) = my(sk=sigma(k)); sumdiv(k, d, (sk % (d*sigma(d))) == 0) == n;
a(n) = my(k=1); while (!isok(k, n), k++); k; \\ Michel Marcus, May 12 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Marius A. Burtea, May 10 2021
STATUS
approved