OFFSET
1,2
COMMENTS
Also numbers k such that rad(k) = rad(omega(k)*bigomega(k)) where rad is A007947. [Charles R Greathouse IV, Jul 18 2011]
LINKS
Arkadiusz Wesolowski, Table of n, a(n) for n = 1..400
MATHEMATICA
Select[Range[12960], PowerMod[PrimeNu[#] PrimeOmega[#], #, #] == 0 &] (* T. D. Noe, Jul 14 2011 *)
PROG
(PARI) for(n=1, 2^14, if(lift(Mod((omega(n)*bigomega(n))^n, n))==0, print1(n, ", "))); \\ Arkadiusz Wesolowski, Jul 18 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
Arkadiusz Wesolowski, Jul 14 2011
STATUS
approved