login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

a(n) = A000720(A348907(n+1)).
2

%I #7 Nov 11 2021 20:05:33

%S 1,2,1,3,2,4,1,3,2,5,4,6,1,3,2,7,5,8,4,6,1,9,3,2,7,5,8,10,4,11,6,1,9,

%T 3,2,12,7,5,8,13,10,14,4,11,6,15,1,9,3,2,12,16,7,5,8,13,10,17,14,18,4,

%U 11,6,15,1,19,9,3,2,20,12,21,16,7,5,8,13,22,10

%N a(n) = A000720(A348907(n+1)).

%C Regarding this sequence as an irregular triangle T(m,j) where the rows m terminate with 1 exhibits row length A338237(m). In such rows m, we have a permutation of the range of natural numbers 1..A338237(m).

%C Records are the natural numbers.

%H Michael De Vlieger, <a href="/A349191/b349191.txt">Table of n, a(n) for n = 1..10237</a> (as an irregular triangle, rows 1 <= n <= 36, flattened)

%H Michael De Vlieger, <a href="/A349191/a349191.png">Log-log scatterplot of a(n)</a> for 1 <= n <= 11636, showing 36 rows if read as an irregular table.

%e Table showing a(n) for the first rows m of this sequence seen as an irregular triangle T(m,j). "New" numbers introduced for prime (n+1) are shown in parentheses:

%e m\j 1 2 3 4 5 6 7 8 9 10 11 A338237(m)

%e ------------------------------------------------------------

%e 1: (1) 1

%e 2: (2) 1 2

%e 3: (3) 2 (4) 1 4

%e 4: 3 2 (5) 4 (6) 1 6

%e 5: 3 2 (7) 5 (8) 4 6 1 8

%e 6: (9) 3 2 7 5 8 (10) 4 (11) 6 1 11

%e ... (End)

%t c = 0; 1 + Reap[Do[Set[a[i], If[PrimeQ[i], i; c++, a[i - c]] ]; Sow[a[i]], {i, 2, 2^24}] ][[-1, -1]]

%Y Cf. A000027, A000040, A000720, A338237, A348907, A349192.

%K nonn,easy

%O 1,2

%A _Michael De Vlieger_, Nov 09 2021