|
%I #15 Jul 20 2019 18:29:08
%S 4,16,18,36,48,64,72,100,144,150,162,180,192,196,256,288,294,324,400,
%T 432,448,450,484,490,576,588,600,648,676,720,768,784,882,900,960,1024,
%U 1134,1152,1156,1176,1200,1210,1296,1350,1444,1458,1584,1600,1620,1728
%N A number m is included if k + m/k divides m for at least one divisor k of m.
%C All terms are even. All even perfect squares are included. If m is included, then 4m is also included.
%C Also the set of all numbers m of the form m = r*s*(r+s)^2*t^2 with positive integers r,s,t; if one additionally requires (r,s) = 1 and r <= s, then there is exactly one such representation for each m. [From _Hagen von Eitzen_, Jul 22 2009]
%H Amiram Eldar, <a href="/A139719/b139719.txt">Table of n, a(n) for n = 1..10000</a>
%e 72 is included because 6 is a divisor of 72 and (6 + 72/6) = 18 divides 72.
%t nkdQ[n_]:=Module[{divs=Most[Divisors[n]]},MemberQ[Divisible[n,#]&/@ (#+n/#&/@ divs),True]]; Select[Range[2000],nkdQ] (* _Harvey P. Dale_, Apr 24 2012 *)
%Y Cf. A139718.
%K nonn
%O 1,1
%A _Leroy Quet_, May 01 2008
%E Extended by _Ray Chandler_, Jul 01 2009
|
