OFFSET
1,2
LINKS
Colin Barker, Table of n, a(n) for n = 1..1000
EXAMPLE
a(3) = 90 because the first 3 hexagonal numbers are [1,6,15], and 3/(1/1+1/6+1/15) = 90/37.
PROG
(PARI)
harmonicmean(v) = #v / sum(k=1, #v, 1/v[k])
s=vector(30); for(n=1, #s, s[n]=numerator(harmonicmean(vector(n, k, k*(2*k-1))))); s
CROSSREFS
KEYWORD
nonn
AUTHOR
Colin Barker, Nov 19 2014
STATUS
approved