login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A307873
The concatenation pkp is the number obtained by placing prime p either side of R_k, the k-th repunit (1, k times); a(n) is the smallest k such that pkp is prime, where p=prime(n), or -1 if no such k exists.
0
-1, 1, -1, 10905, 15, 2, 1, 2, 3, 1, 3, 173, 1, 14, 1, 43, 1, 5, 11, 1, 2, 3, 3, 1, 2, -1, 5, 421, 3, 1, -1, 1, 1, 3, -1, 15, -1, 3, 3, 163, -1, 3, 13, -1, 679, -1, 5, 5, -1, 107, 93, 1, -1, 3, -1, 1, -1, 9, 5, -1, -1, 9, 1089, -1, 3, 7, 3, 15, -1, 27, -1, 1, -1, 27, 17, 25, 1, 15, 3
OFFSET
1,4
COMMENTS
Primes (from above data) for which pkp is composite for all k are 2, 5, 101, 127, 149, 157, 179, 193, 199, 227, 241, 257, 269, 281, 283, 311, 347, 353, 367. In every case the factorization of pkp contains at least one characteristic prime divisor (very different from A306861).
Conjecture: There are an infinite number of -1 terms in this sequence.
FORMULA
If prime(n) is a repunit prime R_k, for some k in A004023 and R_t is the smallest repunit prime such that t > 2*k, then a(n)=R_(t-2*k).
EXAMPLE
2/2k2, 5/5k5, 7/101k101, 11,13/127k127, 11/149k149, for all k, so a(1)=a(3)=a(26)= a(31)=a(35)=-1. For prime(n)=A004023(2)=R_19, a(n)=R_(317-2*19)=R_279.
MAPLE
P(p) := proc (p::prime, N::posint := 5000) local n, k, m0, m; n := length(p); for k from 1 to N do m0 := add(10^i, i = 0 .. k-1); m := p*10^(k+n)+m0*10^n+p; if isprime(m) then return k end if; if `mod`(k, 1000) = 0 then print(k) end if end do end proc; P(p) # substitute a prime p here to run the code, it produces an answer (k) if one exists <=N and terms must be computed one at a time.
CROSSREFS
KEYWORD
sign,base
AUTHOR
STATUS
approved