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 11-20 of 92 results. Next

A001752 Expansion of 1/((1+x)*(1-x)^5).

Original entry on oeis.org

1, 4, 11, 24, 46, 80, 130, 200, 295, 420, 581, 784, 1036, 1344, 1716, 2160, 2685, 3300, 4015, 4840, 5786, 6864, 8086, 9464, 11011, 12740, 14665, 16800, 19160, 21760, 24616, 27744, 31161, 34884, 38931, 43320, 48070, 53200, 58730, 64680, 71071, 77924, 85261
Offset: 0

Views

Author

Keywords

Comments

Define a unit column of a binary matrix to be a column with only one 1. a(n) = number of 3 X n binary matrices with 1 unit column up to row and column permutations (if offset is 1). - Vladeta Jovovic, Sep 09 2000
Generally, number of 3 X n binary matrices with k=0,1,2,... unit columns, up to row and column permutations, is the coefficient of x^k in 1/6*(Z(S_n; 5 + 3*x,5 + 3*x^2, ...) + 3*Z(S_n; 3 + x,5 + 3*x^2,3 + x^3,5 + 3*x^4, ...) + 2*Z(S_n; 2,2,5 + 3*x^3,2,2,5 + 3*x^6, ...)), where Z(S_n; x_1,x_2,...,x_n) is the cycle index of symmetric group S_n of degree n.
First differences of a(n) give number of minimal 3-covers of an unlabeled n-set that cover 4 points of that set uniquely (if offset is 4).
Transform of tetrahedral numbers, binomial(n+3,3), under the Riordan array (1/(1-x^2), x). - Paul Barry, Apr 16 2005
Equals triangle A152205 as an infinite lower triangular matrix * [1, 2, 3, ...]. - Gary W. Adamson, Feb 14 2010
With a leading zero, number of all possible octahedra of any size, formed when intersecting a regular tetrahedron by planes parallel to its sides and dividing its edges into n equal parts. - V.J. Pohjola, Sep 13 2012
With 2 leading zeros and offset 1, the sequence becomes 0,0,1,4,11,... for n=1,2,3,... Call this b(n). Consider the partitions of n into two parts (p,q) with p <= q. Then b(n) is the total volume of the family of rectangular prisms with dimensions p, |q - p| and |q - p|. - Wesley Ivan Hurt, Apr 14 2018
Conjecture: For n > 2, a(n-3) is the absolute value of the coefficient of the term [x^(n-2)] in the characteristic polynomial of an n X n square matrix M(n) defined as the n-th principal submatrix of the array A010751 whose general element is given by M[i,j] = floor((j - i + 1)/2). - Stefano Spezia, Jan 12 2022
Consider the following drawing of the complete graph on n vertices K_n: Vertices 1, 2, ..., n are on a straight line. Any pair of nonconsecutive vertices (i, j) with i < j is connected by a semicircle that goes above the line if i is odd, and below if i is even. With four leading zeros and offset 1, a(n) gives the number of edge crossings of the aforementioned drawing of K_n. - Carlo Francisco E. Adajar, Mar 17 2022

Examples

			There are 4 binary 3 X 2 matrices with 1 unit column up to row and column permutations:
  [0 0] [0 0] [0 1] [0 1]
  [0 0] [0 1] [0 1] [0 1]
  [0 1] [1 1] [1 0] [1 1].
For n=5, the numbers of the octahedra, starting from the smallest size, are Te(5)=35, Te(3)=10, and Te(1)=1, the sum being 46. Te denotes the tetrahedral number A000292. - _V.J. Pohjola_, Sep 13 2012
		

References

  • T. A. Saaty, The Minimum Number of Intersections in Complete Graphs, Proc. Natl. Acad. Sci. USA., 52 (1964), 688-690.

Crossrefs

Cf. A001753 (partial sums), A002623 (first differences), A158454 (signed column k=2), A169792 (binomial transform).

Programs

  • Magma
    [Floor(((n+3)^2-1)*((n+3)^2-3)/48): n in [0..40]]; // Vincenzo Librandi, Aug 15 2011
  • Maple
    A001752:=n->(3*(-1)^n+93+168*n+100*n^2+24*n^3+2*n^4)/96:
    seq(A001752(n), n=0..50); # Wesley Ivan Hurt, Apr 01 2015
  • Mathematica
    a = {1, 4}; Do[AppendTo[a, a[[n - 2]] + (n*(n + 1)*(n + 2))/6], {n, 3, 10}]; a
    (* Number of octahedra *) nnn = 100; Teo[n_] := (n - 1) n (n + 1)/6
    Table[Sum[Teo[n - nn], {nn, 0, n - 1, 2}], {n, 1, nnn}]
    (* V.J. Pohjola, Sep 13 2012 *)
    LinearRecurrence[{4,-5,0,5,-4,1},{1,4,11,24,46,80},50] (* Harvey P. Dale, Feb 07 2019 *)
  • PARI
    a(n)=if(n<0,0,((n+3)^2-1)*((n+3)^2-3)/48-if(n%2,1/16))
    
  • PARI
    a(n)=(n^4+12*n^3+50*n^2+84*n+46)\/48 \\ Charles R Greathouse IV, Sep 13 2012
    

Formula

a(n) = floor(((n+3)^2 - 1)*((n+3)^2 - 3)/48).
G.f.: 1/((1+x)*(1-x)^5).
a(n) - 2*a(n-1) + a(n-2) = A002620(n+2).
a(n) + a(n-1) = A000332(n+4).
a(n) - a(n-2) = A000292(n+1).
a(n) = Sum_{k=0..n} (-1)^(n-k)*C(k+4, 4). - Paul Barry, Jul 01 2003
a(n) = (3*(-1)^n + 93 + 168*n + 100*n^2 + 24*n^3 + 2*n^4)/96. - Cecilia Rossiter (cecilia(AT)noticingnumbers.net), Dec 14 2004
From Paul Barry, Apr 16 2005: (Start)
a(n) = Sum_{k=0..floor(n/2)} binomial(n-2k+3, 3).
a(n) = Sum_{k=0..n} binomial(k+3, 3)*(1-(-1)^(n+k-1))/2. (End)
a(n) = A108561(n+5,n) for n > 0. - Reinhard Zumkeller, Jun 10 2005
From Wesley Ivan Hurt, Apr 01 2015: (Start)
a(n) = 4*a(n-1) - 5*a(n-2) + 5*(n-4) - 4*a(n-5) + a(n-6).
a(n) = Sum_{i=0..n+3} (n+3-i) * floor(i^2/2)/2. (End)
Boas-Buck recurrence: a(n) = (1/n)*Sum_{p=0..n-1} (5 + (-1)^(n-p))*a(p), n >= 1, a(0) = 1. See the Boas-Buck comment in A158454 (here for the unsigned column k = 2 with offset 0). - Wolfdieter Lang, Aug 10 2017
Convolution of A000217 and A004526. - R. J. Mathar, Mar 29 2018
E.g.f.: ((48 + 147*x + 93*x^2 + 18*x^3 + x^4)*cosh(x) + (45 + 147*x + 93*x^2 + 18*x^3 + x^4)*sinh(x))/48. - Stefano Spezia, Jan 12 2022

Extensions

Formulae corrected by Bruno Berselli, Sep 13 2012

A002717 a(n) = floor(n(n+2)(2n+1)/8).

Original entry on oeis.org

0, 1, 5, 13, 27, 48, 78, 118, 170, 235, 315, 411, 525, 658, 812, 988, 1188, 1413, 1665, 1945, 2255, 2596, 2970, 3378, 3822, 4303, 4823, 5383, 5985, 6630, 7320, 8056, 8840, 9673, 10557, 11493, 12483, 13528, 14630, 15790, 17010, 18291, 19635, 21043, 22517, 24058
Offset: 0

