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!)
A045661 a(n) = Product_{d|n} (n/d + d). 1
2, 9, 16, 100, 36, 1225, 64, 2916, 600, 5929, 144, 529984, 196, 18225, 16384, 231200, 324, 3538161, 400, 5143824, 48400, 89401, 576, 1482250000, 6760, 164025, 112896, 26050816, 900, 5679280321, 1024, 50808384, 226576, 442225 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
a(n) is a square when n is not a square or of the form (2m^2)^2 = 4m^4.
If n is prime, then a(n) = (n+1)^2. - Wesley Ivan Hurt, Apr 19 2021
LINKS
FORMULA
a(A006881(n)) = (sopf(A006881(n)) * (A006881(n)+1) )^2. - Wesley Ivan Hurt, May 20 2013
MATHEMATICA
Times@@(#/Divisors[ # ]+Divisors[ # ])& /@ Range[ 48 ]
PROG
(Haskell)
a045661 n = product [n'+d | d <- [1..n], let (n', m) = divMod n d, m == 0]
-- Reinhard Zumkeller, Feb 02 2012, Jan 25 2012
(PARI) a(n)=my(t=1); fordiv(n, d, t*=n/d+d); t \\ Charles R Greathouse IV, Jan 25 2012
(PARI) A045661(n)=my(t=1+#n=divisors(n)); prod(i=1, (t-1)\2, n[i]+n[t-i])^2*if(bittest(t, 0), 1, 2*n[t\2]) \\ M. F. Hasler, Jan 25 2012
CROSSREFS
Sequence in context: A032211 A032136 A032040 * A078707 A085332 A204424
KEYWORD
nonn,nice,easy
AUTHOR
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 September 29 23:48 EDT 2023. Contains 365781 sequences. (Running on oeis4.)