|
%I #18 Oct 30 2014 11:51:53
%S 2,11,101,1013,10037,100019,1000033,10000019,100000037,1000000033,
%T 10000000019,100000000019,1000000000039,10000000000037,
%U 100000000000031,1000000000000037,10000000000000079,100000000000000013,1000000000000000031,10000000000000000051,100000000000000000039
%N Least n-digit prime whose digit sum is also prime.
%C Subsequence of A046704 (primes with digits sum being prime).
%C Some terms of this sequence are also in A003617, the least n-digit primes. - _Michel Marcus_, Oct 30 2014
%H Paolo P. Lava, <a href="/A249447/b249447.txt">Table of n, a(n) for n = 1..100</a>
%e a(1) = 2 because it is the least prime with just one digit.
%e a(2) = 11 because it is the least prime with 2 digits whose sum, 1 + 1 = 2, is a prime.
%e Again, a(7) = 1000033 because it is the least prime with 7 digits whose sum is a prime: 1 + 0 + 0 + 0 + 0 + 3 + 3 = 7.
%p P:=proc(q) local a,b,k,n; for k from 0 to q do
%p for n from 10^k to 10^(k+1)-1 do if isprime(n) then a:=n; b:=0;
%p while a>0 do b:=b+(a mod 10); a:=trunc(a/10); od;
%p if isprime(b) then print(n); break; fi; fi;
%p od; od; end: P(10^3);
%o (PARI) a(n) = {p = nextprime(10^(n-1)); while (!isprime(sumdigits(p)), p = nextprime(p+1)); p;} \\ _Michel Marcus_, Oct 29 2014
%Y Cf. A003617, A046704, A068166, A249448.
%K nonn,base
%O 1,1
%A _Paolo P. Lava_, Oct 29 2014
| 