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

A068650 Primes in which a string of 1's is sandwiched between two 3's.

Original entry on oeis.org

313, 3111111111113, 311111111111113, 3111111111111111111111111111113, 311111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111113
Offset: 1

Views

Author

Amarnath Murthy, Feb 28 2002

Keywords

Comments

The next term (a(6)) has 127 digits. - Harvey P. Dale, Dec 18 2022

Crossrefs

Programs

  • Maple
    a := 3:b := 1:i := 1:for n from 0 to 500 do c := a+10*(10^n-1)/9*b+10^(n+1)*a; if(isprime(c)) then d[i] := c; i := i+1; end if; end do:q := seq(d[j],j=1..i-1);
  • Mathematica
    Select[Table[10*FromDigits[PadRight[{3},n,1]]+3,{n,2,110}],PrimeQ] (* Harvey P. Dale, Dec 18 2022 *)

Extensions

More terms from Sascha Kurz, Mar 19 2002
Edited by Ray Chandler, Nov 04 2014

A082704 Numbers k such that (28*10^(k-1) + 17)/9 is a depression prime.

Original entry on oeis.org

3, 13, 15, 31, 105, 127, 343, 601, 9825
Offset: 1

Views

Author

Patrick De Geest, Apr 13 2003

Keywords

Comments

Prime versus probable prime status and proofs are given in the author's table.
a(10) > 2^16. - Lucas A. Brown, Apr 18 2021
a(10) > 2*10^5. - Tyler Busby, Feb 01 2023

Examples

			k=15 -> (28*10^(15-1) + 17)/9 = 311111111111113.
		

References

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

Crossrefs

Formula

a(n) = A056251(n) + 2.

Extensions

Edited by Ray Chandler, Nov 04 2014

A263182 Smallest k such that k//A002275(n)//k is prime, where // denotes concatenation and A002275(n) is the n-th repunit (R_n).

Original entry on oeis.org

1, 3, 13, 17, 1073, 19, 17, 29, 10000117, 73, 17, 3, 1007, 3, 43, 11, 1000000000000029, 1, 31, 11, 1191, 1, 1143, 31, 10000079, 21, 91, 59, 1019, 3, 67, 117, 10000000000000000000000000000077, 109, 89, 49, 1097, 41, 1053, 43, 10000047, 87, 23, 53, 1149, 83, 57
Offset: 0

Views

Author

Felix Fröhlich, Oct 11 2015

Keywords

Comments

a(n) = 3 if n is in A056251.
From Chai Wah Wu, Nov 05 2019 : (Start)
Theorem: a(2^r*s) >= 10^(2^r-1) for all r >= 0, s > 0.
Proof: Note that if k has m digits, then k//A002275(n)//k = k*(10^(n+m)+1) + A002275(n)*10^m which is a multiple of gcd(A002275(n),10^(n+m)+1).
Next, since 10^(2^r) - 1 = (10^(2^(r-1) -1))*(10^(2^(r-1) + 1)) and 9 does not divide 10^n+1, by induction it is easy to see that 10^(2^w) + 1 is a divisor of A002275(2^r) for 1 <= w < r. Since A002275(2^r) is a divisor of A002275(2^r*s), 10^(2^w) + 1 is also a divisor of A002275(2^r*s).
For 1 <= m < 2^r, let t be the 2-adic valuation of m, i.e. 0 <= t = A007814(m) < r.
Then 10^(2^r*s+m)+1 = 10^(2^t*q)+1 = (10^(2^t))^q + 1 for some odd number q.
Since the sum of two odd powers a^q+b^q is divisible by a+b, this implies that 10^(2^r*s+m)+1 is divisible by 10^(2^t)+1.
This means that for n = 2^r*s and 1 <= m < 2^r, gcd(A002275(n),10^(n+m)+1) >= 10^(2^t)+1 > 1, i.e. k//A002275(n)//k is not prime.
Thus a(2^r*s) must have at least 2^r digits, i.e. a(2^r*s) >= 10^(2^r-1). QED
As a consequence, a(n) >= 10^(A006519(n)-1). This result is still true if some of the digits of k are leading zeros.
(End)

Examples

			R_0 = 0 and the smallest k such that k//0//k is prime is 1, so a(0) = 1.
		

Crossrefs

Programs

  • Mathematica
    Table[k = 1; While[! PrimeQ[f[n, k]], k++]; k, {n, 0, 7}] (* Michael De Vlieger, Oct 13 2015 *)
  • PARI
    a(n) = my(rep=(10^n-1)/9, k=1); while(!ispseudoprime(eval(Str(k, rep, k))), k++); k

Formula

a(A004023(n)-2) = 1. - Chai Wah Wu, Nov 04 2019

Extensions

a(16)-a(46) from Chai Wah Wu, Nov 04 2019
Showing 1-3 of 3 results.