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.

A024495 a(n) = C(n,2) + C(n,5) + ... + C(n, 3*floor(n/3)+2).

Original entry on oeis.org

0, 0, 1, 3, 6, 11, 21, 42, 85, 171, 342, 683, 1365, 2730, 5461, 10923, 21846, 43691, 87381, 174762, 349525, 699051, 1398102, 2796203, 5592405, 11184810, 22369621, 44739243, 89478486, 178956971, 357913941, 715827882, 1431655765, 2863311531, 5726623062
Offset: 0

Views

Author

Keywords

Comments

Trisections give A082365, A132804, A132805. - Paul Curtz, Nov 18 2007
If the offset is changed to 1, this is the maximal number of closed regions bounded by straight lines after n straight line cuts in a plane: a(n) = a(n-1) + n - 3, a(1)=0; a(2)=0; a(3)=1; and so on. - Srikanth K S, Jan 23 2008
M^n * [1,0,0] = [A024493(n), a(n), A024494(n)]; where M = a 3x3 matrix [1,1,0; 0,1,1; 1,0,1]. Sum of terms = 2^n. Example: M^5 * [1,0,0] = [11, 11, 10], sum = 2^5 = 32. - Gary W. Adamson, Mar 13 2009
For n>=1, a(n-1) is the number of generalized compositions of n when there are i^2/2 - 3*i/2 + 1 different types of i, (i=1,2,...). - Milan Janjic, Sep 24 2010
Let M be any endomorphism on any vector space, such that M^3 = 1 (identity). Then (1+M)^n = A024493(n) + A024494(n)*M + a(n)*M^2. - Stanislav Sykora, Jun 10 2012
{A024493, A131708, A024495} is the difference analog of the hyperbolic functions {h_1(x), h_2(x), h_3(x)} of order 3. For the definitions of {h_i(x)} and the difference analog {H_i(n)} see [Erdelyi] and the Shevelev link respectively. - Vladimir Shevelev, Aug 01 2017
This is the p-INVERT of (1,1,1,1,1,...) for p(S) = 1 - S^3; see A291000. - Clark Kimberling, Aug 24 2017

References

  • A. Erdelyi, Higher Transcendental Functions, McGraw-Hill, 1955, Vol. 3, Chapter XVIII.
  • D. E. Knuth, The Art of Computer Programming. Addison-Wesley, Reading, MA, Vol. 1, 2nd. ed., Problem 38, p. 70.

Crossrefs

Sequences of the form 1/((1-x)^m - x^m): A000079 (m=1,2), this sequence (m=3), A000749 (m=4), A049016 (m=5), A192080 (m=6), A049017 (m=7), A290995 (m=8), A306939 (m=9).

Programs

  • Magma
    R:=PowerSeriesRing(Integers(), 30); [0,0] cat Coefficients(R!( x^2/((1-x)^3-x^3) )); // G. C. Greubel, Apr 11 2023
    
  • Maple
    a:= proc(n) option remember; `if`(n=0, 0, 2*a(n-1)+
          [-1, 0, 1, 1, 0, -1, -1][1+(n mod 6)])
        end:
    seq(a(n), n=0..33); # Paul Weisenhorn, May 17 2020
  • Mathematica
    LinearRecurrence[{3,-3,2},{0,0,1},40] (* Harvey P. Dale, Sep 20 2016 *)
  • PARI
    a(n) = sum(k=0,n\3,binomial(n,3*k+2)) /* Michael Somos, Feb 14 2006 */
    
  • PARI
    a(n)=if(n<0, 0, ([1,0,1;1,1,0;0,1,1]^n)[3,1]) /* Michael Somos, Feb 14 2006 */
    
  • SageMath
    def A024495(n): return (2^n - chebyshev_U(n, 1/2) - chebyshev_U(n-1, 1/2))/3
    [A024495(n) for n in range(41)] # G. C. Greubel, Apr 11 2023

Formula

