%I #11 Jan 29 2019 09:24:25
%S 18,36,48,50,72,98,100,144,196,200,242,338,375,392,400,484,578,676,
%T 722,784,800,968,1029,1058,1125,1156,1280,1352,1444,1568,1600,1682,
%U 1922,1936,2116,2312,2704,2738,2888,3087,3136,3200,3362,3364,3375,3645,3698,3844
%N Composite solutions to S(n) = rad(n), where S() is A002034 and rad() is A007947.
%t kemp[n_] := Module[{m = 1}, While[!IntegerQ[m!/n], m++]; m]; rad[n_] := Times @@ (First@# & /@ FactorInteger@ n); aQ[n_] := CompositeQ[n] && kemp[n] == rad[n]; Select[Range[10000], aQ] (* _Amiram Eldar_, Jan 28 2019 *)
%Y Cf. A002034, A007947.
%K nonn
%O 1,1
%A _Jason Earls_, Aug 22 2002
%E More terms from _Amiram Eldar_, Jan 28 2019
|