OFFSET
1,1
LINKS
Michel Marcus, Table of n, a(n) for n = 1..10000
EXAMPLE
MATHEMATICA
pow = Select[Range[1600], # == 1 || PrimePowerQ[#] &]; Select[pow[[Select[Range @ Length[pow], PrimeQ]]], PrimeQ] (* Amiram Eldar, Jun 03 2021 *)
PROG
(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; }
lista(nn) = {my(nb = prime(nn), map = allmps(nb)); forprime (p=1, nn, if( isprime(mapget(map, p)), print1(p, ", ")); ); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Jun 03 2021
STATUS
approved