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.

Showing 1-4 of 4 results.

A305236 Numbers n such that the multiplicative group of integers modulo n is isomorphic to C_m X C_m, m > 1.

Original entry on oeis.org

8, 12, 63, 126, 513, 1026, 2107, 4214, 12625, 25250, 26533, 39609, 53066, 79218, 355023, 710046, 3190833, 4457713, 6381666, 8915426, 19854847, 38463283, 39709694, 76926566, 242138449, 370634743, 484276898, 516465451, 574336561, 701607583, 741269486, 1032930902, 1148673122, 1380336193, 1403215166, 2324581983, 2760672386, 4649163966, 4882890625, 6174434113, 9765781250
Offset: 1

Views

Author

Jianing Song, Jun 19 2018

Keywords

Comments

Note that 24 is only number k such that the multiplicative group of integers modulo k is isomorphic to C_m X C_m X C_m, m > 1.
The number of elements in the multiplicative group of integers modulo a(n) of order d is A007434(d), whenever d is divisible by A002322(a(n)).
The corresponding m (=A002322(a(n))) are 2, 2, 6, 6, 18, 18, 42, 42, 100, 100, 156, 162, 156, 162, 486, 486, 1458, 2028, 1458, 2028, ... Each term in A114874, except for those of the form 2^k, k >= 2, occurs exactly twice in this list.
Numbers k such that A046072(k) = 2 and A316089(k) = 1. - Jianing Song, Sep 15 2018
Except for 8 and 12, these are numbers of the form p^e*((p-1)*p^(e-1) + 1) or 2*p^e*((p-1)*p^(e-1) + 1) where p is an odd prime and (p-1)*p^(e-1) + 1 is prime. - Jianing Song, Apr 13 2019

Examples

			The multiplicative group of integers modulo 63 is isomorphic to C_6 X C_6. There are A007434(1) = 1 element of order 1, A007434(2) = 3 elements of order 2, A007434(3) = 8 elements of order 3, A007434(6) = 24 elements of order 6 modulo 63.
The multiplicative group of integers modulo 513 is isomorphic to C_18 X C_18. There are A007434(1) = 1 element of order 1, A007434(2) = 3 elements of order 2, A007434(3) = 8 elements of order 3, A007434(6) = 24 elements of order 6, A007434(9) = 72 elements of order 9, A007434(18) = 216 elements of order 18 modulo 513.
		

Crossrefs

Cf. A114874.
Odd terms are given by A307527.

