%I #14 Dec 06 2020 18:06:28
%S 4,6,9,10,12,14,15,21,22,25,26,33,34,35,36,38,39,46,49,51,55,56,57,58,
%T 62,65,69,74,77,82,85,86,87,91,93,94,95,106,108,111,115,118,119,121,
%U 122,123,129,133,134,141,142,143,145,146,155,158,159,161,166,169,177,178,183,185,187,194,201
%N Composite numbers k such that A339423(k) divides k.
%H Robert Israel, <a href="/A339424/b339424.txt">Table of n, a(n) for n = 1..10000</a>
%e a(5)=12 = 2*2*3 is a term because 2 + 2*2 = 6 divides 12.
%p filter:= proc(n)
%p local F,T,P,j;
%p if isprime(n) then return false fi;
%p F:= sort(map(t -> t[1]$t[2], ifactors(n)[2]));
%p T:= 0; P:= 1;
%p for j from 1 to nops(F)-1 do
%p P:= P*F[j];
%p T:= T+P;
%p od;
%p n mod T = 0
%p end proc:
%p select(filter, [$4..1000]);
%Y Cf. A002808, A339423.
%Y Union of A001358 and A339425.
%K nonn
%O 1,1
%A _J. M. Bergot_ and _Robert Israel_, Dec 03 2020
|