OFFSET
1,1
COMMENTS
Squares pertaining to A083481.
a(n) == (p*q*r... ) where p,q,r are prime factors of n(n+1).
LINKS
Ivan Neretin, Table of n, a(n) for n = 1..10000
FORMULA
MATHEMATICA
Table[Times @@ ((a = Transpose[FactorInteger[n (n + 1)]])[[1]]^Quotient[a[[2]] + 1, 2]), {n, 54}] (* Ivan Neretin, May 20 2015 *)
PROG
(PARI) a(n)=sqrt(n*(n+1)*core(n*(n+1)))
(Python)
from math import prod
from sympy import factorint
def A083482(n): return n*(n+1)//prod(p**(q>>1) for p, q in factorint(n*(n+1)).items()) # Chai Wah Wu, Mar 20 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy and Meenakshi Srikanth (menakan_s(AT)yahoo.com), May 03 2003
EXTENSIONS
More terms from Benoit Cloitre, May 04 2003
More terms from David Wasserman, Nov 16 2004
STATUS
approved