|
%I #16 Sep 14 2022 07:27:06
%S 2,3,7,11,13,17,19,23,29,37,41,43,47,53,61,67,71,73,79,83,89,97,101,
%T 103,107,109,113,127,131,137,139,149,151,157,163,167,173,181,191,193,
%U 197,199,211,223,227,229,233,239,241,251,257,263,269,271,277,281,283
%N Indices of primes p whose order of primeness A078442(p) is prime.
%C All terms are prime.
%H Alois P. Heinz, <a href="/A333364/b333364.txt">Table of n, a(n) for n = 1..10000</a>
%H N. Fernandez, <a href="http://www.borve.org/primeness/FOP.html">An order of primeness, F(p)</a>
%H N. Fernandez, <a href="/A006450/a006450.html">An order of primeness</a> [cached copy, included with permission of the author]
%F { p in primes : A049076(p) is prime }.
%F a(n) = pi(A333353(n)), with pi = A000720.
%e 11 is a term: prime(11) = 31 -> 11 -> 5 -> 3 -> 2 -> 1, five (a prime number of) steps "->" = pi = A000720.
%p b:= proc(n) option remember;
%p `if`(isprime(n), 1+b(numtheory[pi](n)), 0)
%p end:
%p a:= proc(n) option remember; local p;
%p p:= `if`(n=1, 1, a(n-1));
%p do p:= nextprime(p);
%p if isprime(b(p)+1) then break fi
%p od; p
%p end:
%p seq(a(n), n=1..62);
%t b[n_] := b[n] = If[PrimeQ[n], 1 + b[PrimePi[n]], 0];
%t a[n_] := a[n] = Module[{p}, p = If[n == 1, 1, a[n - 1]];
%t While[True, p = NextPrime[p]; If[PrimeQ[b[p] + 1], Break[]]]; p];
%t Table[a[n], {n, 1, 62}] (* _Jean-François Alcover_, Sep 14 2022, after _Alois P. Heinz_ *)
%Y Cf. A000040, A000720, A049076, A333353.
%K nonn
%O 1,1
%A _Alois P. Heinz_, Mar 16 2020
|
