OFFSET
1,4
COMMENTS
If p is a prime then a(p) = A001142(p)/(p^(p-1)).
EXAMPLE
a(5) = 1*5*10*10*5*1/625= 4.
MAPLE
A001142 := proc(n) local k ; mul(k^(2*k-1-n), k=1..n) ; end: A109873 := proc(n) local a; a := A001142(n) ; while a mod n = 0 and a > 1 do a := a/n ; od; RETURN(a) ; end: seq(A109873(n), n=1..20) ; # R. J. Mathar, Aug 15 2007
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Jul 10 2005
EXTENSIONS
More terms from R. J. Mathar, Aug 15 2007
STATUS
approved