OFFSET
1
LINKS
EXAMPLE
5, being a Pythagorean prime (A002144), when raised to any power, will always give a number of the form 4m+1, therefore a(1) = a(5) = a(25) = a(125) = 1.
7, although itself a prime of the form 4m+3 (A002145), when raised to an even power will always give a number of the form 4m+1, therefore a(49) = a(2401) = 1.
PROG
(PARI) A354030(n) = ((1==n)||((1==(n%4)) && isprimepower(n)));
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, May 15 2022
STATUS
approved