login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A084165
Primes which are 1 mod m, where m is the index of the prime in sequence A002313 (Real primes with corresponding complex primes). The index m can be found in A084166 Primes which are -1 mod m can be found in sequence A084163.
3
5, 13, 17, 37, 89, 97, 181, 2689, 2969, 4621, 7457, 8081, 8161, 36709, 62701, 169489, 169709, 169753, 282809, 770101, 5763577, 9491101, 9491281, 9495121, 42544261, 115195501, 189689041, 189689653, 312315373, 312316409, 2294883817
OFFSET
1,1
COMMENTS
Real primes 2,5,13,17,29,37,... have a unique representation as sum of two squares. Values larger 2 are the primes p with p = 1 mod 4. This is sequence A002313. If p = x^2 + y^2, the corresponding complex prime is x+y*i First complex prime is 1+i with 2 as corresponding real prime, according to reference, page 1-2.
REFERENCES
Handbook of First Complex Prime Numbers, Part1 + 2 Ervand Kogbetliantz and Alice Krikorian, Gordon and Breach, 1971
EXAMPLE
89 is the 11th prime in sequence A002313, 11*8 = 88, so 89 = 1 mod 11
MATHEMATICA
Module[{nn=112*10^6, pr, len}, pr=Select[Prime[Range[nn]], MemberQ[ {1, 2}, Mod[ #, 4]]&]; len=Length[pr]; Select[Thread[{pr, Range[len]}], Mod[ #[[1]], #[[2]]] == 1&]][[All, 1]] (* Harvey P. Dale, Aug 13 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Sven Simon, May 17 2003
STATUS
approved