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

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

A309533 Numbers k such that (144^k + 1)/145 is prime.

Original entry on oeis.org

23, 41, 317, 3371, 45259, 119671
Offset: 1

Views

Author

Paul Bourdelais, Aug 06 2019

Keywords

Comments

The corresponding primes are terms of A059055. - Bernard Schott, Aug 09 2019

Crossrefs

Programs

  • Mathematica
    Do[p=Prime[n]; If[PrimeQ[(144^p + 1)/145], Print[p]], {n, 1, 1000000}]
  • PARI
    is(n)=ispseudoprime((144^n+1)/145)

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

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

A347138 Numbers k such that (100^k + 1)/101 is prime.

Original entry on oeis.org

3, 293, 461, 11867, 90089
Offset: 1

Views

Author

Paul Bourdelais, Aug 19 2021

Keywords

Comments

These are the repunit primes in base -100. It is unusual to represent numbers in a negative base, but it follows the same formulation as any base: numbers are represented as a sum of powers in that base, i.e., a0*1 + a1*b^1 + a2*b^2 + a3*b^3 ... Since the base is negative, the terms will be alternating positive/negative. For repunits the coefficients are all ones so the sum reduces to 1 + b + b^2 + b^3 + ... + b^(k-1) = (b^k-1)/(b-1). Since b is negative and k is an odd prime, the sum equals (|b|^k+1)/(|b|+1). For k=3, the sum is 9901, which is prime. As with all repunits, we only need to PRP test the prime exponents. The factors of repunits base -100 will be of the form p=2*k*m+1 where m must be even, which is common for (negative) bases that are squares.

Examples

			3 is a term since (100^3 + 1)/101 = 9901 is a prime.
		

Crossrefs

Programs

  • Mathematica
    Do[ If[ PrimeQ[ (100^n + 1)/101], Print[n]], {n, 0, 18000}]
  • PARI
    is(n)=isprime((100^n+1)/101)

A348170 Numbers k such that (35^k - 1)/34 is prime.

Original entry on oeis.org

313, 1297, 568453
Offset: 1

Views

Author

Paul Bourdelais, Oct 04 2021

Keywords

Comments

These are the repunit primes in base 35.

Examples

			313 is a term since (35^313 - 1)/34 is a prime. It has 482 digits in base 10.
		

Crossrefs

Programs

  • Mathematica
    Do[ If[ PrimeQ[ (35^n-1)/34], Print[n]], {n, 0, 600000}]
  • PARI
    is(n)=isprime((35^n-1)/34)

A350036 Numbers k such that (81^k + 1)/82 is prime.

Original entry on oeis.org

3, 5, 701, 829, 1031, 1033, 7229, 19463, 370421
Offset: 1

Views

Author

Paul Bourdelais, Dec 09 2021

Keywords

Comments

These are the Repunits in base -81. Since 81=3^4, factors will be of the form p=8nk+1. (Negative) bases that are powers of small numbers appear to have a higher frequency of primes than Repunits in other bases. The best linear fit for this base is currently 0.29918 which is much lower (better) than the conjectured 0.56145948 (see link to conjecture).

Examples

			3 is a term since (81^3 + 1)/82 = 6481 is a prime.
		

Crossrefs

Programs

  • Mathematica
    Do[ If[ PrimeQ[ (81^n+1)/82], Print[n]], {n, 0, 1000000}]
  • PARI
    is(n)=isprime((81^n+1)/82)
Previous Showing 11-18 of 18 results.