This site is supported by donations to The OEIS Foundation. Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A267859 The p-defect p - N(p) of the elliptic curve y^2 = x^3 + x for primes p congruent to 1 modulo 4 (A002144). 4
 2, -6, 2, 10, 2, 10, -14, 10, -6, 10, 18, 2, -6, -14, -22, -14, -22, 26, 18, -14, 2, -30, 26, -30, 2, 26, 18, 10, 34, 26, -22, 18, 10, 34, -14, 34, -38, 2, -6, -30, 34, -14, 42, -38, 10, -22, 42, -38, 26, 2, -46, 10, 34, -38, 50, 26, 50, -46, 2, 10, -30, -54, 18, -38, 50, 34, -22, 10, 50, -54 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,1 COMMENTS See A002172 for a differently signed sequence. The number N(p) of solutions modulo a prime p of the elliptic curve y^2 = x^3 + x (of discriminant -4) is given for all p in A095978. The p-defect a_p = p - N(p) for prime 2 and primes congruent to 3 modulo 4 vanishes. A002144(n) - (a(n)/2)^2  = (2*A002973(n))^2, n >= 1. See the formula for A095978 for primes 1 (mod 4). This sequence gives also the non-vanishing p-defects of the elliptic curve y^2 = x^3 - 4*x. See a comment on A138515 with the Martin and Ono link for the modularity series for these two elliptic curves. - Wolfdieter Lang, May 26 2016 REFERENCES J. H. Silverman, A Friendly Introduction to Number Theory, 3rd ed., Pearson Education, Inc, 2006, p. 398. In the 4th ed., 2014, p. 371. LINKS Seiichi Manyama, Table of n, a(n) for n = 1..1000 FORMULA a(n) = A002144(n) - A095978(m) with A002144(n) = A000040(m), n >= 1. a(n) = A138515((A002144(n) - 1)/4), n >= 1. - Wolfdieter Lang, May 26 2016 EXAMPLE n = 2: p = A002144(2) = 13 = A000040(6), m = 6, a(2) = 13 - A095978(6) = 13 - 19  = -6. n = 2:  -6 = A138515((A002144(2) - 1)/4) = A138515(3) = -6. - Wolfdieter Lang, May 26 2016 MATHEMATICA terms = 100; A002144 = Select[Range[5, 20*terms, 4], PrimeQ]; A095978[n_] := Module[{p, xy, x}, p = Prime[n]; If[n==1 || Mod[p, 4]==3, Return[p]]; xy = {Re[#], Im[#]}& @ FactorInteger[p, GaussianIntegers -> True][[2, 1]]; x = SelectFirst[xy, OddQ]; If[Mod[x, 4]==1, p - 2*x, p + 2*x]]; a[n_] := (p = A002144[[n]]; m = PrimePi[p]; p - A095978[m]); Array[a, terms] (* Jean-François Alcover, Feb 26 2016, after Robert Israel (A095978) *) CROSSREFS Cf. A000040, A002144, A002172, A095978, A138515. Sequence in context: A153190 A010240 A199154 * A002172 A126289 A200044 Adjacent sequences:  A267856 A267857 A267858 * A267860 A267861 A267862 KEYWORD sign AUTHOR Wolfdieter Lang, Feb 06 2016 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified June 17 14:09 EDT 2019. Contains 324185 sequences. (Running on oeis4.)