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-5 of 5 results.

A006549 Numbers k such that k and k+1 are prime powers.

Original entry on oeis.org

1, 2, 3, 4, 7, 8, 16, 31, 127, 256, 8191, 65536, 131071, 524287, 2147483647, 2305843009213693951, 618970019642690137449562111, 162259276829213363391578010288127, 170141183460469231731687303715884105727
Offset: 1

Views

Author

Keywords

Comments

Numbers k such that k + (0, 1) is a prime power pair.
Consecutive prime powers.
k + (0, 2m), m >= 1, being an admissible pattern for prime pairs, since (0, 2m) == (0, 0) (mod 2), has high density.
k + (0, 2m-1), m >= 1, being a non-admissible pattern for prime pairs, since (0, 2m-1) == (0, 1) (mod 2), has low density [the only possible pairs are (2^a - 2m-1, 2^a) or (2^a, 2^a + 2m-1), a >= 0].
Numbers k such that k and k+1 are primes would give only 2, for the prime pair (2, 3).
This sequence corresponds to the least member of each one of the following prime power pairs, ordered by increasing value of least member: (1, 2), (2^3, 3^2), (Fermat primes - 1, Fermat primes), (Mersenne primes, Mersenne primes + 1).
It is not known whether this sequence is infinite, but is conjectured to be since:
(*) 2^3, 3^2 are the only consecutive prime powers with exponents >= 2
(as a consequence of Mihailescu's theorem -- Mihailescu proved Catalan's conjecture in 2002);
(*) Only the first 5 Fermat numbers f_0 to f_4 are known to be prime
(it is conjectured that there might be no others, f_5 to f_32 are all composite);
(*) It has been conjectured that there exist an infinite number of Mersenne primes.
Numbers k such that A003418(k) appears only once in the sequence A003418. This may suggest that k is also characterized by the pairs formed by a 2 whose direct neighbor is a prime number in the sequence A014963. - Eric Desbiaux, Feb 11 2015
The power graph and enhanced power graph of the groups PGL(2,q) have the same clique number iff q>1 is a term of this sequence (Peter Cameron's link). - Bernard Schott, Dec 14 2021

References

  • R. K. Guy, Unsolved Problems in Number Theory, D9.
  • P. Ribenboim, 13 Lect. on Fermat's Last Theorem, p. 236.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • David W. Wilson and Eric Rains (rains(AT)caltech.edu) found a simple proof that in this case of Catalan's conjecture either n or n+1 must be a power of 2 and the other number must be a prime, except for n=8. Using this the sequence is easy to extend.

Crossrefs

Cf. A019434 Fermat primes: primes of form 2^(2^n) + 1, n >= 0.
Cf. A000668 Mersenne primes (of form 2^p - 1 where p is a prime).
Cf. A120431 Numbers n such that n and n+2 are prime powers.
Cf. A164571 Numbers n such that n and n+3 are prime powers.
Cf. A164572 Numbers n such that n and n+4 are prime powers.
Cf. A164573 Numbers n such that n and n+5 are prime powers.
Cf. A164574 Numbers n such that n and n+6 are prime powers.

Programs

  • Haskell
    a006549 n = a006549_list !! (n-1)
    a006549_list = [1,2,3,4,7,8] ++ f (drop 4 a000040_list) where
       f (p:ps) | a010055 (p - 1) == 1 = (p - 1) : f ps
                | a010055 (p + 1) == 1 = p : f ps
                | otherwise            = f ps
    -- Reinhard Zumkeller, Jan 03 2013
    
  • Mathematica
    Do[ a = Length[ FactorInteger[ 2^n - 1 ] ]; b = Length[ FactorInteger[ 2^n ] ]; c = Length[ FactorInteger[ 2^n + 1 ] ]; If[ a == b, Print[ 2^n - 1 ] ]; If[ b == c, Print[ 2^n ] ], {n, 0, 127} ]
    Join[{1},SequencePosition[Boole[PrimePowerQ[Range[600000]]],{1,1}][[All,1]]] (* Requires Mathematica version 10 or later *) (* Generates the first 14 terms of the sequence. Increase Range constant to generate more. *) (* Harvey P. Dale, Apr 12 2020 *)
  • PARI
    is(n)=if(n<5,return(n>0)); isprimepower(n) && isprimepower(n+1) \\ Charles R Greathouse IV, Apr 24 2015

Extensions

More terms from David W. Wilson
Additional comments from Daniel Forgues, Aug 17 2009

A120431 Numbers k such that k and k+2 are prime powers.

Original entry on oeis.org

1, 2, 3, 5, 7, 9, 11, 17, 23, 25, 27, 29, 41, 47, 59, 71, 79, 81, 101, 107, 125, 137, 149, 167, 179, 191, 197, 227, 239, 241, 269, 281, 311, 347, 359, 419, 431, 461, 521, 569, 599, 617, 641, 659, 727, 809, 821, 827, 839, 857, 881, 1019, 1031, 1049, 1061, 1091
Offset: 1

Views

Author

Greg Huber, Jul 13 2006

Keywords

Comments

Twin prime powers, a generalization of the twin primes. The twin primes are a subsequence.
From Daniel Forgues, Aug 17 2009: (Start)
Numbers k such that k + (0, 2) is a prime power pair.
k + (0, 2m), m >= 1, being an admissible pattern for prime pairs has high density.
k + (0, 2m-1), m >= 1, being a non-admissible pattern for prime pairs, has low density [the only possible pairs are (2^a - 2m-1, 2^a) or (2^a, 2^a + 2m-1), a >= 0.] (End)

Examples

			a(5) = 7 since the 5th pair of twin prime powers is (7,9), while the first four pairs are (1,3), (2,4), (3,5) and (5,7).
		

Crossrefs

Programs

  • Magma
    [1] cat [n: n in [2..1200] | IsPrimePower(n) and IsPrimePower(n+2)]; // Vincenzo Librandi, Nov 03 2018
  • Maple
    isppow := proc(n) local pf; pf := ifactors(n)[2]; if nops(pf) = 1 or n =1 then true; else false; fi; end; isA120431 := proc(n) RETURN (isppow(n) and isppow(n+2)); end; for n from 1 to 1500 do if isA120431(n) then printf("%d, ",n); fi; od; # R. J. Mathar, Dec 16 2006
  • Mathematica
    Join[{1}, Select[Range[1100], And@@PrimePowerQ/@{#, # + 2} &]] (* Vincenzo Librandi, Nov 03 2018 *)
  • PARI
    is(n)=if(n<4,return(n>0)); isprimepower(n) && isprimepower(n+2) \\ Charles R Greathouse IV, Apr 24 2015
    

Formula

a(n) = A064076(n-2) for n >= 3. - Georg Fischer, Nov 02 2018

Extensions

More terms from R. J. Mathar, Dec 16 2006

A164572 Numbers k such that k and k+4 are both prime powers.

Original entry on oeis.org

1, 3, 4, 5, 7, 9, 13, 19, 23, 25, 27, 37, 43, 49, 67, 79, 97, 103, 109, 121, 127, 163, 169, 193, 223, 229, 239, 277, 289, 307, 313, 343, 349, 379, 397, 439, 457, 463, 487, 499, 613, 643, 673, 729, 739, 757, 769, 823, 853, 859, 877, 883, 907, 937, 967, 1009, 1087
Offset: 1

Views

Author

Daniel Forgues, Aug 16 2009, Aug 17 2009

Keywords

Comments

Numbers n such that n + (0, 4) is a prime power pair.
A generalization of the cousin primes. The cousin primes are a subsequence.
n + (0, 2m), m >= 1, being an admissible pattern for prime pairs, since (0, 2m) = (0, 0) (mod 2), has high density.
n + (0, 2m-1), m >= 1, being a non-admissible pattern for prime pairs, since (0, 2m-1) = (0, 1) (mod 2), has low density [the only possible pairs are (2^a - 2m-1, 2^a) or (2^a, 2^a + 2m-1), a >= 0.]

Crossrefs

k and (x) are prime powers: A006549 (k+1), A120431 (k+2), A164571 (k+3), this sequence (k+4), A164573 (k+5), A164574 (k+6).

Programs

  • Mathematica
    Select[Range[1000], PrimeNu[#] < 2 && PrimeNu[# + 4] < 2 &] (* Amiram Eldar, Oct 01 2020 *)
  • PARI
    is(n)=if(n==1,return(1)); isprimepower(n) && isprimepower(n+4) \\ Charles R Greathouse IV, Apr 24 2015

A164571 Numbers n such that n and n+3 are prime powers.

Original entry on oeis.org

1, 2, 4, 5, 8, 13, 16, 29, 61, 64, 125, 128, 509, 1021, 4093, 4096, 16381, 32768, 65536, 262144, 1048573, 4194301, 16777213, 268435456, 536870909, 1073741824, 36028797018963968
Offset: 1

Views

Author

Daniel Forgues, Aug 16 2009

Keywords

Comments

Numbers n such that n + (0, 3) is a prime power pair.
n + (0, 2m), m >= 1, being an admissible pattern for prime pairs, since (0, 2m) = (0, 0) (mod 2), has high density.
n + (0, 2m-1), m >= 1, being a non-admissible pattern for prime pairs, since (0, 2m-1) = (0, 1) (mod 2), has low density [the only possible pairs are (2^a - 2m-1, 2^a) or (2^a, 2^a + 2m-1), a >= 0.]
n + (0, 3) being a non-admissible pattern for prime pairs, has only prime power pairs (2^a - 3, 2^a) or (2^a, 2^a + 3), a >= 0.
Numbers n such that n and n+3 are primes would give only 2, for the prime pair (2, 5).
10^18 < a(28) <= 19807040628566084398385987581. - Donovan Johnson, Aug 17 2009

Crossrefs

Cf. A000961.
Cf. A006549 Numbers n such that n and n+1 are prime powers.
Cf. A120431 Numbers n such that n and n+2 are prime powers.
Cf. A164571 Numbers n such that n and n+3 are prime powers.
Cf. A164572 Numbers n such that n and n+4 are prime powers.
Cf. A164573 Numbers n such that n and n+5 are prime powers.
Cf. A164574 Numbers n such that n and n+6 are prime powers.

Programs

  • PARI
    ispp(n) = (n==1) || isprime(n) || (ispower(n, ,&p) && isprime(p));
    isok(n) = ispp(n) && ispp(n+3); \\ Michel Marcus, Aug 31 2013
    
  • PARI
    v=List(); for(n=0, 1e3, if(isprimepower(2^n-3), listput(v, 2^n-3)); if(isprimepower(2^n+3), listput(v, 2^n))); Set(v) \\ Charles R Greathouse IV, Apr 24 2015

Extensions

Edited by Daniel Forgues, Aug 17 2009
a(20)-a(27) from Donovan Johnson, Aug 17 2009

A164574 Numbers k such that k and k+6 are both prime powers.

Original entry on oeis.org

1, 2, 3, 5, 7, 11, 13, 17, 19, 23, 25, 31, 37, 41, 43, 47, 53, 61, 67, 73, 83, 97, 101, 103, 107, 121, 125, 131, 151, 157, 163, 167, 173, 191, 193, 223, 227, 233, 251, 257, 263, 271, 277, 283, 307, 311, 331, 337, 343, 347, 353, 361, 367, 373, 383, 433, 443, 457
Offset: 1

Views

Author

Daniel Forgues, Aug 16 2009

Keywords

Comments

Numbers n such that n + (0, 6) is a prime power pair.
n + (0, 2m), m >= 1, being an admissible pattern for prime pairs, since (0, 2m) = (0, 0) (mod 2), has high density.
n + (0, 2m-1), m >= 1, being a non-admissible pattern for prime pairs, since (0, 2m-1) = (0, 1) (mod 2), has low density [the only possible pairs are (2^a - 2m-1, 2^a) or (2^a, 2^a + 2m-1), a >= 0.]

Crossrefs

k and (x) are prime powers: A006549 (k+1) A120431 (k+2), A164571 (k+3), A164572 (k+4), A164573 (k+5), this sequence (k+6).

Programs

  • Mathematica
    Join[{1},Select[Range[500],AllTrue[{#,#+6},PrimePowerQ]&]] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Sep 30 2018 *)
  • PARI
    is(n)=if(n<4,return(n>0)); isprimepower(n) && isprimepower(n+6) \\ Charles R Greathouse IV, Apr 24 2015

Extensions

Edited by Daniel Forgues, Aug 17 2009
Showing 1-5 of 5 results.