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