OFFSET
0,4
COMMENTS
LINKS
Clark Kimberling, Table of n, a(n) for n = 0..1000
MATHEMATICA
a[0] = 1; a[1] = 1; a[2] = 1;
a[n_] := a[n] = 2*a[n - 1] - a[n - 3] + Sum[a[Floor[n/k]], {k, 2, n}];
Table[a[n], {n, 0, 90}] (* A298406 *)
PROG
(Python)
from functools import lru_cache
@lru_cache(maxsize=None)
def A298406(n):
if n <= 2:
return 1
k1 = n//j
while k1 > 1:
j2 = n//k1 + 1
c += (j2-j)*A298406(k1)
j, k1 = j2, n//j2
return c+n-j+1 # Chai Wah Wu, Mar 31 2021
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Clark Kimberling, Feb 10 2018
STATUS
approved