The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A216436 For g=A201557(n), define a(n) as the prime p|g such that G(g/p) is maximum, where G(k) = sigma(k)/(k*log(log(k))). 2
 17, 17, 19, 19, 19, 19, 19, 23, 23, 31, 31, 31, 7, 7, 7, 7, 7, 37, 37, 41, 43, 43, 43, 43, 43, 43, 2, 43, 2, 43, 2, 47, 47, 47, 47, 47, 47, 2, 47, 2, 47, 2, 47, 2, 53, 61, 61, 61, 61, 61, 61, 5, 5, 5, 67, 71, 11, 11, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,1 COMMENTS The ratio G(A201557(n)/p)/G(A201557(n)) is defined by Nicolas as the Gronwall quotient (A201557 = proper GA1 numbers). LINKS Michel Marcus, Table of n, a(n) for n = 1..1000 G. Caveney, J.-L. Nicolas and J. Sondow, On SA, CA, and GA numbers, Ramanujan J., 29 (2012), 359-384. J.-L. Nicolas, Table of proper GA1 numbers up to 10^60, 2011. J.-L. Nicolas, Computation of GA1 numbers, 2011. FORMULA G(A201557(n)/a(n)) >= G(A201557(n)/q) if prime q|A201557(n). MAPLE # See link "Computation of GA1 numbers". PROG (PARI) findq(i) = {f = factor(i); maxqu = 0.0; qmax = 0; for(iq=1, length(f~), qq = f[iq, 1]; qu = g(i/qq)/g(i); if (qu > maxqu, maxqu = qu; qmax = qq; ) ); return (qmax); } \\ for i in A201557 CROSSREFS Cf. A000203, A201557, A197638. Sequence in context: A102423 A010856 A291515 * A060360 A081702 A345735 Adjacent sequences: A216433 A216434 A216435 * A216437 A216438 A216439 KEYWORD nonn AUTHOR Michel Marcus, Sep 10 2012 EXTENSIONS Definition simplified and formula supplied by Jonathan Sondow, Sep 11 2012 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified August 10 21:39 EDT 2024. Contains 375058 sequences. (Running on oeis4.)