OFFSET
1,2
LINKS
Colin Barker, Table of n, a(n) for n = 1..1000
EXAMPLE
a(3) = 504 because the first 3 octagonal numbers are [1,8,21], and 3/(1/1+1/8+1/21) = 504/197.
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, 3*k^2-2*k)))); s
CROSSREFS
KEYWORD
nonn
AUTHOR
Colin Barker, Nov 21 2014
STATUS
approved