OFFSET
1,6
COMMENTS
Previous name "Largest unitary square divisor of n!" was incorrect. See A374988 for the correct sequence with this name. - Amiram Eldar, Jul 26 2024
FORMULA
a(n) = A056623(n!). - Michel Marcus, Aug 16 2020
a(n) = A056627(n)^2. - Amiram Eldar, Jul 08 2024
MATHEMATICA
f[p_, 1] := 1; f[p_, e_] := If[EvenQ[e], p^e, p^(e-3)]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n!]; Array[a, 27] (* Amiram Eldar, Jul 08 2024 *)
PROG
(PARI)
A055229(n) = { my(c=core(n)); gcd(c, n/c); }; \\ Charles R Greathouse IV, Nov 20 2012
A008833(n) = n/core(n) \\ Michael B. Porter, Oct 17 2009
a(n) = A056623(n!); \\ Michel Marcus, Aug 16 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Aug 08 2000
EXTENSIONS
More terms from Michel Marcus, Aug 16 2020
Incorrect name replaced with a formula by Amiram Eldar, Jul 26 2024
STATUS
approved