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,2

Links

Table of n, a(n) for n=1..70.

Example

        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

Prog

(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

Crossrefs

Cf. A000027, A136119, A137319, A137417, A137418, A136259, A136272, A137292, A136120, A137832, A137838, A137894, A137898, A137901, A137902, A137903.

Sequence in context: A088651 A161019 A137901 * A023837 A240115 A196249

Adjacent sequences: A138092 A138093 A138094 * A138096 A138097 A138098

Keyword

easy,nonn

Author

Ctibor O. Zizka, May 03 2008

Extensions

Corrected and extended by Michel Marcus, Aug 09 2022

Status

approved