a(n) = ( 2^n + 2*cos((n-4)*Pi/3) )/3 = (2^n - A057079(n))/3.
a(n) = 2*a(n-1) + A010892(n-2) = a(n-1) + A024494(n-1). With initial zero, binomial transform of A011655 which is effectively A010892 unsigned. - Henry Bottomley, Jun 04 2001
a(2) = 1, a(3) = 3, a(n+2) = a(n+1) - a(n) + 2^n. - Benoit Cloitre, Sep 04 2002
a(n) = Sum_{k=0..n} 2^k*2*sin(Pi*(n-k)/3 + Pi/3)/sqrt(3) (offset 0). - Paul Barry, May 18 2004
G.f.: x^2/((1-x)^3 - x^3) = x^2 / ( (1-2*x)*(1-x+x^2) ).
a(n) = 3*a(n-1) - 3*a(n-2) + 2*a(n-3). - Paul Curtz, Nov 18 2007
a(n) + A024493(n-1) = A131577(n). - Paul Curtz, Jan 24 2008
From Paul Curtz, May 29 2011: (Start)
a(n) + a(n+3) = 3*2^n = A007283(n).
a(n+6) - a(n) = 21*2^n = A175805(n).
a(n) + a(n+9) = 171*2^n.
a(n+12) - a(n) = 1365*2^n. (End)
a(n) = A113405(n) + A113405(n+1). - Paul Curtz, Jun 05 2011
Start with x(0)=1, y(0)=0, z(0)=0 and set x(n+1) = x(n) + z(n), y(n+1) = y(n) + x(n), z(n+1) = z(n) + y(n). Then a(n) = z(n). - Stanislav Sykora, Jun 10 2012
G.f.: -x^2/( x^3 - 1 + 3*x/Q(0) ) where Q(k) = 1 + k*(x+1) + 3*x - x*(k+1)*(k+4)/Q(k+1) ; (continued fraction). - Sergei N. Gladkovskii, Mar 15 2013
a(n) = 1/18*(-4*(-1)^floor((n - 1)/3) - 6*(-1)^floor(n/3) - 3*(-1)^floor((n + 1)/3) + (-1)^(1 + floor((n + 2)/3)) + 3*2^(n + 1)). - John M. Campbell, Dec 23 2016
a(n) = (1/63)*(-40 + 21*2^n - 42*floor(n/6) + 32*floor((n+3)/6) + 16*floor((n+ 4)/6) - 24*floor((n+5)/6) - 22*floor((n+7)/6) + 21*floor((n+8)/6) + 10*floor((n+9)/6) + 5*floor((n+10)/6) + 3*floor((n+11)/6) + floor((n+ 13)/6)). - John M. Campbell, Dec 24 2016
a(n+m) = a(n)*A024493(m) + A131708(n)*A131708(m) + A024493(n)*a(m). - Vladimir Shevelev, Aug 01 2017
From Kevin Ryde, Sep 24 2020: (Start)
a(n) = (1/3)*2^n - (1/3)*cos((1/3)*Pi*n) - (1/sqrt(3))*sin((1/3)*Pi*n). [Cournot]
a(n) + A111927(n) + A131708(n) = 2^n - 1. [Cournot, page 96 last formula, but misprint should be 2^x - 1 rather than 2^p - 1] (End)
E.g.f.: (exp(2*x) - exp(x/2)*(cos(sqrt(3)*x/2) + sqrt(3)*sin(sqrt(3)*x/2)))/3. - Stefano Spezia, Feb 06 2025

A119363 a(n) = Sum_{k=0..n} C(n,3k)^2.

Original entry on oeis.org

1, 1, 1, 2, 17, 101, 402, 1275, 3921, 14114, 58601, 243695, 950578, 3537847, 13166791, 50514102, 198627921, 782913717, 3054480306, 11824753551, 45823049817, 178682390994, 700285942731, 2747647985241, 10767833451954, 42164261091351, 165225573240651
Offset: 0

Views

Author

Paul Barry, May 16 2006

Keywords

Comments

a(n) - A119364(n) = A119365(n).

Crossrefs

Central coefficients of number triangle A119335.
a(n) = A119335(2n, n).

Programs

  • Mathematica
    Table[Sum[Binomial[n,3k]^2, {k,0,n}], {n,0,30}] (* Vaclav Kotesovec, Mar 12 2019 *)
    Table[HypergeometricPFQ[{1/3 - n/3, 1/3 - n/3, 2/3 - n/3, 2/3 - n/3, -n/3, -n/3}, {1/3, 1/3, 2/3, 2/3, 1}, 1], {n, 0, 30}] (* Vaclav Kotesovec, Mar 12 2019 *)

Formula

From Vaclav Kotesovec, Mar 12 2019: (Start)
Recurrence: (n-2)*(n-1)*n*(637*n^6 - 11466*n^5 + 84364*n^4 - 324394*n^3 + 686227*n^2 - 755060*n + 336132)*a(n) = 3*(n-2)*(n-1)*(1274*n^7 - 23569*n^6 + 180194*n^5 - 733383*n^4 + 1699606*n^3 - 2208294*n^2 + 1449504*n - 351000)*a(n-1) - 3*(n-2)*(3185*n^8 - 63700*n^7 + 539028*n^6 - 2512118*n^5 + 7020469*n^4 - 11971242*n^3 + 12050010*n^2 - 6446736*n + 1362744)*a(n-2) + (14014*n^9 - 315315*n^8 + 3072678*n^7 - 16986046*n^6 + 58535088*n^5 - 129861691*n^4 + 184326992*n^3 - 159830656*n^2 + 75517728*n - 14313456)*a(n-3) + 3*(n-3)*(3185*n^8 - 63700*n^7 + 538391*n^6 - 2501394*n^5 + 6946794*n^4 - 11707256*n^3 + 11530544*n^2 - 5915328*n + 1142208)*a(n-4) + 18*(n-4)*(n-3)*(2*n - 9)*(637*n^6 - 7644*n^5 + 36589*n^4 - 88858*n^3 + 114124*n^2 - 71840*n + 16440)*a(n-5).
a(n) ~ 4^n / (3*sqrt(Pi*n)). (End)

Extensions

Edited by N. J. A. Sloane, Jun 12 2008

A139468 a(n) = Sum{k=0..n} C(n,3k+1)^2.

Original entry on oeis.org

0, 1, 4, 9, 17, 50, 261, 1275, 5028, 17253, 58601, 218042, 876789, 3537847, 13783018, 52301709, 198627921, 767778786, 3010327497, 11824753551, 46200429186, 179787741723, 700285942731, 2738134757118, 10739885115573, 42164261091351, 165467386466802
Offset: 0

Views

Author

N. J. A. Sloane, Jun 12 2008

Keywords

Comments

The recurrence is same as for A119363. - Vaclav Kotesovec, Mar 12 2019

Crossrefs

Programs

  • Mathematica
    Table[Sum[Binomial[n,3k+1]^2,{k,0,n}],{n,0,30}] (* Harvey P. Dale, Sep 08 2018 *)

Formula

a(n) ~ 4^n / (3*sqrt(Pi*n)). - Vaclav Kotesovec, Mar 12 2019
Showing 1-3 of 3 results.