Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).
%I #39 Jan 08 2021 20:51:46
%S 1,2,4,8,16,3,6,12,24,48,96,192,384,5,10,20,40,80,160,320,640,1280,
%T 2560,7,14,28,56,112,224,448,896,1792,9,18,36,72,144,288,576,1152,
%U 2304,4608,9216,11,22,44,88,176,352,704,1408,2816,5632,11264,13,26,52,104,208,416
%N Let S = smallest missing positive number, adjoin S, 2*S, 4*S, 8*S, 16*S, ... to the sequence until reaching a term that has S as a substring; reset S to the smallest missing positive number, repeat.
%C Theorem 1: Every positive numbers appears at least once.
%C Proof from _Keith F. Lynch_, Jan 04 2020:
%C Since no nonzero power of 2 equals a power of 10, i.e., log(10)/log(2) is irrational, any sequence in which each term is the double of the previous term will start with every decimal number infinitely many times. So any S will terminate after a finite number of steps, and the next missing number will be used as S. QED
%C Theorem 2: No term is repeated.
%C Proof:
%C Suppose N is repeated, so there are a pair of chains
%C {S, 2*S, 4*S, ..., N = 2^i*S, ...},
%C {T, 2*T, 4*T, ..., N = 2^j*T, ...},
%C where T occurs after S. There are two cases. If i>=j then T = 2^(i-j)*S, so T was not a missing number. If i<j then S = 2^(j-i)*T, and T would have been a smaller choice for S. QED
%C So this is a permutation of the positive integers.
%C From _Rémy Sigrist_, Jan 23 2020: (Start)
%C The sequence can naturally be seen as an irregular table where:
%C - the n-th row has A331442(n) = 1 + A331619(T(n, 1)) terms, and
%C - T(n, k+1) = 2*T(n, k) for k = 1..A331442(n)-1.
%C (End)
%D Eric Angelini, Posting to Math Fun Mailing List, Jan 04 2020.
%H Rémy Sigrist, <a href="/A331440/b331440.txt">Table of n, a(n) for n = 1..10103</a> (first 168 chains)
%H Rémy Sigrist, <a href="/A331440/a331440.gp.txt">PARI program for A331440</a>
%e The process begins like this:
%e Initially S = 1 is the smallest missing number, so we have:
%e S = 1, 2, 4, 8, 16, stop (because 16 contains S), S = 3, 6, 12, 24, 48, 96, 192, 384, stop, S = 5, 10, 20, 40, 80, 60, 320, 640, 1280, 2560, stop, S = 7, 14, 28, 56, 112, 224, 448, 896, 1792, stop, S = 9, 18, 36, 72, ...
%o (PARI) See Links section.
%Y The inverse permutation is A331441. The lengths of the chains are given in A331442.
%Y Cf. A331619.
%K nonn,base,tabf
%O 1,2
%A _N. J. A. Sloane_, Jan 21 2020