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

A034851 Rows of Losanitsch's triangle T(n, k), n >= 0, 0 <= k <= n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 2, 4, 2, 1, 1, 3, 6, 6, 3, 1, 1, 3, 9, 10, 9, 3, 1, 1, 4, 12, 19, 19, 12, 4, 1, 1, 4, 16, 28, 38, 28, 16, 4, 1, 1, 5, 20, 44, 66, 66, 44, 20, 5, 1, 1, 5, 25, 60, 110, 126, 110, 60, 25, 5, 1, 1, 6, 30, 85, 170, 236, 236, 170, 85, 30, 6, 1, 1, 6, 36, 110, 255
Offset: 0

Views

Author

Keywords

Comments

Sometimes erroneously called "Lossnitsch's triangle". But the author's name is Losanitsch (I have seen the original paper in Chem. Ber.). This is a German version of the Serbian name Lozanic. - N. J. A. Sloane, Jun 29 2008
For n >= 3, a(n-3,k) is the number of series-reduced (or homeomorphically irreducible) trees which become a path P(k+1) on k+1 nodes, k >= 0, when all leaves are omitted (see illustration). Proof by Pólya's enumeration theorem. - Wolfdieter Lang, Jun 08 2001
The number of ways to put beads of two colors in a line, but take symmetry into consideration, so that 011 and 110 are considered the same. - Yong Kong (ykong(AT)nus.edu.sg), Jan 04 2005
Alternating row sums are 1,0,1,0,2,0,4,0,8,0,16,0,... - Gerald McGarvey, Oct 20 2008
The triangle sums, see A180662 for their definitions, link Losanitsch's triangle A034851 with several sequences, see the crossrefs. We observe that the Ze3 and Ze4 sums link Losanitsch's triangle with A005683, i.e., R. K. Guy's Twopins game. - Johannes W. Meijer, Jul 14 2011
T(n-(L-1)k, k) is the number of ways to cover an n-length line by exactly k L-length segments excluding symmetric covers. For L=2 it is corresponds to A102541, for L=3 to A228570 and for L=4 to A228572. - Philipp O. Tsvetkov, Nov 08 2013
Also the number of equivalence classes of ways of placing k 1 X 1 tiles in an n X 1 rectangle under all symmetry operations of the rectangle. - Christopher Hunt Gribble, Feb 16 2014
T(n, k) is the number of non-isomorphic outer planar graphs of order n+3, size n+3+k, and maximum degree k+2. - Christian Barrientos, Oct 18 2018
From Álvar Ibeas, Jun 01 2020: (Start)
T(n, k) is the sum of even-degree coefficients of the Gaussian polynomial [n, k]_q. The area below a NE lattice path between (0,0) and (k, n-k) is even for T(n, k) paths and odd for A034852(n, k) of them.
For a (non-reversible) string of k black and n-k white beads, consider the minimum number of bead transpositions needed to place the black ones to the left and the white ones to the right (in other words, the number of inversions of the permutation obtained by labeling the black beads by integers 1,...,k and the white ones by k+1,...,n, in the same order they take on the string). It is even for T(n, k) strings and odd for A034852(n, k) cases.
(End)
Named after the Serbian chemist, politician and diplomat Simeon Milivoje "Sima" Lozanić (1847-1935). - Amiram Eldar, Jun 10 2021
T(n, k) is the number of caterpillars with a perfect matching, with 2n+2 vertices and diameter 2n-1-k. - Christian Barrientos, Sep 12 2023

Examples

			Triangle begins
  1;
  1,  1;
  1,  1,  1;
  1,  2,  2,  1;
  1,  2,  4,  2,  1;
  1,  3,  6,  6,  3,  1;
  1,  3,  9, 10,  9,  3,  1;
  1,  4, 12, 19, 19, 12,  4,  1;
  1,  4, 16, 28, 38, 28, 16,  4,  1;
  1,  5, 20, 44, 66, 66, 44, 20,  5,  1;
		

Crossrefs

Triangle sums (see the comments): A005418 (Row), A011782 (Related to Row2), A102526 (Related to Kn11, Kn12, Kn13, Kn21, Kn22, Kn23), A005207 (Kn3, Kn4), A005418 (Fi1, Fi2), A102543 (Ca1, Ca2), A192928 (Gi1, Gi2), A005683 (Ze3, Ze4).
Sums of squares of terms in rows equal A211208.

