OFFSET
1,2
COMMENTS
Closed under multiplication.
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..10000
EXAMPLE
Neither 5 or 10 (= 2*5) are included, because the prime factor 5 (of the form 4k+1) occurs just once.
6 = 2*3 is present, as there are no prime factors of 4k+1 present at all, and zero is an even number.
Also 25 (5*5) and 50 (2*5*5) and 75 (3*5*5) and 625 (5*5*5*5) are included, because in all of them, the prime factor 5 (of the form 4k+1) occurs an even number of times.
MATHEMATICA
{1}~Join~Select[Range@ 140, NoneTrue[FactorInteger@ #, And[Mod[First@ #, 4] == 1, OddQ@ Last@ #] &] &] (* Michael De Vlieger, Feb 04 2016, Version 10 *)
PROG
(PARI) isok(n) = {my(f = factor(n)); for (i=1, #f~, if (((f[i, 1] % 4) == 1) && (f[i, 2] % 2), return (0)); ); return (1); } \\ Michel Marcus, Feb 04 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Feb 03 2016
STATUS
approved