|
%I #15 Jul 16 2020 02:38:31
%S 27,105,150,180,231,588,627,650,805,840,897,945,1008,1100,1122,1134,
%T 1581,1755,2079,2106,2625,2660,2800,2958,2964,2967,2970,3055,3125,
%U 3150,3432,3564,3750,3861,4000,4070,4185,4500,4543,4760,4800,5292,5304,5355
%N Composites n such that A001414(n) is odd and divides n.
%H Robert Israel, <a href="/A036346/b036346.txt">Table of n, a(n) for n = 1..10000</a>
%e 5445 = 3*3*5*11*11 -> sum = 3+3+5+11+11 = 33 (=odd) so 33 divides 5445 exactly.
%p filter:= proc(n) local F,s;
%p if isprime(n) then return false fi;
%p F:= ifactors(n)[2];
%p s:= add(t[1]*t[2],t=F);
%p s::odd and ( n mod s = 0)
%p end proc:
%p select(filter, [$1..10000]); # _Robert Israel_, Jul 15 2020
%t With[{s = Map[{#, Total[Times @@@ FactorInteger[#]]} &, Select[Range[4, 6000], CompositeQ]]}, Select[s, Mod[#1, #2] == 0 && OddQ[#2] & @@ # &][[All, 1]] ] (* _Michael De Vlieger_, Jul 15 2020 *)
%Y Cf. A001414, A046346, A036345.
%K nonn
%O 1,1
%A _Patrick De Geest_, Dec 15 1998
%E Name and offset corrected by _Robert Israel_, Jul 15 2020
|
