OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
211 is in this sequence because 4 * 211 + 1 = 845 = 5 * 13^2, 4 * 211 + 2 = 846 = 2 * 3^2 * 47 and 4 * 211 + 3 = 847 = 7 * 11^2.
MAPLE
remove(t->ormap(numtheory:-issqrfree, [4*t+1, 4*t+2, 4*t+3]), [$1..2*10^4]); # Robert Israel, Apr 03 2018
MATHEMATICA
Select[Range[1000], Union[{MoebiusMu[4# + 1], MoebiusMu[4# + 2], MoebiusMu[4# + 3]}] == {0} &] (* Alonso del Arte, May 26 2015 *)
PROG
(Magma) [n: n in [1..20000] | not IsSquarefree(4*n+1) and not IsSquarefree(4*n+2) and not IsSquarefree(4*n+3)];
(PARI) isok(n) = !issquarefree(4*n+1) && !issquarefree(4*n+2) && !issquarefree(4*n+3); \\ Michel Marcus, Apr 04 2018
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Juri-Stepan Gerasimov, May 26 2015
STATUS
approved