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-10 of 22 results. Next

A126856 Numbers n such that (31^n + 1)/32 is prime.

Original entry on oeis.org

109, 461, 1061, 50777
Offset: 1

Views

Author

Alexander Adamchuk, Mar 23 2007

Keywords

Comments

All terms are primes.
a(5) > 10^5. - Robert Price, Jul 12 2013

Crossrefs

Programs

  • Mathematica
    Do[ p=Prime[n]; If[ PrimeQ[ (31^p + 1)/32 ], Print[p] ], {n,1,1100} ]
  • PARI
    is(n)=isprime((31^n+1)/32) \\ Charles R Greathouse IV, Feb 17 2017

Extensions

a(4) from Robert Price, Jul 12 2013

A229145 Numbers k such that (36^k + 1)/37 is prime.

Original entry on oeis.org

31, 191, 257, 367, 3061, 110503, 1145393
Offset: 1

Views

Author

Robert Price, Sep 15 2013

Keywords

Comments

All such numbers k are prime.
Note that a(6) = 110503 corresponds to (36^110503 + 1)/37, which is only a probable prime with 171975 digits.
The primes corresponding to the terms of this sequence have 1 as their last digit and an even number as their next-to-last digit. - Iain Fox, Dec 08 2017

Crossrefs

Programs

  • Mathematica
    Do[ p=Prime[n]; If[ PrimeQ[ (36^p + 1)/37 ], Print[p] ], {n, 1, 9592} ]
  • PARI
    is(n)=isprime((36^n+1)/37) \\ Charles R Greathouse IV, Feb 17 2017

Extensions

a(6) = 110503 (posted by Lelio R. Paula on primenumbers.net) from Paul Bourdelais, Dec 08 2017
a(7) from Paul Bourdelais, Nov 03 2023

A185240 Numbers k such that (35^k + 1)/36 is prime.

Original entry on oeis.org

11, 13, 79, 127, 503, 617, 709, 857, 1499, 3823, 135623, 280979
Offset: 1

Views

Author

Robert Price, Aug 29 2013

Keywords

Comments

All terms are primes. a(11) > 10^5.

Crossrefs

Programs

  • Mathematica
    Do[ p=Prime[n]; If[ PrimeQ[ (35^p + 1)/36 ], Print[p] ], {n, 1, 9592} ]
  • PARI
    is(n)=isprime((35^n+1)/36) \\ Charles R Greathouse IV, Feb 17 2017

Extensions

a(11)=135623 found as probable prime and added by Paul Bourdelais, Jul 05 2018
a(12) from Paul Bourdelais, Sep 13 2021

A229524 Numbers k such that (38^k + 1)/39 is prime.

Original entry on oeis.org

5, 167, 1063, 1597, 2749, 3373, 13691, 83891, 131591
Offset: 1

Views

Author

Robert Price, Sep 25 2013

Keywords

Comments

All terms are primes. a(9) > 10^5.

Crossrefs

Programs

  • Mathematica
    Do[ p=Prime[n]; If[ PrimeQ[ (38^p + 1)/39 ], Print[p] ], {n, 1, 9592} ]
  • PARI
    is(n)=ispseudoprime((38^n+1)/39) \\ Charles R Greathouse IV, Feb 17 2017

Extensions

a(9)=131591 corresponds to a probable prime discovered by Paul Bourdelais, Jul 03 2018

A073936 Numbers k such that 2^k + 1 is the product of two distinct primes.

Original entry on oeis.org

5, 6, 7, 11, 12, 13, 17, 19, 20, 23, 28, 31, 32, 40, 43, 61, 64, 79, 92, 101, 104, 127, 128, 148, 167, 191, 199, 256, 313, 347, 356, 596, 692, 701, 1004, 1228, 1268, 1709, 2617, 3539, 3824, 5807, 10501, 10691, 11279, 12391, 14479, 42737, 83339, 95369, 117239
Offset: 1

