login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A109361 a(n) = Product_{k=1..n} sigma(k)/d(k), where sigma(k) = Sum_{j|k} j and d(k) = Sum_{j|k} 1. Set a(n) = 0 if the corresponding product is not an integer (e.g., for n=2 and n=10). 1
1, 0, 3, 7, 21, 63, 252, 945, 4095, 0, 110565, 515970, 3611790, 21670740, 130024440, 806151528, 7255363752, 47159864388, 471598643880, 3301190507160, 26409524057280, 237685716515520, 2852228598186240, 21391714486396800 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,3
COMMENTS
The product at n = 2 is the noninteger 1.5. The product at n = 10 is the noninteger 18427.5. Jack Brennen's observed that the only values which are not integers occur when n = 2 or 10, for n < 5000. Are all products for n >= 11 integers?
No other nonintegers found up to 200000. - Michel Marcus, Sep 14 2015
No other nonintegers up to 3000000. - Robert Israel, Jan 22 2018
LINKS
FORMULA
Product_{k=1..n} sigma(k)/d(k) = Product_{p=primes} Product_{k>=1} ((p^(k+1)-1)*k/((p^k -1)(k+1)))^floor(n/p^k).
a(n) = A066780(n)/A066843(n) if this is an integer, else 0. - Michel Marcus, Sep 14 2015
EXAMPLE
a(4) = 1 * 3 * 4 * 7 /(1 * 2 * 2 * 3) = 7.
MAPLE
p:= 1: A[1]:= 1:
for n from 2 to 50 do
p:= p * numtheory:-sigma(n)/numtheory:-tau(n);
if p::integer then A[n]:= p else A[n]:= 0 fi
od:
seq(A[n], n=1..50); # Robert Israel, Jan 22 2018
MATHEMATICA
Table[If[IntegerQ[Product[DivisorSigma[1, k]/Length[Divisors[k]], {k, 1, n}]], Product[DivisorSigma[1, k]/Length[Divisors[k]], {k, 1, n}], 0], {n, 1, 30}] (* Stefan Steinerberger, Oct 24 2007 *)
PROG
(PARI) a(n) = my(q = prod(k=1, n, sigma(k)/numdiv(k))); if (denominator(q)==1, q, 0); \\ Michel Marcus, Sep 14 2015
CROSSREFS
Sequence in context: A141495 A151412 A121797 * A052805 A148674 A148675
KEYWORD
nonn
AUTHOR
Leroy Quet, Aug 22 2005
EXTENSIONS
More terms from Stefan Steinerberger, Oct 24 2007
STATUS
approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified February 27 11:52 EST 2024. Contains 370375 sequences. (Running on oeis4.)