OFFSET
1,2
COMMENTS
This sequence represents a class of sequences generated by rules of the form "a(1)=1, and if x is in a then hx+i and jx+k are in a, where h,i,j,k are integers." If m>1, at least one of the numbers b(m)=(a(m)-i)/h and c(m)=(a(m)-k)/j is in the set N of natural numbers. Let d(n) be the n-th b(m) in N, and let e(n) be the n-th c(m) in N. Note that a is a subsequence of both d and e. Examples:
...
...
...
...
...
Representative divisibility properties:
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
David Garth and Adam Gouge, Affinely Self-Generating Sets and Morphisms, Journal of Integer Sequences, Article 07.1.5, 10 (2007) 1-13.
FORMULA
a(1)=1, and if x is in a then 3x-2 and 4x-2 are in a; the terms of a are listed in without repetitions, in increasing order.
EXAMPLE
1 -> 2 -> 4,6 -> 10,14,16,22 ->
MAPLE
N:= 2000: # to get all terms <= N
S:= {}: agenda:= {1}:
while nops(agenda) > 0 do
S:= S union agenda;
agenda:= select(`<=`, map(t -> (3*t-2, 4*t-2), agenda) minus S, N)
od:
sort(convert(S, list)); # Robert Israel, Dec 22 2015
MATHEMATICA
PROG
(Haskell)
import Data.Set (singleton, deleteFindMin, insert)
a191113 n = a191113_list !! (n-1)
a191113_list = 1 : f (singleton 2)
where f s = m : (f $ insert (3*m-2) $ insert (4*m-2) s')
where (m, s') = deleteFindMin s
-- Reinhard Zumkeller, Jun 01 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, May 27 2011
STATUS
approved