Programs

  • Haskell
    a034851 n k = a034851_row n !! k
    a034851_row 0 = [1]
    a034851_row 1 = [1,1]
    a034851_row n = zipWith (-) (zipWith (+) ([0] ++ losa) (losa ++ [0]))
                                ([0] ++ a204293_row (n-2) ++ [0])
       where losa = a034851_row (n-1)
    a034851_tabl = map a034851_row [0..]
    -- Reinhard Zumkeller, Jan 14 2012
  • Maple
    A034851 := proc(n,k) option remember; local t; if k = 0 or k = n then return(1) fi; if n mod 2 = 0 and k mod 2 = 1 then t := binomial(n/2-1,(k-1)/2) else t := 0; fi; A034851(n-1,k-1)+A034851(n-1,k)-t; end: seq(seq(A034851(n, k), k=0..n), n=0..11);
  • Mathematica
    t[n_?EvenQ, k_?OddQ] := Binomial[n, k]/2; t[n_, k_] := (Binomial[n, k] + Binomial[Quotient[n, 2], Quotient[k, 2]])/2; Flatten[Table[t[n, k], {n, 0, 12}, {k, 0, n}]](* Jean-François Alcover, Feb 07 2012, after PARI *)
  • PARI
    {T(n, k) = (1/2) *(binomial(n, k) + binomial(n%2, k%2) * binomial(n\2, k\2))}; /* Michael Somos, Oct 20 1999 */
    

Formula

T(n, k) = (1/2) * (A007318(n, k) + A051159(n, k)).
G.f. for k-th column (if formatted as lower triangular matrix a(n, k)): x^k*Pe(floor((k+1)/2), x^2)/(((1-x)^(k+1))*(1+x)^(floor((k+1)/2))), where Pe(n, x^2) := Sum_{m=0..floor(n/2)} A034839(n, m)*x^(2*m) (row polynomials of Pascal array even numbered columns). - Wolfdieter Lang, May 08 2001
a(n, k) = a(n-1, k-1) + a(n-1, k) - C(n/2-1, (k-1)/2), where the last term is present only if n is even and k is odd (see Sloane link).
T(n, k) = T(n-2, k-2) + T(n-2, k) + C(n-2, k-1), n > 1.
Let P(n, x, y) = Sum_{m=0..n} a(n, m)*x^m*y^(n-m), then for x > 0, y > 0 we have P(n, x, y) = (x+y)*P(n-1, x, y) for n odd and P(n, x, y) = (x+y)*P(n-1, x, y) - x*y*(x^2+y^2)^((n-2)/2) for n even. - Gerald McGarvey, Feb 15 2005
T(n, k) = T(n-1, k-1) + T(n-1, k) - A204293(n-2, k-1), 0 < k <= n and n > 1. - Reinhard Zumkeller, Jan 14 2012
From Christopher Hunt Gribble, Feb 25 2014: (Start)
It appears that:
T(n,k) = C(n,k)/2, n even, k odd;
T(n,k) = (C(n,k) + C(n/2,k/2))/2, n even, k even;
T(n,k) = (C(n,k) + C((n-1)/2,(k-1)/2))/2, n odd, k odd;
T(n,k) = (C(n,k) + C((n-1)/2,k/2))/2, n odd, k even.
(End)

Extensions

More terms from James Sellers, May 04 2000
Name edited by Johannes W. Meijer, Aug 26 2013

A005993 Expansion of (1+x^2)/((1-x)^2*(1-x^2)^2).

Original entry on oeis.org

1, 2, 6, 10, 19, 28, 44, 60, 85, 110, 146, 182, 231, 280, 344, 408, 489, 570, 670, 770, 891, 1012, 1156, 1300, 1469, 1638, 1834, 2030, 2255, 2480, 2736, 2992, 3281, 3570, 3894, 4218, 4579, 4940, 5340, 5740, 6181, 6622, 7106, 7590, 8119, 8648, 9224, 9800
Offset: 0

Views

Author

N. J. A. Sloane, Winston C. Yang (yang(AT)math.wisc.edu)

Keywords

Comments

