%I #74 Jan 08 2024 22:42:15
%S 0,1,-1,1,1,1,1,1,1,1,0,1,1,1,1,1,1,1,-1,1,1,1,1,1,1,1,1,1,0,0,0,1,1,
%T 0,0,1,0,0,1,1,-1,1,-1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,-1,0,1,1,1,1,1,1,1,
%U 1,1,1,1,1,1,1,1,-1,1,-1,1,-1,1,1,-1,0,1,-1,1,0,-1,1
%N Irregular triangle read by rows: coefficients of cyclotomic polynomial Phi_n(x) (exponents in increasing order).
%C We follow Maple in defining Phi_0 to be x; it could equally well be taken to be 1.
%C From _Wolfdieter Lang_, Oct 29 2013: (Start)
%C The length of row n >= 1 of this table is phi(n) + 1 = A000010(n) + 1. Row n = 0 has here length 2.
%C Phi_n(x) is the minimal polynomial of omega_n := exp(i*2*Pi/n) over the rationals. Namely, Phi_n(x) = Product_{k=0..n-1, gcd(k,n)=1} (x - (omega_n)^k). See the Graham et al. reference, 4.50 a, pp. 149, 506.
%C Phi_n(x) = Product_{d|n} (x^d - 1)^(mu(n/d)) with the Moebius function mu(n) = A008683(n), n >= 1. See the Graham et al. reference, 4.50 b, pp. 149, 506.
%C Phi_n(x) = Phi_{rad(n)}(x^(n/rad(n))), n >= 2, with rad(n) = A007947(n), the squarefree kernel of n. Proof from the preceding formula, where only squarefree n/d (A005117) from the set of divisors of n enter, by mapping each factor (numerator or denominator) of the left hand side to one of the right hand side and vice versa.
%C (End)
%C Each row can be considered as the last column of the companion matrix of the cyclotomic polynomial: A000010(n) is the size of such a square matrix, last column has opposite signs and the last term (before last term of each row in A013595) equal to A008683(n). - _Eric Desbiaux_, Dec 14 2015
%D E. R. Berlekamp, Algebraic Coding Theory, McGraw-Hill, 1968; see p. 90.
%D Z. I. Borevich and I. R. Shafarevich, Number Theory. Academic Press, NY, 1966, p. 325.
%D R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics, Addison-Wesley, 1991, p. 137.
%D K. Ireland and M. Rosen, A Classical Introduction to Modern Number Theory, Springer, 1982, p. 194.
%H Robert Israel, <a href="/A013595/b013595.txt">Table of n, a(n) for n = 0..10000</a>
%H Emma Lehmer, <a href="http://dx.doi.org/10.1090/S0002-9904-1936-06309-3">On the magnitude of the coefficients of the cyclotomic polynomial</a>, Bull. Amer. Math. Soc. 42 (1936), 389-392.
%H Rakshith Rajashekar, Marco Di Renzo, K.V.S. Hari, L. Hanzo, <a href="https://eprints.soton.ac.uk/414186/1/Full_duplex_SM_Drones.pdf">A generalised transmit and receive diversity condition for feedback assisted MIMO systems: theory & applications in full-duplex spatial modulation</a>, 2017.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/CyclotomicPolynomial.html">Cyclotomic Polynomial</a>.
%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Cyclotomic_polynomial">Cyclotomic Polynomial</a>.
%F a(n,m) = [x^m] Phi_n(x), n >= 0, 0 <= m <= phi(n), with phi(n) = A000010(n). - _Wolfdieter Lang_, Oct 29 2013
%e Phi_0 = x; Phi_1 = x - 1; Phi_2 = x + 1; Phi_3 = x^2 + x + 1; Phi_4 = x^2 + 1; ...
%e From _Wolfdieter Lang_, Oct 29 2013: (Start)
%e The irregular triangle a(n,m) begins:
%e n\m 0 1 2 3 4 5 6 7 8 9 10 11 12 ...
%e 0: 0 1
%e 1: -1 1
%e 2: 1 1
%e 3: 1 1 1
%e 4: 1 0 1
%e 5: 1 1 1 1 1
%e 6: 1 -1 1
%e 7: 1 1 1 1 1 1 1
%e 8: 1 0 0 0 1
%e 9: 1 0 0 1 0 0 1
%e 10: 1 -1 1 -1 1
%e 11: 1 1 1 1 1 1 1 1 1 1 1
%e 12: 1 0 -1 0 1
%e 13: 1 1 1 1 1 1 1 1 1 1 1 1 1
%e 14: 1 -1 1 -1 1 -1 1
%e 15: 1 -1 0 1 -1 1 0 -1 1
%e ...
%e Phi_15(x) = (x^1 - 1)*((x^3 - 1)^(-1))*((x^5 - 1)^(-1))*(x^15 - 1) because mu(15) = mu(1) = +1 and mu(3) = mu(5) = -1. Hence Phi_15(x) = 1 - x + x^3 - x^4 + x^5 - x^7 + x^8, giving row n = 15.
%e Example for the reduction via the squarefree kernel: Phi_12(x) = Phi_6(x^(12/6)) = Phi_6(x^2). By the formula with the Mobius function Phi_6(x) = Phi_2(x^3)/Phi_2(x) = 1 - x + x^2 and with x -> x^2 this becomes Phi_12(x) = 1 - x^2 + x^4.
%e (End)
%p N:= 100: # to get coefficients up to cyclotomic(N,x)
%p with(numtheory):
%p for n from 0 to N do
%p C:= cyclotomic(n,x);
%p L[n]:= seq(coeff(C,x,i),i=0..degree(C));
%p od:
%p A:= [seq](L[n],n=0..N): # note that A013595(n) = A[n+1]
%p # _Robert Israel_, Apr 17 2014
%t Table[CoefficientList[x^KroneckerDelta[n] Cyclotomic[n, x], x], {n, 0, 15}] // Flatten (* _Peter Luschny_, Dec 27 2016 *)
%o (PARI) row(n) = if (n==0, p=x, p = polcyclo(n)); Vecrev(p); \\ _Michel Marcus_, Dec 14 2015
%Y Cf. A013596, A020500 (row sums, n >= 1), A020513 (alternating row sums).
%Y For record coefficients see A160340, A262404, A262405, A278567.
%Y Column m=1 is A157657.
%K sign,easy,nice,tabf
%O 0,3440
%A _N. J. A. Sloane_
%E Maple program corrected by _Robert Israel_, Apr 17 2014