cp's OEIS Frontend

A241022 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.

%I A241022 #6 Apr 16 2014 17:26:29
%S A241022 13,313,3313,31333,313333,3331333,31333333,333331333,3333133333,
%T A241022 33331333333,333313333333,0,33133333333333,333333313333333,
%U A241022 1333333333333333,33333333331333333,333333133333333333,3333313333333333333,33313333333333333333,333313333333333333333
%N A241022 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 A241022 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 A241022 Michel Lagneau, <a href="/A241022/b241022.txt">Table of n, a(n) for n = 2..150</a>
%p A241022 with(numtheory):nn:=80:T:=array(1..nn):
%p A241022    for n from 2 to nn do:
%p A241022      for i from 1 to n do:
%p A241022      T[i]:=3:
%p A241022      od:
%p A241022        ii:=0:
%p A241022        for j from 1 to n while(ii=0)do:
%p A241022        T[j]:=1:s:=sum('T[i]*10^(n-i)', 'i'=1..n):
%p A241022          if type(s,prime)=true
%p A241022          then
%p A241022          ii:=1: printf(`%d, `,s):
%p A241022          else
%p A241022          T[j]:=3:
%p A241022          fi:
%p A241022        od:
%p A241022           if ii=0
%p A241022            then
%p A241022            printf(`%d, `,0):
%p A241022            else
%p A241022           fi:
%p A241022      od:
%Y A241022 Cf. A241018, A241019, A241020, A241021, A241027.
%K A241022 nonn,base
%O A241022 2,1
%A A241022 _Michel Lagneau_, Apr 15 2014