OFFSET
1,1
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
f:= proc(n) local k; uses numtheory;
for k from 1 to 4 do if sigma(k*n)>=2*k*n then return k*n fi od:
6*n
end proc:
map(f, [$1..100]); # Robert Israel, Feb 10 2019
MATHEMATICA
a[n_] := Do[If[DivisorSigma[1, k*n] >= 2*k*n, Return[k*n]], {k, {1, 2, 3, 4, 6}}];
Table[a[n], {n, 1, 100}] (* Jean-François Alcover, Feb 09 2023 *)
PROG
(PARI) a(n) = forstep(m=n, 6*n, n, if(sigma(m)>=2*m, return(m)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Jeppe Stig Nielsen, Feb 01 2015
STATUS
approved