OFFSET
0,4
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..1000
MAPLE
g:= proc(n) g(n):= `if`(n=0, 1, ilcm(n, g(n-1))) end:
b:= proc(n, i) option remember; local p;
p:= `if`(i<1, 1, ithprime(i));
`if`(n=0 or i<1, 1, max(b(n, i-1),
seq(p^j*b(n-p^j, i-1), j=1..ilog[p](n))))
end:
a:= n->g(n)/b(n, `if`(n<8, 3,
numtheory[pi](ceil(1.328*isqrt(n*ilog(n)))))):
seq(a(n), n=0..40); # Alois P. Heinz, May 22 2013
MATHEMATICA
b[n_, i_] := b[n, i] = Module[{p}, p = If[i<1, 1, Prime[i]]; If[n==0 || i<1, 1, Max[b[n, i-1], Table[p^j*b[n-p^j, i-1], {j, 1, Log[p, n] // Floor }]]]]; a[0]=1; a[n_] := LCM @@ Range[n] / b[n, If[n<8, 3, PrimePi[ Ceiling[ 1.328*Sqrt[n*Log[n] // Floor]]]]]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Feb 27 2016, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, May 10 2013
STATUS
approved