OFFSET
0,2
COMMENTS
a(15) > 2.4*10^12, if it exists.
All the terms are Zeckendorf-Niven numbers (A328208).
EXAMPLE
n | The n iterations
--+---------------------------------------------
1 | 4 -> 2 = Fibonacci(3)
2 | 12 -> 4 -> 2
3 | 24 -> 12 -> 4 -> 2
4 | 180 -> 60 -> 30 -> 10 -> 5 = Fibonacci(5)
5 | 1056 -> 264 -> 66 -> 22 -> 11 -> 11/2
6 | 2592 -> 1296 -> 324 -> 108 -> 27 -> 9 -> 9/2
MATHEMATICA
zeck[n_] := Length[DeleteCases[NestWhileList[# - Fibonacci[Floor[Log[Sqrt[5]*# + 3/2]/Log[GoldenRatio]]] &, n, # > 1 &], 0]]; (* Alonso del Arte at A007895 *)
s[n_] := s[n] = Module[{z = zeck[n]}, If[z == 1, 0, If[!Divisible[n, z], 1, 1 + s[n/z]]]];
seq[len_] := Module[{v = Table[0, {len}], c = 0, k = 1, i}, While[c < len, i = s[k] + 1; If[v[[i]] == 0, c++; v[[i]] = k]; k++]; v]; seq[9]
PROG
(PARI) zeck(n) = if(n<4, n>0, my(k=2, s, t); while(fibonacci(k++)<=n, ); while(k && n, t=fibonacci(k); if(t<=n, n-=t; s++); k--); s) \\ Charles R Greathouse IV at A007895
s(n) = {my(z = zeck(n)); if(z == 1, 0, if(n % z, 1, 1 + s(n/z))); }
lista(len) = {my(v = vector(len), c = 0, k = 1, i); while(c < len, i = s(k) + 1; if(v[i] == 0, c++; v[i] = k); k++); v; }
CROSSREFS
KEYWORD
nonn,base,more,new
AUTHOR
Amiram Eldar, Oct 20 2024
STATUS
approved