A197638 GA1 numbers: composite m with G(m) >= G(m/p) for all prime factors p of m, where G(k) = sigma(k)/(k*log(log(k))) and sigma(k) = sum of divisors of k.

4, 14, 22, 26, 34, 38, 46, 58, 62, 74, 82, 86, 94, 106, 118, 122, 134, 142, 146, 158, 166, 178, 194, 202, 206, 214, 218, 226, 254, 262, 274, 278, 298, 302, 314, 326, 334, 346, 358, 362, 382, 386, 394, 398, 422, 446, 454, 458, 466, 478, 482, 502, 514, 526

Offset

1,1

Comments

The members with exactly two prime divisors counted with multiplicity are 4 and 2*p, for primes p > 5. (See Section 5 of "Robin's theorem, primes, and a new elementary reformulation of the Riemann Hypothesis".)

The smallest member with more than two prime factors is 183783600. Such GA1 numbers are called "proper" - see A201557 and "Table of proper GA1 numbers up to 10^60".

The smallest odd member is 1058462574572984015114271643676625.

See "On SA, CA, and GA numbers".

Links

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from T. D. Noe)

G. Caveney, J.-L. Nicolas, and J. Sondow, Robin's theorem, primes, and a new elementary reformulation of the Riemann Hypothesis, Integers 11 (2011), article A33.

G. Caveney, J.-L. Nicolas and J. Sondow, On SA, CA, and GA numbers, Ramanujan J., 29 (2012), 359-384 and arXiv:1112.6010 [math.NT], 2011-2012.

J.-L. Nicolas, Computation of GA1 numbers, 2011.

J.-L. Nicolas, Table of proper GA1 numbers up to 10^60, 2011.

Example

4 is a member because G(4) > 0 > G(2) = G(4/2).

Maple

See "Computation of GA1 numbers".

Mathematica

g[k_] := g[k] = DivisorSigma[1, k]/(k*Log[Log[k]]); okQ[n_] := Module[{p = Transpose[FactorInteger[n]][[1]]}, i = 1; While[i <= Length[p] && g[n] >= g[n/p[[i]]], i++]; i > Length[p]]; Select[Range[2, 1000], ! PrimeQ[#] && okQ[#] &] (* T. D. Noe, Dec 03 2011 *)

Prog

(PARI) g(k) = sigma(k)/(k*log(log(k)));
isga1(k)=if (isprime(k), return (0)); gk = g(k); f = factor(k); for(i=1, length(f~), if (gk < g(k/f[i, 1]), return(0)); ); 1; \\ Michel Marcus, Sep 09 2012

Crossrefs

Cf. A000203, A197639, A201557, A201558.

Sequence in context: A267768 A185008 A165721 * A051806 A034051 A199984

Adjacent sequences: A197635 A197636 A197637 * A197639 A197640 A197641

Keyword

nonn

Author

Geoffrey Caveney, Jean-Louis Nicolas and Jonathan Sondow, Dec 02 2011

Status

approved