OFFSET
1,1
LINKS
T. D. Noe, Table of n, a(n) for n = 1..1000
FORMULA
x is a term if and only if x = 3*p1*p2*...*pk with primes 2 <= p1 <= p2 <= ... <= pk and 3*p1*p2*...*pi >= p(i+1) for all i < k.
a(n) = 3 * A196149(n). - Reinhard Zumkeller, Sep 28 2011
The number of terms <= x is c*x/log(x) + O(x/(log(x))^2), where c = 0.68514..., and a(n) = C*n*log(n*log(n)) + O(n), where C = 1/c = 1.45954... This follows from the formula just above. - Andreas Weingartner, Jun 30 2021
EXAMPLE
63 is an element because 63 = 3*3*7 and 3 <= 3 and 7 <= 3*3.
MAPLE
N:= 1000: # get all terms <= N
S:= {3}:
New:= {3}:
while New <> {} do
x:= New[1];
New:= subsop(1=NULL, New);
R:= {seq(k*x, k=1..min(x, N/x))} minus S;
S:= S union R;
New:= New union R;
od:
sort(convert(S, list)); # Robert Israel, Aug 27 2015
MATHEMATICA
3 Select[Range[132], Max[#[[2]]/#[[1]] & /@ Partition[Divisors[#], 2, 1]] <= 3 &] (* Michael De Vlieger, Aug 27 2015, after Harvey P. Dale at A196149 *)
PROG
(Haskell)
import Data.List.Ordered (union)
a052287 n = a052287_list !! (n-1)
a052287_list = f [3] where
f (x:xs) = x : f (xs `union` map (x *) [2..x])
-- Reinhard Zumkeller, Jun 25 2015, Sep 28 2011
CROSSREFS
KEYWORD
easy,nice,nonn
AUTHOR
Giuseppe Melfi, Feb 08 2000
EXTENSIONS
More terms from Reinhard Zumkeller, Jun 22 2003
STATUS
approved