OFFSET
1,4
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..395
FORMULA
a(n) = Product_{k=0..n-1} gcd(n,k!) = n^n/A058067(n).
a(p) = 1 for primes p. - Joerg Arndt, Apr 01 2014
a(p^2)= p^(p*(2*p-3)) for primes p. - R. J. Mathar, Apr 02 2014
MAPLE
MATHEMATICA
Table[Times@@Table[GCD[i, n!], {n, 0, i-1}], {i, 30}] (* Harvey P. Dale, Feb 23 2015 *)
PROG
(PARI) a(n) = vecprod(vector(n, k, gcd(n, (k-1)!))); \\ Michel Marcus, Sep 29 2020
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
R. J. Mathar, Apr 01 2014
EXTENSIONS
More terms from Harvey P. Dale, Feb 23 2015
a(20) corrected by Amiram Eldar, Sep 29 2020
STATUS
approved