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

A136117 Pentagonal numbers (A000326) which are the sum of 2 other positive pentagonal numbers.

Original entry on oeis.org

70, 92, 852, 925, 1247, 1426, 1926, 2625, 3577, 5192, 6305, 6501, 7107, 7740, 7957, 8177, 8626, 9560, 10292, 12927, 13207, 14652, 15555, 16172, 18095, 20475, 20827, 21901, 22265, 22632, 23002, 23751, 24130, 28497, 29330, 31032, 33227, 33675
Offset: 1

Views

Author

M. F. Hasler, Dec 15 2007; corrected Dec 25 2007

Keywords

Comments

It is conjectured that every integer and hence every pentagonal number, greater than 33066, hence greater than A000326(149) = 33227, can be represented as the sum of three pentagonal numbers. - Jonathan Vos Post, Dec 18 2007

Examples

			a(1)=70=P(7) is the least pentagonal number which can be written as sum of two other pentagonal numbers, P(7)=P(5)+P(5).

Crossrefs

Cf. A000326, A136112-A136118, A007527.

Programs

  • PARI
    P(n)=n*(3*n-1)>>1 /* a.k.a. A000326 */
isPent(t)=P(sqrtint(t<<1\3)+1)==t
for(i=1,299,for(j=1,(i+1)\sqrt(2),isPent(P(i)-P(j)) && print1(P(i)",") || next(2)))
/* The following is much faster, at the cost of implementing sum2sqr(), cf. A133388*/
A136117next(i)=i=sqrtint(i\3*2)*6+5; until(0, for(j=2,#t=sum2sqr((i+=6)^2+1), t[j]%6==[5,5] && break(2)));i^2\24
A136117vect(n,i)=vector(n,j,i=A136117next(i)) /* 2nd arg =0 by default but allows one to start elsewhere */
A136117(n,i)=until(!n--,i=A136117next(i));i \\ M. F. Hasler, Dec 25 2007

Formula

a(n) = A000326(A136116(n)) = A000326(m)+A136114(m) where m is the index of the n-th nonzero term in A136114 or A136115.

A164579 Integer averages of halves of first cubes of natural numbers (n^3)/2 for some n.

Original entry on oeis.org

6, 56, 81, 198, 480, 578, 950, 1656, 1875, 2646, 3968, 4356, 5670, 7800, 8405, 10406, 13536, 14406, 17238, 21560, 22743, 26550, 32256, 33800, 38726, 46008, 47961, 54150, 63200, 65610, 73206, 84216, 87131, 96278, 109440, 112908, 123750, 139256
Offset: 1

Views

Author

Vladimir Joseph Stephan Orlovsky, Aug 16 2009

Keywords

Comments

Also, integers of the form (1/8)*n*(n+1)^2 for some n. - Zak Seidov, Aug 17 2009

Examples

			1/2, 9/4, 6, 25/2, 45/2, 147/4, 56, 81, ...

Links

Crossrefs

Cf. A050248, A051456, A078617, A078618, A136116, A154293, A164576, A164577, A164578

Programs

  • Mathematica
    s=0;lst={};Do[a=(s+=(n^3)/2)/n;If[Mod[a,1]==0,AppendTo[lst,a]],{n,3*5!}];lst
LinearRecurrence[{1,0,3,-3,0,-3,3,0,1,-1},{6,56,81,198,480,578,950,1656,1875,2646},40] (* Harvey P. Dale, Jul 26 2017 *)
Module[{nn=200,ac},ac=Accumulate[Range[nn]^3/2];Select[#[[1]]/#[[2]]&/@ Thread[{ac,Range[nn]}],IntegerQ]] (* Harvey P. Dale, Jan 28 2020 *)
  • PARI
    forstep(n=3, 150, [4,1,3], print1(n*(n+1)^2>>3, ", ")); \\ Charles R Greathouse IV, Nov 02 2009

Formula

G.f.: ( x*(6+50*x+25*x^2+99*x^3+132*x^4+23*x^5+39*x^6+10*x^7) ) / ( (1+x+x^2)^3*(x-1)^4 ). - R. J. Mathar, Jan 25 2011

A164619 Integers of the form A164577(k)/3.

Original entry on oeis.org

4, 15, 54, 75, 132, 169, 320, 459, 735, 847, 1104, 1250, 1764, 2175, 2904, 3179, 3780, 4107, 5200, 6027, 7425, 7935, 9024, 9604, 11492, 12879, 15162, 15979, 17700, 18605, 21504, 23595, 26979, 28175, 30672, 31974, 36100, 39039, 43740, 45387, 48804
Offset: 1

Views

Author

Vladimir Joseph Stephan Orlovsky, Aug 17 2009

Keywords

Comments

The sequence members are the third of the average of a set of smallest cubes, if integer.

Examples

			A third of the average of the first cube, A164577(1)/3=1/3, is not an integer and does not contribute to the sequence.
A third of the average of the first two cubes, A164577(2)/3=4, is an integer and defines a(1)=4 of the sequence.

Links

Crossrefs

Cf. A050248, A051456, A078617, A078618, A136116, A154293, A164286, A164576, A164577, A164578, A164579.

Programs

  • Mathematica
    s=0;lst={};Do[a=(s+=(n^3)/3)/n;If[Mod[a,1]==0,AppendTo[lst,a]],{n,2*5!}]; lst
LinearRecurrence[{2,-1,-1,2,-1,2,-4,2,2,-4,2,-1,2,-1,-1,2,-1},{4,15,54,75,132,169,320,459,735,847,1104,1250,1764,2175,2904,3179,3780},50] (* Harvey P. Dale, Apr 06 2016 *)
  • PARI
    Vec(x*(x^14 +x^13 +16*x^12 +10*x^11 +47*x^10 -22*x^9 +61*x^8 +10*x^7 +88*x^6 +8*x^5 +43*x^4 -14*x^3 +28*x^2 +7*x +4) / ((x -1)^4*(x +1)^3*(x^2 -x +1)^3*(x^2 +x +1)^2) + O(x^100)) \\ Colin Barker, Oct 27 2014

Formula

a(n) = +2*a(n-1) -a(n-2) -a(n-3) +2*a(n-4) -a(n-5) +2*a(n-6) -4*a(n-7) +2*a(n-8) +2*a(n-9) -4*a(n-10) +2*a(n-11) -a(n-12) +2*a(n-13) -a(n-14) -a(n-15) +2*a(n-16) -a(n-17). - R. J. Mathar, Jan 25 2011
G.f.: x*(x^14 +x^13 +16*x^12 +10*x^11 +47*x^10 -22*x^9 +61*x^8 +10*x^7 +88*x^6 +8*x^5 +43*x^4 -14*x^3 +28*x^2 +7*x +4) / ((x -1)^4*(x +1)^3*(x^2 -x +1)^3*(x^2 +x +1)^2). - Colin Barker, Oct 27 2014

Extensions

Edited by R. J. Mathar, Aug 20 2009
Showing 1-3 of 3 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.