Views

Author

Keywords

Comments

Number of triangles in triangular matchstick arrangement of side n, for n >= 1. Row sums of A085691.
We observe that the sequence is the transform of A006578 by the following transform T: T(u_0,u_1,u_2,u_3,...)=(u_0,u_0+u_1, u_0+u_1+u_2, u_0+u_1+u_2+u_3+u_4,...). In another terms v_p=sum(u_k,k=0..p) and the G.f phi_v of v is given by: phi_v=phi_u/(1-z). - Richard Choulet, Jan 28 2010
Row sums of A220053, for n > 0. - Reinhard Zumkeller, Dec 03 2012
a(n) has the expansion (1*0)+(1*1)+(4*1)+(4*2)+(7*2)+(7*3)+..., where the expansion stops when a(n) has n+1 number of terms. The expansion starts at (1*0), and progresses by alternating addition of 1 to the second number and 3 to the first number. - Arlu Genesis A. Padilla, Jun 04 2014
Taking the absolute values of each n-th difference and excluding the first n terms of each mentioned sequence, A002717 has the first difference A006578 (see formula of Michael Somos dated Jun 09 2014), the second difference A032766 (see 'partial sum' crossref), the third difference A000034, the fourth difference A000012, and the fifth to n-th difference A000004. - Arlu Genesis A. Padilla, Jun 12 2014

Examples

			f(3)=13 because the following figure contains 13 triangles if horizontal bars are added:
....... /\
...... /\/\
..... /\/\/\
G.f. = x + 5*x^2 + 13*x^3 + 27*x^4 + 48*x^5 + 78*x^6 + 118*x^7 + 170*x^8 + ...
		

References

  • J. H. Conway and R. K. Guy, The Book of Numbers, p. 83.
  • 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

Cf. A000292 number of triangles with same orientation as largest triangle, A002623 number of triangles pointing in opposite direction to largest triangle, A085691 number of triangles of side k in arrangement of side n.
Bisections: A135712 (odd part), A135713 (even part).

Programs

  • Magma
    [Floor(n*(n+2)*(2*n+1)/8): n in [0..50]]; // Wesley Ivan Hurt, Jun 04 2014
  • Maple
    A002717:=n->floor(n*(n+2)*(2*n+1)/8); seq(A002717(n), n=0..100);
  • Mathematica
    Table[Floor[n(n+2)(2n+1)/8],{n,0,50}] (* or *) LinearRecurrence[{3,-2,-2,3,-1},{0,1,5,13,27},50] (* Harvey P. Dale, Jan 20 2013 *)
  • PARI
    {a(n) = n * (n+2) * (2*n+1) \ 8};
    

Formula

a(n) = (1/16)*[2n(2n+1)(n+2)+cos(Pi*n)-1]. - Justin C. Bozonier (justinb67(AT)excite.com), Dec 05 2000
a(m+1)-2a(m)+2a(m-2)-a(m-3) = 3. - Len Smiley, Oct 08 2001
a(n) = (2n(2n+1)(n+2)+(-1)^n-1)/16. - Wesley Petty (Wesley.Petty(AT)mail.tamucc.edu), Oct 25 2003
a(n) = A000292(n-1) + A002623(n-2). - Hugo Pfoertner, Mar 06 2004
a(n) = Sum_{k=0..n} (-1)^(n-k)*k*binomial(k+1,2).
G.f.: x(1+2x)/((1+x)(1-x)^4). - Simon Plouffe in his 1992 dissertation (with a different offset).
a(0)=0, a(1)=1, a(2)=5, a(3)=13, a(4)=27, a(n)=3*a(n-1)-2*a(n-2)-2*a(n-3)+ 3*a(n-4)- a(n-5). - Harvey P. Dale, Jan 20 2013
a(n) = a(n-1) + A016777(floor(0.5*n))*floor(0.5+0.5*n). - Arlu Genesis A. Padilla, Jun 04 2014
a(-n) = - A045947(n). a(n) = a(n-1) + A006578(n). - Michael Somos, Jun 09 2014
a(n) = Sum_{i=1..n} T(n-i+1)+T(n-2*i+1), where T(n)=n*(n+1)/2=A000217(n) if n>0 and 0 if n<=0. So we have a(n+2)-a(n)=(n+2)^2+(n+1)*(n+2)/2. - Maurice Mischler, Sep 08 2014
E.g.f.: (x*(2*x^2 + 11*x + 9)*cosh(x) + (2*x^3 + 11*x^2 + 9*x - 1)*sinh(x))/8. - Stefano Spezia, Jul 19 2022

A028657 Triangle read by rows: T(n,k) = number of n-node graphs with k nodes in distinguished bipartite block, k = 0..n.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 7, 4, 1, 1, 5, 13, 13, 5, 1, 1, 6, 22, 36, 22, 6, 1, 1, 7, 34, 87, 87, 34, 7, 1, 1, 8, 50, 190, 317, 190, 50, 8, 1, 1, 9, 70, 386, 1053, 1053, 386, 70, 9, 1, 1, 10, 95, 734, 3250, 5624, 3250, 734, 95, 10, 1, 1, 11, 125, 1324, 9343, 28576, 28576, 9343, 1324, 125, 11, 1
Offset: 0

Views

Author

Vladeta Jovovic, Jun 16 2000

Keywords

Comments

Also, row n gives the number of unlabeled bicolored graphs having k nodes of one color and n-k nodes of the other color; the color classes are not interchangeable.
Also the number of principal transversal matroids (also known as fundamental transversal matroids) of size n and rank k (originally enumerated by Brylawski). - Gordon F. Royle, Oct 30 2007
This sequence is also obtained if we read the array A(m,n) = number of inequivalent m X n binary matrices by antidiagonals, where equivalence means permutations of rows or columns (m>=0, n>=0) [Kerber]. - N. J. A. Sloane, Sep 01 2013

Examples

			The triangle T(n,k) begins:
  1;
  1,  1;
  1,  2,  1;
  1,  3,  3,  1;
  1,  4,  7,  4,  1;
  1,  5, 13, 13,  5,  1;
  1,  6, 22, 36, 22,  6,  1;
  ...
For example, there are 36 graphs on 6 nodes with a distinguished bipartite block with 3 nodes.
The array A(m,n) (m>=0, n>=0) (see Comments) begins:
  1 1  1    1     1      1        1         1           1 ...
  1 2  3    4     5      6        7         8           9 ...
  1 3  7   13    22     34       50        70          95 ...
  1 4 13   36    87    190      386       734        1324 ...
  1 5 22   87   317   1053     3250      9343       25207 ...
  1 6 34  190  1053   5624    28576    136758      613894 ...
  1 7 50  386  3250  28576   251610   2141733    17256831 ...
  1 8 70  734  9343 136758  2141733  33642660   508147108 ...
  1 9 95 1324 25207 613894 17256831 508147108 14685630688 ...
... - _N. J. A. Sloane_, Sep 01 2013
		

References

  • R. W. Robinson, Numerical implementation of graph counting algorithms, AGRC Grant, Math. Dept., Univ. Newcastle, Australia, 1976.

Crossrefs

