%I #28 Nov 05 2023 10:27:09
%S 1,4,9,16,36,64,81,144,256,324,576,729,1024,1296,2304,2916,4096,5184,
%T 6561,9216,11664,16384,20736,26244,36864,46656,59049,65536,82944,
%U 104976,147456,186624,236196,262144,331776,419904,531441,589824,746496,944784
%N Numbers of the form 4^i * 9^j, with i, j >= 0.
%C Numbers of the form 2^(2*i) * 3^(2*j)) or 3-smooth squares: intersection of A003586 and A000290; A001221(a(n)) <= 2; A001222(a(n)) is even; A006530(a(n)) <= 3. - _Reinhard Zumkeller_, May 16 2015
%C Closed under multiplication. - _Klaus Purath_, Oct 06 2023
%H Reinhard Zumkeller, <a href="/A025620/b025620.txt">Table of n, a(n) for n = 1..10000</a>
%F Sum_{n>=1} 1/a(n) = (4*9)/((4-1)*(9-1)) = 3/2. - _Amiram Eldar_, Sep 24 2020
%F a(n) ~ exp(sqrt(8*log(2)*log(3)*n)) / 6. - _Vaclav Kotesovec_, Sep 24 2020
%t n = 10^6; Flatten[Table[4^i*9^j, {i, 0, Log[4, n]}, {j, 0, Log[9, n/4^i]}]] // Sort (* _Amiram Eldar_, Sep 24 2020 *)
%o (Haskell)
%o import Data.Set (singleton, deleteFindMin, insert)
%o a025620 n = a025620_list !! (n-1)
%o a025620_list = f $ singleton 1 where
%o f s = y : f (insert (4 * y) $ insert (9 * y) s')
%o where (y, s') = deleteFindMin s
%o -- _Reinhard Zumkeller_, May 16 2015
%o (PARI) list(lim)=my(v=List(), N); for(n=0, logint(lim\=1, 9), N=9^n; while(N<=lim, listput(v, N); N<<=2)); Set(v) \\ _Charles R Greathouse IV_, Jan 10 2018
%Y Cf. A003586, A000290, A001221, A001222, A006530, subsequence of A036667.
%K easy,nonn
%O 1,2
%A _David W. Wilson_
|