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