A328558 Squares (A000290) that are not arithmetic numbers (A003601).
4, 9, 16, 25, 36, 64, 81, 100, 121, 144, 196, 225, 256, 289, 324, 400, 441, 484, 529, 576, 625, 676, 729, 784, 841, 900, 1024, 1089, 1156, 1225, 1296, 1444, 1521, 1600, 1681, 1764, 1936, 2025, 2116, 2209, 2304, 2401, 2500, 2601, 2704, 2809, 2916, 3025, 3136
Offset: 1
Keywords
Programs
-
Magma
[m: m in [1..10^5] | not IsIntegral(SumOfDivisors(m) / NumberOfDivisors(m)) and IsIntegral(Sqrt(m))]
-
Mathematica
Select[Range[56]^2, !Divisible[DivisorSigma[1, #], DivisorSigma[0, #]] &] (* Amiram Eldar, Oct 20 2019 *)
-
PARI
isA328558(n) = (issquare(n)&&(sigma(n)%numdiv(n))); \\ Antti Karttunen, Oct 19 2019
Comments