Views

Author

Labos Elemer, Aug 13 2002

Keywords

Comments

Original name: "2^n + 1 is squarefree and has exactly 2 prime factors."
From Giuseppe Coppoletta, May 08 2017: (Start)
As 3 divides 2^a(n) + 1 for any odd term a(n), all odd terms are prime and exactly the Wagstaff primes (A000978), at the exclusion of 3 (which gives 2^3 + 1 = 3^2 not squarefree).
For the even terms, let a(n) = d * 2^j with d odd integer and j > 0. If d > 1, as (2^2^j)^q + 1 divides 2^a(n) + 1 for any odd prime q dividing d, then d must be prime.
So the even terms are all given by the following two class:
a) (d = 1) a(n) = 2^j such that Fj is a semiprime Fermat number. Up to now, only j = 5, 6, 7, 8 are known to give a Fermat semiprime, giving the even terms 32, 64, 128 and 256. We are also assured that 2^j is not a term for j = 9..19 because Fj is not a semiprime for those value of j (see Wagstaf's link). F20 is the first composite Fermat number which could give another even term (it would be 2^20 = 1048576). However, it seems highly unlikely that other Fermat semiprimes could exist.
b) (d = p odd prime) a(n) = p * 2^j with j such that Fj is a Fermat prime and p a prime verifying ((Fj - 1)^p + 1)/Fj is a prime.
Exemplifying that, we have:
for j = 1 this gives only the even term a(2) = 2 * 3 = 6 (see Jack Brennen's result in ref),
for j = 2 we have all the terms of type 2^2 * A057182.
for j = 3 the even terms are of type 2^3 * A127317.
For j = 4 at least up to 200000, there is only the term a(41) = 2^4 * 239 = 3824 (see comment in A127317).
All terms after a(50) refer to probabilistic primality tests for 2^a(n) + 1 (see Caldwell's link for the list of the largest certified Wagstaff primes).
After a(56), from the above, the primes 267017, 269987, 374321, 986191, 4031399 and the even value 4101572 are also terms, but still remains the (remote) possibility of some gaps in between. In addition, 13347311 and 13372531 are also terms, but possibly much further in the numbering (see comments in A000978).
(End)
Intersection of A092559 and A066263. - Eric Chen, Jun 13 2018

Examples

			11 is a member because 1 + 2^11 = 2049 = 3 * 683.
9 is not a term because 1 + 2^9 = 513 = 3^3 * 19
		

Crossrefs

Cf. A000005, A000051, A046798, A092559, A000978. Different from A066263.

Programs

  • Mathematica
    Do[ If[ Length[ Divisors[1 + 2^n]] == 4, Print[n]], {n, 1, 200}]
    (* Second program: *)
    Select[Range@ 200, DivisorSigma[0, 2^# + 1] == 4 &] (* Michael De Vlieger, May 09 2017 *)
  • Sage
    [n for n in xsrange(3,200) if sigma(2^n+1,0)==4]
    # Second program (faster):
    
  • Sage
    v=[]; N=2000
    for n in xsrange(4,N):
        j=valuation(n,2)
        if j<5:
            Fj=2^2^j+1; p=ZZ(n/2^j); q=ZZ((2^n+1)/Fj)
            if p.is_prime() and q.is_prime(proof=false): v.append(n)
        elif j<9 and n.is_power_of(2): v.append(n)
    print(v) # Giuseppe Coppoletta, May 11 2017

Formula

Solutions to A000005[A000051(x)]=4 or A046798[x]=4

Extensions

Edited by Robert G. Wilson v, Aug 19 2002
a(28)-a(51) by Giuseppe Coppoletta, May 02 2017
Name reworded by Jon E. Schoenfield, Jun 15 2018

A229663 Numbers n such that (40^n + 1)/41 is prime.

Original entry on oeis.org

53, 67, 1217, 5867, 6143, 11681, 29959
Offset: 1

Views

Author

Robert Price, Sep 27 2013

Keywords

Comments

All terms are primes.
a(8) > 10^5.

Crossrefs

Programs

  • Mathematica
    Do[ p=Prime[n]; If[ PrimeQ[ (40^p + 1)/41 ], Print[p] ], {n, 1, 9592} ]
  • PARI
    is(n)=ispseudoprime((40^n+1)/41) \\ Charles R Greathouse IV, Feb 17 2017

A230036 Numbers n such that (39^n + 1)/40 is prime.

Original entry on oeis.org

3, 13, 149, 15377
Offset: 1

Views

Author

Robert Price, Oct 05 2013

Keywords

Comments

All terms are primes.
a(5) > 10^5.

Crossrefs

Programs

  • Mathematica
    Do[ p=Prime[n]; If[ PrimeQ[ (39^p + 1)/40 ], Print[p] ], {n, 1, 9592} ]
  • PARI
    is(n)=ispseudoprime((39^n+1)/40) \\ Charles R Greathouse IV, Feb 17 2017

A066263 Numbers k such that 2^k + 1 has just two distinct prime factors.

Original entry on oeis.org

5, 6, 7, 9, 10, 11, 12, 13, 17, 19, 20, 23, 28, 31, 32, 40, 43, 61, 64, 79, 92, 101, 104, 127, 128, 148, 167, 191, 199, 256, 313, 347, 356, 596, 692, 701, 1004, 1228, 1268, 1709, 2617, 3539, 3824, 5807, 10501, 10691, 11279, 12391, 14479, 42737, 83339, 95369
Offset: 1

Views

Author

Benoit Cloitre, Dec 31 2001

Keywords

Comments

From Giuseppe Coppoletta, May 16 2017: (Start)
All terms after a(52) refer to probabilistic primality tests for 2^a(n) + 1 (see Caldwell's link for the list of the largest certified Wagstaff primes). After a(58), 267017, 269987, 374321, 986191, 4031399 and 4101572 are also terms, but there still remains the remote possibility of some gaps in between. In addition, 13347311 and 13372531 are also terms, but possibly much farther in the numbering (see comments in A000978).
For the relation with Fermat numbers and for the primality of odd terms, see comments in A073936. The terms 9 and 10 give a value of 2^n + 1 which is not squarefree, so they are not in A073936. For the rest, the actually known terms of the two sequences coincide. In order to verify if any other term could be found hereafter that is not in A073936, all we have to do is to examine the terms for which 2^n + 1 is not squarefree. Considering that 3 divides 2^a(n) + 1 for any odd term a(n) and using Zsigmondy's and Mihăilescu-Catalan's theorems (see links), one can verify that any nonsquarefree term greater than 9 has to be of the form a(n) = 2^j * Fj, where Fj is the Fermat prime 2^2^j + 1. So basically we have to see if ((Fj-1)^Fj + 1)/(Fj)^2 is a prime or the power of a prime for any Fermat prime Fj. The case j = 1 gives the term a(n) = 10 because (4^5 + 1)/5^2 = 41 is a prime, while for j = 2, (16^17 + 1)/17^2 = 354689 * 2879347902817 is composite. Similarly (256^257 + 1)/257^2 is neither a prime nor the power of a prime, so there is no contribution from the cases j = 2, 3 (see also comments in A127317).
For j = 4 and for any possible other Fermat prime which could be found later, the question is still open, in the sense that it is not actually known if n = 16 * F4 = 1048592 is a term or not. That seems very unlikely, but in order to decide that question for j = 4, one would have to check if (2^1048592 + 1)/65537^2 is a prime or the power of a prime. As this number has 315649 digits, I wonder if it is possible to handle it with the existing primality tests.
(End)

Examples

			3 and 4 are not terms because 2^3 + 1 and 2^4 + 1 have only a single prime factor (counted without multiplicity).
6 and 10 are terms because 2^6 + 1 = 5 * 13 and 2^10 + 1 = 5^2 * 41 have two distinct prime factors.
		

Crossrefs

Programs

  • Mathematica
    f[n_] := First[ Transpose[ FactorInteger[2^n + 1]]]; Select[ Range[100], Length[f[ # ]] == 2 & ]
    Select[Range[1300],PrimeNu[2^#+1]==2&] (* Harvey P. Dale, Nov 28 2014 *)
  • PARI
    isok(k) = #factor(2^k+1)~ == 2; \\ Michel Marcus, Nov 14 2017

Formula

A001221(2^a(n) + 1) = 2.

Extensions

Edited by Robert G. Wilson v, Jan 03 2002
a(40)-a(52) by Giuseppe Coppoletta, May 02 2017

A231604 Numbers n such that (42^n + 1)/43 is prime.

Original entry on oeis.org

3, 709, 1637, 17911, 127609, 172663
Offset: 1

Views

Author

Robert Price, Nov 11 2013

Keywords

Comments

The first 5 terms are primes.

Crossrefs

Programs

  • Mathematica
    Do[ p=Prime[n]; If[ PrimeQ[ (42^p + 1)/43 ], Print[p] ], {n, 1, 9592} ]
  • PARI
    is(n)=ispseudoprime((42^n+1)/43) \\ Charles R Greathouse IV, Feb 20 2017

Extensions

a(5)=127609 corresponds to a probable prime discovered by Paul Bourdelais, Jul 02 2018
a(6)=172663 corresponds to a probable prime discovered by Paul Bourdelais, Jul 29 2019

A127317 Numbers n such that (256^n + 1)/257 is prime.

Original entry on oeis.org

5, 13, 23029
Offset: 1

Views

Author

Alexander Adamchuk, Mar 29 2007

Keywords

Comments

All terms are primes. Largest currently known prime of the form (2^n + 1)/257 is (256^23029 + 1)/257 found by Donovan Johnson 03/2005. The only currently known prime of the form (2^n + 1)/65537 is (65536^239 + 1)/65537.
From Giuseppe Coppoletta, May 18 2017: (Start)
In general, for any j > 1, if (2^(n*2^j) + 1)/Fj is a prime (where Fj = 2^2^j + 1 is the corresponding Fermat number), then n needs to be prime, as for any odd proper factor q of n, 2^(q*2^j) + 1 is another factor of the numerator. The same for j = 0, apart for the particular value n = 3^2.
For the case j = 4, I checked it again, and (65536^p + 1)/65537 indeed is not a prime at least for 239 < p < 12500, i.e. (2^n + 1)/65537 is not a prime at least up to n = 200000. Any higher upper bound available?
One can also remark that 65536 = 2^16 and 239 = 2^8 - 2^4 - 1. Is there any special reason (see Brennen's link) for that?
I checked also that (2^(p*2^j) + 1)/Fj is never a proper power (in particular it is not a prime power) for j = 0..4 and for any prime p, at least for any exponent p*2^j < 200000.
We can even conjecture that ((Fj-1)^p + 1)/Fj is always squarefree for any odd prime p and for any Fermat number Fj with j >= 0. Note that this is not true if p is not restricted to be a prime, even if p and Fj are coprime, as shown by the following counterexample relative to the case j = 1, f1 = 5: 4^91 + 1 == 0 mod 1093^2. Remark that any such counterexample has to be a Wieferich prime (A001220), but not every Wieferich prime gives a counterexample, as shown by the second known Wieferich prime (3511), which cannot match here because it belongs to A072936.
(End)

Crossrefs

Cf. A000978 = numbers n such that (2^n + 1)/3 is prime. Cf. A057182 = numbers n such that (16^n + 1)/17 is a prime.

Programs

  • Mathematica
    Do[n=8*Prime[k];f=2^n+1;If[PrimeQ[f/257],Print[{n,n/8}]],{k,1,2570}]
Showing 1-10 of 22 results. Next