cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A087788 3-Carmichael numbers: Carmichael numbers equal to the product of 3 primes: k = p*q*r, where p < q < r are primes such that a^(k-1) == 1 (mod k) if a is prime to k.

Original entry on oeis.org

561, 1105, 1729, 2465, 2821, 6601, 8911, 10585, 15841, 29341, 46657, 52633, 115921, 162401, 252601, 294409, 314821, 334153, 399001, 410041, 488881, 512461, 530881, 1024651, 1152271, 1193221, 1461241, 1615681, 1857241, 1909001, 2508013
Offset: 1

Views

Author

Miklos Kristof, Oct 07 2003

Keywords

Comments

It is interesting that most of the numbers have the last digit 1. For example 530881, 3581761, 7207201, etc.
Granville & Pomerance conjecture that there are ~ c x^(1/3)/(log x)^3 terms of this sequence up to x. Heath-Brown proves that, for any e > 0, there are O(x^(7/20 + e)) terms of this sequence up to x. - Charles R Greathouse IV, Nov 19 2012

Examples

			a(6)=6601=7*23*41: 7-1|6601-1, 23-1|6601-1, 41-1|6601-1, i.e., 6|6600, 22|6600, 40|6600.
		

References

  • O. Ore, Number Theory and Its History, McGraw-Hill, 1948, Reprinted by Dover Publications, 1988, Chapter 14.

Crossrefs

Intersection of A002997 and A007304.
Cf. A162290.

Programs

  • PARI
    list(lim)=my(v=List());forprime(p=3,(lim)^(1/3), forprime(q=p+1, sqrt(lim\p),forprime(r=q+1,lim\(p*q),if((q*r-1)%(p-1)||(p*r-1)%(q-1)||(p*q-1)%(r-1),,listput(v,p*q*r)))));vecsort(Vec(v)) \\ Charles R Greathouse IV, Nov 19 2012

Formula

k is composite and squarefree and for p prime, p|k => p-1|k-1. A composite odd number k is a Carmichael number if and only if k is squarefree and p-1 divides k-1 for every prime p dividing k (Korselt, 1899) k = p*q*r, p-1|k-1, q-1|k-1, r-1|k-1.

Extensions

Minor edit to definition by N. J. A. Sloane, Sep 14 2009