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.
%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