A138095 Limiting sequence when we start with positive integers (A000027) and at step n >= 1 add: to the term at position n + 1*a(n) the value 1, to the term at position n + 2*a(n) the value 1, ..., to the term at position n + n*a(n) the value 1.
1, 3, 3, 4, 6, 7, 7, 10, 10, 10, 12, 14, 14, 15, 15, 17, 18, 19, 20, 23, 22, 22, 25, 24, 25, 27, 29, 30, 32, 32, 31, 32, 34, 35, 39, 36, 38, 39, 41, 42, 43, 43, 45, 46, 46, 46, 48, 52, 51, 53, 51, 52, 55, 55, 56, 59, 57, 60, 63, 62, 62, 64, 63, 65, 66, 68, 69, 70, 72, 71
Offset: 1
Examples
n | 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
--------+-----------------------------------------------------------
Step 1 | 3
Step 2 | 6 9
Step 3 | 7 10 13
Step 4 | 10 14 17 21
Step 5 | 12 18
Step 6 | 14 22
Step 7 | 15
Step 8 | 19
Step 9 | 20
Step 10 | 23
--------+-----------------------------------------------------------
a(n) | 1 3 3 4 6 7 7 10 10 10 12 14 14 15 15 17 18 19 20 23
Crossrefs
Programs
-
PARI
lista(nn) = my(va = [1..nn]); for (n=1, nn, for (k=1, n, my(j = n+k*va[n]); if (j <= #va, va[j]++); )); va; \\ Michel Marcus, Aug 09 2022
Extensions
Corrected and extended by Michel Marcus, Aug 09 2022