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.

Previous Showing 41-50 of 207 results. Next

A031507 a(n) = smallest k>0 such that the elliptic curve y^2 = x^3 + k has rank n, or -1 if no such k exists.

Original entry on oeis.org

1, 2, 15, 113, 2089, 66265, 1358556
Offset: 0

Views

Author

Keywords

Comments

See A031508 for the smallest negative k. - Artur Jasinski, Nov 21 2011
See A060950 for the rank of y^2 = x^3 + n. - Jonathan Sondow, Sep 10 2013
Gebel, Pethö, & Zimmer: "One experimental observation derived from the tables is that the rank r of Mordell's curves grows according to r = O(log |k|/|log log |k||^(2/3))." Hence this fit suggests a(n) >> exp(n (log n)^(1/3)) where >> is the Vinogradov symbol. - Charles R Greathouse IV, Sep 10 2013
The curves for k and -27*k are isogenous (as Noam Elkies points out---see Womack), so they have the same rank. - Jonathan Sondow, Sep 10 2013
Womack (2003) gives further upper bounds: a(7) <= 47550317, a(8) <= 1632201497, a(9) <= 185418133372, a(10) <= 68513487607153. - M. F. Hasler, Jul 01 2024
The three questions for arbitrary k, positive k, and negative k are not very far from each other because the curves for k and -27k are related by a 3-isogeny and therefore have the same rank. It would be most natural to ask for the minimal |k| for k of either sign [see A373795]. - Noam D. Elkies, Jul 02 2024
a(16) <= 1160221354461565256631205207888 (Elkies, ANTS-XVI, 2024). The same article also establishes the existence of a value of k which has rank >= 17. - N. J. A. Sloane, Jul 05 2024

Examples

			a(12) <= 27*A031508(12) <= 27*6533891544658786928 = 176415071705787247056 (from Quer 1987 and Womack). - _Jonathan Sondow_, Sep 10 2013
		

References

  • Noam D. Elkies, Rank of an elliptic curve and 3-rank of a quadratic field via the Burgess bounds, 2024 Algorithmic Number Theory Symposium, ANTS-XVI, MIT, July 2024.

Crossrefs

Programs

  • PARI
    {A031507(n)=for(k=1, oo, ellrank(ellinit([0, k]))[1]==n && return(k))} \\ Use ellanalyticrank() for PARI version < 2.14. - M. F. Hasler, Jul 01 2024

Formula

a(n) <= 27*A031508(n) and A031508(n) <= 27*a(n). - Jonathan Sondow, Sep 10 2013

Extensions

Definition clarified by Jonathan Sondow, Oct 26 2013
Escape clause added to definition by N. J. A. Sloane, Jun 29 2024, because, as John Cremona reminds me, it is not known if k always exists.

A060951 Rank of elliptic curve y^2 = x^3 - n.

Original entry on oeis.org

0, 1, 0, 1, 0, 0, 1, 0, 0, 0, 2, 0, 1, 0, 1, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 2, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 2, 1, 0, 0, 1, 1, 1, 0, 2, 1, 1, 1, 1, 0, 2, 1, 1, 1, 1, 1, 1, 1, 2, 0, 1, 0, 1, 1, 2, 0, 0, 0, 1, 1, 0, 1, 1, 2, 0, 0, 1, 0, 1, 0, 2, 1, 1, 0, 1, 0, 2
Offset: 1

Views

Author

N. J. A. Sloane, May 10 2001

Keywords

Comments

The curves for n and -27*n are isogenous (as Noam Elkies points out--see Womack), so they have the same rank. - Jonathan Sondow, Sep 10 2013

Examples

			a(1) = A060950(27) = a(729) = 0. - _Jonathan Sondow_, Sep 10 2013
		

Crossrefs

