OFFSET
1,1
FORMULA
a(n) = denominator(n*sum(1/prime(k), k=1..n)). [Gary Detlefs, Dec 06 2011]
MAPLE
seq(lcm(product(ithprime(i), i=1..j), j)/j, j=1..40);
PROG
(PARI) a(n) = denominator(n*sum(k=1, n, 1/prime(k))); \\ Michel Marcus, May 26 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Oct 17 2002
EXTENSIONS
Corrected and extended by Sascha Kurz, Jan 03 2003
STATUS
approved