OFFSET
1,1
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
remove(k-> issqr(add(numtheory[pi](i[1])*i[2], i=ifactors(k)[2])), [$1..99])[]; # Alois P. Heinz, Jan 22 2019
MATHEMATICA
Select[Range[100], !IntegerQ[Sqrt[Sum[PrimePi[f[[1]]]*f[[2]], {f, FactorInteger[#]}]]]&]
PROG
(PARI) isok(n) = {my(f=factor(n)); !issquare(sum(k=1, #f~, primepi(f[k, 1])*f[k, 2])); } \\ Michel Marcus, Jan 18 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Jan 17 2019
STATUS
approved