Cf. A081120 (number of integral solutions to Mordell's equation y^2 = x^3 - n).

Programs

  • PARI
    {a(n) = if( n<1, 0, length( ellgenerators( ellinit( [ 0, 0, 0, 0, -n], 1))))} /* Michael Somos, Mar 17 2011 */
    
  • PARI
    apply( {A060951(n)=ellrank(ellinit([0,-n]))[1]}, [1..99]) \\ For version < 2.14, use ellanalyticrank(...). - M. F. Hasler, Jul 01 2024

Formula

a(n) = A060950(27*n) and A060950(n) = a(27*n), so a(n) = a(729*n). - Jonathan Sondow, Sep 10 2013

Extensions

Corrected Apr 08 2005 at the suggestion of James R. Buddenhagen. There were errors caused by the fact that Mishima lists each curve of rank two twice, once for each generator.

A002153 Numbers k for which the rank of the elliptic curve y^2 = x^3 + k is 1.

Original entry on oeis.org

2, 3, 5, 8, 9, 10, 11, 12, 18, 19, 22, 26, 28, 30, 31, 33, 35, 36, 38, 39, 40, 41, 44, 46, 47, 48, 50, 52, 54, 55, 56, 58, 61, 62, 66, 67, 68, 69, 71, 72, 74, 76, 77, 80, 82, 83, 91, 92, 94, 97, 98, 99, 100, 102, 103, 105, 107, 108, 110, 111, 112, 117, 118, 119
Offset: 1

Views

Author

Keywords

References

  • 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

  • Magma
    for k in[1..200] do if Rank(EllipticCurve([0,0,0,0,k])) eq 1 then print k; end if; end for; // Vaclav Kotesovec, Jul 07 2019

Extensions

Corrected and extended by James R. Buddenhagen, Feb 18 2005

A008881 a(n) = Product_{j=0..5} floor((n+j)/6).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 1, 2, 4, 8, 16, 32, 64, 96, 144, 216, 324, 486, 729, 972, 1296, 1728, 2304, 3072, 4096, 5120, 6400, 8000, 10000, 12500, 15625, 18750, 22500, 27000, 32400, 38880, 46656, 54432, 63504, 74088, 86436, 100842, 117649, 134456, 153664, 175616, 200704
Offset: 0

Views

Author

Keywords

Comments

For n >= 6, a(n) is the maximal product of 6 positive integers with sum n. - Wesley Ivan Hurt, Jun 29 2022
The maximal product of k positive variables when their sum is equal to s is obtained when each term = s/k; hence, a(6m) = m^6 (A001014). - Bernard Schott, Jul 28 2022

Crossrefs

Maximal product of k positive integers with sum n, for k = 2..10: A002620 (k=2), A006501 (k=3), A008233 (k=4), A008382 (k=5), this sequence (k=6), A009641 (k=7), A009694 (k=8), A009714 (k=9), A354600 (k=10).
Cf. A001014 (6th power), A008588 (multiples of 6), A013664.

Programs

  • GAP
    List([0..50], n-> Product([0..5], j-> Int((n+j)/6))); # G. C. Greubel, Sep 13 2019
  • Magma
    [(&*[Floor((n+j)/6): j in [0..5]]): n in [0..50]]; // G. C. Greubel, Sep 13 2019
    
  • Maple
    seq( mul( floor((n+i)/6), i=0..5 ), n=0..80);
  • Mathematica
    Product[Floor[(Range[51]+j-2)/6], {j,6}] (* G. C. Greubel, Sep 13 2019 *)
  • PARI
    vector(50, n, prod(j=0,5, (n+j)\6) ) \\ G. C. Greubel, Sep 13 2019
    
  • Sage
    [product(floor((n+j)/6) for j in (0..5)) for n in (0..50)] # G. C. Greubel, Sep 13 2019
    

Formula

Sum_{n>=6} 1/a(n) = 1 + zeta(6). - Amiram Eldar, Jan 10 2023

A017490 a(n) = (11*n + 8)^6.

Original entry on oeis.org

262144, 47045881, 729000000, 4750104241, 19770609664, 62523502209, 164206490176, 377149515625, 782757789696, 1500730351849, 2699554153024, 4608273662721, 7529536000000, 11853911588401, 18075490334784, 26808753332089, 38806720086016, 54980371265625
Offset: 0

Views

Author

N. J. A. Sloane, Dec 11 1996

Keywords

Crossrefs

Powers of the form (11*n+8)^m: A017485 (m=1), A017486 (m=2), A017487 (m=3), A017488 (m=4), A017489 (m=5), this sequence (m=6), A017491 (m=7), A017492 (m=8), A017493 (m=9), A017494 (m=10), A017495 (m=11), A017496 (m=12).

Programs

  • GAP
    List([0..20], n-> (11*n+8)^6); # G. C. Greubel, Sep 22 2019
  • Magma
    [(11*n+8)^6: n in [0..20]]; // Vincenzo Librandi, Sep 04 2011
    
  • Maple
    A017490:=n->(11*n+8)^6; seq(A017490(n), n=0..20); # Wesley Ivan Hurt, May 21 2014
  • Mathematica
    (11*Range[0,20]+8)^6 (* or *) LinearRecurrence[{7,-21,35,-35,21,-7,1}, {262144, 47045881, 729000000, 4750104241, 19770609664, 62523502209, 164206490176}, 20] (* Harvey P. Dale, Nov 08 2013 *)
  • Maxima
    makelist( (11*n+8)^6, n, 0, 20); /* Martin Ettl, Oct 21 2012 */
    
  • PARI
    vector(20, n, (11*n-3)^6) \\ G. C. Greubel, Sep 22 2019
    
  • Sage
    [(11*n+8)^6 for n in (0..20)] # G. C. Greubel, Sep 22 2019
    

Formula

a(0)=262144, a(1)=47045881, a(2)=729000000, a(3)=4750104241, a(4)=19770609664, a(5)=62523502209, a(6)=164206490176, a(n) = 7*a(n-1) - 21*a(n-2) + 35*a(n-3) - 35*a(n-4) + 21*a(n-5) - 7*a(n-6) + a(n-7). - Harvey P. Dale, Nov 08 2013
a(n) = A001014(A017485(n)). - Wesley Ivan Hurt, May 21 2014
From G. C. Greubel, Sep 22 2019: (Start)
G.f.: (262144 +45210873*x +405183857*x^2 +625892702*x^3 +191449182*x^4 +7524433*x^5 +729*x^6)/(1-x)^7.
E.g.f.: (262144 +46783737*x +317585191*x^2 +450663290*x^3 +206511305*x^4 + 34303863*x^5 +1771561*x^6)*exp(x). (End)

A017502 a(n) = (11*n + 9)^6.

Original entry on oeis.org

531441, 64000000, 887503681, 5489031744, 22164361129, 68719476736, 177978515625, 404567235136, 832972004929, 1586874322944, 2839760855281, 4826809000000, 7858047974841, 12332795428864, 18755369578009
Offset: 0

Views

Author

Keywords

Crossrefs

Powers of the form (11*n+9)^m: A017497 (m=1), A017498 (m=2), A017499 (m=3), A017500 (m=4), A017501 (m=5), this sequence (m=6), A017503 (m=7), A017504 (m=8), A017505 (m=9), A017506 (m=10), A017607 (m=11), A017508 (m=12).
Subsequence of A001014.

Programs

  • GAP
    List([0..20], n-> (11*n+9)^6); # G. C. Greubel, Oct 28 2019
  • Magma
    [(11*n+9)^6: n in [0..20]]; // G. C. Greubel, Oct 28 2019
    
  • Maple
    seq((11*n+9)^6, n=0..20); # G. C. Greubel, Oct 28 2019
  • Mathematica
    (11Range[0,20]+9)^6 (* or *) LinearRecurrence[{7,-21,35,-35,21,-7,1}, {531441,64000000,887503681,5489031744,22164361129,68719476736, 177978515625}, 20] (* Harvey P. Dale, Dec 06 2018 *)
  • Maxima
    makelist((11*n+9)^6, n, 0, 30); /* Martin Ettl, Oct 21 2012 */
    
  • PARI
    vector(21, n, (11*n-2)^6) \\ G. C. Greubel, Oct 28 2019
    
  • Sage
    [(11*n+9)^6 for n in (0..20)] # G. C. Greubel, Oct 28 2019
    

Formula

From G. C. Greubel, Oct 28 2019: (Start)
G.f.: (531441 + 60279913*x + 450663942*x^2 + 601905542*x^3 + 157316657*x^4 + 4826361*x^5 + 64*x^6)/(1-x)^7.
E.g.f.: (531441 + 63468559*x + 380017561*x^2 + 502998210*x^3 + 219907820*x^4 + 35270169*x^5 + 1771561*x^6)*exp(x). (End)

A273429 Number of ordered ways to write n as x^6 + y^2 + z^2 + w^2, where x,y,z,w are nonnegative integers with y <= z <= w.

Original entry on oeis.org

1, 2, 2, 2, 2, 2, 2, 1, 1, 3, 3, 2, 2, 2, 2, 1, 1, 3, 4, 3, 2, 2, 2, 1, 1, 3, 4, 4, 2, 2, 3, 1, 1, 3, 4, 3, 3, 3, 3, 2, 1, 4, 4, 2, 2, 3, 3, 1, 1, 3, 5, 5, 3, 3, 5, 3, 1, 3, 3, 3, 2, 2, 4, 2, 2, 5, 7, 5, 4, 5, 4, 1, 3, 6, 6, 6, 4, 4, 4, 1, 2
Offset: 0

Views

Author

Zhi-Wei Sun, May 22 2016

Keywords

Comments

The author proved in arXiv:1604.06723 that for each c = 1, 4 any natural number can be written as c*x^6 + y^2 + z^2 + w^2 with x,y,z,w nonnegative integers. Thus a(n) > 0 for all n = 0,1,2,....
We note that a(n) = 1 for the following values of n not divisible by 2^6: 7, 8, 15, 16, 23, 24, 31, 32, 40, 47, 48, 56, 71, 79, 92, 112, 143, 176, 191, 240, 304, 368, 560, 624, 688, 752, 1072, 1136, 1456, 1520, 1840, 1904, 2608, 2672, 3760, 3824, 6512, 6896.
For more conjectural refinements of Lagrange's four-square theorem, one may consult the author's preprint arXiv:1604.06723.

Examples

			a(7) = 1 since 7 = 1^6 + 1^2 + 1^2 + 2^2 with 1 = 1 < 2.
a(8) = 1 since 8 = 0^6 + 0^2 + 2^2 + 2^2 with 0 < 2 = 2.
a(15) = 1 since 15 = 1^6 + 1^2 + 2^2 + 3^2 with 1 < 2 < 3.
a(16) = 1 since 16 = 0^6 + 0^2 + 0^2 + 4^2 with 0 = 0 < 4.
a(56) = 1 since 56 = 0^6 + 2^2 + 4^2 + 6^2 with 2 < 4 < 6.
a(71) = 1 since 71 = 1^6 + 3^2 + 5^2 + 6^2 with 3 < 5 < 6.
a(79) = 1 since 79 = 1^6 + 2^2 + 5^2 + 7^2 with 2 < 5 < 7.
a(92) = 1 since 92 = 1^6 + 1^2 + 3^2 + 9^2 with 1 < 3 < 9.
a(143) = 1 since 143 = 1^6 + 5^2 + 6^2 + 9^2 with 5 < 6 < 9.
a(191) = 1 since 191 = 1^6 + 3^2 + 9^2 + 10^2 with 3 < 9 < 10.
a(624) = 1 since 624 = 2^6 + 4^2 + 12^2 + 20^2 with 4 < 12 < 20.
a(2672) = 1 since 2672 = 2^6 + 4^2 + 36^2 + 36^2 with 4 < 36 = 36.
a(3760) = 1 since 3760 = 0^6 + 4^2 + 12^2 + 60^2 with 4 < 12 < 60.
a(3824) = 1 since 3824 = 2^6 + 4^2 + 12^2 + 60^2 with 4 < 12 < 60.
a(6512) = 1 since 6512 = 2^6 + 12^2 + 52^2 + 60^2 with 12 < 52 < 60.
a(6896) = 1 since 6896 = 2^6 + 36^2 + 44^2 + 60^2 with 36 < 44 < 60.
		

Crossrefs

Programs

  • Mathematica
    SQ[n_]:=SQ[n]=IntegerQ[Sqrt[n]]
    Do[r=0;Do[If[SQ[n-x^6-y^2-z^2],r=r+1],{x,0,n^(1/6)},{y,0,Sqrt[(n-x^6)/3]},{z,y,Sqrt[(n-x^6-y^2)/2]}];Print[n," ",r];Continue,{n,0,80}]

A101104 a(1)=1, a(2)=12, a(3)=23, and a(n)=24 for n>=4.

Original entry on oeis.org

1, 12, 23, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24
Offset: 1

Views

Author

Cecilia Rossiter, Dec 15 2004

Keywords

Comments

Original name: The first summation of row 4 of Euler's triangle - a row that will recursively accumulate to the power of 4.

Crossrefs

For other sequences based upon MagicNKZ(n,k,z):
..... | n = 1 | n = 2 | n = 3 | n = 4 | n = 5 | n = 6 | n = 7
---------------------------------------------------------------------------
z = 0 | A000007 | A019590 | .......MagicNKZ(n,k,0) = A008292(n,k+1) .......
z = 1 | A000012 | A040000 | A101101 | thisSeq | A101100 | ....... | .......
z = 2 | A000027 | A005408 | A008458 | A101103 | A101095 | ....... | .......
z = 3 | A000217 | A000290 | A003215 | A005914 | A101096 | ....... | .......
z = 4 | A000292 | A000330 | A000578 | A005917 | A101098 | ....... | .......
z = 5 | A000332 | A002415 | A000537 | A000583 | A022521 | ....... | A255181
Cf. A101095 for an expanded table and more about MagicNKZ.

Programs

  • Mathematica
    MagicNKZ = Sum[(-1)^j*Binomial[n+1-z, j]*(k-j+1)^n, {j, 0, k+1}];Table[MagicNKZ, {n, 4, 4}, {z, 1, 1}, {k, 0, 34}]
    Join[{1, 12, 23},LinearRecurrence[{1},{24},56]] (* Ray Chandler, Sep 23 2015 *)

Formula

a(k) = MagicNKZ(4,k,1) where MagicNKZ(n,k,z) = Sum_{j=0..k+1} (-1)^j*binomial(n+1-z,j)*(k-j+1)^n (cf. A101095). That is, a(k) = Sum_{j=0..k+1} (-1)^j*binomial(4, j)*(k-j+1)^4.
a(1)=1, a(2)=12, a(3)=23, and a(n)=24 for n>=4. - Joerg Arndt, Nov 30 2014
G.f.: x*(1+11*x+11*x^2+x^3)/(1-x). - Colin Barker, Apr 16 2012

Extensions

New name from Joerg Arndt, Nov 30 2014
Original Formula edited and Crossrefs table added by Danny Rorabaugh, Apr 22 2015

A102833 Numbers n for which rank of the elliptic curve y^2=x^3+n is 3.

Original entry on oeis.org

113, 141, 316, 346, 359, 427, 443, 506, 537, 568, 659, 681, 730, 745, 873, 892, 899, 940, 997, 1016, 1025, 1090, 1149, 1157, 1171, 1213, 1304, 1305, 1342, 1367, 1373, 1478, 1522, 1639, 1646, 1737, 1753, 1772, 1811, 1841, 1897, 1907, 1954, 2024, 2143
Offset: 1

Views

Author

James R. Buddenhagen, Feb 18 2005. Entry revised by N. J. A. Sloane, Jun 10 2012

Keywords

Crossrefs

Programs

  • Magma
    for k in[1..2000] do if Rank(EllipticCurve([0,0,0,0,k])) eq 3 then print k; end if; end for; // Vaclav Kotesovec, Jul 07 2019

Extensions

More terms from T. D. Noe, Jul 24 2007

A123866 a(n) = n^6 - 1.

Original entry on oeis.org

0, 63, 728, 4095, 15624, 46655, 117648, 262143, 531440, 999999, 1771560, 2985983, 4826808, 7529535, 11390624, 16777215, 24137568, 34012223, 47045880, 63999999, 85766120, 113379903, 148035888, 191102975, 244140624, 308915775, 387420488
Offset: 1

Views

Author

Reinhard Zumkeller, Oct 16 2006

Keywords

Comments

a(n) mod 7 = 0 iff n mod 7 > 0: a(A008589(n))=6; a(A047304(n)) = 0; a(n) mod 7 = 6*(1-A082784(n)).
a(n) = A005563(n-1)*A059826(n) = A068601(n)*A001093(n). - Reinhard Zumkeller, Feb 02 2007

Crossrefs

Programs

Formula

G.f.: x^2*(63 + 287*x + 322*x^2 + 42*x^3 + 7*x^4 - x^5)/(1-x)^7. - Colin Barker, May 08 2012
a(n) = 7*a(n-1) - 21*a(n-2) + 35*a(n-3) - 35*a(n-4) + 21*a(n-5) - 7*a(n-6) + a(n-7); a(1)=0, a(2)=63, a(3)=728, a(4)=4095, a(5)=15624, a(6)=46655, a(7)=117648. - Harvey P. Dale, Nov 18 2012
Sum_{n>=2} 1/a(n) = 11/12 - Pi*sqrt(3)*tanh(Pi*sqrt(3)/2)/6. - Vaclav Kotesovec, Feb 14 2015
E.g.f.: 1 + (-1 + x + 31*x^2 + 90*x^3 + 65*x^4 + 15*x^5 + x^6)*exp(x). - G. C. Greubel, Aug 08 2019
Product_{n>=2} (1 + 1/a(n)) = 6*Pi^2*sech(sqrt(3)*Pi/2)^2. - Amiram Eldar, Jan 20 2021
Previous Showing 41-50 of 207 results. Next