OFFSET
2,8
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 2..10000
MATHEMATICA
stol[n_] := stol[n] = If[n == 1, {}, If[n != Round[Round[n/GoldenRatio]*GoldenRatio], Join[stol[Floor[n/GoldenRatio^2] + 1], {0}], Join[stol[Round[n/GoldenRatio]], {1}]]];
a[n_] := IntegerExponent[FromDigits[stol[n]], 10]; Array[a, 100, 2]
PROG
(PARI) stol(n) = {my(phi=quadgen(5)); if(n==1, [], if(n != round(round(n/phi)*phi), concat(stol(floor(n/phi^2) + 1), [0]), concat(stol(round(n/phi)), [1]))); }
a(n) = valuation(fromdigits(stol(n)), 10);
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Amiram Eldar, Jul 07 2023
STATUS
approved