A059898 Duplicate of A014752.
31, 43, 109, 127, 157, 223, 229, 277, 283, 307, 397, 433, 439, 457, 499, 601, 643, 691
Offset: 1
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.
QuadPrimes2[2, 0, 13, 10000] (* see A106856 *)
list(lim)=my(v=List([2,13]),t); for(y=1,sqrtint(lim\13), for(x=1,sqrtint((lim-13*y^2)\2), if(isprime(t=2*x^2+13*y^2), listput(v,t)))); Set(v) \\ Charles R Greathouse IV, Feb 07 2017
[ p: p in PrimesUpTo(433) | exists(t){x : x in ResidueClassRing(p) | x^3 eq 2} ]; // Klaus Brockhaus, Dec 02 2008
f[p_] := Block[{k = 2}, While[k < p && Mod[k^3, p] != 2, k++ ]; If[k == p, 0, 1]]; Select[ Prime[ Range[100]], f[ # ] == 1 &] (* Robert G. Wilson v, Jul 26 2004 *)
select(p->ispower(Mod(2,p),3),primes(100)) \\ Charles R Greathouse IV, Apr 28 2015
[ p: p in PrimesUpTo(4597) | r eq 1 and Order(R!2) eq q where q,r is Quotrem(p,3) where R is ResidueClassRing(p) ]; // Klaus Brockhaus, Dec 02 2008
Reap[For[p = 2, p < 10^4, p = NextPrime[p], If[MultiplicativeOrder[2, p] == (p-1)/3, Sow[p]]]][[2, 1]] (* Jean-François Alcover, Dec 10 2015 *)
forprime(p=3,10^4,if(znorder(Mod(2,p))==(p-1)/3,print1(p,", "))); \\ Joerg Arndt, May 17 2013
[ p: p in PrimesUpTo(1000) | #[ x: x in ResidueClassRing(p) | x^3 eq 2 ] eq 1 ]; // Klaus Brockhaus, Apr 11 2009
Select[Prime[Range[150]],MemberQ[{0,2},Mod[#,3]]&] (* Harvey P. Dale, Jun 14 2011 *)
is(n)=isprime(n) && n%3!=1 \\ Charles R Greathouse IV, Apr 20 2015
from itertools import count, islice from sympy import isprime def A045309_gen(): # generator of terms yield from (2,3) yield from filter(isprime, count(5,6)) A045309_list = list(islice(A045309_gen(),48)) # Chai Wah Wu, Apr 28 2025
For n = 2, the first two primes 2 and 3 each have three cube roots mod 307: 2 has cube roots 52, 270, 292 and 3 has cube roots 192, 194, 228. - _Robert Israel_, Aug 02 2016
Primes:= [2]: pp:= 7: for n from 1 to 12 do for p from pp by 6 while not(isprime(p) and andmap(t -> t &^ ((p-1)/3) mod p = 1, Primes)) do od: A[n]:= p; pp:= p; Primes:= [op(Primes), nextprime(Primes[-1])]; od: seq(A[i],i=1..12); # Robert Israel, Aug 02 2016
(* This naive program being very slow, limit is set to 8 terms *) lim=8; np[] := While[p=NextPrime[p]; Mod[p,3]!=1]; crQ[n_, p_] := Reduce[ 0A002225={}; While[Length[A002225] < lim, If[And @@ (crQ[#,p]& /@ pp), AppendTo[pp, NextPrime[ Last[pp]]]; Print[p]; AppendTo[A002225,p], np[] ] ]; A002225 (* Jean-François Alcover, Sep 09 2011 *)
p6 = Select[6*Range[0, 400]+1, PrimeQ]; Select[p6, (Reduce[3 == k^3+m*#, {k, m}, Integers] =!= False)&] (* Jean-François Alcover, Feb 20 2014 *)
forprime(p=1, 9999, p%6==1&&ispower(Mod(3, p), 3)&&print1(p", ")) \\ M. F. Hasler, Feb 18 2014
is_A014753(p)={p%6==1&&ispower(Mod(3, p), 3)&&isprime(p)} \\ M. F. Hasler, Feb 18 2014
a(3) = 57, since 109 is the third term of A014752, 57, 58 and 103 are the solutions mod 109 of x^3 = 2 and 57 is the least one.
a(3) = 58, since 109 is the third term of A014752 and 58 is the second solution mod 109 of x^3 = 2.
a(3) = 103, since 109 is the third term of A014752 and 103 is the third solution mod 109 of x^3 = 2.
Comments