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.

Previous Showing 11-19 of 19 results.

A231865 Numbers n such that (43^n + 1)/44 is prime.

Original entry on oeis.org

5, 7, 19, 251, 277, 383, 503, 3019, 4517, 9967, 29573
Offset: 1

Views

Author

Robert Price, Nov 14 2013

Keywords

Comments

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

Crossrefs

Programs

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

A127727 Primes of the form p^e - p^(e-1) + p^(e-2) - ... + (-1)^e, where p is prime.

Original entry on oeis.org

2, 3, 5, 7, 11, 43, 61, 157, 521, 547, 683, 2731, 4423, 6163, 13421, 19183, 22651, 26407, 37057, 43691, 113233, 121453, 143263, 174763, 208393, 292141, 375157, 398581, 412807, 527803, 590593, 843643, 981091, 1041421, 1193557, 1246573
Offset: 1

Views

Author

T. D. Noe, Jan 25 2007

Keywords

Comments

These primes are important in studying k-imperfect numbers (A127724), see Iannucci-link. Except for the cases p^e = 3 and 8, which yield primes 2 and 5, e is an even number such that e+1 is prime. In fact, except for those two cases, all the primes are of the form (1+p^q)/(1+p), where q is an odd prime; that is, repunit primes with negative prime base.

Examples

			From _David A. Corneth_, Oct 28 2017: (Start)
For (p, e) = (3, 1) we have the prime 3^1 - 3^0 = 2.
For (p, e) = (2, 3) we have the prime 2^3 - 2^2 + 2^1 - 2^0 = 5.
The examples above are the cases mentioned in the comments not of the form (1+p^q)/(1+p). A prime of that form is below;
For (p, e) = (2, 4) we have the prime 2^4 - 2^3 + 2^2 - 2^1 + 2^0 = 11 = (1+p^(e+1)) / (1+p) = 33/3. (End)
		

Crossrefs

Programs

  • PARI
    upto(n) = {my(res = List([2,5])); forprime(p = 2, sqrtnint(n, 2), forprime(q = 3, logint(n * (1+p), p), r = (1+p^q)/(1+p); if(isprime(r), listput(res, r)))); listsort(res, 1); res} \\ David A. Corneth, Oct 28 2017

A235683 Numbers n such that (46^n + 1)/47 is prime.

Original entry on oeis.org

7, 23, 59, 71, 107, 223, 331, 2207, 6841, 94841
Offset: 1

Views

Author

Robert Price, Jan 13 2014

Keywords

Comments

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

Crossrefs

Programs

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

A237052 Numbers n such that (49^n + 1)/50 is prime.

Original entry on oeis.org

7, 19, 37, 83, 1481, 12527, 20149
Offset: 1

Views

Author

Robert Price, Feb 02 2014

Keywords

Comments

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

Crossrefs

Programs

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

Extensions

Typo in description corrected by Ray Chandler, Feb 20 2017

A126589 Numbers n>1 such that prime of the form (n^k-1)/(n-1) does not exist for k>2; or A128164(n) = 0.

Original entry on oeis.org

4, 9, 16, 25, 32, 36, 49, 64, 81, 100, 121, 125, 144, 169, 196, 216, 225, 243, 256, 289, 324, 343, 361, 400, 441, 484, 529, 576, 625, 676, 729, 784, 841, 900, 961, 1000, 1024, 1089, 1156, 1225, 1296, 1369, 1444, 1521, 1600, 1681, 1728, 1764, 1849, 1936, 2025
Offset: 1

Views

Author

Alexander Adamchuk, Mar 13 2007

Keywords

Comments

Appears to be the union of the perfect squares k^2 (for k>1) and the prime powers p^k (for k>1) with some exceptions, such as 2^3, 3^3, 2^7, etc.
The perfect powers except those of the form n^(p^m) where p and (n^(p^(m+1))-1)/(n^(p^m)-1) are primes, p>2 and m>=1. - Max Alekseyev, Mar 09 2009

Examples

			A128164 begins with offset 2: {3, 3, 0, 3, 3, 5, 3, 0, 19, 17, 3, 5, 3, 3, 0, 3, ...}. Thus a(1) = 4, a(2) = 9, a(3) = 16.
		

Crossrefs

Extensions

Extended by Max Alekseyev, Mar 09 2009

A236167 Numbers k such that (47^k + 1)/48 is prime.

Original entry on oeis.org

5, 19, 23, 79, 1783, 7681
Offset: 1

Views

Author

Robert Price, Jan 19 2014

Keywords

Comments

a(7) > 10^5.

Crossrefs

Programs

  • Mathematica
    Do[ p=Prime[n]; If[ PrimeQ[ (47^p + 1)/48 ], Print[p] ], {n, 1, 9592} ]
  • PARI
    is(n)=ispseudoprime((47^n+1)/48) \\ Charles R Greathouse IV, Jun 06 2017
    
  • Python
    from sympy import isprime
    def afind(startat=0, limit=10**9):
      pow47 = 47**startat
      for k in range(startat, limit+1):
        q, r = divmod(pow47+1, 48)
        if r == 0 and isprime(q): print(k, end=", ")
        pow47 *= 47
    afind(limit=300) # Michael S. Branicky, May 19 2021

A343589 Smallest prime of the form n^k-(n-1) or 0 if no such prime exists.

Original entry on oeis.org

3, 7, 13, 3121, 31, 43, 549755813881, 73, 991, 1321, 248821, 157, 2731, 211, 241, 34271896307617, 307, 6841, 13107199999999999999981, 421, 463, 141050039560662968926081, 331753, 601, 17551, 7625597484961, 757, 1816075630094014572464024421543167816955354437761
Offset: 2

Views

Author

Blake Branstool, Apr 20 2021

Keywords

Comments

All values up to n=70 have been found and proved to be primes. n=71 has k=3019 and gives a probable prime.
See A113516, which gives the k values and is the main entry for these primes, for more extensively researched information. - Peter Munn, Nov 20 2021

Examples

			For n=2 and k=2, 2^2-(2-1)=3 thus a(2)=3. k is 2 as well for n=3,4.
For n=5 the first k to result in a prime is 5, 5^5-(5-1)=3121 thus a(5)=3121.
		

Crossrefs

A113516 gives the k values.

Programs

  • PARI
    a(n) = my(k=1, p); while (!isprime(p=n^k-(n-1)), k++); p; \\ Michel Marcus, Nov 17 2021

Extensions

Name revised by Peter Munn, Nov 16 2021

A185230 Numbers n such that (33^n + 1)/34 is prime.

Original entry on oeis.org

5, 67, 157, 12211, 313553
Offset: 1

Views

Author

Robert Price, Aug 29 2013

Keywords

Comments

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

Crossrefs

Programs

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

Extensions

a(5) from Paul Bourdelais, Feb 26 2021

A236530 Numbers n such that (48^n + 1)/49 is prime.

Original entry on oeis.org

5, 17, 131, 84589
Offset: 1

Views

Author

Robert Price, Jan 27 2014

Keywords

Comments

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

Crossrefs

Programs

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

Extensions

Incorrect first term deleted by Robert Price, Feb 21 2014
Previous Showing 11-19 of 19 results.