OFFSET
0,3
LINKS
Robert Israel, Table of n, a(n) for n = 0..350
EXAMPLE
b(4) = b'(9/3) + b'(8/3) + b'(5/3) + b'(0/3) = b(3) + ((1/3)*b(2) + (2/3)*b(3)) + ((1/3)*b(1) + (2/3)*b(2)) + b(0) = 65/6.
MAPLE
bp:= proc(x) local t; option remember; if x::integer then b(x)
else t:= frac(x); t*b(ceil(x))+(1-t)*b(floor(x))
fi
end proc:
b:= proc(m) local k; option remember; add(bp(((m-1)^2-k^2)/(m-1)), k=0..m-1); end proc:
b(0):= 1: b(1):= 1:
map(numer@b, [$0..20]); # Robert Israel, May 20 2019
MATHEMATICA
bp[x_] := bp[x] = Module[{t}, If[IntegerQ[x], b[x],
t = FractionalPart[x]; t*b[Ceiling[x]] + (1-t)*b[Floor[x]]]];
b[m_] := b[m] = Sum[bp[((m-1)^2 - k^2)/(m-1)], {k, 0, m-1}];
b[0] = 1; b[1] = 1;
Map[Numerator@b@#&, Range[0, 20]] (* Jean-François Alcover, Jan 12 2024, after Robert Israel *)
CROSSREFS
KEYWORD
nonn,easy,frac
AUTHOR
Michele Dondi (blazar(AT)lcm.mi.infn.it), Jun 10 2002
EXTENSIONS
More terms from David Wasserman
Edited by N. J. A. Sloane, May 20 2019
STATUS
approved