OFFSET
1,4
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..5000
MATHEMATICA
b[j_]:= b[j]= Sum[KroneckerDelta[j, Binomial[m+2, 3]], {m, 0, 15}];
Table[A025120[n], {n, 100}] (* G. C. Greubel, Sep 14 2022 *)
PROG
(Magma)
A023533:= func< n | Binomial(Floor((6*n-1)^(1/3)) +2, 3) ne n select 0 else 1 >;
[A025120(n): n in [1..100]]; // G. C. Greubel, Sep 14 2022
(SageMath)
@CachedFunction
def b(j): return sum(bool(j==binomial(m+2, 3)) for m in (0..13))
@CachedFunction
def A025120(n): return sum(floor((n-j+2)*golden_ratio)*b(j) for j in (((n+4)//2)..n+1))
[A025120(n) for n in (1..100)] # G. C. Greubel, Sep 14 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved