|
%I #13 Jun 20 2018 22:16:28
%S 13,31,97,163,181,283,409,499,709,787,811,877,1087,1399,1423,1609,
%T 1777,1801,1879,2347,2677,2719,3457,3517,3919,4273,4483,5701,6043,
%U 6121,6211,6481,6691,7573,8941,9733,9739,10069,10093,10159,10243,10789,11161,11251,11689,12799,12907
%N Primes p such that p+10, p+100 and p+1000 are all prime.
%C Number of terms < 10^k: 0, 3, 12, 37, 159, 789, 3960, 21708, 129910, ..., . - _Robert G. Wilson v_, Jun 20 2018
%H Robert Israel, <a href="/A271575/b271575.txt">Table of n, a(n) for n = 1..10000</a>
%e p=13; p+10=23 (is prime), p+100=113 (is prime), p+1000=1013 (is prime).
%p select(t -> isprime(t+1000) and isprime(t+100) and isprime(t+10) and isprime(t), [seq(i,i=7..20000, 6)]); # _Robert Israel_, Jun 20 2018
%t Select[Prime[Range[10000]], PrimeQ[# + 10] && PrimeQ[# + 100] && PrimeQ[# + 1000] &] (* _Robert Price_, Apr 10 2016 *)
%o (PARI) lista(nn) = forprime(p=2, nn, if (isprime(p+10) && isprime(p+100) && isprime(p+1000), print1(p, ", "))); \\ _Michel Marcus_, Apr 10 2016
%Y Cf. A000040, A002385, A015916, A023203.
%K nonn
%O 1,1
%A _Emre APARI_, Apr 10 2016
|
