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.

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

Original entry on oeis.org

0, 1, 3, 5, 8, 11, 15, 19, 24, 29, 35, 41, 48, 55, 63, 71, 80, 89, 99, 109, 120, 131, 143, 155, 168, 181, 195, 209, 224, 239, 255, 271, 288, 305, 323, 341, 360, 379, 399, 419, 440, 461, 483, 505, 528, 551, 575, 599, 624, 649, 675, 701, 728, 755, 783, 811, 840
Offset: 1

Views

Author

Keywords

Comments

a(n+1) is the number of 2 X n binary matrices with no zero rows or columns, up to row and column permutation.
[ (4th elementary symmetric function of S(n))/(3rd elementary symmetric function of S(n)) ], where S(n) = {first n+3 odd positive integers}.
First differences are 1, 2, 2, 3, 3, 4, 4, 5, 5, ... .
Let M_n denotes the n X n matrix m(i,j) = 1 if i =j; m(i,j) = 1 if (i+j) is odd; m(i,j) = 0 if i+j is even, then a(n) = -det M_(n+1) - Benoit Cloitre, Jun 19 2002
a(n) is the number of squares with corners on an n X n grid, distinct up to translation. See also A002415, A108279.
Starting (1, 3, 5, 8, 11, ...), = row sums of triangle A135841. - Gary W. Adamson, Dec 01 2007
Number of solutions to x+y >= n-1 in integers x,y with 1 <= x <= y <= n-1. - Franz Vrabec, Feb 22 2008
Let A be the Hessenberg matrix of order n, defined by: A[1,j]=1, A[i,i]:=-1, A[i,i-1]=-1, and A[i,j]=0 otherwise. Then, for n>=5, a(n-4)=-coeff(charpoly(A,x),x^2). - Milan Janjic, Jan 26 2010
Equals row sums of a triangle with alternate columns of (1,2,3,...) and (1,1,1,...). - Gary W. Adamson, May 21 2010
Conjecture: if a(n) = p#(primorial)-1 for some prime number p, then q=(n+1) is also a prime number where p#=floor(q^2/4). Tested up to n=10^100000 no counterexamples are found. It seems that the subsequence is very scattered. So far the triples (p,q,a(q-1)) are {(2,3,1), (3,5,5), (5,11,29), (7,29,209), (17,1429,510509)}. - David Morales Marciel, Oct 02 2015
Numbers of an Ulam spiral starting at 0 in which the shape of the spiral is exactly a rectangle. E.g., a(4)=5 the Ulam spiral is including at that moment only the elements 0,1,2,3,4,5 and the shape is a rectangle. The area is always a(n)+1. E.g., for a(4) the area of the rectangle is 2(rows) X 3(columns) = 6 = a(4) + 1. - David Morales Marciel, Apr 05 2016
Numbers of different quadratic forms (quadrics) in the real projective space P^n(R). - Serkan Sonel, Aug 26 2020
a(n+1) is the number of one-dimensional subspaces of (F_3)^n, counted up to coordinate permutation. E.g.: For n=4, there are five one-dimensional subspaces in (F_3)^3 up to coordinate permutation: [1 2 2] [0 2 2] [1 0 2] [0 0 2] [1 1 1]. This example suggests a bijection (which has to be adjusted for the all-ones matrix) with the binary matrices of the first comment. - Álvar Ibeas, Sep 21 2021

Examples

			There are five 2 X 3 binary matrices with no zero rows or columns up to row and column permutation:
   [1 0 0]  [1 0 0]  [1 1 0]  [1 1 0]  [1 1 1]
   [0 1 1]  [1 1 1]  [0 1 1]  [1 1 1]  [1 1 1].
		

References

  • O. Giering, Vorlesungen über höhere Geometrie, Vieweg, Braunschweig, 1982. See p. 59.

Crossrefs

Cf. A014616, A135841, A034856, A005744 (partial sums), A008619 (1st differences).
A row or column of the array A196416 (possibly with 1 subtracted from it).
Cf. A008619.
Second column of A232206.

