OFFSET
1,1
PROG
(PARI) a(n) = my(phi = (1+sqrt(5))/2); sum(k=1, n, floor(k*phi)*floor((n-k+1)*phi^2)); \\ Michel Marcus, Jun 10 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Missing a(13) = 1735 inserted by Sean A. Irvine, Jun 09 2019
STATUS
approved