OFFSET
1,1
COMMENTS
Inspired by A249270.
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1000
James Grime and Brady Haran, 2.920050977316, Numberphile video, Nov 26 2020.
MAPLE
b:= proc(n) option remember; `if`(n=1, Pi,
(f-> (frac(f)+1)*floor(f))(b(n-1)))
end:
a:= n-> floor(b(n)):
seq(a(n), n=1..50); # Alois P. Heinz, Dec 03 2020
MATHEMATICA
Block[{a = {Pi}, $MaxExtraPrecision = 10^3}, Do[AppendTo[a, (FractionalPart[#] + 1) Floor[#]] &@ a[[-1]], 41]; Floor /@ a] (* Michael De Vlieger, Dec 04 2020 *)
PROG
(NARS2000) {(⌊{(⌊⍵)×1+1|⍵}⍣⍵)○1x}¨0, ⍳100
(PARI) lista(nn) = {localprec(500); my(vx = vector(nn)); vx[1] = Pi; for (n=2, nn, vx[n] = (frac(vx[n-1])+1)*floor(vx[n-1]); ); apply(floor, vx); } \\ Michel Marcus, Dec 03 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Michael Turniansky, Dec 03 2020
STATUS
approved