OFFSET
1,1
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..1000
EXAMPLE
12*rad(12)+1 = 12*rad(3*2^2)+1 = 12*3*2+1 = 72+1 = 73, therefore 73 is a term.
a(33) = 10369 = 10368 + 1: A078310(1728) = (2*3)*(2^6*3^3) = 10368.
MATHEMATICA
powQ[n_] := n == 1 || AllTrue[FactorInteger[n][[;; , 2]], # > 1 &]; Select[Prime[Range[2400]], powQ[# - 1] &] (* Amiram Eldar, Jul 31 2022 *)
PROG
(Haskell)
a078324 n = a078324_list !! (n-1)
a078324_list = filter ((== 1) . a010051') a224866_list
-- Reinhard Zumkeller, Jul 23 2013
(PARI) is(n) = isprime(n) && ispowerful(n-1); \\ Amiram Eldar, Jul 31 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Nov 23 2002
EXTENSIONS
Missing terms 10369, 16001, 17497 and 18433 inserted by Reinhard Zumkeller, Jul 23 2013
STATUS
approved