|
%I #44 Sep 08 2022 08:45:39
%S 5,11,13,19,23,29,41,43,53,61,71,79,83,89,109,113,131,139,149,163,173,
%T 179,181,193,211,223,233,239,251,263,293,313,331,349,379,383,401,421,
%U 431,439,443,449,461,491,503,523,569,599,601,613,641,643,659,673,683
%N Primes p such that p+18 is also prime.
%C Both p and p+18 have the same digital root (A010888). - _Zak Seidov_, Sep 14 2015
%C No term belongs to A030432. - _Michel Marcus_, Sep 14 2015
%C No term belongs to A045437. - _Bruno Berselli_, Sep 14 2015
%H Vincenzo Librandi, <a href="/A153418/b153418.txt">Table of n, a(n) for n = 1..1000</a>
%e 5 is in sequence because 5+18=23 is also prime;
%e 11 is in sequence because 11+18=29 is also prime.
%t lst={};d=18;Do[p=Prime[n];If[PrimeQ[p+d],AppendTo[lst,p]],{n,6!}];lst
%t Select[Prime[Range[150]], PrimeQ[(# + 18)]&] (* _Vincenzo Librandi_, Apr 14 2013 *)
%o (Magma) [p: p in PrimesUpTo(1000) | IsPrime(p+18)]; // _Vincenzo Librandi_, Apr 14 2013
%o (PARI) list(n)=forprime(p=1,n,if(isprime(p+18),print1(p", "))) \\ _Anders Hellström_, Sep 13 2015
%o (Sage) [p for p in primes(700) if is_prime(p+18)]; # _Bruno Berselli_, Sep 14 2015
%Y Cf. A007953, A010888, A045437, A049488, A030432, A153417.
%Y A031936 is a subsequence. - _Zak Seidov_, Sep 13 2015
%K nonn,easy
%O 1,1
%A _Vladimir Joseph Stephan Orlovsky_, Dec 25 2008
%E Definition improved by _Bruno Berselli_, Oct 31 2012
|
