A159236 -id:A159236 - cp's OEIS frontend

A306920 a(n) is the smallest prime > 10 where a string of exactly n zeros can be inserted somewhere into the decimal expansion such that the resulting number is also prime.

11, 19, 17, 13, 13, 23, 17, 17, 31, 13, 23, 41, 137, 61, 23, 13, 13, 67, 53, 89, 19, 107, 17, 29, 61, 263, 31, 37, 127, 53, 269, 199, 137, 23, 31, 89, 61, 13, 43, 163, 53, 131, 109, 19, 79, 283, 109, 19, 269, 223, 97, 97, 223, 89, 13, 79, 67, 107, 17, 389, 197

Offset: 1

Felix Fröhlich, Mar 16 2019

For many small n, if the decimal expansion of a(n) contains the digit 0, then a(n+1) is a(n) with one zero digit removed. However, this is not true in general. The counterexamples' indices in this sequence are given by A344860.

			For n = 13: If a string of 13 zeros is inserted between the digits 1 and 3 in 137, the resulting number is 1000000000000037, which is prime. Since 137 is the smallest prime where such a string of 13 zeros can be inserted to get another prime, a(13) = 137.

Felix Fröhlich, Table of n, a(n) for n = 1..2300

Column 1 of A306926.

Cf. A159236, A164329, A215417, A290174, A344860.

