|
|
A007821
|
|
Primes p such that pi(p) is not prime.
|
|
80
|
|
|
2, 7, 13, 19, 23, 29, 37, 43, 47, 53, 61, 71, 73, 79, 89, 97, 101, 103, 107, 113, 131, 137, 139, 149, 151, 163, 167, 173, 181, 193, 197, 199, 223, 227, 229, 233, 239, 251, 257, 263, 269, 271, 281, 293, 307, 311, 313, 317, 337, 347, 349, 359, 373
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
The complement of A006450 (primes with prime index) within the primes A000040.
|
|
REFERENCES
|
C. Kimberling, Fractal sequences and interspersions, Ars Combinatoria, vol. 45 p 157 1997.
|
|
LINKS
|
|
|
FORMULA
|
G.f. over nonprime powers: Sum_{k >= 1} prime(k)*x^k-prime(prime(k))*x^prime(k). - Benedict W. J. Irwin, Jun 11 2016
|
|
MAPLE
|
|
|
MATHEMATICA
|
With[{nn=100}, Pick[Prime[Range[nn]], Table[If[PrimeQ[n], 0, 1], {n, nn}], 1]] (* Harvey P. Dale, Aug 14 2020 *)
|
|
PROG
|
(Haskell)
a007821 = a000040 . a018252
a007821_list = map a000040 a018252_list
(PARI) forprime(p=2, 1e3, if(!isprime(primepi(p)), print1(p, ", "))) \\ Felix Fröhlich, Aug 16 2014
|
|
CROSSREFS
|
Cf. A049076, A049078, A049079, A049080, A049081, A058322, A058324, A058325, A058326, A058327, A058328, A093046, A006450.
Let A = primes A000040, B = nonprimes A018252. The 2-level compounds are AA = A006450, AB = A007821, BA = A078782, BB = A102615. The 3-level compounds AAA, AAB, ..., BBB are A038580, A049078, A270792, A102617, A270794, A270795, A270796, A102616.
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|