OFFSET
1,1
LINKS
Charlie Neder, Table of n, a(n) for n = 1..1000
PROG
(Python)
from math import gcd
for n in range(1, 1001):
a_n = 1
res = []
for addend in range(1, n//2+1): res.append(addend*(n - addend))
for dividend in range(1, n+1):
if not n%dividend: res.append(dividend + n//dividend)
if dividend*dividend >= n: break
for i in res: a_n *= i // gcd(a_n, i)
print(n, a_n)
# Charlie Neder, Oct 09 2018
CROSSREFS
KEYWORD
nonn,nice
AUTHOR
Andrew Weimholt, Jul 05 2010
STATUS
approved
