A022796 - OEIS

%I #24 May 15 2021 06:18:34

%S 1,2,3,4,6,8,10,13,16,19,23,27,31,36,41,46,52,58,64,71,78,85,93,101,

%T 109,117,126,135,144,154,164,174,185,196,207,219,231,243,256,269,282,

%U 296,310,324,339,354,369,384,400,416,432,449,466,483,501,519

%N Place where n-th 1 occurs in A023134.

%C This is column 1 of the transposable interspersion A283944.

%H Clark Kimberling, <a href="/A022796/b022796.txt">Table of n, a(n) for n = 1..1001</a>

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

%F a(n) <= floor(binomial(n, 2)/Pi) + n. - _David A. Corneth_, Mar 26 2017

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

%o (PARI) for(n=1, 100, print1(n + sum(k=1, n, floor((n - k)/Pi)),", ")); \\ _Indranil Ghosh_, Mar 26 2017

%o (Python)

%o import math

%o from mpmath import *

%o mp.dps = 100

%o 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

%Y Cf. A000796, A283944.

%K nonn

%O 1,2

%A _Clark Kimberling_

%E Updated by _Clark Kimberling_, Mar 27 2017

%E Offset changed to 1 by _Sean A. Irvine_, May 21 2019