Row sums give A049312.
A246106 is a very similar array.
Diagonals of the array A(m,n) give A002724, A002725, A002728.
Rows (or columns) give A002623, A002727, A006148, A052264.
A(n,k) = A353585(2, n, k).

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0, {0}, `if`(i<1, {},
          {seq(map(p-> p+j*x^i, b(n-i*j, i-1) )[], j=0..n/i)}))
        end:
    g:= proc(n, k) option remember; add(add(2^add(add(igcd(i, j)*
          coeff(s, x, i)* coeff(t, x, j), j=1..degree(t)),
          i=1..degree(s))/mul(i^coeff(s, x, i)*coeff(s, x, i)!,
          i=1..degree(s))/mul(i^coeff(t, x, i)*coeff(t, x, i)!,
          i=1..degree(t)), t=b(n+k$2)), s=b(n$2))
        end:
    A:= (n, k)-> g(min(n, k), abs(n-k)):
    seq(seq(A(n, d-n), n=0..d), d=0..14); # Alois P. Heinz, Aug 01 2014
  • Mathematica
    b[n_, i_] := b[n, i] = If[n == 0, {0}, If[i<1, {}, Union[ Flatten[ Table[ Function[ {p}, p + j*x^i] /@ b[n - i*j, i-1], {j, 0, n/i}]]]]];
    g[n_, k_] := g[n, k] = Sum[ Sum[ 2^Sum[ Sum[GCD[i, j] * Coefficient[s, x, i] * Coefficient[t, x, j], {j, 1, Exponent[t, x]}], {i, 1, Exponent[s, x]}] / Product[i^Coefficient[s, x, i] * Coefficient[s, x, i]!, {i, 1, Exponent[s, x]}] / Product[i^Coefficient[t, x, i] * Coefficient[t, x, i]!, {i, 1, Exponent[t, x]}], {t, b[n+k, n+k]}], {s, b[n, n]}];
    A[n_, k_] := g[Min[n, k], Abs[n-k]];
    Table[Table[A[n, d-n], {n, 0, d}], {d, 0, 14}] // Flatten (* Jean-François Alcover, Jan 28 2015, after Alois P. Heinz *)
  • PARI
    permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
    K(q, t)={sum(j=1, #q, gcd(t, q[j]))}
    A(n, m)={my(s=0); forpart(q=m, s+=permcount(q)*polcoef(exp(sum(t=1, n, 2^K(q, t)/t*x^t) + O(x*x^n)), n)); s/m!}
    { for(r=0, 10, for(k=0, r, print1(A(r-k,k), ", ")); print) } \\ Andrew Howroyd, Mar 25 2020
    
  • PARI
    \\ G(k,x) gives k-th column as rational function (see Jovovic link).
    permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
    Fix(q,x)={my(v=divisors(lcm(Vec(q))), u=apply(t->2^sum(j=1, #q, gcd(t, q[j])), v)); 1/prod(i=1, #v, my(t=v[i]); (1-x^t)^(sum(j=1, i, my(d=t/v[j]); if(!frac(d), moebius(d)*u[j]))/t))}
    G(m,x)={my(s=0); forpart(q=m, s+=permcount(q)*Fix(q,x)); s/m!}
    T(n,k)={my(m=max(k, n-k)); polcoef(G(n-m, x + O(x*x^m)), m)} \\ Andrew Howroyd, Mar 26 2020
    
  • PARI
    A028657(n,k)=A353585(2, n, k) \\ M. F. Hasler, May 01 2022

Formula

A(m,n) = Sum_{p in P(m), q in P(n)} 2^Sum_{i in p, j in q} gcd(i,j) / (N(p) N(q)) where P(m) are the partition of m (see e.g., A036036), N(p) = Product_{distinct parts x in p} x^m(x)*m(x)!, m(x) = multiplicity of x in p. [corrected by Anders Kaseorg, Oct 04 2024]

A016061 a(n) = n*(n+1)*(4*n+5)/6.

Original entry on oeis.org

0, 3, 13, 34, 70, 125, 203, 308, 444, 615, 825, 1078, 1378, 1729, 2135, 2600, 3128, 3723, 4389, 5130, 5950, 6853, 7843, 8924, 10100, 11375, 12753, 14238, 15834, 17545, 19375, 21328, 23408, 25619, 27965, 30450, 33078, 35853, 38779, 41860
Offset: 0

Views

Author

Keywords

Comments

Number of ZnS molecules in cluster of n layers in zinc blende crystal.
(Zinc sulfide crystallizes in two different forms: wurtzite and zinc blende, the latter is also spelled zincblende.) - Jonathan Vos Post, Jan 22 2013
The Kn4 triangle sums of the Connell-Pol triangle A159797 lead to the sequence given above. For the definitions of the Kn4 and other triangle sums see A180662. - Johannes W. Meijer, May 20 2011
If one generated primitive Pythagorean triangles (2n+1, 2n+3) the collective sum of their perimeters for each n is four times the numbers listed in this sequence. - J. M. Bergot, Jul 18 2011
a(n) is the number of 3-tuples (w,x,y) having all terms in {0,...,n} and nA000292(n)+A000292(n+1)=n^3. - Clark Kimberling, Jun 04 2012
Degrees of the Hilbert polynomials for B_3 and C_3, per p. 13 of Gashi et al. - Jonathan Vos Post, Dec 14 2013
Number of solutions to a + b = c + d when 0 < a <= k, 0 <= b, c, d <= k, k = 0, 1, 2, 3.... Taken from Step 1 2007 problem #1(i) using 4 digit balanced numbers. - Bobby Milazzo, Mar 09 2013
From J. M. Bergot, Jun 18 2013: (Start)
Consider the lower half, including the main diagonal, of the array in A144216 as a triangle. The rows begin:
0;
1, 2;
3, 4, 6;
6, 7, 9, 12, ...
The sum of the terms in row(n) is a(n). (End)
This sequence is related to A008865 by a(n) = n*A008865(n+1) - Sum_{i=1..n} A008865(i) for n>0. - Bruno Berselli, Aug 06 2015

References

  • P. Jena and S. N. Behera, Clusters and Nanostructured Materials, Nova Science Publishers, 1996.

Crossrefs

Bisection of A002623.
Row sums of triangle A120070.

Programs

  • Magma
    I:=[0,3,13,34]; [n le 4 select I[n] else 4*Self(n-1)-6*Self(n-2)+4*Self(n-3)-Self(n-4): n in [1..40]]; // Vincenzo Librandi, Jul 25 2013
  • Maple
    A016061 := proc(n)
        n*(n+1)*(4*n+5)/6 ;
    end proc: # R. J. Mathar, Sep 26 2013
  • Mathematica
    CoefficientList[Series[x (3 + x) / (1 - x)^4, {x, 0, 40}], x] (* Vincenzo Librandi, Jul 25 2013 *)
    Table[n(n+1)(4*n+5)/6, {n,0,100}] (* Wesley Ivan Hurt, Sep 25 2013 *)
  • PARI
    v=vector(40,i,t(i)); s=0; forstep(i=2,40,2,s+=v[i]; print1(s","))
    

Formula

G.f.: x*(3+x)/(1-x)^4. - Paul Barry, Feb 27 2003
Partial sums of A014105. - Jon Perry, Jul 23 2003
a(n) = Sum_{i=0..n-1} 2*i^2 + i. - Jani Nurminen (slinky(AT)iki.fi), May 14 2006
a(n) = 2*n^3/3 +3*n^2/2 + 5*n/6. - Jonathan Vos Post, Dec 14 2013
a(n) = (4*n+5)/(2*n+1)*A000330(n). - Alexander R. Povolotsky, Mar 09 2013
a(n) = 4*a(n-1) -6*a(n-2) +4*a(n-3) -a(n-4). - Bobby Milazzo, Mar 10 2013
Sum_{n>=1} 1/a(n) = 12*Pi/5 + 72*log(2)/5 - 426/25. - Amiram Eldar, Jan 04 2022
E.g.f.: exp(x)*x*(18 + 21*x + 4*x^2)/6. - Stefano Spezia, Jul 31 2022

A058187 Expansion of (1+x)/(1-x^2)^4: duplicated tetrahedral numbers.

Original entry on oeis.org

1, 1, 4, 4, 10, 10, 20, 20, 35, 35, 56, 56, 84, 84, 120, 120, 165, 165, 220, 220, 286, 286, 364, 364, 455, 455, 560, 560, 680, 680, 816, 816, 969, 969, 1140, 1140, 1330, 1330, 1540, 1540, 1771, 1771, 2024, 2024, 2300, 2300, 2600, 2600, 2925, 2925, 3276, 3276
Offset: 0

Views

Author

Henry Bottomley, Nov 20 2000

Keywords

Comments

For n >= i, i = 6,7, a(n - i) is the number of incongruent two-color bracelets of n beads, i of which are black (cf. A005513, A032280), having a diameter of symmetry. The latter means the following: if we imagine (0,1)-beads as points (with the corresponding labels) dividing a circumference of a bracelet into n identical parts, then a diameter of symmetry is a diameter (connecting two beads or not) such that a 180-degree turn of one of two sets of points around it (obtained by splitting the circumference by this diameter) leads to the coincidence of the two sets (including their labels). - Vladimir Shevelev, May 03 2011
From Johannes W. Meijer, May 20 2011: (Start)
The Kn11, Kn12, Kn13, Fi1 and Ze1 triangle sums, see A180662 for their definitions, of the Connell-Pol triangle A159797 are linear sums of shifted versions of the duplicated tetrahedral numbers, e.g., Fi1(n) = a(n-1) + 5*a(n-2) + a(n-3) + 5*a(n-4).
The Kn11, Kn12, Kn13, Kn21, Kn22, Kn23, Fi1, Fi2, Ze1 and Ze2 triangle sums of the Connell sequence A001614 as a triangle are also linear sums of shifted versions of the sequence given above. (End)
The number of quadruples of integers [x, u, v, w] that satisfy x > u > v > w >= 0, n + 5 = x + u. - Michael Somos, Feb 09 2015
Also, this sequence is the fourth column in the triangle of the coefficients of the sum of two consecutive Fibonacci polynomials F(n+1, x) and F(n, x) (n>=0) in ascending powers of x. - Mohammad K. Azarian, Jul 18 2018

Crossrefs

Cf. A057884. Sum of 2 consecutive terms gives A006918, whose sum of 2 consecutive terms gives A002623, whose sum of 2 consecutive terms gives A000292, which is this sequence without the duplication. Continuing to sum 2 consecutive terms gives A000330, A005900, A001845, A008412 successively.

Programs

  • Haskell
    a058187 n = a058187_list !! n
    a058187_list = 1 : f 1 1 [1] where
       f x y zs = z : f (x + y) (1 - y) (z:zs) where
         z = sum $ zipWith (*) [1..x] [x,x-1..1]
    -- Reinhard Zumkeller, Dec 21 2011
    
  • Maple
    A058187:= proc(n) option remember; A058187(n):= binomial(floor(n/2)+3, 3) end: seq(A058187(n), n=0..51); # Johannes W. Meijer, May 20 2011
  • Mathematica
    a[n_]:= Length @ FindInstance[{x>u, u>v, v>w, w>=0, x+u==n+5}, {x, u, v, w}, Integers, 10^9]; (* Michael Somos, Feb 09 2015 *)
    With[{tetra=Binomial[Range[30]+2,3]},Riffle[tetra,tetra]] (* Harvey P. Dale, Mar 22 2015 *)
  • PARI
    {a(n) = binomial(n\2+3, 3)}; /* Michael Somos, Jun 07 2005 */
    
  • Sage
    [binomial((n//2)+3, 3) for n in (0..60)] # G. C. Greubel, Feb 18 2022

Formula

a(n) = A006918(n+1) - a(n-1).
a(2*n) = a(2*n+1) = A000292(n) = (n+1)*(n+2)*(n+3)/6.
a(n) = (2*n^3 + 21*n^2 + 67*n + 63)/96 + (n^2 + 7*n + 11)(-1)^n/32. - Paul Barry, Aug 19 2003
a(n) = A108299(n-3,n)*(-1)^floor(n/2) for n > 2. - Reinhard Zumkeller, Jun 01 2005
Euler transform of finite sequence [1, 3]. - Michael Somos, Jun 07 2005
G.f.: 1 / ((1 - x) * (1 - x^2)^3) = 1 / ((1 + x)^3 * (1 - x)^4). a(n) = -a(-7-n) for all n in Z.
a(n) = binomial(floor(n/2) + 3, 3). - Vladimir Shevelev, May 03 2011
a(-n) = -a(n-7); a(n) = A000292(A008619(n)). - Guenther Schrack, Sep 13 2018
Sum_{n>=0} 1/a(n) = 3. - Amiram Eldar, Aug 18 2022

A059260 Triangle read by rows giving coefficient T(i,j) of x^i y^j in 1/(1-y-x*y-x^2) = 1/((1+x)(1-x-y)) for (i,j) = (0,0), (1,0), (0,1), (2,0), (1,1), (0,2), ...

Original entry on oeis.org

1, 0, 1, 1, 1, 1, 0, 2, 2, 1, 1, 2, 4, 3, 1, 0, 3, 6, 7, 4, 1, 1, 3, 9, 13, 11, 5, 1, 0, 4, 12, 22, 24, 16, 6, 1, 1, 4, 16, 34, 46, 40, 22, 7, 1, 0, 5, 20, 50, 80, 86, 62, 29, 8, 1, 1, 5, 25, 70, 130, 166, 148, 91, 37, 9, 1, 0, 6, 30, 95, 200, 296, 314, 239, 128, 46, 10, 1
Offset: 0

Views

Author

N. J. A. Sloane, Jan 23 2001

Keywords

Comments

Coefficients of the (left, normalized) shifted cyclotomic polynomial. Or, coefficients of the basic n-th q-series for q=-2. Indeed, let Y_n(x) = Sum_{k=0..n} x^k, having as roots all the n-th roots of unity except for 0; then coefficients in x of (-1)^n Y_n(-x-1) give exactly the n-th row of A059260 and a practical way to compute it. - Olivier Gérard, Jul 30 2002
The maximum in the (2n)-th row is T(n,n), which is A026641; also T(n,n) ~ (2/3)*binomial(2n,n). The maximum in the (2n-1)-th row is T(n-1,n), which is A014300 (but T does not have the same definition as in A026637); also T(n-1,n) ~ (1/3)*binomial(2n,n). Here is a generalization of the formula given in A026641: T(i,j) = Sum_{k=0..j} binomial(i+k-x,j-k)*binomial(j-k+x,k) for all x real (the proof is easy by induction on i+j using T(i,j) = T(i-1,j) + T(i,j-1)). - Claude Morin, May 21 2002
The second greatest term in the (2n)-th row is T(n-1,n+1), which is A014301; the second greatest term in the (2n+1)-th row is T(n+1,n) = 2*T(n-1,n+1), which is 2*A014301. - Claude Morin
Diagonal sums give A008346. - Paul Barry, Sep 23 2004
Riordan array (1/(1-x^2), x/(1-x)). As a product of Riordan arrays, factors into the product of (1/(1+x),x) and (1/(1-x),1/(1-x)) (binomial matrix). - Paul Barry, Oct 25 2004
Signed version is A239473 with relations to partial sums of sequences. - Tom Copeland, Mar 24 2014
From Robert Coquereaux, Oct 01 2014: (Start)
Columns of the triangle (cf. Example below) give alternate partial sums along nw-se diagonals of the Pascal triangle, i.e., sequences A000035, A004526, A002620 (or A087811), A002623 (or A173196), A001752, A001753, A001769, A001779, A001780, A001781, A001786, A001808, etc.
The dimension of the space of closed currents (distributional forms) of degree p on Gr(n), the Grassmann algebra with n generators, equivalently, the dimension of the space of Gr(n)-valued symmetric multilinear forms with vanishing graded divergence, is V(n,p) = 2^n T(p,n-1) - (-1)^p.
If p is odd V(n,p) is also the dimension of the cyclic cohomology group of order p of the Z2 graded algebra Gr(n).
If p is even the dimension of this cohomology group is V(n,p)+1.
Cf. A193844. (End)
From Peter Bala, Feb 07 2024: (Start)
The following remarks assume the row indexing starts at n = 1.
The sequence of row polynomials R(n,x), beginning R(1,x) = 1, R(2,x) = x, R(3,x) = 1 + x + x^2 , ..., is a strong divisibility sequence of polynomials in the ring Z[x]; that is, for all positive integers n and m, poly_gcd( R(n,x), R(m,x)) = R(gcd(n, m), x) - apply Norfleet (2005), Theorem 3. Consequently, the polynomial sequence {R(n,x): n >= 1} is a divisibility sequence; that is, if n divides m then R(n,x) divides R(m,x) in Z[x]. (End)
From Miquel A. Fiol, Oct 04 2024: (Start)
For j>=1, T(i,j) is the independence number of the (i-j)-supertoken graph FF_(i-j)(S_j) of the star graph S_j with j points.
(Given a graph G on n vertices and an integer k>=1, the k-supertoken (or reduced k-th power) FF_k(G) of G has vertices representing configurations of k indistinguishable tokens in the (not necessarily different) vertices of G, with two configurations being adjacent if one can be obtained from the other by moving one token along an edge. See an example below.)
Following the suggestion of Peter Munn, the k-supertoken graph FF_k(S_j) can also be defined as follows: Consider the Lattice graph L(k,j), whose vertices are the k^j j-vectors with elements in the set {0,..,k-1}, two being adjacent if they differ in just one coordinate by one unity. Then, FF_k(S_j) is the subgraph of L(k+1,j) induced by the vertices at distance at most k from (0,..,0). (End)

Examples

			Triangle begins
  1;
  0,  1;
  1,  1,  1;
  0,  2,  2,  1;
  1,  2,  4,  3,  1;
  0,  3,  6,  7,  4,  1;
  1,  3,  9, 13, 11,  5,  1;
  0,  4, 12, 22, 24, 16,  6,  1;
  1,  4, 16, 34, 46, 40, 22,  7,  1;
  0,  5, 20, 50, 80, 86, 62, 29,  8,  1;
Sequences obtained with _Miquel A. Fiol_'s Sep 30 2024 formula of A(n,c1,c2) for other values of (c1,c2). (In the table, rows are indexed by c1=0..6 and columns by c2=0..6):
A000007  A000012  A000027  A025747  A000292* A000332* A000389*
A059841  A008619  A087811* A002623  A001752  A001753  A001769
A193356  A008794* A005993  A005994  -------  -------  -------
-------  -------  -------  A005995  A018210  -------  A052267
-------  -------  -------  -------  A018211  A018212  -------
-------  -------  -------  -------  -------  A018213  A018214
-------  -------  -------  -------  -------  -------  A062136
*requires offset adjustment.
The 2-supertoken FF_2(S_3) of the star graph S_3 with central vertex 1 and peripheral vertices 2,3,4. (The vertex `ij' of FF_2(S_3) represents the configuration of one token in `ì' and the other token in `j'). The T(5,3)=7 independent vertices are 22, 24, 44, 23, 11, 34, and 33.
     22--12---24---14---44
          | \    / |
         23   11   34
            \  |  /
              13
               |
              33
		

Crossrefs

Cf. A059259. Row sums give A001045.
Seen as a square array read by antidiagonals this is the coefficient of x^k in expansion of 1/((1-x^2)*(1-x)^n) with rows A002620, A002623, A001752, A001753, A001769, A001779, A001780, A001781, A001786, A001808 etc. (allowing for signs). A058393 would then effectively provide the table for nonpositive n. - Henry Bottomley, Jun 25 2001

Programs

  • Maple
    read transforms; 1/(1-y-x*y-x^2); SERIES2(%,x,y,12); SERIES2TOLIST(%,x,y,12);
  • Mathematica
    t[n_, k_] := Sum[ (-1)^(n-j)*Binomial[j, k], {j, 0, n}]; Flatten[ Table[t[n, k], {n, 0, 12}, {k, 0, n}]] (* Jean-François Alcover, Oct 20 2011, after Paul Barry *)
  • PARI
    T(n, k) = sum(j=0, n, (-1)^(n - j)*binomial(j, k));
    for(n=0, 12, for(k=0, n, print1(T(n, k),", ");); print();) \\ Indranil Ghosh, Apr 11 2017
    
  • Python
    from sympy import binomial
    def T(n, k): return sum((-1)**(n - j)*binomial(j, k) for j in range(n + 1))
    for n in range(13): print([T(n, k) for k in range(n + 1)]) # Indranil Ghosh, Apr 11 2017
  • Sage
    def A059260_row(n):
        @cached_function
        def prec(n, k):
            if k==n: return 1
            if k==0: return 0
            return -prec(n-1,k-1)-sum(prec(n,k+i-1) for i in (2..n-k+1))
        return [(-1)^(n-k+1)*prec(n+1, n-k+1) for k in (1..n)]
    for n in (1..9): print(A059260_row(n)) # Peter Luschny, Mar 16 2016
    

Formula

G.f.: 1/(1-y-x*y-x^2) = 1 + y + x^2 + xy + y^2 + 2x^2y + 2xy^2 + y^3 + ...
E.g.f: (exp(-t)+(x+1)*exp((x+1)*t))/(x+2). - Tom Copeland, Mar 19 2014
O.g.f. (n-th row): ((-1)^n+(x+1)^(n+1))/(x+2). - Tom Copeland, Mar 19 2014
T(i, 0) = 1 if i is even or 0 if i is odd, T(0, i) = 1 and otherwise T(i, j) = T(i-1, j) + T(i, j-1); also T(i, j) = Sum_{m=j..i+j} (-1)^(i+j+m)*binomial(m, j). - Robert FERREOL, May 17 2002
T(i, j) ~ (i+j)/(2*i+j)*binomial(i+j, j); more precisely, abs(T(i, j)/binomial(i+j, j) - (i+j)/(2*i+j) )<=1/(4*(i+j)-2); the proof is by induction on i+j using the formula 2*T(i, j) = binomial(i+j, j)+T(i, j-1). - Claude Morin, May 21 2002
T(n, k) = Sum_{j=0..n} (-1)^(n-j)binomial(j, k). - Paul Barry, Aug 25 2004
T(n, k) = Sum_{j=0..n-k} binomial(n-j, j)*binomial(j, n-k-j). - Paul Barry, Jul 25 2005
Equals A097807 * A007318. - Gary W. Adamson, Feb 21 2007
Equals A128173 * A007318 as infinite lower triangular matrices. - Gary W. Adamson, Feb 17 2007
Equals A130595*A097805*A007318 = (inverse Pascal matrix)*(padded Pascal matrix)*(Pascal matrix) = A130595*A200139. Inverse is A097808 = A130595*(padded A130595)*A007318. - Tom Copeland, Nov 14 2016
T(i, j) = binomial(i+j, j)-T(i-1, j). - Laszlo Major, Apr 11 2017
Recurrence for row polynomials (with row indexing starting at n = 1): R(n,x) = x*R(n-1,x) + (x + 1)*R(n-2,x) with R(1,x) = 1 and R(2,x) = x. - Peter Bala, Feb 07 2024
From Miquel A. Fiol, Sep 30 2024: (Start)
The triangle can be seen as a slice of a 3-dimensional table that links it to well-known sequences as follows.
The j-th column of the triangle, T(i,j) for i >= j, equals A(n,c1,c2) = Sum_{k=0..floor(n/2)} binomial(c1+2*k-1,2*k)*binomial(c2+n-2*k-1,n-2*k) when c1=1, c2=j, and n=i-j.
This gives T(i,j) = Sum_{k=0..floor((i-j)/2)} binomial(i-2*k-1, j-1). For other values of (c1,c2), see the example below. (End)

Extensions

Formula corrected by Philippe Deléham, Jan 11 2014

A002727 Number of 3 X n binary matrices up to row and column permutations.

Original entry on oeis.org

1, 4, 13, 36, 87, 190, 386, 734, 1324, 2284, 3790, 6080, 9473, 14378, 21323, 30974, 44159, 61898, 85440, 116286, 156240, 207446, 272432, 354162, 456097, 582238, 737205, 926298, 1155567, 1431892, 1763074, 2157904, 2626276, 3179278, 3829294, 4590118, 5477081
Offset: 0

Views

Author

Keywords

Comments

Also, number of unlabeled bipartite graphs with three left vertices and n right vertices. - Yavuz Oruc, Jan 22 2018

Examples

			G.f. = 1 + 4*x + 13*x^2 + 36*x^3 + 87*x^4 + 190*x^5 + 386*x^6 + 734*x^7 + ...
		

References

  • A. Kerber, Experimentelle Mathematik, Séminaire Lotharingien de Combinatoire. Institut de Recherche Math. Avancée, Université Louis Pasteur, Strasbourg, Actes 19 (1988), 77-83.
  • 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

A row of the array A(m,n) described in A028657. - N. J. A. Sloane, Sep 01 2013

Programs

  • Magma
    I:=[1,4,13,36,87,190,386,734,1324,2284,3790,6080,9473, 14378]; [n le 14 select I[n] else 4*Self(n-1)-4*Self(n-2)-2*Self(n-3)+2*Self(n-4)+4*Self(n-5)+3*Self(n-6)-12*Self(n-7)+ 3*Self(n-8)+4*Self(n-9)+2*Self(n-10)-2*Self(n-11)-4*Self(n-12)+4*Self(n-13)-Self(n-14): n in [1..50]]; // Vincenzo Librandi, Oct 13 2015
    
  • Mathematica
    CoefficientList[Series[(x^6+x^4+2x^3+x^2+1)/((1-x)^4(1-x^2)^2(1-x^3)^2),{x,0,40}],x] (* or *) LinearRecurrence[{4,-4,-2,2,4,3,-12,3,4,2,-2,-4,4,-1},{1,4,13,36,87,190,386,734,1324,2284,3790,6080,9473,14378},41] (* Harvey P. Dale, Nov 10 2011 *)
    Table[Which[
    Mod[n, 3] == 0,
    1/6 (1/27 (54 + 45 n + 12 n^2 + n^3) + 1/320 (4 + n) *(225 + 15 (-1)^n + 352 n + 172 n^2 + 32 n^3 + 2 n^4) + Binomial[7 + n, 7]),
    Mod[n, 3] == 1,
    1/6 (1/27 (50 + 45 n + 12 n^2 + n^3) + 1/320 (4 + n) *(225 + 15 (-1)^n + 352 n + 172 n^2 + 32 n^3 + 2 n^4) + Binomial[7 + n, 7]),
    Mod[n, 3] == 2,
    1/6 (1/27 (28 + 39 n + 12 n^2 + n^3) + 1/320 (4 + n) *(225 + 15 (-1)^n + 352 n + 172 n^2 + 32 n^3 + 2 n^4) + Binomial[7 + n, 7])
    ], {n, 0, 100}] (* Yavuz Oruc, Jan 22 2018 *)
  • PARI
    {a(n) = (6*n^7 + 168*n^6 + 2121*n^5 + 15540*n^4 + 70084*n^3 + 190512*n^2 + n*[284544, 281709, 277824, 281709, 284544, 274989][n%6+1]) \ 181440 + 1}; /* Michael Somos, Aug 22 2016 */
    
  • PARI
    x='x+O('x^99); Vec((1+x^2+2*x^3+x^4+x^6)/((1-x)^2*((1-x)*(1-x^2)*(1-x^3))^2)) \\ Altug Alkan, Mar 03 2018
    
  • PARI
    Vec(G(3, x) + O(x^40)) \\ G defined in A028657. - Andrew Howroyd, Feb 28 2023

Formula

G.f.: (x^6+x^4+2*x^3+x^2+1)/((1-x)^4*(1-x^2)^2*(1-x^3)^2). - Vladeta Jovovic, Feb 04 2000.
a(0)=1, a(1)=4, a(2)=13, a(3)=36, a(4)=87, a(5)=190, a(6)=386, a(7)=734, a(8)=1324, a(9)=2284, a(10)=3790, a(11)=6080, a(12)=9473, a(13)=14378. For n>13, a(n)=4*a(n-1)-4*a(n-2)-2*a(n-3)+2*a(n-4)+4*a(n-5)+3*a(n-6)- 12*a(n-7)+ 3*a(n-8)+4*a(n-9)+2*a(n-10)-2*a(n-11)-4*a(n-12)+4*a(n-13)-a(n-14). - Harvey P. Dale, Nov 10 2011
a(n) = -a(-8 - n) for all n in Z. - Michael Somos, Aug 22 2016
From Yavuz Oruc, Jan 22 2018: (Start)
If n == 0 (mod 3) then a(n)=(1/6)*(binomial(n+7,7) + (3(n+4)(2n^4 + 32n^3 + 172n^2 + 352n + 15(-1)^n + 225))/960 + (2(n^3 + 12n^2 + 45n + 54))/54).
If n == 1 (mod 3) then a(n)=(1/6)*(binomial(n+7,7) + (3(n+4)(2n^4 + 32n^3 + 172n^2 + 352n + 15(-1)^n + 225))/960 + (2(n^3 + 12n^2 + 45n + 50))/54).
If n == 2 (mod 3) then a(n)=(1/6)*(binomial(n+7,7) + (3(n+4)(2n^4 + 32n^3 + 172n^2 + 352n + 15(-1)^n + 225))/960 + (2(n^3 + 12n^2 + 39n + 28))/54). (End)

Extensions

More terms from Vladeta Jovovic, Feb 04 2000
Definition corrected by Max Alekseyev, Feb 05 2010

A014125 Bisection of A001400.

Original entry on oeis.org

1, 3, 6, 11, 18, 27, 39, 54, 72, 94, 120, 150, 185, 225, 270, 321, 378, 441, 511, 588, 672, 764, 864, 972, 1089, 1215, 1350, 1495, 1650, 1815, 1991, 2178, 2376, 2586, 2808, 3042, 3289, 3549, 3822, 4109, 4410, 4725, 5055, 5400, 5760, 6136, 6528, 6936, 7361, 7803
Offset: 0

Views

Author

Keywords

Comments

Also Schoenheim bound L_1(n,5,4).
Degrees of polynomials defined by p(n) = (x^(n+1)*p(n-1)p(n-3) + p(n-2)^2)/p(n-4), p(-4)=p(-3)=p(-2)=p(-1)=1. - Michael Somos, Jul 21 2004
Degrees of polynomial tau-functions of q-discrete Painlevé I, which generate sequence A095708 when q=2 (up to an offset of 3). - Andrew Hone, Jul 29 2004
Because of the Laurent phenomenon for the general q-discrete Painlevé I tau-function recurrence p(n) = (a*x^(n+1)*p(n-1)*p(n-3) + b*p(n-2)^2)/p(n-4), p(n) for n > -1 will always be a polynomial in x and a Laurent polynomial in a,b and the initial data p(-4),p(-3),p(-2),p(-1). - Andrew Hone, Jul 29 2004
Create the sequence 0,0,0,0,0,6,18,36,66,108,... so that the sum of three consecutive terms b(n) + b(n+1) + b(n+2) = A007531(n), with b(0)=0; then a(n) = b(n+5)/6. - J. M. Bergot, Jul 30 2013
Number of partitions of n into three kinds of part 1 and one kind of part 3. - Joerg Arndt, Sep 28 2015
First differences are A001840(k) starting with k=2; second differences are A086161(k) starting with k=1. - Bob Selcoe, Sep 28 2015
Maximum Wiener index of all maximal planar graphs with n+2 vertices. The extremal graphs are cubes of paths. - Allan Bickle, Jul 09 2022
Maximum Wiener index of all maximal 3-degenerate graphs with n+2 vertices. (A maximal 3-degenerate graph can be constructed from a 3-clique by iteratively adding a new 3-leaf (vertex of degree 3) adjacent to three existing vertices.) The extremal graphs are cubes of paths, so the bound also applies to 3-trees. - Allan Bickle, Sep 18 2022

Examples

			Polynomials: p(0)=x+1, p(1)=x^3+x^2+1, p(2)=x^6+x^5+x^3+x^2+2x+1, ...
a(12)=185:  A000217(13)=91 + a(9)=94 == 91+55+28+10+1 = 185. - _Bob Selcoe_, Sep 27 2015
a(3)=11: the 11 partitions of 3 are {1a,1a,1a}, {1a,1a,1b}, {1a,1a,1c}, {1a,1b,1b}, {1a,1b,1c}, {1a,1c,1c}, {1b,1b,1b}, {1b,1b,1c}, {1b,1c,1c}, {1c,1c,1c}, {3}. - _Bob Selcoe_, Oct 04 2015
		

References

  • W. H. Mills and R. C. Mullin, Coverings and packings, pp. 371-399 of Jeffrey H. Dinitz and D. R. Stinson, editors, Contemporary Design Theory, Wiley, 1992. See Eq. 1.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • L. Smiley, Hidden Hexagons (preprint).

Crossrefs

A column of A036838.
Maximum Wiener index of all maximal k-degenerate graphs for k=1..6: A000292, A002623, A014125, A122046, A122047, A175724.

Programs

  • Magma
    [n^3/18+n^2/2+4*n/3+1+(((n+1) mod 3)-1)/9 : n in [0..50]]; // Wesley Ivan Hurt, Apr 14 2015
    
  • Magma
    I:=[1,3,6,11,18,27]; [n le 6 select I[n] else 3*Self(n-1) -3*Self(n-2) +2*Self(n-3)-3*Self(n-4)+3*Self(n-5)-Self(n-6): n in [1..50]]; // Vincenzo Librandi, Apr 15 2015
    
  • Maple
    L := proc(v,k,t,l) local i,t1; t1 := l; for i from v-t+1 to v do t1 := ceil(t1*i/(i-(v-k))); od: t1; end; # gives Schoenheim bound L_l(v,k,t). Current sequence is L_1(n,n-3,n-4,1).
  • Mathematica
    CoefficientList[Series[1/((1 - x)^3*(1 - x^3)), {x, 0, 50}], x] (* Wesley Ivan Hurt, Apr 14 2015 *)
  • PARI
    a(n)=if(n<-5,-a(-6-n),polcoeff(1/(1-x)^3/(1-x^3)+x^n*O(x),n)) /* Michael Somos, Jul 21 2004 */
    
  • PARI
    my(x='x+O('x^50)); Vec(1/((1-x)^3*(1-x^3))) \\ Altug Alkan, Oct 16 2015
    
  • PARI
    a(n)=(n^3 + 9*n^2 + 24*n + 19)\/18 \\ Charles R Greathouse IV, Jun 29 2020
    
  • Sage
    [(binomial(n+4,3) - ((n+4)//3))/3 for n in (0..50)] # G. C. Greubel, Apr 28 2019

Formula

G.f.: 1/((1-x)^3*(1-x^3)).
a(n) = -a(-6-n) = 3*a(n-1) -3*a(n-2) +2*a(n-3) -3*a(n-4) +3*a(n-5) -a(n-6).
The simplest recurrence is fourth order: a(n) = a(n-1) + a(n-3) - a(n-4) + n + 1, which gives the g.f.: 1/((1-x)^3*(1-x^3)), with a(-4) = a(-3) = a(-2) = a(-1) = 0.
a(n) = n^3/18 + n^2/2 + 4*n/3 + 1 + 2/(9*sqrt(3))*sin(2*Pi*n/3). - Andrew Hone, Jul 29 2004
a(n) = n^3/18 + n^2/2 + 4*n/3 + 1 + (((n+1) mod 3) - 1)/9. - same formula, simplified by Gerald Hillier, Apr 14 2015
a(n) = (2*A000027(n+1) + 3*A000292(n+1) + A049347(n-1) + 1 + 3*A000217(n+1))/9. - R. J. Mathar, Nov 16 2007
From Johannes W. Meijer, May 20 2011: (Start)
a(n) = A144677(n) + A144677(n-1) + A144677(n-2).
a(n) = A190717(n-4) + 2*A190717(n-3) + 3*A190717(n-2) + 2*A190717(n-1) + A190717(n). (End)
3*a(n) = binomial(n+4,3) - floor((n+4)/3). - Bruno Berselli, Nov 08 2013
a(n) = A000217(n+1) + a(n-3) = Sum_{j>=0, n>=3*j} (n-3*j+1)*(n-3*j+2)/2. - Bob Selcoe, Sep 27 2015
a(n) = round(((2*n+5)^3 + 3*(2*n+5)^2 - 9*(2*n+5))/144). - Giacomo Guglieri, Jun 28 2020
a(n) = floor(((n+2)^3 + 3*(n+2)^2)/18). - Allan Bickle, Aug 01 2020
a(n) = Sum_{j=0..n} (n-j+1)*floor((j+3)/3). - G. C. Greubel, Oct 18 2021
E.g.f.: exp(x) + exp(x)*x*(34 + 12*x + x^2)/18 + 2*exp(-x/2)*sin(sqrt(3)*x/2)/(9*sqrt(3)). - Stefano Spezia, Apr 05 2023

Extensions

More terms from James Sellers, Dec 24 1999

A108561 Triangle read by rows: T(n,0)=1, T(n,n)=(-1)^n, T(n+1,k)=T(n,k-1)+T(n,k) for 0 < k < n.

Original entry on oeis.org

1, 1, -1, 1, 0, 1, 1, 1, 1, -1, 1, 2, 2, 0, 1, 1, 3, 4, 2, 1, -1, 1, 4, 7, 6, 3, 0, 1, 1, 5, 11, 13, 9, 3, 1, -1, 1, 6, 16, 24, 22, 12, 4, 0, 1, 1, 7, 22, 40, 46, 34, 16, 4, 1, -1, 1, 8, 29, 62, 86, 80, 50, 20, 5, 0, 1, 1, 9, 37, 91, 148, 166, 130, 70, 25, 5, 1, -1, 1, 10, 46, 128, 239, 314, 296, 200, 95, 30, 6, 0, 1, 1, 11, 56, 174, 367
Offset: 0

Views

Author

Reinhard Zumkeller, Jun 10 2005

Keywords

Comments

Sum_{k=0..n} T(n,k) = A078008(n);
Sum_{k=0..n} abs(T(n,k)) = A052953(n-1) for n > 0;
T(n,1) = n - 2 for n > 1;
T(n,2) = A000124(n-3) for n > 2;
T(n,3) = A003600(n-4) for n > 4;
T(n,n-6) = A001753(n-6) for n > 6;
T(n,n-5) = A001752(n-5) for n > 5;
T(n,n-4) = A002623(n-4) for n > 4;
T(n,n-3) = A002620(n-1) for n > 3;
T(n,n-2) = A008619(n-2) for n > 2;
T(n,n-1) = n mod 2 for n > 0;
T(2*n,n) = A072547(n+1).
Sum_{k=0..n} T(n,k)*x^k = A232015(n), A078008(n), A000012(n), A040000(n), A001045(n+2), A140725(n+1) for x = 2, 1, 0, -1, -2, -3 respectively. - Philippe Deléham, Nov 17 2013, Nov 19 2013
(1,a^n) Pascal triangle with a = -1. - Philippe Deléham, Dec 27 2013
T(n,k) = A112465(n,n-k). - Reinhard Zumkeller, Jan 03 2014

Examples

			From _Philippe Deléham_, Nov 17 2013: (Start)
Triangle begins:
  1;
  1, -1;
  1,  0,  1;
  1,  1,  1, -1;
  1,  2,  2,  0,  1;
  1,  3,  4,  2,  1, -1;
  1,  4,  7,  6,  3,  0,  1; (End)
		

Crossrefs

Cf. A007318 (a=1), A008949(a=2), A164844(a=10).
Similar to the triangles A035317, A059259, A080242, A112555.
Cf. A072547 (central terms).

Programs

  • GAP
    Flat(List([0..13],n->List([0..n],k->Sum([0..k],i->Binomial(n,i)*(-2)^(k-i))))); # Muniru A Asiru, Feb 19 2018
  • Haskell
    a108561 n k = a108561_tabl !! n !! k
    a108561_row n = a108561_tabl !! n
    a108561_tabl = map reverse a112465_tabl
    -- Reinhard Zumkeller, Jan 03 2014
    
  • Maple
    A108561 := (n, k) -> add(binomial(n, i)*(-2)^(k-i), i = 0..k):
    seq(seq(A108561(n,k), k = 0..n), n = 0..12); # Peter Bala, Feb 18 2018
  • Mathematica
    Clear[t]; t[n_, 0] = 1; t[n_, n_] := t[n, n] = (-1)^Mod[n, 2]; t[n_, k_] := t[n, k] = t[n-1, k] + t[n-1, k-1]; Table[t[n, k], {n, 0, 13}, {k, 0, n}] // Flatten (* Jean-François Alcover, Mar 06 2013 *)
  • Sage
    def A108561_row(n):
        @cached_function
        def prec(n, k):
            if k==n: return 1
            if k==0: return 0
            return -prec(n-1,k-1)-sum(prec(n,k+i-1) for i in (2..n-k+1))
        return [(-1)^k*prec(n, k) for k in (1..n-1)]+[(-1)^(n+1)]
    for n in (1..12): print(A108561_row(n)) # Peter Luschny, Mar 16 2016
    

Formula

G.f.: (1-y*x)/(1-x-(y+y^2)*x). - Philippe Deléham, Nov 17 2013
T(n,k) = T(n-1,k) + T(n-2,k-1) + T(n-2,k-2), T(0,0)=T(1,0)=1, T(1,1)=-1, T(n,k)=0 if k < 0 or if k > n. - Philippe Deléham, Nov 17 2013
From Peter Bala, Feb 18 2018: (Start)
T(n,k) = Sum_{i = 0..k} binomial(n,i)*(-2)^(k-i), 0 <= k <= n.
The n-th row polynomial is the n-th degree Taylor polynomial of the rational function (1 + x)^n/(1 + 2*x) about 0. For example, for n = 4, (1 + x)^4/(1 + 2*x) = 1 + 2*x + 2*x^2 + x^4 + O(x^5). (End)

Extensions

Definition corrected by Philippe Deléham, Dec 26 2013

A236679 Number T(n,k) of equivalence classes of ways of placing k 2 X 2 tiles in an n X n square under all symmetry operations of the square; irregular triangle T(n,k), n>=2, 0<=k<=floor(n/2)^2, read by rows.

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 4, 2, 1, 1, 3, 13, 20, 14, 1, 6, 37, 138, 277, 273, 143, 39, 7, 1, 1, 6, 75, 505, 2154, 5335, 7855, 6472, 2756, 459, 1, 10, 147, 1547, 10855, 50021, 153311, 311552, 416825, 361426, 200996, 71654, 16419, 2363, 211, 11, 1, 1, 10, 246, 3759, 39926, 291171
Offset: 2

Views

Author

Keywords

Comments

Changing the offset from 2 to 1, this is also the sequence: "Triangle read by rows: T(n,k) is the number of nonequivalent ways to place k non-attacking kings on an n X n board." (For if each king is represented by a 2 X 2 tile with the king in the upper left corner, the kings do not attack each other.) For example, with offset 1, T(4,3) = 20 because there are 20 nonequivalent ways to place 3 kings on a 4 X 4 chessboard so that no king threatens any other. - Heinrich Ludwig and N. J. A. Sloane, Dec 21 2016
It appears that rows 2n and 2n-1 both contain n^2 + 1 entries. Rotations and reflections of placements are not counted. If they are to be counted, see A193580. - Heinrich Ludwig, Dec 11 2016

Examples

			T(4,2) = 4 because the number of equivalence classes of ways of placing 2 2 X 2 square tiles in a 4 X 4 square under all symmetry operations of the square is 4. The portrayal of an example from each equivalence class is:
._______        _______        _______        _______
| . | . |      | . |___|      | . |   |      |_______|
|___|___|      |___| . |      |___|___|      | . | . |
|       |      |   |___|      |   | . |      |___|___|
|_______|      |_______|      |___|___|      |_______|
The first 6 rows of T(n,k) are:
.\ k  0    1    2    3    4    5    6    7    8    9
n
2     1    1
3     1    1
4     1    3    4    2    1
5     1    3   13   20   14
6     1    6   37  138  277  273  143   39    7    1
7     1    6   75  505 2154 5335 7855 6472 2756  459
		

Crossrefs

Row sums give A275869.
Diagonal T(n,n) is A279117.
Cf. A193580.

Formula

It appears that:
T(n,0) = 1, n>= 2
T(n,1) = (floor((n-2)/2)+1)*(floor((n-2)/2+2))/2, n >= 2
T(c+2*2,2) = A131474(c+1)*(2-1) + A000217(c+1)*floor(2^2/4) + A014409(c+2), 0 <= c < 2, c even
T(c+2*2,2) = A131474(c+1)*(2-1) + A000217(c+1)*floor((2-1)(2-3)/4) + A014409(c+2), 0 <= c < 2, c odd
T(c+2*2,3) = (c+1)(c+2)/2(2*A002623(c-1)*floor((2-c-1)/2) + A131941(c+1)*floor((2-c)/2)) + S(c+1,3c+2,3), 0 <= c < 2 where
S(c+1,3c+2,3) =
A054252(2,3), c = 0
A236679(5,3), c = 1

Extensions

More terms from Heinrich Ludwig, Dec 11 2016 (The former entry A279118 from Heinrich Ludwig was merged into this entry by N. J. A. Sloane, Dec 21 2016)
Previous Showing 11-20 of 92 results. Next