login
Smallest prime numbers p of length n having a decimal expansion x(1)x(2)... x(n) such that there exists an index j where x(j) = 1 and x(i) = 3 for i<>j, or 0 if no such prime exists.
3

%I #6 Apr 16 2014 17:26:29

%S 13,313,3313,31333,313333,3331333,31333333,333331333,3333133333,

%T 33331333333,333313333333,0,33133333333333,333333313333333,

%U 1333333333333333,33333333331333333,333333133333333333,3333313333333333333,33313333333333333333,333313333333333333333

%N Smallest prime numbers p of length n having a decimal expansion x(1)x(2)... x(n) such that there exists an index j where x(j) = 1 and x(i) = 3 for i<>j, or 0 if no such prime exists.

%C The corresponding index of the decimal digit 1 are 1, 2, 3, 2, 2, 4, 2, 6, 5, 5, 5, 0, 3, 8, 1, 11, 7, 6, 4, 5,...(A241019).

%H Michel Lagneau, <a href="/A241022/b241022.txt">Table of n, a(n) for n = 2..150</a>

%p with(numtheory):nn:=80:T:=array(1..nn):

%p for n from 2 to nn do:

%p for i from 1 to n do:

%p T[i]:=3:

%p od:

%p ii:=0:

%p for j from 1 to n while(ii=0)do:

%p T[j]:=1:s:=sum('T[i]*10^(n-i)', 'i'=1..n):

%p if type(s,prime)=true

%p then

%p ii:=1: printf(`%d, `,s):

%p else

%p T[j]:=3:

%p fi:

%p od:

%p if ii=0

%p then

%p printf(`%d, `,0):

%p else

%p fi:

%p od:

%Y Cf. A241018, A241019, A241020, A241021, A241027.

%K nonn,base

%O 2,1

%A _Michel Lagneau_, Apr 15 2014