%I #17 Mar 11 2022 07:33:49
%S 2999,3989,4799,4889,5879,5897,5987,6599,6689,6779,6869,6959,6977,
%T 7499,7589,7877,7949,8597,8669,8849,8867,9479,9497,9587,9677,9749,
%U 9767,9839,9857,9929,12899,13799,13997,14699,14879,14897,14969,15797,15887,15959
%N Prime quartet leaders: largest number of a prime quartet.
%C A prime quartet is a set of four different prime numbers such that the fourth number is a 1-digit number which is the sum of the digits of the third number, the third number is the sum of the digits of the second number and the second number is the sum of the digits of the first number.
%C Different from A106766.
%C Comment from _Joshua Zucker_, Apr 24 2007, on the difference between this sequence and A106766: The digit sum must be the largest member of a prime trio, so the first number where the sequences differ must be with digit sum 47 and thus have at least 6 digits - so until then you get all the primes with 4 or 5 digits that have digit sum 29.
%C a(2322)=389999 is the first value different from A106766, where A106766(2322)=390359. See also A106778 = primes with digit sum = 47: A106778(1)=389999. - _Martin Fuller_ and _Ray Chandler_, Apr 24 2007
%C The sequence of prime quintet leaders is probably too large for the OEIS; its first term is the 334-digit prime 5*10^333-10^330-10^328-1 with sum of digits a(1) = 2999. - _Charles R Greathouse IV_, Mar 11 2022
%H Harvey P. Dale, <a href="/A119892/b119892.txt">Table of n, a(n) for n = 1..2000</a>
%H L. Stevens, <a href="http://www.lucstevens.com/Primeensembles.htm">Prime ensembles</a>
%e 2999 is in the sequence because it is the largest number of the prime quartet (2999,29,11,2).
%t pqQ[n_]:=Module[{p1=NestList[Total[IntegerDigits[#]]&,n,3]},AllTrue[ Take[ p1,3],#>9&]&&AllTrue[p1,PrimeQ]]; Select[Range[16000],pqQ] (* Requires Mathematica version 10 or later *) (* _Harvey P. Dale_, Oct 02 2020 *)
%o (PARI) DigitSum(n,b=10)=local(x);x=0;while(n,x+=n%b;n\=b);x
%o PrimeEnsemble(n,b=10)=local(x);x=1;while(ispseudoprime(n),if(n<b,return(x));n=DigitSum(n,b);x++);0
%o forprime(p=2, 16000, if(PrimeEnsemble(p)>=4, print1(p", "))); \\ _Martin Fuller_
%Y Cf. A119889, A119890, A119891.
%K base,nonn
%O 1,1
%A Luc Stevens (lms022(AT)yahoo.com), May 27 2006
|