A022796 Place where n-th 1 occurs in A023134.

1, 2, 3, 4, 6, 8, 10, 13, 16, 19, 23, 27, 31, 36, 41, 46, 52, 58, 64, 71, 78, 85, 93, 101, 109, 117, 126, 135, 144, 154, 164, 174, 185, 196, 207, 219, 231, 243, 256, 269, 282, 296, 310, 324, 339, 354, 369, 384, 400, 416, 432, 449, 466, 483, 501, 519

Offset

1,2

Comments

This is column 1 of the transposable interspersion A283944.

Links

Clark Kimberling, Table of n, a(n) for n = 1..1001

Formula

a(n) = n + Sum_({k=1..n} floor((n-k)/Pi.

a(n) <= floor(binomial(n, 2)/Pi) + n. - David A. Corneth, Mar 26 2017

Mathematica

Table[n + Sum[Floor[(n - k)/Pi], {k, 1, n}], {n, 1, 100}]

Prog

(PARI) for(n=1, 100, print1(n + sum(k=1, n, floor((n - k)/Pi)), ", ")); \\ Indranil Ghosh, Mar 26 2017
(Python)
import math
from mpmath import *
mp.dps = 100
print([n + sum([int(math.floor((n - k)/pi)) for k in range(1, n + 1)]) for n in range(1, 101)]) # Indranil Ghosh, May 28 2019

Crossrefs

Cf. A000796, A283944.

Sequence in context: A280706 A025695 A025694 * A008748 A089649 A049700

Adjacent sequences: A022793 A022794 A022795 * A022797 A022798 A022799

Keyword

nonn

Author

Clark Kimberling

Extensions

Updated by Clark Kimberling, Mar 27 2017

Offset changed to 1 by Sean A. Irvine, May 21 2019

Status

approved