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.

Showing 1-2 of 2 results.

A063980 Pillai primes: primes p such that there exists an integer m such that m! + 1 == 0 (mod p) and p != 1 (mod m).

Original entry on oeis.org

23, 29, 59, 61, 67, 71, 79, 83, 109, 137, 139, 149, 193, 227, 233, 239, 251, 257, 269, 271, 277, 293, 307, 311, 317, 359, 379, 383, 389, 397, 401, 419, 431, 449, 461, 463, 467, 479, 499, 503, 521, 557, 563, 569, 571, 577, 593, 599, 601, 607
Offset: 1

Views

Author

R. K. Guy, Sep 08 2001

Keywords

Comments

Hardy & Subbarao prove that this sequence is infinite. An upper bound can be extracted from their proof: a(n) < e^e^...^e^O(n log n) with e appearing n times. This tetrational bound could be improved with results on the disjointness of the factorizations of numbers of the form k! + 1. - Charles R Greathouse IV, Sep 15 2015
Named after the Indian mathematician Subbayya Sivasankaranarayana Pillai (1901-1950). - Amiram Eldar, Jun 16 2021

Crossrefs

Smallest m is given in A063828, largest in A211411.

Programs

  • Mathematica
    ok[p_] := (r = False; Do[If[Mod[m! + 1, p] == 0 && Mod[p, m] != 1, r = True; Break[]], {m, 2, p}]; r); Select[Prime /@ Range[111], ok] (* Jean-François Alcover, Apr 22 2011 *)
    nn=1000; fact=1+Rest[FoldList[Times,1,Range[nn]]]; t={}; Do[p=Prime[i]; m=2; While[m
  • PARI
    is(p)=my(t=Mod(5040,p)); for(m=8, p-2, t*=m; if(t==-1 && p%m!=1, return(isprime(p)))); 0 \\ Charles R Greathouse IV, Feb 10 2013

Extensions

More terms from David W. Wilson, Sep 08 2001

A063828 Smallest m associated with n-th Pillai prime (A063980).

Original entry on oeis.org

14, 18, 15, 8, 18, 9, 23, 13, 86, 16, 16, 50, 102, 61, 64, 210, 97, 31, 9, 93, 40, 45, 63, 220, 91, 122, 35, 85, 198, 93, 128, 316, 366, 74, 300, 151, 290, 15, 400, 282, 22, 188, 167, 191, 360, 426, 274, 271, 456, 278, 229, 324, 135, 498, 189
Offset: 1

Views

Author

N. J. A. Sloane, Sep 23 2001

Keywords

Examples

			14! + 1 == 0 (mod 23), while 23 != 1 (mod 14), where 23=A063980(1) so a(1)=14.
		

Crossrefs

Programs

  • Mathematica
    nn=1000; fact=1+Rest[FoldList[Times,1,Range[nn]]]; t={}; Do[p=Prime[i]; m=2; While[m
  • PARI
    first(p)=my(t=Mod(5040, p)); for(m=8, p, t*=m; if(t==-1 && p%m!=1, return(m))); 0
    Pillai(p)=my(t=Mod(5040, p)); for(m=8, p-2, t*=m; if(t==-1 && p%m!=1, return(1))); 0
    apply(first, select(Pillai, primes(300))) \\ Charles R Greathouse IV, Feb 10 2013

Extensions

More terms from Vladeta Jovovic, Sep 27 2001
Showing 1-2 of 2 results.