Programs

  • PARI
    for(n=1,10^7,if(#znstar(n)[2]==2 && znstar(n)[2][1]==znstar(n)[2][2], print1(n, ", "))) \\ Jianing Song, Sep 15 2018
    
  • PARI
    the_first_entries(nn) = my(u=[]); for(n=2, sqrt(nn), my(v=factor(n), d=#v[, 1], p=v[d, 1], e=v[d, 2]); if(isprime(n+1) && p!=2 && n==(p-1)*p^e, u=concat(u, [(n+1)*p^(e+1)]))); t=concat([8, 12], concat(u, 2*u)); t=vecsort(select(i->(iJianing Song, Apr 13 2019

Formula

A302257(a(n)) = A258615(a(n))/2.

Extensions

Missing a(40) inserted by Jianing Song, Apr 20 2019

A316089 Number of ordered (k_1, k_2, ..., k_r) such that (Z/nZ)* is isomorphic to C_{k_1} X C_{k_2} X ... X C_{k_r}. Here (Z/nZ)* is the multiplicative group of integers modulo n, r = rank((Z/nZ)*).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 4, 2, 1, 1, 4, 3, 1, 2, 1, 2, 4, 1, 1, 3, 1, 1, 2, 4, 1, 1, 4, 3, 2, 1, 1, 3, 1, 1, 1, 2, 2, 2, 1, 2, 2, 4, 1, 3, 1, 1, 4, 2, 4, 4, 1, 3, 1, 1, 1, 3, 2, 1, 4
Offset: 1

Views

Author

Jianing Song, Jun 27 2018

Keywords

Comments

The rank of a finitely generated group rank(G) is defined as the size of the minimal generating sets (the generating sets with the least elements) of G. By convention a(1) = a(2) = 1, since the multiplicative group of integers modulo 1 or 2 has rank 0, and the empty direct product yields the trivial group.
This sequence is related to the number of minimal generating sets {x_1, x_2, ..., x_r} of (Z/nZ)* such that Product_{i=1..r} x_i^e_i == 1 (mod n) implies that x_i^e_i == 1 (mod n) for 1 <= i <= r. See A302257.
For n >= 3, also the number of ordered (k_1, k_2, ..., k_r) such that (Z/nZ)* is isomorphic to C_{k_1} X C_{k_2} X ... X C_{k_r} and gcd(k_1, k_2, ..., k_r) > 1.
a(n) = 1 iff there exists an integer m such that (Z/nZ)* = (C_m)^r, that is, n is a term in A033948 or A305236 or n = 24.
First occurrence of k, or -1 if k never occurs: 1, 15, 40, 35, -1, 160, -1, 143, 104, -1, -1, 480, -1, -1, -1, ... It's easy to show that prime number >= 5 never occurs.
In general, let (Z/nZ)* = C_(Product_{j=1..k} p_j^e_1j) X C_(Product_{j=1..k} p_j^e_2j) X ... X C_(Product_{j=1..k} p_j^e_rj) where p_1, p_2, ..., p_k are different primes, 0 <= e_1j <= e_2j <= ... <= e_rj for 1 <= j <= k (for n = 1 or 2, (Z/nZ)* is the trivial group, thus r = k = 0), then a(n) is the number of matrices C = (c_ij) such that (c_1j, c_2j, ..., c_rj) is a permutation of (e_1j, e_2j, ..., e_rj) for 1 <= j <= k. - Jianing Song, Feb 04 2019

Examples

			For n = 35, rank((Z/35Z)*) = 2, and (Z/35Z)* = C_2 X C_12 = C_12 X C_2 = C_4 X C_6 = C_6 X C_4, so a(35) = 4. Also, (Z/35Z)* = C_(2^1*3^0) X C_(2^2*3^1), e_11 = 1, e_12 = 0, e_21 = 2, e_22 = 1, so b_11 = b_12 = b_21 = b_22 = 1, then a(35) = (2!)^2 = 4.
For n = 105, rank((Z/105Z)*) = 3, and (Z/105Z)* = C_2 X C_2 X C_12 (along with its two permutations) = C_2 X C_4 X C_6 (along with its five permutations), so a(105) = 9. Also, (Z/35Z)* = C_(2^1*3^0) X C_(2^1*3^0) X C_(2^2*3^1), e_11 = 1, e_12 = 0, e_21 = 1, e_22 = 0, e_31 = 2, e_32 = 1, so b_11 = b_12 = b_21 = b_22 = 2, b_31 = b_32 = 1, then a(105) = (3!)^2/(2!*2!) = 9.
		

Crossrefs

Programs

  • PARI
    zp(g)={sum(i=1, #g, my(f=factor(g[i])); sum(j=1, #f~, x^f[j,1]*(y^f[j,2]-1)))}
    a(n)={my(g=znstar(n).cyc); my(p=zp(g), r=#g); prod(i=1, poldegree(p), my(u=Vec(r + polcoeff(p,i))); r!/prod(j=1, #u, u[j]!))} \\ Andrew Howroyd, Jun 30 2018

Formula

Let (Z/nZ)* = C_(Product_{j=1..k} p_j^e_1j) X C_(Product_{j=1..k} p_j^e_2j) X ... X C_(Product_{j=1..k} p_j^e_rj) where p_1, p_2, ..., p_k are different primes, 0 <= e_1j <= e_2j <= ... <= e_rj for 1 <= j <= k. Define b_ij = max{1 <= t <= r | e_tj = e_ij} - min{1 <= t <= r | e_tj = e_ij} + 1, then a(n) = (r!)^k/Product_{i=1..r, j=1..k} (b_ij)!^(1/b_ij).

A323021 A way to generate a minimal generating set of (Z/nZ)* whose elements are multiplicatively independent for every n: row n lists the elements of S(n), where S(n) is defined in the comment section. Here (Z/nZ)* is the multiplicative group of integers modulo n.

Original entry on oeis.org

2, 3, 2, 5, 3, 5, 7, 2, 7, 2, 5, 7, 2, 3, 7, 11, 5, 15, 3, 11, 2, 11, 17, 8, 10, 13, 5, 7, 13, 17, 2, 15, 2, 15, 17, 2, 7, 11, 3, 5, 31, 13, 23, 3, 22, 31, 19, 29, 2, 21, 14, 28, 17, 21, 31, 6, 29, 31, 3, 13, 23, 11, 37, 5, 5, 17, 31, 37, 3, 27, 35, 37, 27, 41
Offset: 3

Views

Author

Jianing Song, Jan 11 2019

Keywords

Comments

For any finite abelian multiplicative group G and a generating set (not necessarily minimal) S = {x_1, x_2, ..., x_m}, we say the elements in S are multiplicatively independent if Product_{i=1..m} (x_i)^(e_i) = o implies that (x_i)^(e_i) = o for i = 1..m", where o is the group identity. For example, let G = (Z/16Z)*, the elements in {3, 15} are multiplicatively independent, while the elements in {3, 5} are not because 3^2 * 5^2 == 1 (mod 16) but 3^2 == 5^2 == 9 (mod 16). - Jianing Song, Mar 16 2019
S(n) is defined as: S(p^e) = {the least primitive root modulo p^e} if p is an odd prime; S(2) = the empty set, S(4) = {3}, S(2^e) = {5, 2^e - 1} for e >= 3. If gcd(m_1, m_2) = 1, then S((m_1)*(m_2)) = {1 <= x <= (m_1)*(m_2) : x belongs to S(m_1) and x == 1 (mod m_2), or x belongs to S(m_2) and x == 1 (mod x_1)}. The reason we choose S(2^e) = {5, 2^e - 1} is that for n >= 3, 1, 5, 5^2, ..., 5^(2^(e-2)-1) are all elements in (Z/(2^e)Z)* that are congruent to 1 modulo 4 and their additive inverses are all elements that are congruent to 3 modulo 4.
Offset 3, because the generating set of (Z/1Z)* or (Z/2Z)* is the empty set. The length of the n-th row is A046072(n) for n >= 3.
We now show that S(n) is a generating set of (Z/nZ)* whose elements are multiplicatively independent by induction. According to A302257, we only need to show that for every n, S(n) = {x_1, x_2, ..., x_r} satisfies that a one-to-one mapping can be set between all products of the form Product_{i=1..r} (x_i)^(e_i) and the elements in (Z/nZ)*, where 0 <= e_i < ord(x_i,n) for 1 <= i <= r. Suppose gcd(m_1, m_2) = 1. Let S(m_1) = {x_1, x_2, ..., x_s}, S(m_2) = {y_1, y_2, ..., y_t}, then S((m_1)*(m_2)) = {x'1, x'_2, ..., x'_s, y'_1, y'_2, ..., y'_t} where x'_i == x_i (mod m_1), x'_i == 1 (mod m_2), y'_i == 1 (mod m_1), y'_i == y_i (mod m_2). For every a coprime to n, solving a == Product{i=1..s} (x'i)^(e_i) * Product{i=1..t} (y'i)^(f_i) (mod (m_1)*(m_2)) is equivalent to solving a == Product{i=1..s} (x'i)^(e_i) (mod m_1) and a == Product{i=1..s} (y'i)^(f_i) (mod m_2) separately, which are immediately reduced to a == Product{i=1..s} (x_i)^(e_i) (mod m_1) and a == Product_{i=1..s} (y_i)^(f_i) (mod m_2). Since the elements in S(m_1) are multiplicatively independent, the elements in S(m_2) are multiplicatively independent as well, hence the elements in S((m_1)*(m_2)) must also be multiplicatively independent. For example, S(5) = {2}, S(7) = {3}, so S(35) = {22, 31}. For every a coprime to 35, solving a == 22^e * 31^f (mod 35) is equivalent to solving a == 2^e * 1^f (mod 35) and a == 1^e * 3^f (mod 7) separately, which are immediately reduced to a == 2^e (mod 5) and a == 3^f (mod 7). See A302257 for more information. [Revised by Jianing Song, Mar 16 2019]
Also, S(n) can be used to show that if gcd(m_1, m_2) = 1, then every Dirichlet character modulo (m_1)*(m_2) can be written as the product of a character modulo m_1 and a character modulo m_2. Define S(m_1), S(m_2) and S((m_1)*(m_2)) as above. Let {Chi(n, k)} be a Dirichlet character modulo n. {Chi((m_1)*(m_2), k)} is wholly determined by Chi((m_1)*(m_2), x'_1), Chi((m_1)*(m_2), x'_2), ..., Chi((m_1)*(m_2), x'_r), Chi((m_1)*(m_2), y'_1), Chi((m_1)*(m_2), y'_2), ..., Chi((m_1)*(m_2), y'_s). We have Chi(m_1, t)*Chi(m_2, t) = Chi((m_1)*(m_2), t) for every t belongs to S((m_1)*(m_2)). These equations then reduce to Chi(m_1, x_i) = Chi((m_1)*(m_2), x'_i) and Chi(m_2, y_i) = Chi((m_1)*(m_2), y'_i), so {Chi(m_1, k)} and {Chi(m_2, k)} are determined. It's easy to see that the product of a character modulo m_1 and a character modulo m_2 is a character modulo (m_1)*(m_2), so all characters modulo (m_1)*(m_2) can be constructed using the characters modulo m_1 and the characters modulo m_2.

Examples

			Table begins
(1)  Empty set;
(2)  Empty set;
(3)  {2};
(4)  {3};
(5)  {2};
(6)  {5};
(7)  {3};
(8)  {5, 7};
(9)  {2};
(10) {7};
...
Example shows that S(560) = {241, 337, 351, 421}: (Start)
Since 560 = 16*5*7, we have S(560) = {1 <= x <= 560: x == 2 (mod 5) and x == 1 (mod 16*7), or x == 3 (mod 7) and x == 1 (mod 16*5), or x == 5 (mod 16) and x == 1 (mod 5*7), or x == 15 (mod 16) and x == 1 (mod 5*7)} = {241, 337, 351, 421}.
To find the index of a number, say, 403, with respect to the base (241, 337, 351, 421) modulo 560, suppose 241^a * 337^b * 351^c * 421^d == 403 (mod 560). Then we have:
- Modulo 16: 15^c * 5^d == 3 (mod 16) => c = 1, d = 3;
- Modulo 5: 2^b == 3 (mod 5) => b = 3;
- Modulo 7: 3^a == 4 (mod 7) => a = 4.
So the index of 403 with respect to the base (241, 337, 351, 421) modulo 560 is (4, 3, 1, 3). Note that the bases here are ordered. So the index of 403 with respect to the base (337, 241, 421, 351) modulo 560 is (3, 4, 3, 1). (End)
		

Crossrefs

Programs

  • PARI
    a(n) = my(v=vector(#znstar(n)[2]), f=factor(n)); if(n%2||n%8==4, for(i=1, #v, v[i]=lift(chinese(znprimroot(f[i,1]^f[i,2]), Mod(1, n/f[i,1]^f[i,2]))))); if(n%4==2, for(i=1, #v, v[i]=lift(chinese(znprimroot(f[i+1,1]^f[i+1,2]), Mod(1, n/f[i+1,1]^f[i+1,2]))))); if(n%8==0, v[1]=lift(chinese(Mod(5, 2^f[1,2]), Mod(1, n/2^f[1,2]))); v[2]=lift(chinese(Mod(-1, 2^f[1,2]), Mod(1, n/2^f[1,2]))); for(i=3, #v, v[i]=lift(chinese(znprimroot(f[i-1,1]^f[i-1,2]), Mod(1, n/f[i-1,1]^f[i-1,2]))))); v=vecsort(v); v

Extensions

Name simplified by Jianing Song, Mar 16 2019

A327791 Number of ways, up to order, of decomposing (Z/nZ)* as the internal direct product of r cyclic subgroups, where r = rank((Z/nZ)*).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 3, 1, 1, 4, 4, 1, 1, 1, 4, 6, 1, 1, 28, 1, 1, 1, 6, 1, 4, 1, 4, 6, 1, 8, 6, 1, 1, 8, 48, 1, 6, 1, 6, 8, 1, 1, 48, 1, 1, 4, 8, 1, 1, 8, 84, 6, 1, 1, 48, 1, 1, 36, 4, 24, 6, 1, 4, 6, 8, 1, 84, 1, 1, 8, 6, 12, 8, 1, 192, 1, 1, 1, 84, 16, 1, 8, 84, 1, 8
Offset: 1

Views

Author

Jianing Song, Sep 25 2019

Keywords

Comments

In general, let G be a finite abelian group and r = rank(G), let s be the number of arrays (k_1, k_2, ..., k_r) such that G = C_(k_1) X C_(k_2) X ... X C_(k_r). We can see from A302257 that:
(b) Let M = Product_{i=1..r} phi(k_i), then the number of sets {H_1, H_2, ..., H_r} such that G = H_1 X H_2 X ... X H_r is (|Aut(G)|*s/r!)/M. Given {H_1, H_2, ..., H_r} such that G = H_1 X H_2 X ... X H_r, there are phi(|H_i|) ways to choose a generator of H_i, that is, there are Product_{i=1..r} phi(|H_i|) sets {a_1, a_2, ..., a_r} such that = H_i. From (a), the total number of sets {a_1, a_2, ..., a_r} as {H_1, H_2, ..., H_r} runs through all possible choices is |Aut(G)|*s/r!. Note that Product_{i=1..r} phi(|H_i|) is always equal to M regardless of the choice of {H_1, H_2, ..., H_r}. That is to say, each set of cyclic subgroups {H_1, H_2, ..., H_r} corresponds to M sets of generators {a_1, a_2, ..., a_r}. Thus, the number of sets {H_1, H_2, ..., H_r} is (|Aut(G)|*s/r!)/M.

Examples

			n = 8: (Z/8Z)* = {1, 3} X {1, 5} = {1, 3} X {1, 7} = {1, 5} X {1, 7}, so a(8) = 3;
n = 16: (Z/16Z)* = {1, 3, 9, 11} X {1, 7} = {1, 3, 9, 11} X {1, 15} = {1, 5, 9, 13} X {1, 7} = {1, 5, 9, 13} X {1, 15}, so a(16) = 4;
n = 35: (Z/35Z)* = {1, 2, 4, 8, 9, 11, 16, 18, 22, 23, 29, 32} X {1, 6} = {1, 2, 4, 8, 9, 11, 16, 18, 22, 23, 29, 32} X {1, 34} = {1, 3, 4, 9, 11, 12, 13, 16, 17, 27, 29, 33} X {1, 6} = {1, 3, 4, 9, 11, 12, 13, 16, 17, 27, 29, 33} X {1, 34} = {1, 8, 22, 29} X {1, 11, 16, 19, 24, 34} = {1, 8, 22, 29} X {1, 6, 11, 16, 26, 31} = {1, 13, 27, 29} X {1, 11, 16, 19, 24, 34} = {1, 13, 27, 29} X {1, 6, 11, 16, 26, 31}, so a(35) = 8.
		

Crossrefs

Programs

Formula

a(n) = A302257(n)/A327790(n).
Showing 1-4 of 4 results.