OFFSET
1,1
COMMENTS
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..10000
H. Rich, Hofstadter-Conway $10,000 Sequence, Programming Archives, 2008.
Eric Weisstein's World of Mathematics, Hofstadter-Conway $10,000 Sequence.
FORMULA
a(n) = Sum of denominators of Egyptian fraction expansion of A004001(n) - n/2 .
EXAMPLE
MATHEMATICA
HC[n_]:= HC[n]= If[n<3, 1, HC[HC[n-1]] +HC[n-HC[n-1]]]; (*HC=A004001*)
f[x_]:= 4*x -3*Floor[x];
A143232[n_]:= f[HC[n] -n/2];
Table[A143232[n], {n, 100}] (* G. C. Greubel, Sep 10 2024 *)
PROG
(J) a004001 =: ((] +&:$: -) $:@:<:)@.(2&<) M.
a143232 =: (+ 3 * 1&|)@:(a004001 - -:)"0
(Magma)
A004001:=[n le 2 select 1 else Self(Self(n-1))+ Self(n-Self(n-1)):n in [1..125]];
f:= func< x | 4*x - 3*Floor(x) >;
[A143232(n): n in [1..100]]; // G. C. Greubel, Sep 10 2024
(SageMath)
@CachedFunction
def b(n): # b = A004001
if n<3: return 1
else: return b(b(n-1)) + b(n-b(n-1))
def f(x): return x + 3 * (x - floor(x))
def A143232(n): return f(b(n) - n/2)
[A143232(n) for n in range(1, 101)] # G. C. Greubel, Sep 10 2024
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Dan Bron (j (at) bron.us), Jul 31 2008
STATUS
approved