OFFSET
1,1
LINKS
Harry J. Smith, Table of n, a(n) for n = 1..1000
FORMULA
a(n) = 2*A053813(n). - Amiram Eldar, Jul 07 2022
MATHEMATICA
Do[ d = Drop[ Divisors[ n ], -1 ]; l = Length[ d ]; ev = 0; k = 2; While[ k <= l, If[ EvenQ[ d[ [ k ] ] ], ev = ev + d[ [ k ] ] ]; k++ ]; If[ IntegerQ[ n/ev ], Print[ n ] ], {n, 2, 10^3, 2} ]
Select[Range[550], Divisible[#, Total[Select[Most[Divisors[#]], EvenQ]]]&]//Quiet (* Harvey P. Dale, Mar 11 2011 *)
PROG
(PARI) { n=0; for (m=1, 10^9, d=divisors(m); s=0; for (i=2, numdiv(m) - 1, if (d[i]%2 == 0, s += d[i])); if (s > 0 && m%s == 0, write("b066190.txt", n++, " ", m); if (n==1000, return)) ) } \\ Harry J. Smith, Feb 05 2010
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Dec 15 2001
EXTENSIONS
a(56) from Harry J. Smith, Feb 05 2010
STATUS
approved