OFFSET
1,2
COMMENTS
Numbers n such that the sum of exponents of primes == 1 (mod 4) in the prime factorization of n is not congruent to n mod 2. - Robert Israel, Jan 22 2017
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= proc(n) local F, t;
F:= select(t -> t[1] mod 4 = 1, ifactors(n)[2]);
(add(t[2], t=F) - n) mod 2 = 1;
end proc:
select(filter, [$1..1000]); # Robert Israel, Jan 22 2017
MATHEMATICA
Select[Range[150], Mod[DivisorSigma[1, #^2], 4]==1&] (* Harvey P. Dale, Apr 07 2012 *)
PROG
(PARI) isok(n) = (sigma(n^2) % 4) == 1; \\ Michel Marcus, Jan 22 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Apr 20 2003
STATUS
approved