OFFSET
1,2
COMMENTS
a(n) is divisible by n, for all n >= 1.
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..445
MATHEMATICA
f[n_] := f[n] = Sum[Binomial[n, k] GCD[n, k] f[k - 1], {k, n}]; f[0] = 1; Table[f@n/n, {n, 22}] (* Robert G. Wilson v, May 27 2009 *)
PROG
(SageMath)
@CachedFunction
def b(n): return 1 if (n==0) else sum(binomial(n, j)*gcd(n, j)*b(j-1) for j in range(1, n+1)) # b = A159666
def A159667(n): return b(n)/n
[A159667(n) for n in range(1, 40)] # G. C. Greubel, Sep 27 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Apr 19 2009
EXTENSIONS
More terms from Robert G. Wilson v, May 27 2009
STATUS
approved