%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