OFFSET
0,3
LINKS
Andrew Howroyd, Table of n, a(n) for n = 0..10000
PROG
(PARI) \\ here b(n)=A053615(n)
b(n)={my(t=(sqrtint(4*n) + 1)\2); t-abs(t^2-n)}
seq(n)={my(f=vector(sqrtint(n)+1), a=vector(n+1)); for(i=0, n, my(k=f[1+b(i)]); a[1+i]=k; if(k<#f, f[1+k]++)); a} \\ Andrew Howroyd, Nov 12 2024
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Ralf Stephan, May 15 2005
STATUS
approved