OFFSET
1,1
COMMENTS
All terms of A107978 are in the sequence.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
Select[4 Range@ 150 + 1, Mod[#, 4] == 3 &[FactorInteger[#][[-1, 1]]] &] (* Michael De Vlieger, Feb 07 2016 *)
PROG
(PARI) isok(n) = (n != 1) && ((n % 4) == 1) && ((vecmax(factor(n)[, 1]) % 4) == 3); \\ Michel Marcus, Feb 07 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Shevelev, Apr 13 2010
EXTENSIONS
Corrected and extended by Michel Marcus, Feb 07 2016
STATUS
approved