Alkane (or paraffin) numbers l(6,n).
Dimension of the space of homogeneous degree n polynomials in (x1, y1, x2, y2) invariant under permutation of variables x1<->y1, x2<->y2.
Also multidigraphs with loops on 2 nodes with n arcs (see A138107). - Vladeta Jovovic, Dec 27 1999
Euler transform of finite sequence [2,3,0,-1]. - Michael Somos, Mar 17 2004
a(n-2) is the number of plane partitions with trace 2. - Michael Somos, Mar 17 2004
With offset 4, a(n) is the number of bracelets with n beads, 3 of which are red, 1 of which is blue. For odd n, a(n) = C(n-1,3)/2. For even n, a(n) = C(n-1,3)/2 +(n-2)/4. For n >= 6, with K = (n-1)(n-2)/((n-5)(n-4)), for odd n, a(n) = K*a(n-2). For even n, a(n) = K*a(n-2) -(n-2)/(n-5). - Washington Bomfim, Aug 05 2008
Equals (1,2,3,4,...) convolved with (1,0,3,0,5,...). - Gary W. Adamson, Feb 16 2009
Equals row sums of triangle A177878.
Equals (1/2)*((1, 4, 10, 20, 35, 56, ...) + (1, 0, 2 0, 3, 0, 4, ...)).
From Ctibor O. Zizka, Nov 21 2014: (Start)
With offset 4, a(n) is the number of different patterns of the 2-color 4-partition of n.
P(n)_(k;t) gives the number of different patterns of the t-color, k-partition of n.
P(n;i;j) = Sum(r=1..m) c_(i,j)*v_r*F_r(X_1,...,X_i).
m partition number of i.
c_(i,j) number of different coloring patterns on the r-th form (X_1,...,X_i) of i-partition with j-colors.
v_r number of i-partitions of n of the r-th form (X_1,...,X_i).
F_r(X_1,...,X_i) number of different patterns of the r-th form i-partition of n.
Some simple results:
P(1)(k;t)=1, P(2)(k;t)=2, P(3)(k;t)=4, P(4)(k;t)=11, etc.
P(n;1;1) = P(n;n;n) = 1 for all n;
P(n;2;2) = floor(n/2) (A004526);
P(n;3;2) = (n*n - 2*n + n mod 2)/4 (A002620).
This sequence is a(n) = P(n;4;2).
2-coloring of 4-partition is (A,B,A,B) or (B,A,B,A).
Each 4-partition of n has one of the form (X_1,X_1,X_1,X_1),(X_1,X_1,X_1,X_2), (X_1,X_1,X_2,X_2),(X_1,X_1,X_2,X_3),(X_1,X_2,X_3,X_4).
The number of forms is m=5 which is the partition number of k=4.
Partition form (X_1,X_1,X_1,X_1) gives 1 pattern ((X_1A,X_1B,X_1A,X_1B), (X_1,X_1,X_1,X_2) gives 2 patterns, (X_1,X_1,X_2,X_2) gives 4 patterns, (X_1,X_1,X_2,X_3) gives 6 patterns and (X_1,X_2,X_3,X_4) gives 12 patterns.
Thus a(n) = P(n;4;2) = 1*1*v_1 + 1*2*v_2 + 1*4*v_3 + 1*6*v_4 + 1*12*v_5 where v_r is the number of different 4-partitions of the r-th form (X_1,X_2,X_3,X_4) for a given n.
Example:
The 4-partitions of 8 are (2,2,2,2), (1,1,1,5), (1,1,3,3), (1,1,2,4), and (1,2,2,3):
(2,2,2,2) 1 pattern
(1,1,1,5), (1,1,5,1) 2 patterns
(1,1,3,3), (1,3,3,1), (3,1,1,3), (1,3,1,3) 4 patterns
(1,1,2,4), (1,1,4,2), (1,2,1,4), (1,2,4,1), (1,4,1,2), (2,1,1,4) 6 patterns
(2,2,1,3), (2,2,3,1), (2,1,2,3), (2,1,3,2), (2,3,2,1), (1,2,2,3) 6 patterns
Thus a(8) = P(8,4,2) = 1 + 2 + 4 + 6 + 6 = 19. (End)
a(n) = length of run n+2 of consecutive 1's in A254338. - Reinhard Zumkeller, Feb 27 2015
Take a chessboard of (n+2) X (n+2) unit squares in which the a1 square is black. a(n) is the number of composite squares having black unit squares on their vertices. - Ivan N. Ianakiev, Jul 19 2018
a(n) is the number of 1423-avoiding odd Grassmannian permutations of size n+2. Avoiding any of the patterns 2314 or 3412 gives the same sequence. - Juan B. Gil, Mar 09 2023

Examples

			a(2) = 6, since ( x1*y1, x2*y2, x1*x1+y1*y1, x2*x2+y2*y2, x1*x2+y1*y2, x1*y2+x2*y1 ) are a basis for homogeneous quadratic invariant polynomials.
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • L. Smith, Polynomial Invariants of Finite Groups, A K Peters, 1995, p. 96.

Crossrefs

Cf. A177878.
Partial sums of A008794 (without 0). - Bruno Berselli, Aug 30 2013

Programs

  • Haskell
    Following Gary W. Adamson.
    import Data.List (inits, intersperse)
    a005993 n = a005994_list !! n
    a005993_list = map (sum . zipWith (*) (intersperse 0 [1, 3 ..]) . reverse) $
                       tail $ inits [1..]
    -- Reinhard Zumkeller, Feb 27 2015
    
  • Magma
    I:=[1,2,6,10,19,28]; [n le 6 select I[n] else 2*Self(n-1)+Self(n-2)-4*Self(n-3)+Self(n-4)+2*Self(n-5)-Self(n-6): n in [1..60]]; // Vincenzo Librandi, Jul 19 2015
    
  • Maple
    g := proc(n) local i; add(floor(i/2)^2,i=1..n+1) end: # Joseph S. Riel (joer(AT)k-online.com), Mar 22 2002
    a:= n-> (Matrix([[1, 0$3, -1, -2]]).Matrix(6, (i,j)-> if (i=j-1) then 1 elif j=1 then [2, 1, -4, 1, 2, -1][i] else 0 fi)^n)[1,1]; seq (a(n), n=0..44); # Alois P. Heinz, Jul 31 2008
  • Mathematica
    CoefficientList[Series[(1+x^2)/((1-x)^2*(1-x^2)^2),{x,0,44}],x]  (* Jean-François Alcover, Apr 08 2011 *)
    LinearRecurrence[{2,1,-4,1,2,-1},{1,2,6,10,19,28},50] (* Harvey P. Dale, Feb 20 2012 *)
  • PARI
    a(n)=polcoeff((1+x^2)/(1-x)^2/(1-x^2)^2+x*O(x^n),n)
    
  • PARI
    a(n) = (binomial(n+3, n) + (1-n%2)*binomial((n+2)/2, n>>1))/2 \\ Washington Bomfim, Aug 05 2008
    
  • PARI
    a = vector(50); a[1]=1; a[2]=2;
    for(n=3, 50, a[n] = ((n+2)*a[n-2]+2*a[n-1]-n)/(n-2)); a \\ Gerry Martens, Jun 03 2018
    
  • Sage
    def A005993():
        a, b, to_be = 0, 0, True
        while True:
            yield (a*(a*(2*a+9)+13)+b*(b+1)*(2*b+1)+6)//6
            if to_be: b += 1
            else: a += 1
            to_be = not to_be
    a = A005993()
    [next(a) for  in range(48)] # _Peter Luschny, May 04 2016

Formula

l(c, r) = 1/2 C(c+r-3, r) + 1/2 d(c, r), where d(c, r) is C((c + r - 3)/2, r/2) if c is odd and r is even, 0 if c is even and r is odd, C((c + r - 4)/2, r/2) if c is even and r is even, C((c + r - 4)/2, (r - 1)/2) if c is odd and r is odd.
G.f.: (1+x^2)/((1-x)^2*(1-x^2)^2) = (1+x^2)/((1+x)^2*(x-1)^4) = (1/(1-x)^4 +1/(1-x^2)^2)/2.
a(2n) = (n+1)(2n^2+4n+3)/3, a(2n+1) = (n+1)(n+2)(2n+3)/3. a(-4-n) = -a(n).
From Yosu Yurramendi, Sep 12 2008: (Start)
a(n+1) = a(n) + A008794(n+3) with a(1)=1.
a(n) = A027656(n) + 2*A006918(n).
a(n+2) = a(n) + A000982(n+2) with a(1)=1, a(2)=2. (End)
a(n) = 2*a(n-1) + a(n-2) - 4*a(n-3) + a(n-4) + 2*a(n-5) - a(n-6). - Jaume Oliver Lafont, Dec 05 2008
a(n) = (n^3 + 6*n^2 + 11*n + 6)/12 + ((n+2)/4)[n even] (the bracket means that the second term is added if and only if n is even). - Benoit Jubin, Mar 31 2012
a(n) = (1/12)*n*(n+1)*(n+2) + (1/4)*(n+1)*(1/2)*(1-(-1)^n), with offset 1. - Yosu Yurramendi, Jun 20 2013
a(n) = Sum_{i=0..n+1} ceiling(i/2) * round(i/2) = Sum_{i=0..n+2} floor(i/2)^2. - Bruno Berselli, Aug 30 2013
a(n) = (n + 2)*(3*(-1)^n + 2*n^2 + 8*n + 9)/24. - Ilya Gutkovskiy, May 04 2016
Recurrence formula: a(n) = ((n+2)*a(n-2)+2*a(n-1)-n)/(n-2), a(1)=1, a(2)=2. - Gerry Martens, Jun 10 2018
E.g.f.: exp(-x)*(6 - 3*x + exp(2*x)*(18 + 39*x + 18*x^2 + 2*x^3))/24. - Stefano Spezia, Feb 23 2020
a(n) = Sum_{j=0..n/2} binomial(c+2*j-1,2*j)*binomial(c+n-2*j-1,n-2*j) where c=2. For other values of c we have: A008619 (c=1), A005995 (c=3), A018211 (c=4), A018213 (c=5), A062136 (c=6). - Miquel A. Fiol, Sep 24 2024

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

A062135 Odd-numbered columns of Losanitsch triangle A034851 formatted as triangle with an additional first column.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 2, 2, 1, 0, 2, 6, 3, 1, 0, 3, 10, 12, 4, 1, 0, 3, 19, 28, 20, 5, 1, 0, 4, 28, 66, 60, 30, 6, 1, 0, 4, 44, 126, 170, 110, 42, 7, 1, 0, 5, 60, 236, 396, 365, 182, 56, 8, 1, 0, 5, 85, 396, 868, 1001, 693, 280, 72, 9, 1
Offset: 0

Views

Author

Wolfdieter Lang, Jun 19 2001

Keywords

Comments

Because the sequence of column m=2*k, k >= 1, of A034851 is the partial sum sequence of the one of column m=2*k-1 the present triangle is essentially Losanitsch's triangle A034851.
Row sums give A051450 with A051450(0) := 1. Column sequences (without leading zeros) are for m=0..6: A000007, A008619, A005993, A005995, A018211, A018213, A062136.

Examples

			Triangle begins:
  {1};
  {0,1};
  {0,1,1};
  {0,2,2,1};
  ...
Pe(4,x^2)=1+6*x^2+x^4.
		

Crossrefs

Programs

  • Mathematica
    t[n_?EvenQ, k_?OddQ] := Binomial[n, k]/2; t[n_, k_] := (Binomial[n, k] + Binomial[Quotient[n, 2], Quotient[k, 2]])/2; Flatten[Table[t[n - 1 + m, n - m], {n, 0, 12}, {m, 0, n}]] (* Michael De Vlieger, Sep 28 2024, after Jean-François Alcover at A034851  *)

Formula

T(n, m) = A034851(n-1+m, n-m), n >= m >= 0; A034851(n-1, n) := 0, n >= 1, A034851(-1, 0) := 1.
T(n, m) = 0 if n= 1; T(n, m) = T(n-1, m)+sum(T(k, m-1), k=m-1..n-1) if n+m even and T(n, m) = T(n-1, m)+sum(T(k, m-1), k=m-1..n-1)-binomial((n+m-3)/2, m-1) if n+m odd, n >= m >= 1.
G.f. for column m: x^m*Pe(m, x^2)/(((1-x)^(2*m))*(1+x)^m), m >= 0, with Pe(m, x^2)= sum(A034839(m, k)*x^(2*k), k=0..floor(n/2)), the row polynomial of array A034839 (even-indexed entries of the rows of Pascal's triangle).

Extensions

More terms from Michael De Vlieger, Sep 28 2024

A141783 Number of bracelets (turn over necklaces) with n beads: 1 blue, 12 green, and r = n - 13 red.

Original entry on oeis.org

1, 7, 49, 231, 924, 3108, 9324, 25236, 63090, 147070, 323554, 676270, 1352540, 2600612, 4829708, 8692788, 15212379, 25949469, 43249115, 70562765, 112900424, 177412664, 274183208, 417232088, 625848132, 926250780, 1353751140
Offset: 13

Views

Author

Washington Bomfim, Aug 17 2008

Keywords

Crossrefs

Programs

  • Maple
    A141783:=n->(1/2)*(binomial(n - 1, 12) + binomial((n - 2 + (n mod 2))/2, 6)); seq(A141783(n), n=13..50); # Wesley Ivan Hurt, Jan 30 2014
  • Mathematica
    Table[(1/2) (Binomial[n - 1, 12] + Binomial[(n - 2 + Mod[n, 2])/2, 6]), {n, 13, 50}] (* Wesley Ivan Hurt, Jan 30 2014 *)

Formula

a(n) = 1/2*(binomial(n-1,12) + binomial((n-2+n mod 2)/2, 6)).
a(n) = (1/(2*12!))*(n+2)*(n+4)*(n+6)*(n+8)*(n+10)*(n+12)*((n+1)*(n+3)*(n+5)*(n+7)*(n+9)*(n+11) + 1*3*5*7*9*11) - (1/15)*(1/2^10)*(n^5+(65/2)*n^4+400*n^3+(4615/2)*n^2+6154*n+(11895/2))*(1/2)*(1-(-1)^n) [Yosu Yurramendi, Jun 24 2013]

Extensions

Revised by Washington Bomfim, Jul 24 2012
Showing 1-5 of 5 results.