Programs

  • GAP
    a:=[0,1,3,5];; for n in [5..65] do a[n]:=2*a[n-1]-2*a[n-3]+a[n-4]; od; a; # Muniru A Asiru, Oct 23 2018
    
  • Haskell
    a024206 n = (n - 1) * (n + 3) `div` 4
    a024206_list = scanl (+) 0 $ tail a008619_list
    -- Reinhard Zumkeller, Dec 18 2013
    
  • Magma
    [(2*n^2+4*n-7-(-1)^n)/8 : n in [1..100]]; // Wesley Ivan Hurt, Jul 22 2014
    
  • Maple
    A024206:=n->(2*n^2+4*n-7-(-1)^n)/8: seq(A024206(n), n=1..100);
  • Mathematica
    f[x_, y_] := Floor[ Abs[ y/x - x/y]]; Table[ Floor[ f[2, n^2 + 2 n - 2] /2], {n, 57}] (* Robert G. Wilson v, Aug 11 2010 *)
    LinearRecurrence[{2,0,-2,1},{0,1,3,5},60] (* Harvey P. Dale, Jun 14 2013 *)
    Rest[CoefficientList[Series[x^2 (1 + x - x^2)/((1 - x^2) (1 - x)^2), {x, 0, 70}], x]] (* Vincenzo Librandi, Oct 02 2015 *)
  • PARI
    a(n)=(n-1)*(n+3)\4 \\ Charles R Greathouse IV, Jun 26 2013
    
  • PARI
    x='x+O('x^99); concat(0, Vec(x^2*(1+x-x^2)/ ((1-x^2)*(1-x)^2))) \\ Altug Alkan, Apr 05 2016
    
  • Python
    def A024206(n): return (n+1)**2//4 - 1 # Ya-Ping Lu, Jan 01 2024

Formula

G.f.: x^2*(1+x-x^2)/((1-x^2)*(1-x)^2) = x^2*(1+x-x^2) / ( (1+x)*(1-x)^3 ).
a(n+1) = A002623(n) - A002623(n-1) - 1.
a(n) = A002620(n+1) - 1 = A014616(n-2) + 1.
a(n+1) = A002620(n) + n, n >= 0. - Philippe Deléham, Feb 27 2004
a(0)=0, a(n) = floor(a(n-1) + sqrt(a(n-1)) + 1) for n > 0. - Gerald McGarvey, Jul 30 2004
a(n) = floor((n+1)^2/4) - 1. - Franz Vrabec, Feb 22 2008
a(n) = A005744(n-1) - A005744(n-2). - R. J. Mathar, Nov 04 2008
a(n) = a(n-1) + [side length of the least square > a(n-1) ], that is a(n) = a(n-1) + ceiling(sqrt(a(n-1) + 1)). - Ctibor O. Zizka, Oct 06 2009
For a(1)=0, a(2)=1, a(n) = 2*a(n-1) - a(n-2) + 1 if n is odd; a(n) = 2*a(n-1) - a(n-2) if n is even. - Vincenzo Librandi, Dec 23 2010
a(n) = A181971(n, n-1) for n > 0. - Reinhard Zumkeller, Jul 09 2012
a(n) = 2*a(n-1) - 2*a(n-3) + a(n-4); a(1)=0, a(2)=1, a(3)=3, a(4)=5. - Harvey P. Dale, Jun 14 2013
a(n) = floor( (n-1)*(n+3)/4 ). - Wesley Ivan Hurt, Jun 23 2013
a(n) = (2*n^2 + 4*n - 7 - (-1)^n)/8. - Wesley Ivan Hurt, Jul 22 2014
a(n) = a(-n-2) = n-1 + floor( (n-1)^2/4 ). - Bruno Berselli, Feb 03 2015
a(n) = (1/4)*(n+3)^2 - (1/8)*(1 + (-1)^n) - 1. - Serkan Sonel, Aug 26 2020
a(n) + a(n+1) = A034856(n). - R. J. Mathar, Mar 13 2021
a(2*n) = n^2 + n - 1, a(2*n+1) = n^2 + 2*n. - Greg Dresden and Zijie He, Jun 28 2022
Sum_{n>=2} 1/a(n) = 7/4 + tan(sqrt(5)*Pi/2)*Pi/sqrt(5). - Amiram Eldar, Dec 10 2022
E.g.f.: (4 + (x^2 + 3*x - 4)*cosh(x) + (x^2 + 3*x - 3)*sinh(x))/4. - Stefano Spezia, Aug 06 2024

