A135354 a(0)=1, a(n) = largest divisor of n! that is coprime to a(n-1).

1, 1, 2, 3, 8, 15, 16, 315, 128, 2835, 256, 155925, 1024, 6081075, 2048, 638512875, 32768, 10854718875, 65536, 1856156927625, 262144, 194896477400625, 524288, 49308808782358125, 4194304

Offset

0,3

Links

Table of n, a(n) for n=0..24.

Formula

a(2n) = the largest power of 2 that divides (2n)!. a(2n+1) = the largest odd divisor of (2n+1)! = (2n+1)!/a(2n).

Mathematica

a = {1}; For[n = 1, n < 25, n++, AppendTo[a, Select[Divisors[n! ], GCD[a[[ -1]], # ] == 1 &][[ -1]]]]; a (* Stefan Steinerberger, Dec 10 2007 *)
ldnf[{n_, a_}]:={n+1, Max[Select[Divisors[(n+1)!], CoprimeQ[#, a]&]]}; Transpose[ NestList[ldnf, {0, 1}, 30]][[2]] (* Harvey P. Dale, Jan 21 2016 *)

Crossrefs

Cf. A060818, A049606.

Sequence in context: A265694 A238778 A004731 * A240974 A356879 A293688

Adjacent sequences: A135351 A135352 A135353 * A135355 A135356 A135357

Keyword

nonn

Author

Leroy Quet, Dec 07 2007

Extensions

More terms from Stefan Steinerberger, Dec 10 2007

Status

approved