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.

A385220 Primes p such that multiplicative order of 3 modulo p is odd.

Original entry on oeis.org

2, 11, 13, 23, 47, 59, 71, 83, 107, 109, 131, 167, 179, 181, 191, 227, 229, 239, 251, 263, 277, 311, 313, 347, 359, 383, 419, 421, 431, 433, 443, 467, 479, 491, 503, 541, 563, 587, 599, 601, 647, 659, 683, 709, 719, 733, 743, 757, 827, 829, 839, 863, 887, 911, 947, 971, 983
Offset: 1

Views

Author

Jianing Song, Jun 22 2025

Keywords

Comments

The multiplicative order of 3 modulo a(n) is A385226(n).
Without 2, contained in primes congruent to 1 or 11 modulo 12 (primes p such that 3 is a quadratic residue modulo p; A097933), and contains primes congruent to 11 modulo 12 (A068231).
Conjecture: this sequence has density 1/3 among the primes.

Crossrefs

A068231 < this sequence < A045317 < A040101 < A097933 (ignoring terms 2, 3), where Ax < Ay means that Ax is a subsequence of Ay.
Complement of A301916 in {primes} \ {3}.
Cf. A385226 (the actual multiplicative orders).
Cf. other bases: A014663 (base 2), this sequence (base 3), A385221 (base 4), A385192 (base 5), A163183 (base -2), A385223 (base -3), A385224 (base -4), A385225 (base -5).

Programs

  • Mathematica
    Select[Prime[Range[200]], OddQ[MultiplicativeOrder[3, #]] &] (* Paolo Xausa, Jun 28 2025 *)
  • PARI
    isA385220(p) = isprime(p) && (p!=3) && znorder(Mod(3,p))%2

A040101 Primes p such that x^4 = 3 has a solution mod p.

Original entry on oeis.org

2, 3, 11, 13, 23, 47, 59, 71, 83, 107, 109, 131, 167, 179, 181, 191, 193, 227, 229, 239, 251, 263, 277, 311, 313, 347, 359, 383, 419, 421, 431, 433, 443, 467, 479, 491, 503, 541, 563, 577, 587, 599, 601, 647, 659
Offset: 1

Views

Author

Keywords

Comments

Union of 2, 3, A068231 (primes congruent to 11 modulo 12), primes p == 1 (mod 4) such that 3^((p-1)/4) == 1 (mod p). - Jianing Song, Jun 22 2025

Crossrefs

A subsequence of A038874.
A068231 < A385220 < A045317 < this sequence < A097933 (ignoring terms 2, 3), where Ax < Ay means that Ax is a subsequence of Ay.

Programs

  • Magma
    [p: p in PrimesUpTo(800) | exists(t){x : x in ResidueClassRing(p) | x^4 eq 3}]; // Vincenzo Librandi, Sep 11 2012
    
  • Mathematica
    ok [p_]:=Reduce[Mod[x^4- 3, p] == 0, x, Integers] =!= False;  Select[Prime[Range[200]], ok] (* Vincenzo Librandi, Sep 11 2012 *)
  • PARI
    isA040101(p) = isprime(p) && (p==2 || p==3 || p%12==11 || (p%4==1 && Mod(3, p)^((p-1)/4) == 1)) \\ Jianing Song, Jun 22 2025

A045318 Primes p such that x^8 = 3 has no solution mod p.

Original entry on oeis.org

5, 7, 17, 19, 29, 31, 37, 41, 43, 53, 61, 67, 73, 79, 89, 97, 101, 103, 113, 127, 137, 139, 149, 151, 157, 163, 173, 193, 197, 199, 211, 223, 233, 241, 257, 269, 271, 281, 283, 293, 307, 317, 331, 337, 349, 353, 367, 373, 379, 389, 397
Offset: 1

Views

Author

Keywords

Comments

Complement of A045317 relative to A000040. - Vincenzo Librandi, Sep 19 2012

Crossrefs

Programs

  • Magma
    [p: p in PrimesUpTo(500) | not exists{x : x in ResidueClassRing(p) | x^8 eq 3} ]; // Vincenzo Librandi, Sep 19 2012
  • Mathematica
    ok[p_]:= Reduce[Mod[x^8 - 3, p] == 0, x, Integers] == False; Select[Prime[Range[200]], ok] (* Vincenzo Librandi, Sep 19 2012 *)

A320481 Primes in A301916 but not in A045318.

Original entry on oeis.org

2, 769, 1297, 6529, 7057, 8017, 8737, 12097, 12289, 13297, 13441, 14929, 15073, 15361, 15937, 16273, 18913, 19441, 20593, 21601, 21649, 22273, 22369, 23857, 25633, 26017, 26449, 26497, 27793, 28513, 30529, 31249, 34369, 34849, 36913, 37057, 37441, 37633, 38833, 38977, 39409
Offset: 1

Views

Author

N. J. A. Sloane, Oct 17 2018

Keywords

Comments

Is there a simpler characterization of these primes?
Answer from Don Reble, Oct 25 2018: (Start)
Let POT(x) be the largest power of 2 which divides x (A006519).
Apart from the initial 2, this sequence consists of those primes P such that
2 <= POT(the order of 3 modulo P) <= POT(P-1)/8.
The condition "2 <=" ensures that P divides some 3^k+1, and the condition "<= POT(P-1)/8" is so that 3 has an eighth root modulo P. A062117 is the order of 3 modulo prime(n). (End)
Comments from Richard Bumby, Nov 12 2018: (Start)
When considering methods for finding square roots mod p one is led to filtering the nonzero elements by the power of 2 dividing the multiplicative order of the element. The lowest level -- elements of odd order -- have easily computed square roots, and the square roots of other elements can be found if you can discover at least one element at a higher level.
To say that "x^8 = 3 has no solution mod p" is to say that 3 is in one of the top three levels and that there are more than 3 levels (so that 8 divides p-1).
To say that primes "divide numbers of the form 3^k + 1" is to say that -1 is a power of 3 mod p, or that 3 is not at the lowest level. If there are only four levels (9 mod 16), these statements are equivalent. Otherwise, the two statements are different. An interesting case has 3 at the second level, so that (-3) has odd order allowing cube roots of unity to be found quickly.
I was told that Odoni had some results on finding the number of primes with k levels for which a given number (e.g., 3) is at level j, but I never tracked down a reference. If the asymptotic behavior is what one would expect, A045318 and A301916 are really far from being "almost the same", except in the trivial sense of "zero density". (End)

References

  • Georg Fischer, email to N. J. A. Sloane, Oct 16 2018.

Crossrefs

Programs

  • Mathematica
    Select[Prime@Range@4200,PowerModList[3,1/8,#]!={}&&IntegerQ@MultiplicativeOrder[3,#,-1]&] (* Giorgos Kalogeropoulos, Feb 23 2022 *)

Extensions

More terms from Michel Marcus, Oct 17 2018
Showing 1-4 of 4 results.