login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A085787 Generalized heptagonal numbers: n*(5*n-3)/2, n=0, +-1, +-2 +-3, ... 65

%I

%S 0,1,4,7,13,18,27,34,46,55,70,81,99,112,133,148,172,189,216,235,265,

%T 286,319,342,378,403,442,469,511,540,585,616,664,697,748,783,837,874,

%U 931,970,1030,1071,1134,1177,1243,1288,1357,1404,1476,1525,1600,1651,1729

%N Generalized heptagonal numbers: n*(5*n-3)/2, n=0, +-1, +-2 +-3, ...

%C Zero together with the partial sums of A080512. - _Omar E. Pol_, Sep 10 2011

%C Second heptagonal numbers (A147875) and positive terms of A000566 interleaved. - _Omar E. Pol_, Aug 04 2012

%C These numbers appear in a theta function identity. See the Hardy-Wright reference, Theorem 355 on p. 284. See the g.f. of A113429. - _Wolfdieter Lang_, Oct 28 2016

%C Characteristic function is A133100. - _Michael Somos_, Jan 30 2017

%C 40*a(n) + 9 is a square. - _Bruno Berselli_, Apr 18 2018

%D G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, Fifth ed., Clarendon Press, Oxford, 2003, p. 284.

%H Vincenzo Librandi, <a href="/A085787/b085787.txt">Table of n, a(n) for n = 0..10000</a>

%H <a href="/index/Rec#order_05">Index entries for linear recurrences with constant coefficients</a>, signature (1,2,-2,-1,1)

%F a(n) = A000217(n) + A000217(floor(n/2)).

%F a(2*n-1) = A000566(n).

%F a(2*n) = A147875(n). - _Bruno Berselli_, Apr 18 2018

%F G.f.: x * (1 + 3*x + x^2) / ((1 - x) * (1 - x^2)^2). a(n) = a(-1-n) for all n in Z. - _Michael Somos_, Oct 17 2006

%F a(n) = 5*n*(n + 1)/8 - 1/16 + (-1)^n*(2*n + 1)/16. - _R. J. Mathar_, Jun 29 2009

%F a(n) = (A000217(n) + A001082(n))/2 = (A001318(n) + A118277(n))/2. - _Omar E. Pol_, Jan 11 2013

%F a(n) = A002378(n) - A001318(n). - _Omar E. Pol_, Oct 23 2013

%F Sum_{n>=1} 1/a(n) = 10/9 + (2*sqrt(1 - 2/sqrt(5))*Pi)/3. - _Vaclav Kotesovec_, Oct 05 2016

%e From the first formula: a(5) = A000217(5) + A000217(2) = 15 + 3 = 18.

%t Select[Table[(n*(n+1)/2-1)/5,{n,500}],IntegerQ] (* _Vladimir Joseph Stephan Orlovsky_, Feb 06 2012 *)

%o (PARI) t(n)=n*(n+1)/2

%o for(i=0,40,print1(t(i)+t(floor(i/2)), ", "))

%o (PARI) {a(n) = (5*(-n\2)^2 - (-n\2)*3*(-1)^n) / 2}; /* _Michael Somos_, Oct 17 2006 */

%o (MAGMA) [5*n*(n+1)/8-1/16+(-1)^n*(2*n+1)/16: n in [0..60]]; // _Vincenzo Librandi_, Sep 11 2011

%o (Haskell)

%o a085787 n = a085787_list !! n

%o a085787_list = scanl (+) 0 a080512_list

%o -- _Reinhard Zumkeller_, Apr 06 2015

%Y Column 3 of A195152.

%Y Cf. A001318: A000217(n) - A000217(floor(n/2)).

%Y Generalized k-gonal numbers, k>=5: A001318, A000217, this sequence, A001082, A118277, A074377, A195160, A195162, A195313, A195818.

%Y Cf. A080512, A113429, A133100.

%K nonn,easy

%O 0,3

%A _Jon Perry_, Jul 23 2003

%E New name from _T. D. Noe_, Apr 21 2006

%E Formula in sequence name added by _Omar E. Pol_, May 28 2012

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

License Agreements, Terms of Use, Privacy Policy .

Last modified May 23 12:12 EDT 2018. Contains 304467 sequences. (Running on oeis4.)