insert(n, len, pos) = my(d=digits(n), v=[], w=[]); for(y=1, pos, v=concat(v, d[y])); v=concat(v, vector(len)); for(z=pos+1, #d, v=concat(v, d[z])); subst(Pol(v), x, 10)
a(n) = forprime(p=10, , for(k=1, #digits(p)-1, my(zins=insert(p, n, k)); if(ispseudoprime(zins), return(p))))

A217045 Primes that remain prime when a single "4" digit is inserted between any two adjacent decimal digits.

Original entry on oeis.org

19, 37, 43, 61, 67, 73, 97, 109, 199, 211, 223, 241, 349, 409, 421, 457, 463, 541, 571, 751, 757, 823, 991, 1033, 1087, 1321, 1423, 1447, 1543, 2749, 3361, 3469, 3499, 3847, 4111, 4273, 4483, 5059, 5437, 5443, 5449, 6373, 6709, 6793, 7687, 8089, 8221, 8443

Offset: 1

Paolo P. Lava, Sep 25 2012

			87697 is prime and also 876947, 876497, 874697 and 847697.

Paolo P. Lava, Table of n, a(n) for n = 1..141

Cf. A050674, A050711-A050719, A069246, A159236, A215417, A215419-A215421, A217044, A217046, A217047, A217062-A217065

with(numtheory);
A217045:=proc(q,x)
local a,b,c,i,n,ok;
for n from 5 to q do
a:=ithprime(n); b:=0;
while a>0 do b:=b+1; a:=trunc(a/10); od; a:=ithprime(n); ok:=1;
  for i from 1 to b-1 do
    c:=a+9*10^i*trunc(a/10^i)+10^i*x;
    if not isprime(c) then ok:=0; break; fi; od;
  if ok=1 then print(ithprime(n)); fi;
od; end:
A217045(100000,4)

Select[Prime[Range[5,1500]],AllTrue[Table[FromDigits[Insert[ IntegerDigits[ #],4,n]],{n,2,IntegerLength[#]}],PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Dec 04 2017 *)

is(n)=my(v=concat([""], digits(n))); for(i=2, #v-1, v[1]=Str(v[1], v[i]); v[i]=4; if(i>2, v[i-1]=""); if(!isprime(eval(concat(v))), return(0))); isprime(n) \\ Charles R Greathouse IV, Sep 26 2012

A217046 Primes that remain prime when a single "6" digit is inserted between any two adjacent decimal digits.

Original entry on oeis.org

13, 17, 23, 29, 37, 41, 43, 47, 53, 59, 61, 71, 79, 83, 97, 101, 109, 113, 137, 157, 163, 167, 263, 277, 293, 307, 313, 317, 331, 397, 421, 443, 457, 463, 569, 607, 653, 659, 661, 673, 691, 739, 769, 787, 809, 823, 829, 863, 881, 977, 997, 1063, 1087, 1453

Offset: 1

Paolo P. Lava, Sep 25 2012

			185917 is prime and also 1859167, 1859617, 1856917, 1865917 and 1685917.

Harvey P. Dale, Table of n, a(n) for n = 1..500 (First 262 terms from Paolo P. Lava)

Cf. A050674, A050711-A050719, A069246, A159236, A215417, A215419-A215421, A217044, A217045, A217047, A217062-A217065

with(numtheory);
A217044:=proc(q,x)
local a,b,c,i,n,ok;
for n from 5 to q do
a:=ithprime(n); b:=0;
while a>0 do b:=b+1; a:=trunc(a/10); od; a:=ithprime(n); ok:=1;
  for i from 1 to b-1 do
    c:=a+9*10^i*trunc(a/10^i)+10^i*x;
    if not isprime(c) then ok:=0; break; fi; od;
  if ok=1 then print(ithprime(n)); fi;
od; end:
A217044(100000,6)

Select[Prime[Range[5,1200]],And@@PrimeQ[FromDigits/@Table[ Insert[ IntegerDigits[ #],6,i],{i,2,IntegerLength[#]}]]&] (* Harvey P. Dale, Oct 09 2012 *)

is(n)=my(v=concat([""], digits(n))); for(i=2, #v-1, v[1]=Str(v[1], v[i]); v[i]=6; if(i>2, v[i-1]=""); if(!isprime(eval(concat(v))), return(0))); isprime(n) \\ Charles R Greathouse IV, Sep 26 2012

A217063 Primes that remain prime when a single "3" digit is inserted between any two adjacent decimal digits.

Original entry on oeis.org

11, 17, 19, 23, 29, 31, 37, 41, 43, 61, 73, 79, 89, 97, 101, 103, 127, 167, 173, 181, 211, 233, 239, 251, 271, 283, 307, 331, 359, 373, 439, 491, 509, 523, 547, 599, 673, 709, 733, 769, 877, 887, 937, 941, 991, 1033, 1229, 1381, 1619, 1721, 1759, 1789, 1901

Offset: 1

Paolo P. Lava, Sep 26 2012

			212881 is prime and also 2128831, 2128381, 2123881, 213288 and 2312881.

Bruno Berselli, Table of n, a(n) for n = 1..1000 (first 275 terms from Paolo Lava)

Cf. A050674, A050711-A050719, A069246, A159236, A215417, A215419-A215421, A217044-A217047, A217062-A217065.

[p: p in PrimesInInterval(11, 2000) | forall{m: t in [1..#Intseq(p)-1] | IsPrime(m) where m is (Floor(p/10^t)*10+3)*10^t+p mod 10^t}]; // Bruno Berselli, Sep 26 2012

with(numtheory);
A217063:=proc(q,x)
local a,b,c,i,n,ok;
for n from 5 to q do
  a:=ithprime(n); b:=0; while a>0 do b:=b+1; a:=trunc(a/10); od; a:=ithprime(n); ok:=1;
    for i from 1 to b-1 do
      c:=a+9*10^i*trunc(a/10^i)+10^i*x;  if not isprime(c) then ok:=0; break; fi; od;
    if ok=1 then print(ithprime(n)); fi; od; end:
A217063(1000000,3);

is(n)=my(v=concat([""], digits(n))); for(i=2, #v-1, v[1]=Str(v[1], v[i]); v[i]=3; if(i>2, v[i-1]=""); if(!isprime(eval(concat(v))), return(0))); isprime(n) \\ Charles R Greathouse IV, Sep 26 2012

from sympy import isprime, primerange
def ok(p):
    if p < 10: return False
    s = str(p)
    return all(isprime(int(s[:i] + "3" + s[i:])) for i in range(1, len(s)))
def aupto(limit): return [p for p in primerange(1, limit+1) if ok(p)]
print(aupto(1901)) # Michael S. Branicky, Nov 17 2021

A217064 Primes that remain prime when a single "5" digit is inserted between any two adjacent decimal digits.

Original entry on oeis.org

11, 17, 47, 71, 83, 89, 149, 167, 179, 251, 257, 293, 347, 359, 383, 419, 461, 467, 491, 557, 563, 569, 653, 773, 911, 1193, 1217, 1277, 1451, 1559, 1667, 1823, 1901, 2243, 2309, 2357, 2579, 2657, 2999, 3527, 3533, 4289, 5051, 5351, 5501, 5843, 6089, 6551, 6581

Offset: 1

Paolo P. Lava, Sep 26 2012

			290183 is prime and also 2901853, 2901583, 2905183, 2950183 and 2590183.

Harvey P. Dale, Table of n, a(n) for n = 1..500 (First 140 terms from Paolo P. Lava)

Cf. A050674, A050711-A050719, A069246, A159236, A215417, A215419-A215421, A217044-A217047, A217062, A217063, A217065

with(numtheory);
A217064:=proc(q,x)
local a,b,c,i,n,ok;
for n from 5 to q do
  a:=ithprime(n); b:=0; while a>0 do b:=b+1; a:=trunc(a/10); od; a:=ithprime(n); ok:=1;
    for i from 1 to b-1 do
      c:=a+9*10^i*trunc(a/10^i)+10^i*x;  if not isprime(c) then ok:=0; break; fi; od;
    if ok=1 then print(ithprime(n)); fi; od; end:
A217064(1000000,5);

Select[Prime[Range[5,1000]],AllTrue[FromDigits/@Table[ Insert[ IntegerDigits[ #],5,n],{n,2,IntegerLength[#]}],PrimeQ]&] (* Harvey P. Dale, Feb 20 2022 *)

is(n)=my(v=concat([""], digits(n))); for(i=2, #v-1, v[1]=Str(v[1], v[i]); v[i]=5; if(i>2, v[i-1]=""); if(!isprime(eval(concat(v))), return(0))); isprime(n) \\ Charles R Greathouse IV, Sep 26 2012

A304246 Numbers that yield a prime whenever a '1' is inserted between any two digits.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 13, 21, 31, 33, 37, 49, 63, 67, 69, 79, 81, 91, 99, 103, 109, 117, 123, 151, 163, 181, 193, 211, 213, 231, 241, 279, 309, 319, 363, 367, 391, 411, 427, 429, 453, 457, 459, 501, 513, 519, 547, 571, 601, 613, 621, 631, 697, 703, 709, 721, 729, 777, 787, 801, 811, 817, 879, 903, 951, 981, 987

Offset: 1

M. F. Hasler, Jun 01 2018

The single-digit terms voidly satisfy the condition: no '1' can be inserted anywhere, so all possible insertions yield a prime.
Motivated by sequence A164329 which is the analog for inserting 0.
Compare to A068673 where 1 is prefixed or appended, and to A068679 where 1 is prefixed, appended or inserted anywhere - which is therefore the intersection between this sequence and A068673.
See also A050711 where 1 is inserted between all adjacent digits. - R. J. Mathar, Feb 28 2020

			21 is in the sequence, because if '1' is inserted between "any" pair consecutive digits (the only possibility being to insert it between the first and second digit, which yields 211), the result is always prime. The definition does not require the term itself to be prime.
103 is in the sequence because inserting 1 between the first and second, or between the second and third digit, would yield 1103 or 1013, respectively, which are both prime.

Cf. A164329 (prime when 0 is inserted anywhere), A216169 (subset of composite terms), A215417 (subset of primes), A159236 (prime when 0 is inserted between all digits).
Cf. A068679 (1 is prefixed, appended or inserted anywhere), A069246 (primes among these), A068673 (1 is prefixed, or appended).
Cf. A158594 (3 is prefixed, appended or inserted anywhere), A215419 (primes among these).
Cf. A069832 (7 is prefixed, appended or inserted anywhere), A215420 (primes among these), A068677 (7 is prefixed or appended).
Cf. A069833 (9 is prefixed, appended or inserted anywhere), A215421 (primes among these).
Cf. A158232 (13 is prefixed or appended).
Cf. A304243 (2 is prefixed or prime(k+2) is inserted after the k-th digit), A304244 (prime(k) is inserted after the k-th digit), A304245 (prime(k+1) is inserted after the k-th digit, k > 1, or '2' after the first digit).

[0] cat [k:k in [1..1000]| forall{i:i in [1..#Intseq(k)-1]| IsPrime(Seqint(Reverse(v[1..i] cat [1] cat v[i+1..#v]))) where v is Reverse(Intseq(k)) }]; // Marius A. Burtea, Feb 09 2020

filter:= proc(n) local j,t;
  for j from 1 to ilog10(n) do
    if not isprime(10*n-9*(n mod 10^j)+10^j) then return false fi
  od;
  true
end proc:
select(filter, [$0..1000]); # Robert Israel, Jun 01 2018

is(n)=!for(k=1,logint(n+!n,10),isprime(10*n-n%10^k*9+10^k)||return)

A304247 Numbers which yield a prime whenever a '2' is inserted between any single pair of adjacent digits.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 17, 23, 27, 29, 41, 51, 53, 77, 81, 83, 87, 89, 99, 101, 113, 123, 129, 131, 137, 149, 183, 207, 221, 243, 251, 297, 303, 321, 329, 357, 359, 399, 401, 417, 419, 429, 441, 443, 453, 461, 471, 473, 527, 533, 581, 597, 611, 621

Offset: 1

M. F. Hasler, Jun 01 2018

Motivated by existing sequences defined in an analog way for other digits to be inserted, e.g., A164329 for the digit 0, cf. cross-references.
For single-digit terms, the condition is voidly satisfied: nothing can be inserted.
See also A050712 where 2 is inserted between each pair of adjacent digits. - R. J. Mathar, Feb 28 2020

			123 is in the sequence because it yields a prime when a '2' is inserted after the first or after the second digit, which yields the prime 1223 in both cases. The term itself does not need to be prime.

Chai Wah Wu, Table of n, a(n) for n = 1..10000

Cf. A164329 (prime when 0 is inserted anywhere), A216169 (subset of composite terms), A215417 (subset of primes), A159236 (0 is inserted between all digits).
Cf. A068679 (1 is prefixed, appended or inserted anywhere), A069246 (primes among these), A068673 (1 is prefixed, or appended), A304246 (1 is inserted anywhere).
Cf. A158594 (3 is prefixed, appended or inserted anywhere), A215419 (primes among these).
Cf. A069832 (7 is prefixed, appended or inserted anywhere), A215420 (primes among these), A068677 (7 is prefixed or appended).
Cf. A069833 (9 is prefixed, appended or inserted anywhere), A215421 (primes among these).
Cf. A158232 (13 is prefixed or appended).
Cf. A304243 (2 is prefixed or prime(k+2) is inserted after the k-th digit), A304244 (prime(k) is inserted after the k-th digit), A304245 (prime(k+1) is inserted after the k-th digit, k > 1, or '2' after the first digit).

filter:= proc(n) local j,t;
  for j from 1 to ilog10(n) do
    if not isprime(10*n-9*(n mod 10^j)+2*10^j) then return false fi
  od;
  true
end proc:
select(filter, [$0..10000]); # Robert Israel, Jun 01 2018

is(n,p=2,L=logint(n+!n,10)+1,d,P)=!for(k=1,L-1,isprime((d=divrem(n,P=10^(L-k)))[2]+(10*d[1]+p)*P)||return)

A304248 Numbers that yield a prime whenever a '3' is inserted between any two digits.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 11, 17, 19, 23, 29, 31, 37, 41, 43, 49, 61, 73, 79, 89, 97, 101, 103, 121, 127, 167, 173, 181, 209, 211, 233, 239, 247, 251, 271, 283, 299, 307, 331, 343, 359, 361, 373, 391, 437, 439, 473, 491, 497, 509, 523, 533, 547, 551, 599

Offset: 1

M. F. Hasler, Jun 01 2018

Motivated by existing sequences defined in a similar way for other digits (e.g., A164329 for digit 0), subsequence A158594 = intersection of this and A068674 ('3' is prefixed or appended), and others: cf. cross-references.

			121 is in the sequence because it yields a prime when a digit 3 is inserted after the first or after the second digit, which yields the prime 1321 or 1231, respectively. The term itself does not need to be prime.
The single-digit numbers 0..9 are in the sequence because they satisfy the condition voidly: nothing can be inserted, so no insertion yields a nonprime, so all possible insertions always yield a prime.

Chai Wah Wu, Table of n, a(n) for n = 1..10000

Cf. A164329 (prime when 0 is inserted anywhere), A216169 (subset of composite terms), A215417 (subset of primes), A159236 (prime when 0 is inserted between all digits).
Cf. A068679 (1 is prefixed, appended or inserted anywhere), A069246 (primes among these), A068673 (1 is prefixed, or appended), A304246 (1 is inserted anywhere).
Cf. A304247 (2 is inserted anywhere).
Cf. A158594 (3 is prefixed, appended or inserted anywhere), A215419 (primes among these), A068674 (3 is prefixed or appended).
Cf. A069832 (7 is prefixed, appended or inserted anywhere), A215420 (primes among these), A068677 (7 is prefixed or appended).
Cf. A069833 (9 is prefixed, appended or inserted anywhere), A215421 (primes among these).
Cf. A158232 (13 is prefixed or appended).
Cf. A304243 (2 is prefixed or prime(k+2) is inserted after the k-th digit), A304244 (prime(k) is inserted after the k-th digit), A304245 (prime(k+1) is inserted after the k-th digit, k > 1, or '2' after the first digit).

[0] cat [k:k in [1..600]| forall{i:i in [1..#Intseq(k)-1]| IsPrime(Seqint(Reverse(v[1..i] cat [3] cat v[i+1..#v]))) where v is Reverse(Intseq(k))}]; // Marius A. Burtea, Feb 09 2020

Select[Range[0,600],AllTrue[FromDigits/@Table[Insert[IntegerDigits[#],3,n],{n,2,IntegerLength[ #]}],PrimeQ]&] (* Harvey P. Dale, Nov 06 2022 *)

is(n, p=3, L=logint(n+!n, 10)+1, d, P)=!for(k=1, L-1, isprime((d=divrem(n, P=10^(L-k)))[2]+(10*d[1]+p)*P)||return)

A119680 Prime numbers obtained by inserting a 0 between each pair of adjacent digits of a prime number > 10.

Original entry on oeis.org

101, 103, 107, 109, 307, 401, 503, 509, 601, 607, 701, 709, 809, 907, 10007, 10009, 10103, 10301, 10501, 10607, 10709, 10903, 10909, 20101, 20507, 20707, 20903, 30103, 30307, 30509, 30703, 30803, 30809, 40009, 40507, 40709, 50707, 50909, 60103, 60107, 60509

Offset: 1

Roger L. Bagula, Jun 11 2006

From Rémy Sigrist, Oct 08 2017: (Start)
See A159236 for the original prime numbers.
The least prime numbers > 10 remaining prime during exactly k iterations of the operation of inserting a 0 between each pair of adjacent digits are, for small values of k:
k prime
- -----
0 23
1 11
2 19
3 17
4 220333
5 8677267
(End)

			The first four terms arise from 11 -> 101, 13 -> 103, 17 -> 107, 19 -> 109.
23 -> 203 is not prime, so 203 is not a term.

Cf. A159236.

a = Table[Table[Mod[Floor[Prime[m]/10^n], 10], {n, 4, 0, -1}], {m, 5, 200}]; Dimensions[a] b = Table[Sum[(If[Mod[n - 1, 2] == 0, a[[m, 1 + Floor[(n - 1)/2]]], 0])*10^(9 - n), {n, 1, 9}], {m, 1, 195}]; c = Flatten[Table[If[PrimeQ[b[[m]]], b[[m]], {}], {m, 1, 195}]]

forprime (p=10, 599, if (isprime(pp=fromdigits(digits(p), 100)), print1 (pp ", "))) \\ Rémy Sigrist, Oct 08 2017

from itertools import count, islice
from sympy import isprime, nextprime
def ok(n):
    return n > 10 and isprime(n) and isprime(int("0".join(list(str(n)))))
def agen():
    p = 11
    while True:
        t = int("0".join(list(str(p))))
        if isprime(t): yield t
        p = nextprime(p)
print(list(islice(agen(), 50))) # Michael S. Branicky, Jul 11 2022

Name edited by Rémy Sigrist, Oct 08 2017

a(39)-a(41) from Michael S. Branicky, Jul 11 2022

cp's OEIS Frontend

A306920 a(n) is the smallest prime > 10 where a string of exactly n zeros can be inserted somewhere into the decimal expansion such that the resulting number is also prime.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

A217045 Primes that remain prime when a single "4" digit is inserted between any two adjacent decimal digits.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

A217046 Primes that remain prime when a single "6" digit is inserted between any two adjacent decimal digits.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

A217063 Primes that remain prime when a single "3" digit is inserted between any two adjacent decimal digits.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Magma

Maple

PARI

Python

A217064 Primes that remain prime when a single "5" digit is inserted between any two adjacent decimal digits.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

A304246 Numbers that yield a prime whenever a '1' is inserted between any two digits.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Magma

Maple

PARI

A304247 Numbers which yield a prime whenever a '2' is inserted between any single pair of adjacent digits.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

PARI

A304248 Numbers that yield a prime whenever a '3' is inserted between any two digits.

Views