cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A344948 Primes whose position among the powers of primes (A000961) is also prime.

This page as a plain text file.
%I A344948 #11 Jun 06 2021 10:25:38
%S A344948 2,3,5,19,29,47,67,73,101,113,137,167,193,199,239,263,313,349,389,419,
%T A344948 431,449,461,487,571,599,641,701,719,751,797,823,857,887,911,977,991,
%U A344948 1019,1097,1193,1223,1231,1277,1301,1307,1399,1439,1481,1511,1531,1571,1601
%N A344948 Primes whose position among the powers of primes (A000961) is also prime.
%H A344948 Michel Marcus, <a href="/A344948/b344948.txt">Table of n, a(n) for n = 1..10000</a>
%e A344948 The position of 2 in A000961 is 2; so 2 is a term.
%e A344948 The position of 19 in A000961 is 13; so 19 is a term.
%t A344948 pow = Select[Range[1600], # == 1 || PrimePowerQ[#] &]; Select[pow[[Select[Range @ Length[pow], PrimeQ]]], PrimeQ] (* _Amiram Eldar_, Jun 03 2021 *)
%o A344948 (PARI) allmps(nn) = {my(map = Map()); mapput(map, 1, 1); my(nb=1); for (n=2, nn, if (isprimepower(n), nb++; mapput(map, n, nb));); map;}
%o A344948 lista(nn) = {my(nb = prime(nn), map = allmps(nb)); forprime (p=1, nn, if( isprime(mapget(map, p)), print1(p, ", ")););}
%Y A344948 Cf. A000040, A000961, A024620.
%K A344948 nonn
%O A344948 1,1
%A A344948 _Michel Marcus_, Jun 03 2021