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
KEYWORD
nonn
AUTHOR
EXTENSIONS
Updated by Clark Kimberling, Mar 27 2017
Offset changed to 1 by Sean A. Irvine, May 21 2019
STATUS
approved