OFFSET
0,3
EXAMPLE
a(6) = lcm(b(0),b(1),b(2),b(3),b(4),b(5),b(6)) = lcm(1,1,2,3,8,5,16) = 240. 6!! = 6*4*2 = 48 divides 240.
MAPLE
A130480 := proc(nmax) local a, n, df, k, f, p, l, b; a := [1, 1] ; while nops(a)< nmax do n := nops(a) ; df := ifactors(doublefactorial(n))[2] ; l := ilcm(op(a)) ; k := 1 ; for f in df do p := op(1, f)^op(2, f) ; if l mod p <> 0 then k := k*p ; fi; od: a := [op(a), k] ; od; b := [] ; for i from 1 to nops(a) do b := [op(b), ilcm(op(1..i, a))] ; od: RETURN(b) ; end: A130480(30) ; # R. J. Mathar, Oct 16 2007
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, May 29 2007
EXTENSIONS
More terms from R. J. Mathar, Oct 16 2007
STATUS
approved