A385169 Primes p == 3 (mod 4) such that the multiplicative order of 2+-i modulo p in Gaussian integers (A385165) is odd.

331, 571, 599, 691, 839, 971, 1051, 1171, 1291, 1451, 1571, 1879, 2131, 2411, 2971, 3251, 3331, 3491, 3571, 3691, 3851, 4051, 4091, 4211, 4651, 4679, 4691, 4919, 4931, 5051, 5171, 5479, 5531, 5651, 5839, 5851, 5879, 6011, 6599, 6679, 6691, 7079, 7211, 7331, 7691, 8011, 8039, 8171, 8731, 8839, 9011, 9371, 9811

Offset

1,1

Comments

Primes p == 3 (mod 4) are precisely the rational primes in the ring of Gaussian integers.

Let ord(a,m) be the multiplicative order of a modulo m. (Of course if a and m are integers, it doesn't matter if the base ring is Z or Z[i]). For a prime p == 3 (mod 4), we have that ord(2+-i,p) is divisible by ord(5,p), and that ord(2+-i,p) divides (p+1) * ord(5,p). What's more, ord(2+-i,p) divides (p^2-1)/2 if and only if 5 is a quadratic residue of integers modulo p. (See A385165).

As a result, if ord(2+-i,p) is not divisible by 8, then ord(5,p) is odd:

- Of course this is true if ord(2+-i,p) is odd.

- If ord(2+-i,p) == 2 (mod 4) and ord(5,p) is even, then ord(2+-i,p)/ord(5,p) is odd, and so ord(2+-i,p) divides ((p+1)/4) * ord(5,p), then ord(5,p) is odd. This implies that ord(2+-i,p) is odd, a contradiction.

- If ord(2+-i,p) == 4 (mod 8) and ord(5,p) is even (we have ord(5,p) == 2 (mod 4) since p == 3 (mod 4)), then ord(2+-i,p)/ord(5,p) == 2 (mod 4), and so ord(2+-i,p) divides ((p+1)/2) * ord(5,p), then ord(5,p) is odd. This implies that ord(2+-i,p) == 2 (mod 4), a contradiction.

From the above paragraph, this sequence is also primes p == 3 (mod 4) such that ord(2+-i,p)/ord(5,p) is odd.

Links

Jianing Song, Table of n, a(n) for n = 1..10000

Example

8731 is a term since (2+-i)^635253 == 1 (mod 8731), and 635253 is odd.
8839 is a term since (2+-i)^57447 == 1 (mod 8839), and 57447 is odd.
9011 is a term since (2+-i)^2029953 == 1 (mod 9011), and 2029953 is odd.

Prog

(PARI) ord(p) = my(d = divisors((p+1)*znorder(Mod(5, p)))); for(i=1, #d, if(Mod([2, -1; 1, 2], p)^d[i] == 1, return(d[i]))) \\ for a prime p == 3 (mod 4), returns ord(2+-i, p)
isA385169(p) = isprime(p) && p%4==3 && ord(p)%2

Crossrefs

Cf. A385165, A385179, A385192, A385217 (the actual multiplicative orders).

A385188 < this sequence < A385180 < A385167 < intersection of A122869 and A385168, where Ax < Ay means that Ax is a subsequence of Ay.

Sequence in context: A142763 A020373 A385180 * A142601 A210534 A210511

Adjacent sequences: A385166 A385167 A385168 * A385170 A385171 A385172

Keyword

nonn,easy

Author

Jianing Song, Jun 20 2025

Status

approved