%I #12 Aug 24 2024 10:35:10
%S 6,12,18,24,36,48,54,72,108,144,162,192,216,324,432,576,648,768,1296,
%T 1458,1536,1728,2304,2916,4608,5184,5832,6912,11664,12288,13122,13824,
%U 20736,26244,36864,39366,41472,46656,52488,78732,104976,110592,157464,186624,196608,314928,331776,373248,419904,589824,786432
%N 3-smooth numbers 2^i*3^j where i and j are regular 3-smooth numbers.
%e 54 is in the list because it is 2^1*3^3 and 1 and 3 are 3-smooth numbers.
%p interface(prettyprint=0) ; isA003586 := proc(n) local ifa,f,bas ; ifa := op(2,ifactors(n)) ; for f from 1 to nops(ifa) do bas := op(1,op(f,ifa)) ; if bas > 3 then RETURN(false) ; fi ; od ; RETURN(true) ; end: isA120942 := proc(n) local ifa,f,bas,expo ; ifa := op(2,ifactors(n)) ; if nops(ifa) <> 2 then RETURN(false) ; fi ; for f from 1 to nops(ifa) do bas := op(1,op(f,ifa)) ; if bas > 3 then RETURN(false) ; else expo := op(2,op(f,ifa)) ; if not isA003586(expo) then RETURN(false) ; fi ; fi ; od ; RETURN(true) ; end: for n from 2 to 40000 do if isA120942(n) then printf("%d,",n) ; fi ; od; # _R. J. Mathar_, Aug 23 2006
%t maxTerm = 10^6;
%t maxExponent = Log2[maxTerm] // Floor;
%t ee = Select[Range[maxExponent], FactorInteger[#][[-1, 1]] <= 3&];
%t ff = Flatten[Table[{{2, i}, {3, j}}, {i, ee}, {j, ee}], 1];
%t okQ[n_] := MemberQ[ff, FactorInteger[n]];
%t Select[Range[maxTerm], okQ] (* _Jean-François Alcover_, Aug 24 2024 *)
%Y Cf. A003586.
%K nonn
%O 1,1
%A _Ben Paul Thurston_, Aug 18 2006
%E Corrected and extended by _Ray Chandler_ and _R. J. Mathar_, Aug 23 2006
%E More terms from _Jean-François Alcover_, Aug 24 2024