login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Primes and their indices such that when their respective SOD's are both prime, the SOD of the index is the nextprime of the prime SOD.
5

%I #14 Mar 17 2021 23:44:16

%S 227,313,1723,2003,2311,2609,3329,3701,4007,4027,4801,4931,5107,7253,

%T 10457,12143,12163,12211,12433,12547,13063,14143,14341,14831,15139,

%U 15373,16091,17027,17047,19403,20047,20261,21059,21149,22157,23053,23293,23431,24229,24623

%N Primes and their indices such that when their respective SOD's are both prime, the SOD of the index is the nextprime of the prime SOD.

%C "SOD" = "sum of digits".

%F Individually sum the digits of both prime index and associated prime. If both sums are prime and the sum of the index is the next prime of the prime SOD, add to sequence.

%e 313 is a term, 65 is the index of the prime 313. The SOD(65) = 11. The SOD(313) = 7. 11 is the next prime after 7.

%t Select[Range[25000], PrimeQ[#] && PrimeQ[Plus @@ IntegerDigits[#]] &&

%t NextPrime[Plus @@ IntegerDigits[#]] == Plus @@ IntegerDigits[PrimePi[#]] &] (* _Tanya Khovanova_, Mar 16 2021 *)

%Y Cf. A000720, A007953, A117459, A117460, A117461, A117462, A117463.

%K easy,nonn,base

%O 1,1

%A _Enoch Haga_, Mar 18 2006

%E Offset changed to 1 by _Alois P. Heinz_, Mar 16 2021