OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..136 from R. J. Mathar)
MAPLE
n :=1 :
for k from 1 to 50000 do
if modp (k, A184389(k)) = 0 then
printf("%d %d\n", n, k) ;
n := n+1 ;
end if;
end do: # R. J. Mathar, Oct 04 2014
MATHEMATICA
t[n_] := n*(n + 1)/2; Select[Range[10^4], Divisible[#, t[DivisorSigma[0, #]]] &] (* Amiram Eldar, Jan 17 2021 *)
dsiQ[n_]:=With[{d=DivisorSigma[0, n]}, IntegerQ[n/((d(d+1))/2)]]; Select[Range[10000], dsiQ] (* Harvey P. Dale, Aug 20 2023 *)
PROG
(PARI) lista(nn) = {for (n=1, nn, if (2*n % (numdiv(n)*(numdiv(n)+1)) == 0, print1(n, ", ")); ); } \\ Michel Marcus, Jun 02 2013
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Ctibor O. Zizka, May 30 2009
EXTENSIONS
Corrected by Michel Marcus, Jun 02 2013
STATUS
approved