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.

A183177 Numbers n such that (10^(2n+1)+15*10^n-1)/3 is prime.

Original entry on oeis.org

1, 7, 85, 94, 273, 356, 1077, 1797, 6758, 30232
Offset: 1

Views

Author

Ray Chandler, Dec 28 2010

Keywords

Comments

a(11) > 10^5. - Robert Price, Apr 21 2016

References

  • C. Caldwell and H. Dubner, "Journal of Recreational Mathematics", Volume 28, No. 1, 1996-97, pp. 1-9.

Crossrefs

Programs

  • Mathematica
    Do[If[PrimeQ[(10^(2n + 1) + 15*10^n - 1)/3], Print[n]], {n, 3000}]
  • PARI
    is(n)=ispseudoprime((10^(2*n+1)+15*10^n-1)/3) \\ Charles R Greathouse IV, Jun 13 2017

Formula

a(n) = (A077792(n)-1)/2.

Extensions

a(10) from Robert Price, Apr 21 2016

A332138 a(n) = (10^(2*n+1)-1)/3 + 5*10^n.

Original entry on oeis.org

8, 383, 33833, 3338333, 333383333, 33333833333, 3333338333333, 333333383333333, 33333333833333333, 3333333338333333333, 333333333383333333333, 33333333333833333333333, 3333333333338333333333333, 333333333333383333333333333, 33333333333333833333333333333, 3333333333333338333333333333333
Offset: 0

Views

Author

M. F. Hasler, Feb 09 2020

Keywords

Comments

See A183177 = {1, 7, 85, 94, 273, 356, ...} for the indices of primes.

Crossrefs

Cf. (A077792-1)/2 = A183177: indices of primes.
Cf. A002275 (repunits R_n = (10^n-1)/9), A002277 (3*R_n), A011557 (10^n).
Cf. A138148 (cyclops numbers with binary digits), A002113 (palindromes).
Cf. A332118 .. A332178, A181965 (variants with different repeated digit 1, ..., 9).
Cf. A332130 .. A332139 (variants with different middle digit 0, ..., 9).

Programs

  • Maple
    A332138 := n -> (10^(2*n+1)-1)/3+5*10^n;
  • Mathematica
    Array[ (10^(2 # + 1)-1)/3 + 5*10^# &, 15, 0]
  • PARI
    apply( {A332138(n)=10^(n*2+1)\3+5*10^n}, [0..15])
    
  • Python
    def A332138(n): return 10**(n*2+1)//3+5*10**n

Formula

a(n) = 3*A138148(n) + 8*10^n = A002277(2n+1) + 5*10^n.
G.f.: (8 - 505*x + 200*x^2)/((1 - x)(1 - 10*x)(1 - 100*x)).
a(n) = 111*a(n-1) - 1110*a(n-2) + 1000*a(n-3) for n > 2.
Showing 1-2 of 2 results.