OFFSET
1,1
COMMENTS
n is in the sequence if and only if A007947(n) is in the sequence. - Robert Israel, Aug 06 2018
LINKS
Robert Israel, Table of n, a(n) for n = 1..2000
MAPLE
isA013593 := proc(n)
numtheory[cyclotomic](n, x) ;
{coeffs(%, x)} ;
map(abs, %) ;
if % minus {1, 2, 3, 4} = {} then
false;
else
true;
end if;
end proc:
for n from 1 do
if isA013593(n) then
print(n);
end if;
end do: # R. J. Mathar, Nov 28 2016
MATHEMATICA
S[ n_ ] := For[ j=1; t=0, j<n, j++, t=Cases[ CoefficientList[ Cyclotomic[ j, x ], x ], k_ /; Abs[ k ]>4 ]; If[ Length[ t ]!=0, Print[ j ] ] ]; S[ 4845 ]
Select[Range[8000], Max[Abs[CoefficientList[Cyclotomic[#, x], x]]]>4&] (* Harvey P. Dale, Feb 27 2015 *)
PROG
(PARI) isok(n) = #select(x->abs(x) >= 5, Vec(polcyclo(n))) > 0; \\ Michel Marcus, Aug 07 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter T. Wang (peterw(AT)cco.caltech.edu)
EXTENSIONS
More terms from Don Reble, Dec 09 2001
Definition clarified by Harvey P. Dale, Feb 27 2015
STATUS
approved