Extensions

Corrected and extended by Vladeta Jovovic, Jun 02 2000

A032132 Number of series-reduced dyslexic planted planar trees with n leaves.

Original entry on oeis.org

1, 1, 2, 6, 17, 57, 191, 684, 2482, 9275, 35127, 135156, 525545, 2064329, 8173895, 32600082, 130823306, 527888023, 2140454687, 8716907165, 35638352814, 146221542191, 601870210193, 2484682879348, 10285116277096, 42679973961811, 177514171393035, 739881841810694, 3089914920914855, 12927860306782626
Offset: 1

Views

Author

Keywords

Comments

Apparently, beginning with a(3), number of non-equivalent canonical forms of separation coordinates on the hyperspheres. Cf. Schöbel and Veselov for this and other interpretations. - Tom Copeland, Nov 21 2017
From Petros Hadjicostas, Jan 17 2018: (Start)
Let A(x) = Sum_{n>=1} a(n)*x^n. For a derivation of the formula BIK(A(x)) = (1/2)*(A(x)/(1-A(x)) + (A(x) + A(x^2))/(1-A(x^2))), see the comments for sequence A001224 and the weblink below containing Bower's theory of transforms.
We clarify the comment by T. Copeland above. Consider the material in Section 3 of Devadoss and Read (2001). According to their terminology, let b(m,n) be "the number of A-clusters having m cells and n outside edges not counting the root edge." Let B(x,y) = Sum_{m>=0, n>=0} b(m,n)*x^m*y^n. (See p. 78 in their paper, where they use the notations a_{m,n} and A(x,y) rather than b(m,n) and B(x,y), respectively, that we use here.)
On p. 79 (Eq. (3.1)) of their paper, they prove that B(x,y) = y + (x/2)*(B(x,y)^2/(1-B(x,y)) + (1 + B(x,y))*B(x^2, y^2)/(1-B(x^2,y^2))). Unfortunately, the factor x in the previous formula is left out (i.e., it is a typo), and the same typo is reproduced in Schöbel and Veselov (2014, 2015).
Using Table 2 (p. 92) from Devadoss and Read (2001) (and the material on p. 79), we get that B(x,y) = y+ x*y^2 + (x^2 + x)*y^3 + (2*x^3 + 3*x^2 + x)*y^4 + (3*x^4 + 8*x^3 + 5*x^2 + x)*y^5 + ...
We claim that a(n) = Sum_{m>=0} b(m,n) and A(y) = Sum_{n>=1} a(n)*y^n = B(x=1, y). To prove these claims, note that, for x=1, the above series becomes B(x=1,y) = y + y^2 + 2*y^3 + 6*y^4 + 17*y^5 + ..., while the functional equation above becomes B(1, y) = y + (1/2)*(B(1,y)^2/(1-B(1,y)) + (1 + B(1,y))*B(1,y^2)/(1-B(1,y^2))), which is equivalent to 2*B(1,y) = y + (1/2)*(B(1,y)/(1-B(1,y)) + (B(1,y) + B(1,y^2))/(1-B(1,y^2))). The latter formula is the one given in the formula section below (derived from Bower's theory) with x replaced with y and A(x) replaced with B(1,y). This proves that B(x=1, y) = A(y), from which we can easily get that a(n) = Sum_{m>=0} b(m,n).
Note that b(m=0, n) = 0 for n <> 1, but b(m=0, n=1) = 1; b(m,n) = 0 when m >= n >= 1; and b(m=1, n) = 1 for n>=2. Also, b(m,m+1) = A001190(m+1) for m>=1, which are the Wedderburn-Etherington numbers, and apparently b(m=2, n) = A024206(n-1) for n>=2 (conjecture).
In Section 6 of their paper, Schöbel and Veselov (2014, 2015) prove that b(m,n) is the "number of non-equivalent faces of [the Stasheff polytope] K_n of codimension m-1." Apparently then, for n>=2 and k>=0, b(n-k,n+1) is the "number of canonical forms for separation coordinates of [hypersphere] S^n" with k "independent continuous parameters". For k=0 and n>=2, b(n,n+1) = A001190(n+1) = "number of canonical forms for separation coordinates" of hypersphere S^n with 0 continuous parameters.
It turns out that for k, the number of continuous parameters of S^n, we have 0 <= k <= n-1 (see pp. 1269-1270 in Shobel and Veselov (2015)). Hence, for n>=2, Sum_{k=0..n-1} b(n-k, n+1) = Sum_{m=1..n} b(m, n+1) = Sum_{m=0..n} b(m, n+1) = a(n+1) (see above). As a result, for n>=2, a(n+1) is the "total number of [non-equivalent] canonical forms for separation coordinates on [hypersphere] S^n", which is the comment made by T. Copeland above.
(End)
For an explanation on the meaning of clusters of types A, B, and C see Section 3 (pp. 78-81) in Devadoos and Read (2001). See also the comments for sequence A232206. - Petros Hadjicostas, Mar 02 2018

Crossrefs

Programs

  • Mathematica
    BIK[p_] := (1/(1-p) + (1+p)/(1-p /. x -> x^2))/2;
    seq[n_] := Module[{p=x}, For[i=2, i <= n, i++, p += x^i Coefficient[BIK[p] + x O[x]^i // Normal, x, i]]; CoefficientList[p/x, x]];
    seq[30] (* Jean-François Alcover, Nov 22 2018, after Andrew Howroyd *)
  • PARI
    BIK(p)={(1/(1-p) + (1+p)/subst(1-p, x, x^2))/2}
    seq(n)={my(p=x); for(i=2, n, p+=x^i*polcoeff(BIK(p) + O(x*x^i), i)); Vecrev(p/x)} \\ Andrew Howroyd, Aug 30 2018

Formula

Doubles (index 2+) under "BIK" (reversible, indistinct, unlabeled) transform.
G.f.: If A(x) = Sum_{n>=1} a(n)*x^n, then 2*A(x) = x + BIK(A(x)) = x + (1/2)*(A(x)/(1-A(x)) + (A(x) + A(x^2))/(1-A(x^2))). - Petros Hadjicostas, Jan 17 2018

Extensions

a(25)-a(30) from Petros Hadjicostas, Jan 17 2018

A295380 Number of canonical forms for separation coordinates on hyperspheres S_n, ordered by increasing number of independent continuous parameters.

Original entry on oeis.org

1, 1, 1, 2, 3, 1, 3, 8, 5, 1, 6, 20, 22, 8, 1, 11, 49, 73, 46, 11, 1, 23, 119, 233, 206, 87, 15, 1, 46, 288, 689, 807, 485, 147, 19, 1, 98, 696, 1988, 2891, 2320, 1021, 236, 24, 1, 207, 1681, 5561, 9737, 9800, 5795, 1960, 356, 29, 1, 451, 4062, 15322, 31350, 38216, 28586, 13088, 3525, 520, 35, 1, 983, 9821, 41558, 97552, 139901, 127465, 74280, 27224, 5989, 730, 41, 1
Offset: 1

Views

Author

Tom Copeland, Nov 21 2017

Keywords

Comments

Table 1 of the Schöbel and Veselov paper with initial 1 added. Reverse of Table 2 of the Devadoss and Read paper.
Apparently A032132 contains the row sums.
From Petros Hadjicostas, Jan 28 2018: (Start)
In this triangle, which is read by rows, for 0 <= k <= n-1 and n>=1, let T(n,k) be the number of inequivalent canonical forms for separation coordinates of the hypersphere S^n with k independent continuous parameters. It is the mirror image of sequence A232206, that is, T(n, k) = A232206(n+1, n-k) for 0 <= k <= n-1 and n>=1. (Triangular array A232206(N, K) is defined for N >= 2 and 1 <= K <= N-1.)
If B(x,y) = Sum_{n,k>=0} T(n,k)*x^n*y^k (with T(0,0) = 1, T(0,k) = 0 for k>=1, and T(n,k) = 0 for 1 <= n <= k), then B(x,y) = 1 + (x/2)*(B(x,y)^2/(1-x*y*B(x,y)) + (1 + x*y*B(x,y))*B(x^2,y^2)/(1-x^2*y^2*B(x^2,y^2))). This can be derived from the bivariate g.f. of A232206. See the comments for that sequence.
Let S(n) := Sum_{k>=0} T(n,k). The g.f. of S(n) is B(x, y=1). If we let y=1 in the above functional equation, we get x*B(x,1) = x + (1/2)*((x*B(x,1))^2/(1-x*B(x,1)) + (1 + x*B(x,1))*x^2*B(x^2,1)/(1-x^2*B(x^2,1))). After some algebra, we get 2*x*B(x,1) = x + (1/2)(x*B(x,1)/(1-x*B(x,1)) + (x*B(x,1) + x^2*B(x^2,1))/(1-x^2*B(x,1))), i.e., 2*x*B(x,1) = x + BIK(x*B(x,1)), where we have the "BIK" (reversible, indistinct, unlabeled) transform of C. G. Bower. This proves that S(n) = A032132(n+1) for n>=0, which is Copeland's claim above.
Note that for the second column we have T(n,k=2) = A048739(n-2) for 2 <= n < = 10, but T(11,2) = 4062 <> 4059 = A048739(9). In any case, they have different g.f.s (see the formula section below).
(End)

Examples

			From _Petros Hadjicostas_, Jan 27 2018: (Start)
Triangle T(n,k) begins:
n\k      0     1     2     3     4     5     6    7   8  9
----------------------------------------------------------------
(S^1)    1,
(S^2)    1,    1,
(S^3)    2,    3,    1,
(S^4)    3,    8,    5,    1,
(S^5)    6,   20,   22,    8,    1,
(S^6)   11,   49,   73,   46,   11,    1,
(S^7)   23,  119,  233,  206,   87,   15,    1,
(S^8)   46,  288,  689,  807,  485,  147,   19,   1,
(S^9)   98,  696, 1988, 2891, 2320, 1021,  236,  24,  1,
(S^10) 207, 1681, 5561, 9737, 9800, 5795, 1960, 356, 29, 1,
...
(End)
		

Crossrefs

Formula

From Petros Hadjicostas, Jan 28 2018: (Start)
G.f.: If B(x,y) = Sum_{n,k>=0} T(n,k)*x^n*y^k (with T(0,0) = 1, T(0,k) = 0 for k>=1, and T(n,k) = 0 for 1 <= n <= k), then B(x,y) = 1 + (x/2)*(B(x,y)^2/(1-x*y*B(x,y)) + (1 + x*y*B(x,y))*B(x^2,y^2)/(1-x^2*y^2*B(x^2,y^2))).
If c(N,K) = A232206(N,K) and C(x,y) = Sum_{N,K>=0} c(N,K)*x^N*y^K (with c(1,0) = 1 and c(N,K) = 0 for 0 <= N <= K), then C(x,y) = x*B(x*y, 1/y) and B(x,y) = C(x*y, 1/y)/(x*y).
Setting y=0 in the above functional equation, we get x*B(x,0) = x + (1/2)*((x*B(x,0))^2 + x^2*B(x^2,0)), which is the functional equation for the g.f. of the first column. This proves that T(n,k=0) = A001190(n+1) for n>=0 (assuming T(0,0) = 1).
The g.f. of the second column is B_1(x,0) = Sum_{n>=0} T(n,2)*x^n = lim_{y->0} (B(x,y)-B(x,0))/y, where B(x,0) = 1 + x + x^2 + ... is the g.f. of the first column. We get B_1(x,0) = x*B(x,0)*(B(x,0) - 1)/(1 - x*B(x,0)).
(End)

Extensions

Typo for T(11,3)=15322 corrected by Petros Hadjicostas, Jan 28 2018
Showing 1-3 of 3 results.