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-2 of 2 results.

A002338 x such that p = (x^2 + 27*y^2)/4, where p is the n-th prime of the form 3k+1.

Original entry on oeis.org

1, 5, 7, 4, 11, 8, 1, 5, 7, 17, 19, 13, 2, 20, 23, 19, 14, 25, 7, 23, 11, 13, 28, 22, 17, 29, 26, 32, 16, 35, 1, 5, 37, 35, 13, 29, 34, 31, 19, 2, 28, 10, 23, 25, 32, 43, 29, 1, 31, 11, 26, 49, 47, 17, 43, 40, 49, 37, 8, 53, 44, 50, 16, 41, 29, 49, 31, 56, 5, 7, 35, 13, 59, 47, 19, 52, 61, 41, 61, 10, 43, 14, 53, 59, 64, 65, 62, 55, 22, 65, 35, 67, 7
Offset: 1

Views

Author

Keywords

Comments

A123489 is a signed version. - Michael Somos, Aug 27 2012

References

  • A. J. C. Cunningham, Quadratic Partitions. Hodgson, London, 1904, p. 1.
  • B. Engquist and Wilfried Schmid, Mathematics Unlimited - 2001 and Beyond, Chapter on Error-correcting codes and curves over finite fields, see pp. 1118-1119. [From Neven Juric, Oct 16 2008.]
  • D. H. Lehmer, Guide to Tables in the Theory of Numbers. Bulletin No. 105, National Research Council, Washington, DC, 1941, p. 55.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Mathematica
    Reap[For[p = 2, p<2000, p = NextPrime[p], For[x = 1, x <= Floor[2Sqrt[p]], x++, px = 4p - x^2; If[Mod[px, 27] == 0, If[IntegerQ[Sqrt[px/27]], Sow[x] ]]]]][[2, 1]] (* Jean-François Alcover, Sep 06 2018, after Ruperto Corso *)
  • PARI
    forprime(p=2,10000,for(x=1,floor(2*sqrt(p)),px=4*p-x^2;if(px%27==0,if(issquare(px/27,&y),print1(x","))))) /* Ruperto Corso, Dec 14 2011 */

Extensions

Corrected and extended by Ruperto Corso, Dec 14 2011

A123489 a(n) = Sum_{k=0..p-1} Kronecker(4k^3+1, p) where p is the n-th prime of the form 6k+1.

Original entry on oeis.org

1, -5, 7, 4, -11, -8, 1, -5, 7, -17, 19, 13, -2, -20, -23, 19, -14, 25, 7, -23, -11, 13, 28, 22, -17, -29, -26, -32, 16, -35, 1, -5, 37, -35, 13, -29, 34, 31, 19, -2, 28, 10, -23, 25, -32, 43, -29, 1, 31, -11, -26, 49, -47, -17, 43, 40, 49, 37, -8, -53, -44, -50, 16, -41, -29, 49, 31, -56, -5, 7, -35, 13, -59, -47
Offset: 1

Views

Author

Michael Somos, Sep 30 2006

Keywords

Comments

Given a prime p == 1 (mod 6), the sum x is the unique solution to 4*p = x^2 + 27*y^2 where x == 1 (mod 3) and y is a positive integer.
Given a prime p == 1 (mod 6), the number of solutions to u^3 + v^3 == 1 (mod p) is p - 2 + x. - Michael Somos, Jul 11 2022

Examples

			If p = 37, then 4*37 = (-11)^2 + 27*(1)^2 where -11 = Sum_{k=0..36} Kronecker(4k^3+1, 37) and 37 is the 5th prime of the form 6k+1 so a(5) = -11.
If p = 13, then 4*p = x^2 + 26*y^2 where x = -5, y = 1. The solutions to u^3 + v^3 == 1 (mod p) is {(0,1), (1,0), (0,3), (3,0), (0,9), (9,0)} with cardinality = 6 = 13 - 2 + (-5). - _Michael Somos_, Jul 11 2022
		

References

  • L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 3, p. 55.

Crossrefs

Cf. A002338 is the unsigned version, A002339 is the y.
Cf. A002476 (primes of the form 6k+1).

Programs

  • PARI
    {a(n) = my(p, c); if(n<1, 0, c=0; p=0; while(c
    				
Showing 1-2 of 2 results.