cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A327818 Square array read by ascending antidiagonals: T(n,m) is the number of irreducible factors in the factorization of the m-th cyclotomic polynomial over GF(k), k = A246655(n) (counted with multiplicity).

This page as a plain text file.
%I A327818 #52 Jul 21 2024 12:49:12
%S A327818 1,1,1,1,1,1,1,1,2,2,1,1,2,1,1,1,1,1,2,1,1,1,1,2,2,2,2,2,1,1,1,1,4,2,
%T A327818 1,4,1,1,2,2,1,1,2,2,1,1,1,1,2,1,2,1,4,6,1,1,1,2,1,2,1,6,2,2,1,1,1,1,
%U A327818 2,2,4,2,6,2,1,2,2,2,1,1,1,2,1,1,2,4,2,4,2,2,1
%N A327818 Square array read by ascending antidiagonals: T(n,m) is the number of irreducible factors in the factorization of the m-th cyclotomic polynomial over GF(k), k = A246655(n) (counted with multiplicity).
%C A327818 T(n,m) = 1 if and only if (let k = A246655(n)):
%C A327818 (a) gcd(m,k) = 1, and k is a primitive root modulo m;
%C A327818 (b) k is a power of 2, m == 2 (mod 4), and k is a primitive root modulo m/2.
%C A327818 As a result, T(n,m) > 1 if at least one of the following holds:
%C A327818 (i) There is no primitive root modulo m, that is, m is in A033949;
%C A327818 (ii) k is a square number and m > 2.
%C A327818 If p = A246655(n) is prime, then T(n,m) is also the number of irreducible factors in the factorization of the ideal (p) in Z[zeta_m], zeta_m = exp(2*Pi*i/m). Actually, if the m-th cyclotomic polynomial factors as Product_{i=1..T(n,m)} F_i(x) over GF(p), then the factorization of (p) in Z[zeta_m] is (p) = Product_{i=1..T(n,m)} (p,F_i(zeta_m)). As a result, p remains inert in Q(zeta_m) <=> T(n,m) = 1. See Page 47-48, Proposition 8.3 and Page 61-62, Proposition 10.3 of the Neukirch link for a proof. - _Jianing Song_, Sep 13 2022
%H A327818 Amiram Eldar, <a href="/A327818/b327818.txt">Table of n, a(n) for n = 1..10011</a> (antidiagonals n = 1..141, flattened)
%H A327818 Jürgen Neukirch, <a href="http://www.math.toronto.edu/~ila/Neukirch_Algebraic_number_theory.pdf">Algebraic_number_theory</a>.
%H A327818 Jianing Song, <a href="/A327818/a327818_2.pdf">My notes on the factorization of cyclotomic polynomials over finite fields</a>.
%H A327818 Hongfeng Wu, Li Zhu, Rongquan Feng, and Siman Yang, <a href="https://doi.org/10.1007/s10623-016-0224-5">Explicit factorizations of cyclotomic polynomials over finite fields</a>, Des. Codes Cryptogr. 83, 197-217 (2017).
%F A327818 Let m = p^e*s, where p is the prime factor of k = A246655(n), gcd(p,s) = 1, then T(n,m) = phi(m)/ord(k,s), where phi = A000010, ord(k,s) is the multiplicative order of k modulo s. Proof:
%F A327818 (a) First consider the case e = 0. Let F be the algebraic closure of GF(k), then Phi_s(x) = 0 has solutions in F, where Phi_s(x) is the s-th cyclotomic polynomial. Let a be any of the solution.
%F A327818     In F, a belongs to GF(k^d) <=> a^(k^d-1) = 1 <=> s|(k^d-1) (note that in F, s is the smallest integer such that a^s = 1). As a result, a belongs to GF(k^ord(k,s)) but not to GF(k^d) for any d < ord(k,s), that is to say, a has algebraic degree ord(k,s) over GF(k). Because a is any of the roots of Phi_s(x) in F, every irreducible factor of Phi_s(x) over GF(k) is of degree ord(k,s), so the number of irreducible factors is phi(s)/ord(k,s);
%F A327818 (b) If e > 0, we can see from the Moebius inversion formula such that Phi_m(x) = (Phi_s(x))^phi(p^e), that the number of irreducible factors is phi(m)/ord(k,s).
%F A327818 The Introduction part in Page 2 and Lemma 1 in Page 3 of Hongfeng Wu's paper (see Links section) also mentions part (a) of this result.
%e A327818 Table starts
%e A327818   n  A246655(n)  m=1  2  3  4  5  6  7  8  9  10
%e A327818   1       2       1   1  1  2  1  1  2  4  1   1
%e A327818   2       3       1   1  2  1  1  2  1  2  6   1
%e A327818   3       4       1   1  2  2  2  2  2  4  2   2
%e A327818   4       5       1   1  1  2  4  1  1  2  1   4
%e A327818   5       7       1   1  2  1  1  2  6  2  2   1
%e A327818   6       8       1   1  1  2  1  1  6  4  3   1
%e A327818   7       9       1   1  2  2  2  2  2  4  6   2
%e A327818   8      11       1   1  1  1  4  1  2  2  1   4
%t A327818 T[max_] := Module[{pp = Select[Range[max], PrimePowerQ], m, t}, m = Length[pp]; t[i_, j_] := Module[{p = FactorInteger[pp[[j]]][[1, 1]]}, EulerPhi[i] / MultiplicativeOrder[pp[[j]], i/p^IntegerExponent[i, p]]]; Table[t[j, i - j + 1], {i, 1, m}, {j, 1, i}] // Flatten]; T[24] (* _Amiram Eldar_, Jul 21 2024 *)
%o A327818 (PARI) f(k,m) = if(isprimepower(k), my(p=factor(k)[1,1], s=m/p^valuation(m,p)); eulerphi(m)/znorder(Mod(k,s)))
%o A327818 A246655(n) = my(i=0); for(t=1, oo, if(isprimepower(t), i++); if(i==n, return(t)))
%o A327818 a(n,m) = f(A246655(n),m)
%Y A327818 Cf. A246655, A000010, A033949.
%Y A327818 Rows 1..7 give: A318622, A327812, A327813, A327814, A327815, A327816, A327817.
%K A327818 nonn,easy,tabl
%O A327818 1,9
%A A327818 _Jianing Song_, Sep 26 2019