A354223 The initial sequence of the positive integers gradually becomes interspersed with the copies of its terms and so is transformed into this sequence. The method is described in the Comments section.
1, 2, 3, 1, 3, 4, 2, 1, 2, 4, 1, 4, 5, 3, 2, 1, 2, 3, 1, 3, 5, 2, 1, 2, 5, 1, 5, 6, 4, 3, 2, 1, 2, 3, 1, 3, 4, 2, 1, 2, 4, 1, 4, 6, 3, 2, 1, 2, 3, 1, 3, 6, 2, 1, 2, 6, 1, 6, 7, 5, 4, 3, 2, 1, 2, 3, 1, 3, 4, 2, 1, 2, 4, 1, 4, 5, 3, 2, 1, 2, 3, 1, 3, 5, 2, 1, 2, 5, 1, 5, 7, 4, 3, 2, 1, 2, 3, 1, 3, 4, 2, 1, 2, 4, 1, 4, 7, 3, 2, 1, 2, 3, 1, 3, 7, 2, 1, 2, 7, 1, 7
Offset: 1
Keywords
Examples
The initial line of the stamps: (1) (2) (3) (4) (5) ... Stamp (1) leaps over two others as is governed by stamp (2), leaving its imprint 1 in front: 1 (2) (3) (1) (4) (5) ... Stamp (2) leaps over three others as stamp (3) dictates. Its imprint 2 remains on the paper: 1 2 (3) (1) (4) (2) (5) ... Stamp (3) leaps over only the stamp (1), leaving its imprint 3 behind: 1 2 3 (1) (3) (4) (2) (5) ... Now stamp (1) leaps again, over three others into the gap between stamps (2) and (5). Its print 1 remains: 1 2 3 1 (3) (4) (2) (1) (5) ... (And so on.)
Links
- Thomas Scheuerle, Table of n, a(n) for n = 1..6000
- Thomas Scheuerle, First differences of a(n) as scatter plot from 1 to 10000.
Programs
-
MATLAB
function a = A354223( max_n ) a = [1:max_n]; for n = 1:max_n m = a(n); j = a(n+1); a = [a(1:n+j) m a(n+j+1:end)]; end a = a(1:max_n); end % Thomas Scheuerle, May 20 2022
Formula
From Thomas Scheuerle, May 19 2022: (Start)
a(2^(m+1) + n - 1) - a(n) < 2, if n < 2^(m+1) - m. This evaluates to zero for all n with m+1 exceptions. For these exceptions it will become 1. The exceptions are n = {A132045(m), A132045(m)+2^m, A132045(m)+2^m+2^(m-1), A132045(m)+2^m+2^(m-1)+2^(m-2), ..., A132045(m)+2^(m+1)-1}.
(1/n)*Sum_{k=1..n} a(k) <= 3. It appears that this arithmetic mean converges to 3 if n approaches infinity.
a(Sum_{k=1..n} (1 + A195986(k))) = 2. (End)
Comments