OFFSET
1,2
PROG
(PARI) ispp1(n) = isprimepower(n) || (n==1); \\ A000961
lista(nn) = {my(s=1); for (n=1, nn, if (ispp1(n), s*= n; print1(s, ", ")); ); } \\ Michel Marcus, Mar 26 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Den Roussel (DenRoussel(AT)webtv.net)
EXTENSIONS
Offset 1 and more terms from Michel Marcus, Mar 26 2020
STATUS
approved