OFFSET
1,3
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..1000
Steven Finch, Cilleruelo's LCM Constants, 2013. [Cached copy, with permission of the author]
Eric Weisstein's MathWorld, Dirichlet's theorem
FORMULA
A(a) = (a/phi(a))*Sum_{j=1..a, gcd(j,a)=1} 1/j.
log(lcm_{k=0..n} a*k+b) ~ A(a)*n for gcd(a,b)=1.
EXAMPLE
Sequence A(a) begins 1, 2, 9/4, 8/3, 125/48, 18/5, 343/120, 352/105, 1863/560, ...
MATHEMATICA
A[a_] := (a/EulerPhi[a])*Sum[If[GCD[j, a] == 1, 1/j, 0], {j, 1, a}]; Array[A, 40] // Denominator
PROG
(PARI) a(n)={denominator(n*sum(j=1, n, if(gcd(j, n)==1, 1/j))/eulerphi(n))} \\ Andrew Howroyd, Mar 16 2018
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Jean-François Alcover, Oct 23 2014
STATUS
approved