OFFSET
1,5
COMMENTS
a(n) >= A000010(n)-1 since if 2<=k<n and (k,n)=1, the resultant is 1. - corrected by Robert Israel, Jul 24 2016
For n>1 a(n) = number of roots of the n-th polynomial in A275345, equal to 1. - Mats Granvik, Jul 24 2016
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
T. M. Apostol, Resultants of Cyclotomic Polynomials, Proc. Amer. Math. Soc. 24, 457-462, 1970.
T. M. Apostol, The Resultant of the Cyclotomic Polynomials Fm(ax) and Fn(bx), Math. Comput. 29, 1-6, 1975.
Eric Weisstein's World of Mathematics, Cyclotomic polynomials.
FORMULA
a(n) = n - A073093(n).
a(n) = n - A001222(n) - 1. - Michel Marcus, Jul 24 2016
MAPLE
seq(n -numtheory:-bigomega(n)-1, n=1..1000); # Robert Israel, Jul 25 2016
MATHEMATICA
Table[n - PrimeOmega@ n - 1, {n, 73}] (* Michael De Vlieger, Jul 26 2016 *)
PROG
(PARI) a(n)=sum(k=1, n, if(1-polresultant(polcyclo(n), polcyclo(k)), 0, 1))
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Oct 13 2002
EXTENSIONS
a(30)=2 and a(31)=6 merged into a(30)=26 by Mats Granvik, Jul 24 2016
STATUS
approved