This site is supported by donations to The OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A299788 a(n) = denominator of Product_{d|n} (sigma(d)/d) where sigma(k) = the sum of the divisors of k (A000203). 1
 1, 2, 3, 8, 5, 1, 7, 64, 27, 25, 11, 3, 13, 49, 25, 1024, 17, 27, 19, 1000, 441, 121, 23, 16, 125, 169, 729, 7, 29, 625, 31, 32768, 121, 289, 1225, 5832, 37, 361, 1521, 6400, 41, 2401, 43, 1331, 3375, 529, 47, 3072, 343, 31250, 289, 17576, 53, 6561, 3025, 49 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,2 LINKS FORMULA a(n) = 1 for n = 1, 6, ... (no other n <= 5*10^6). a(n) = n for n = primes (A000040). a(n) = n * sqrt(n) for n in A280076 (union of 1 and squares of primes (A001248)). EXAMPLE For n=4; Product_{d|4} (sigma(d)/d) = (sigma(1)/1)*(sigma(2)/2)*(sigma(4)/4) = (1/1)*(3/2)*(7/4) = 21/8; a(4) = 8. MATHEMATICA Array[Denominator@ Product[DivisorSigma[1, d]/d, {d, Divisors@ #}] &, 51] PROG (MAGMA) [Denominator(&*[&+[c: c in Divisors(d)] / d: d in Divisors(n)]): n in [1..100]] (PARI) a(n) = my(p=1); fordiv(n, d, p *= sigma(d)/d); denominator(p); \\ Michel Marcus, Jan 21 2019 CROSSREFS Cf. A000203, A318491, A318492, A322673 (numerator). Sequence in context: A229870 A202651 A084110 * A244915 A244668 A192646 Adjacent sequences:  A299785 A299786 A299787 * A299789 A299790 A299791 KEYWORD nonn,frac AUTHOR Jaroslav Krizek, Jan 21 2019 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified May 21 09:15 EDT 2019. Contains 323441 sequences. (Running on oeis4.)