login
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.
7

%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