A076928 a(1) = 1, a(n+1)= a(n)*(n+1) divided by the largest prime divisor of n+1.

1, 1, 1, 2, 2, 4, 4, 16, 48, 96, 96, 384, 384, 768, 2304, 18432, 18432, 110592, 110592, 442368, 1327104, 2654208, 2654208, 21233664, 106168320, 212336640, 1911029760, 7644119040, 7644119040, 45864714240, 45864714240, 733835427840

Offset

1,4

Links

Table of n, a(n) for n=1..32.

Maple

a[1] := 1:for n from 2 to 100 do q := ifactors(n)[2]:a[n] := n*a[n-1]/q[nops(q)][1]: od:seq(a[j], j=1..100);

Mathematica

nxt[{a_, b_}]:={a(b+1)/FactorInteger[b+1][[-1, 1]], b+1}; Transpose[ NestList[ nxt[#]&, {1, 1}, 40]][[1]] (* Harvey P. Dale, Aug 11 2011 *)

Crossrefs

Cf. A076929.

Sequence in context: A086117 A095240 A110428 * A286737 A025557 A285909

Adjacent sequences: A076925 A076926 A076927 * A076929 A076930 A076931

Keyword

nonn

Author

Amarnath Murthy, Oct 18 2002

Extensions

Corrected and extended by Sascha Kurz, Jan 21 2003

Status

approved