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
Keywords
Links
- Clark Kimberling, Table of n, a(n) for n = 1..1001
Programs
-
Mathematica
Table[n + Sum[Floor[(n - k)/Pi], {k, 1, n}], {n, 1, 100}] -
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
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
Extensions
Updated by Clark Kimberling, Mar 27 2017
Offset changed to 1 by Sean A. Irvine, May 21 2019
Comments