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-10 of 126 results. Next

A173403 Inverse binomial transform of A002416.

Original entry on oeis.org

1, 1, 13, 469, 63577, 33231721, 68519123173, 562469619451069, 18442242396353040817, 2417685638793025070212561, 1267626422541873052658376446653, 2658442047546208031914776455678477989, 22300713297142388711251601783864453690641417
Offset: 0

Views

Author

Brian Drake, Feb 17 2010

Keywords

Comments

a(n) is the number of n X n matrices of 0's and 1's with the property that there is no index k such that both the k-th column and the k-th row consist of all zeros.
a(n) is the number of binary relations on n labeled vertices with no vertex of indegree and outdegree = 0. - Geoffrey Critzer, Oct 02 2012

References

  • E. A. Bender and S. G. Williamson, Foundations of Combinatorics with Applications, Dover, 2005, exercise 4.1.6.

Crossrefs

Programs

  • Maple
    N:=8: seq( sum(binomial(n,i)*2^((n-i)^2)*(-1)^(i), i=0..n), n=0..N);
  • Mathematica
    Table[Sum[(-1)^k Binomial[n,k] 2^(n-k)^2,{k,0,n}],{n,0,20}]  (* Geoffrey Critzer, Oct 02 2012 *)

Formula

a(n) = Sum_{k=0..n} (-1)^k*binomial(n,k)*2^((n-k)^2).
a(n) ~ 2^(n^2). - Vaclav Kotesovec, Oct 30 2017

A006125 a(n) = 2^(n*(n-1)/2).

Original entry on oeis.org

1, 1, 2, 8, 64, 1024, 32768, 2097152, 268435456, 68719476736, 35184372088832, 36028797018963968, 73786976294838206464, 302231454903657293676544, 2475880078570760549798248448, 40564819207303340847894502572032, 1329227995784915872903807060280344576
Offset: 0

Views

Author

Keywords

Comments

Number of graphs on n labeled nodes; also number of outcomes of labeled n-team round-robin tournaments.
Number of perfect matchings of order n Aztec diamond. [see Speyer]
Number of Gelfand-Zeitlin patterns with bottom row [1,2,3,...,n]. [Zeilberger]
For n >= 1 a(n) is the size of the Sylow 2-subgroup of the Chevalley group A_n(2) (sequence A002884). - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 30 2001
From James Propp: (Start)
a(n) is the number of ways to tile the region
o-----o
|.....|
o--o.....o--o
|...........|
o--o...........o--o
|.................|
o--o.................o--o
|.......................|
|.......................|
|.......................|
o--o.................o--o
|.................|
o--o...........o--o
|...........|
o--o.....o--o
|.....|
o-----o
(top-to-bottom distance = 2n) with dominoes like either of
o--o o-----o
|..| or |.....|
|..| o-----o
|..|
o--o
(End)
The number of domino tilings in A006253, A004003, A006125 is the number of perfect matchings in the relevant graphs. There are results of Jockusch and Ciucu that if a planar graph has a rotational symmetry then the number of perfect matchings is a square or twice a square - this applies to these 3 sequences. - Dan Fux (dan.fux(AT)OpenGaia.com or danfux(AT)OpenGaia.com), Apr 12 2001
Let M_n denotes the n X n matrix with M_n(i,j)=binomial(2i,j); then det(M_n)=a(n+1). - Benoit Cloitre, Apr 21 2002
Smallest power of 2 which can be expressed as the product of n distinct numbers (powers of 2), e.g., a(4) = 1024 = 2*4*8*16. Also smallest number which can be expressed as the product of n distinct powers. - Amarnath Murthy, Nov 10 2002
The number of binary relations that are both reflexive and symmetric on an n-element set. - Justin Witt (justinmwitt(AT)gmail.com), Jul 12 2005
The number of symmetric binary relations on an (n-1)-element set. - Peter Kagey, Feb 13 2021
To win a game, you must flip n+1 heads in a row, where n is the total number of tails flipped so far. Then the probability of winning for the first time after n tails is A005329 / A006125. The probability of having won before n+1 tails is A114604 / A006125. - Joshua Zucker, Dec 14 2005
a(n) = A126883(n-1)+1. - Zerinvary Lajos, Jun 12 2007
Equals right border of triangle A158474 (unsigned). - Gary W. Adamson, Mar 20 2009
a(n-1) is the number of simple labeled graphs on n nodes such that every node has even degree. - Geoffrey Critzer, Oct 21 2011
a(n+1) is the number of symmetric binary matrices of size n X n. - Nathan J. Russell, Aug 30 2014
Let T_n be the n X n matrix with T_n(i,j) = binomial(2i + j - 3, j-1); then det(T_n) = a(n). - Tony Foster III, Aug 30 2018
k^(n*(n-1)/2) is the determinant of n X n matrix T_(i,j) = binomial(k*i + j - 3, j-1), in this case k=2. - Tony Foster III, May 12 2019
Let B_n be the n+1 X n+1 matrix with B_n(i, j) = Sum_{m=max(0, j-i)..min(j, n-i)} (binomial(i, j-m) * binomial(n-i, m) * (-1)^m), 0<=i,j<=n. Then det B_n = a(n+1). Also, deleting the first row and any column from B_n results in a matrix with determinant a(n). The matrices B_n have the following property: B_n * [x^n, x^(n-1) * y, x^(n-2) * y^2, ..., y^n]^T = [(x-y)^n, (x-y)^(n-1) * (x+y), (x-y)^(n-2) * (x+y)^2, ..., (x+y)^n]^T. - Nicolas Nagel, Jul 02 2019
a(n) is the number of positive definite (-1,1)-matrices of size n X n. - Eric W. Weisstein, Jan 03 2021
a(n) is the number of binary relations on a labeled n-set that are both total and antisymmetric. - José E. Solsona, Feb 05 2023

Examples

			From _Gus Wiseman_, Feb 11 2021: (Start)
This sequence counts labeled graphs on n vertices. For example, the a(0) = 1 through a(2) = 8 graph edge sets are:
  {}  {}  {}    {}
          {12}  {12}
                {13}
                {23}
                {12,13}
                {12,23}
                {13,23}
                {12,13,23}
This sequence also counts labeled graphs with loops on n - 1 vertices. For example, the a(1) = 1 through a(3) = 8 edge sets are the following. A loop is represented as an edge with two equal vertices.
  {}  {}    {}
      {11}  {11}
            {12}
            {22}
            {11,12}
            {11,22}
            {12,22}
            {11,12,22}
(End)
		

References

  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 547 (Fig. 9.7), 573.
  • G. Everest, A. van der Poorten, I. Shparlinski, and T. Ward, Recurrence Sequences, Amer. Math. Soc., 2003; p. 178.
  • J. L. Gross and J. Yellen, eds., Handbook of Graph Theory, CRC Press, 2004; p. 517.
  • F. Harary, Graph Theory. Addison-Wesley, Reading, MA, 1969, p. 178.
  • F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973, p. 3, Eq. (1.1.2).
  • J. Propp, Enumeration of matchings: problems and progress, in: New perspectives in geometric combinatorics, L. Billera et al., eds., Mathematical Sciences Research Institute series, vol. 38, Cambridge University Press, 1999.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A000568 for the unlabeled analog, A053763, A006253, A004003.
Cf. A001187 (connected labeled graphs).
Cf. A158474. - Gary W. Adamson, Mar 20 2009
Cf. A136652 (log). - Paul D. Hanna, Dec 04 2009
The unlabeled version is A000088, or A002494 without isolated vertices.
The directed version is A002416.
The covering case is A006129.
The version for hypergraphs is A058891, or A016031 without singletons.
Row sums of A143543.
The case of connected edge set is A287689.

Programs

Formula

Sequence is given by the Hankel transform of A001003 (Schroeder's numbers) = 1, 1, 3, 11, 45, 197, 903, ...; example: det([1, 1, 3, 11; 1, 3, 11, 45; 3, 11, 45, 197; 11, 45, 197, 903]) = 2^6 = 64. - Philippe Deléham, Mar 02 2004
a(n) = 2^floor(n^2/2)/2^floor(n/2). - Paul Barry, Oct 04 2004
G.f. satisfies: A(x) = 1 + x*A(2x). - Paul D. Hanna, Dec 04 2009
a(n) = 2 * a(n-1)^2 / a(n-2). - Michael Somos, Dec 30 2012
G.f.: G(0)/x - 1/x, where G(k) = 1 + 2^(k-1)*x/(1 - 1/(1 + 1/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jul 26 2013
E.g.f. satisfies A'(x) = A(2x). - Geoffrey Critzer, Sep 07 2013
Sum_{n>=1} 1/a(n) = A299998. - Amiram Eldar, Oct 27 2020
a(n) = s_lambda(1,1,...,1) where s is the Schur polynomial in n variables and lambda is the partition (n,n-1,n-2,...,1). - Leonid Bedratyuk, Feb 06 2022
a(n) = Product_{1 <= j <= i <= n-1} (i + j)/(2*i - 2*j + 1). Cf. A007685. - Peter Bala, Oct 25 2024

Extensions

More terms from Vladeta Jovovic, Apr 09 2000

A053763 a(n) = 2^(n^2 - n).

Original entry on oeis.org

1, 1, 4, 64, 4096, 1048576, 1073741824, 4398046511104, 72057594037927936, 4722366482869645213696, 1237940039285380274899124224, 1298074214633706907132624082305024, 5444517870735015415413993718908291383296, 91343852333181432387730302044767688728495783936
Offset: 0

Views

Author

Stephen G Penrice, Mar 29 2000

Keywords

Comments

Nilpotent n X n matrices over GF(2). Also number of simple digraphs (without self-loops) on n labeled nodes (see also A002416).
For n >= 1 a(n) is the size of the Sylow 2-subgroup of the Chevalley group A_n(4) (sequence A053291). - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 30 2001
(-1)^ceiling(n/2) * resultant of the Chebyshev polynomial of first kind of degree n and Chebyshev polynomial of first kind of degree (n+1) (cf. A039991). - Benoit Cloitre, Jan 26 2003
The number of reflexive binary relations on an n-element set. - Justin Witt (justinmwitt(AT)gmail.com), Jul 12 2005
From Rick L. Shepherd, Dec 24 2008: (Start)
Number of gift exchange scenarios where, for each person k of n people,
i) k gives gifts to g(k) of the others, where 0 <= g(k) <= n-1,
ii) k gives no more than one gift to any specific person,
iii) k gives no single gift to two or more people and
iv) there is no other person j such that j and k jointly give a single gift.
(In other words -- but less precisely -- each person k either gives no gifts or gives exactly one gift per person to 1 <= g(k) <= n-1 others.) (End)
In general, sequences of the form m^((n^2 - n)/2) enumerate the graphs with n labeled nodes with m types of edge. a(n) therefore is the number of labeled graphs with n nodes with 4 types of edge. To clarify the comment from Benoit Cloitre, dated Jan 26 2003, in this context: simple digraphs (without self-loops) have four types of edge. These types of edges are as follows: the absent edge, the directed edge from A -> B, the directed edge from B -> A and the bidirectional edge, A <-> B. - Mark Stander, Apr 11 2019

Examples

			a(2)=4 because there are four 2 x 2 nilpotent matrices over GF(2):{{0,0},{0,0}},{{0,1},{0,0}},{{0,0},{1,0}},{{1,1,},{1,1}} where 1+1=0. - _Geoffrey Critzer_, Oct 05 2012
		

References

  • J. L. Gross and J. Yellen, eds., Handbook of Graph Theory, CRC Press, 2004; p. 521.
  • F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973, p. 5, Eq. (1.1.5).

Crossrefs

Programs

Formula

Sequence given by the Hankel transform (see A001906 for definition) of A059231 = {1, 1, 5, 29, 185, 1257, 8925, 65445, 491825, ...}; example: det([1, 1, 5, 29; 1, 5, 29, 185; 5, 29, 185, 1257; 29, 185, 1257, 8925]) = 4^6 = 4096. - Philippe Deléham, Aug 20 2005
a(n) = 4^binomial(n, n-2). - Zerinvary Lajos, Jun 16 2007
a(n) = Sum_{i=0..n^2-n} binomial(n^2-n, i). - Rick L. Shepherd, Dec 24 2008
G.f. A(x) satisfies: A(x) = 1 + x * A(4*x). - Ilya Gutkovskiy, Jun 04 2020
Sum_{n>=1} 1/a(n) = A319016. - Amiram Eldar, Oct 27 2020
Sum_{n>=0} a(n)*u^n/A002884(n) = Product_{r>=1} 1/(1-u/q^r). - Geoffrey Critzer, Oct 28 2021

A000595 Number of binary relations on n unlabeled points.

Original entry on oeis.org

1, 2, 10, 104, 3044, 291968, 96928992, 112282908928, 458297100061728, 6666621572153927936, 349390545493499839161856, 66603421985078180758538636288, 46557456482586989066031126651104256, 120168591267113007604119117625289606148096, 1152050155760474157553893461743236772303142428672
Offset: 0

Views

Author

Keywords

Comments

Number of orbits under the action of permutation group S(n) on n X n {0,1} matrices. The action is defined by f.M(i,j)=M(f(i),f(j)).
Equivalently, the number of digraphs on n unlabeled nodes with loops allowed but no more than one arc with the same start and end node. - Andrew Howroyd, Oct 22 2017

Examples

			From _Gus Wiseman_, Jun 17 2019: (Start)
Non-isomorphic representatives of the a(2) = 10 relations:
  {}
  {1->1}
  {1->2}
  {1->1, 1->2}
  {1->1, 2->1}
  {1->1, 2->2}
  {1->2, 2->1}
  {1->1, 1->2, 2->1}
  {1->1, 1->2, 2->2}
  {1->1, 1->2, 2->1, 2->2}
(End)
		

References

  • F. Bergeron, G. Labelle and P. Leroux, Combinatorial Species and Tree-Like Structures, Cambridge, 1998, p. 76 (2.2.30)
  • M. D. McIlroy, Calculation of numbers of structures of relations on finite sets, Massachusetts Institute of Technology, Research Laboratory of Electronics, Quarterly Progress Reports, No. 17, Sept. 15, 1955, pp. 14-22.
  • 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

Programs

  • GAP
    NSeq := function ( n ) return Sum(List(ConjugacyClasses(SymmetricGroup(n)), c -> (2^Length(Orbits(Group(Representative(c)), CartesianProduct([1..n],[1..n]), OnTuples))) * Size(c)))/Factorial(n); end; # Dan Hoey, May 04 2001
    
  • Mathematica
    Join[{1,2}, Table[CycleIndex[Join[PairGroup[SymmetricGroup[n],Ordered], Permutations[Range[n^2-n+1,n^2]],2],s] /. Table[s[i]->2, {i,1,n^2-n}], {n,2,7}]] (* Geoffrey Critzer, Nov 02 2011 *)
    permcount[v_] := Module[{m=1, s=0, k=0, t}, For[i=1, i <= Length[v], i++, t = v[[i]]; k = If[i>1 && t == v[[i-1]], k+1, 1]; m *= t*k; s += t]; s!/m];
    edges[v_] := Sum[2*GCD[v[[i]], v[[j]]], {i, 2, Length[v]}, {j, 1, i - 1}] + Total[v];
    a[n_] := (s=0; Do[s += permcount[p]*2^edges[p], {p, IntegerPartitions[n]}]; s/n!);
    Table[a[n], {n, 0, 15}] (* Jean-François Alcover, Jul 08 2018, after Andrew Howroyd *)
    dinorm[m_]:=If[m=={},{},If[Union@@m!=Range[Max@@Flatten[m]],dinorm[m/.Apply[Rule,Table[{(Union@@m)[[i]],i},{i,Length[Union@@m]}],{1}]],First[Sort[dinorm[m,1]]]]];
    dinorm[m_,aft_]:=If[Length[Union@@m]<=aft,{m},With[{mx=Table[Count[m,i,{2}],{i,Select[Union@@m,#1>=aft&]}]},Union@@(dinorm[#1,aft+1]&)/@Union[Table[Map[Sort,m/.{par+aft-1->aft,aft->par+aft-1},{0}],{par,First/@Position[mx,Max[mx]]}]]]];
    Table[Length[Union[dinorm/@Subsets[Tuples[Range[n],2]]]],{n,0,3}] (* Gus Wiseman, Jun 17 2019 *)
  • 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}
    edges(v) = {sum(i=2, #v, sum(j=1, i-1, 2*gcd(v[i],v[j]))) + sum(i=1, #v, v[i])}
    a(n) = {my(s=0); forpart(p=n, s+=permcount(p)*2^edges(p)); s/n!} \\ Andrew Howroyd, Oct 22 2017
    
  • Python
    from itertools import product
    from math import prod, factorial, gcd
    from fractions import Fraction
    from sympy.utilities.iterables import partitions
    def A000595(n): return int(sum(Fraction(1<Chai Wah Wu, Jul 02 2024

Formula

a(n) = sum {1*s_1+2*s_2+...=n} (fixA[s_1, s_2, ...] / (1^s_1*s_1!*2^s_2*s_2!*...)) where fixA[s_1, s_2, ...] = 2^sum {i, j>=1} (gcd(i, j)*s_i*s_j). - Christian G. Bower, Jan 05 2004
a(n) ~ 2^(n^2)/n! [McIlroy, 1955]. - Vaclav Kotesovec, Dec 19 2016

Extensions

More terms from Vladeta Jovovic, Feb 07 2000
Still more terms from Dan Hoey, May 04 2001

A060690 a(n) = binomial(2^n + n - 1, n).

Original entry on oeis.org

1, 2, 10, 120, 3876, 376992, 119877472, 131254487936, 509850594887712, 7145544812472168960, 364974894538906616240640, 68409601066028072105113098240, 47312269462735023248040155132636160, 121317088003402776955124829814219234385920
Offset: 0

Views

Author

Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 19 2001

Keywords

Comments

Also the number of n X n (0,1) matrices modulo rows permutation (by symmetry this is the same as the number of (0,1) matrices modulo columns permutation), i.e., the number of equivalence classes where two matrices A and B are equivalent if one of them is the result of permuting the rows of the other. The total number of (0,1) matrices is in sequence A002416.
Row sums of A220886. - Geoffrey Critzer, Nov 20 2014

Crossrefs

Sequences of the form binomial(2^n +p*n +q, n): A136556 (0,-1), A014070 (0,0), A136505 (0,1), A136506 (0,2), this sequence (1,-1), A132683 (1,0), A132684 (1,1), A132685 (2,0), A132686 (2,1), A132687 (3,-1), A132688 (3,0), A132689 (3,1).
Main diagonal of A092056.
Central terms of A137153.

Programs

  • Magma
    [Binomial(2^n +n-1, n): n in [0..20]]; // G. C. Greubel, Mar 14 2021
    
  • Maple
    with(combinat): for n from 0 to 20 do printf(`%d,`,binomial(2^n+n-1, n)) od:
  • Mathematica
    Table[Binomial[2^n+n-1,n],{n,0,20}] (* Harvey P. Dale, Apr 19 2012 *)
  • PARI
    a(n)=binomial(2^n+n-1,n)
    
  • PARI
    {a(n)=polcoeff(sum(k=0,n,(-log(1-2^k*x +x*O(x^n)))^k/k!),n)} \\ Paul D. Hanna, Dec 29 2007
    
  • PARI
    a(n) = sum(k=0, n, stirling(n,k,1)*(2^n+n-1)^k/n!); \\ Paul D. Hanna, Nov 20 2014
    
  • Python
    from math import comb
    def A060690(n): return comb((1<Chai Wah Wu, Jul 05 2024
  • Sage
    [binomial(2^n +n-1, n) for n in (0..20)] # G. C. Greubel, Mar 14 2021
    

Formula

a(n) = [x^n] 1/(1-x)^(2^n).
a(n) = (1/n!)*Sum_{k=0..n} ( (-1)^(n-k)*Stirling1(n, k)*2^(k*n) ). - Vladeta Jovovic, May 28 2004
a(n) = Sum_{k=0..n} (-1)^(n-k)*binomial(2^n+n,k) - Vladeta Jovovic, Jan 21 2008
a(n) = Sum_{k=0..n} Stirling1(n,k)*(2^n+n-1)^k/n!. - Vladeta Jovovic, Jan 21 2008
G.f.: A(x) = Sum_{n>=0} [ -log(1 - 2^n*x)]^n / n!. More generally, Sum_{n>=0} [ -log(1 - q^n*x)]^n/n! = Sum_{n>=0} C(q^n+n-1,n)*x^n ; also Sum_{n>=0} log(1 + q^n*x)^n/n! = Sum_{n>=0} C(q^n,n)*x^n. - Paul D. Hanna, Dec 29 2007
a(n) ~ 2^(n^2) / n!. - Vaclav Kotesovec, Jul 02 2016
a(n) = A163767(2^n). - Alois P. Heinz, Jun 12 2024

Extensions

More terms from James Sellers, Apr 20 2001
Edited by N. J. A. Sloane, Mar 17 2008

A064062 Generalized Catalan numbers C(2; n).

Original entry on oeis.org

1, 1, 3, 13, 67, 381, 2307, 14589, 95235, 636925, 4341763, 30056445, 210731011, 1493303293, 10678370307, 76957679613, 558403682307, 4075996839933, 29909606989827, 220510631755773, 1632599134961667, 12133359132082173
Offset: 0

Views

Author

Wolfdieter Lang, Sep 13 2001

Keywords

Comments

a(n+1) = Y_{n}(n+1) = Z_{n}, n >= 0, in the Derrida et al. 1992 reference (see A064094) for alpha=2, beta=1 (or alpha=1, beta=2).
a(n) = number of Dyck n-paths (A000108) in which each upstep (U) not at ground level is colored red (R) or blue (B). For example, a(3)=3 counts URDD, UBDD, UDUD (D=downstep). - David Callan, Mar 30 2007
The Hankel transform of this sequence is A002416. - Philippe Deléham, Nov 19 2007
The sequence a(n)/2^n, with g.f. 1/(1-xc(x)/2), has Hankel transform 1/2^n. - Paul Barry, Apr 14 2008
The REVERT transform of the odd numbers [1,3,5,7,9,...] is [1, -3, 13, -67, 381, -2307, 14589, -95235, 636925, ...] - N. J. A. Sloane, May 26 2017

Crossrefs

Generalized Catalan numbers C(m; n): A000012 (m = 0), A000108 (m = 1), A064063 (m = 3) and A064087 - A064093 (m = 4 thru 10); A064310 (m = -1), A064311 (m = -2) and A064325 - A064333 (m = -3 thru -11).

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 30);
    Coefficients(R!( (3 - Sqrt(1-8*x))/(2*(1+x)) )); // G. C. Greubel, Sep 27 2024
  • Maple
    1, seq(simplify(hypergeom([1-n,n],[-n],2)), n=1..100); # Robert Israel, Nov 30 2014
  • Mathematica
    a[0]=1; a[1]=1; a[n_]/;n>=2 := a[n] = a[n-1] + Sum[(a[k] + a[k-1])a[n-k],{k,n-1}]; Table[a[n],{n,0,10}] (* David Callan, Aug 27 2009 *)
    a[n_] := 2*Sum[ (-1)^j*2^(n-j-1)*Binomial[2*(n-j-1), n-j-1]/(n-j), {j, 0, n-1}] + (-1)^n; Table[a[n], {n, 0, 21}] (* Jean-François Alcover, Jul 03 2013 *)
  • PARI
    {a(n)=polcoeff((3-sqrt(1-8*x+x*O(x^n)))/(2+2*x),n)}
    
  • PARI
    {a(n)=local(A=1+x); for(i=1, n, A=1+A^4*intformal(1/(A^2+x*O(x^n)))); polcoeff(A, n)} \\ Paul D. Hanna, Dec 24 2013
    for(n=0, 25, print1(a(n), ", "))
    
  • PARI
    {a(n)=polcoeff(1/(1 - serreverse(x-2*x^2 +x^2*O(x^n))),n)}
    for(n=0,30,print1(a(n),", ")) \\ Paul D. Hanna, Nov 30 2014
    
  • Sage
    def a(n):
        if n==0: return 1
        return hypergeometric([1-n, n], [-n], 2).simplify()
    [a(n) for n in range(22)] # Peter Luschny, Dec 01 2014
    

Formula

G.f.: (1 + 2*x*C(2*x)) / (1+x) = 1/(1 - x*C(2*x)) with C(x) g.f. of Catalan numbers A000108.
a(n) = A062992(n-1) = Sum_{m = 0..n-1} (n-m)*binomial(n-1+m, m)*(2^m)/n, n >= 1, a(0) = 1.
a(n) = Sum_{k = 0..n} A059365(n, k)*2^(n-k). - Philippe Deléham, Jan 19 2004
G.f.: 1/(1-x/(1-2x/(1-2x/(1-2x/(1-.... = 1/(1-x-2x^2/(1-4x-4x^2/(1-4x-4x^2/(1-.... (continued fractions). - Paul Barry, Jan 30 2009
a(n) = (32/Pi)*Integral_{x = 0..1} (8*x)^(n-1)*sqrt(x*(1-x)) / (8*x+1). - Groux Roland, Dec 12 2010
a(n+2) = 8^(n+2)*( c(n+2)-c(1)*c(n+1) - Sum_{i=0..n-1} 8^(-i-2)*c(n-i)*a(i+2) ) with c(n) = Catalan(n+2)/2^(2*n+1). - Groux Roland, Dec 12 2010
a(n) = the upper left term in M^n, M = the production matrix:
1, 1
2, 2, 1
4, 4, 2, 1
8, 8, 4, 2, 1
... - Gary W. Adamson, Jul 08 2011
D-finite with recurrence: n*a(n) + (12-7n)*a(n-1) + 4*(3-2n)*a(n-2) = 0. - R. J. Mathar, Nov 16 2011 (This follows easily from the generating function. - Robert Israel, Nov 30 2014)
G.f. satisfies: A(x) = 1 + A(x)^4 * Integral 1/A(x)^2 dx. - Paul D. Hanna, Dec 24 2013
G.f. satisfies: Integral 1/A(x)^2 dx = x - x^2*G(x), where G(x) is the o.g.f. of A000257, the number of rooted bicubic maps. - Paul D. Hanna, Dec 24 2013
G.f. A(x) satisfies: A(x - 2*x^2) = 1/(1-x). - Paul D. Hanna, Nov 30 2014
a(n) = hypergeometric([1-n, n], [-n], 2) for n > 0. - Peter Luschny, Nov 30 2014
G.f.: (3 - sqrt(1-8*x))/(2*(x+1)). - Robert Israel, Nov 30 2014
a(n) ~ 2^(3*n+1) / (9*sqrt(Pi)*n^(3/2)). - Vaclav Kotesovec, Dec 22 2014
O.g.f. A(x) = 1 + series reversion of (x*(1 - x)/(1 + x)^2). Logarithmically differentiating (A(x) - 1)/x gives 3 + 17*x + 111*x^2 + ..., essentially a g.f for A119259. - Peter Bala, Oct 01 2015
From Peter Bala, Jan 06 2022: (Start)
exp( Sum_{n >= 1} a(n)*x^n/n ) = 1 + x + 2*x^2 + 6*x^3 + 23*x^4 + ... is a g.f. for A022558.
The Gauss congruences a(n*p^k) == a(n^p^(k-1)) (mod p^k) hold for prime p and positive integers n and k. (End)

A117401 Triangle T(n,k) = 2^(k*(n-k)), read by rows.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 4, 4, 1, 1, 8, 16, 8, 1, 1, 16, 64, 64, 16, 1, 1, 32, 256, 512, 256, 32, 1, 1, 64, 1024, 4096, 4096, 1024, 64, 1, 1, 128, 4096, 32768, 65536, 32768, 4096, 128, 1, 1, 256, 16384, 262144, 1048576, 1048576, 262144, 16384, 256, 1
Offset: 0

Views

Author

Paul D. Hanna, Mar 12 2006

Keywords

Comments

Matrix power T^m satisfies: [T^m](n,k) = [T^m](n-k,0)*T(n,k) for all m and so the triangle has an invariant character.

Examples

			A(x,y) = 1/(1-xy) + x/(1-2xy) + x^2/(1-4xy) + x^3/(1-8xy) + ...
Triangle begins:
  1;
  1,   1;
  1,   2,     1;
  1,   4,     4,      1;
  1,   8,    16,      8,       1;
  1,  16,    64,     64,      16,       1;
  1,  32,   256,    512,     256,      32,      1;
  1,  64,  1024,   4096,    4096,    1024,     64,     1;
  1, 128,  4096,  32768,   65536,   32768,   4096,   128,   1;
  1, 256, 16384, 262144, 1048576, 1048576, 262144, 16384, 256, 1;
		

Crossrefs

Cf. A117402 (row sums), A117403 (antidiagonal sums), A002416 (central terms).
Cf. this sequence (m=0), A118180 (m=1), A118185 (m=2), A118190 (m=3), A158116 (m=4), A176642 (m=6), A158117 (m=8), A176627 (m=10), A176639 (m=13), A156581 (m=15).

Programs

  • Magma
    A117401:= func< n, k, m | (m+2)^(k*(n-k)) >;
    [A117401(n, k, 0): k in [0..n], n in [0..12]]; // G. C. Greubel, Jun 28 2021
    
  • Mathematica
    Table[2^((n-k)k),{n,0,10},{k,0,n}]//Flatten (* Harvey P. Dale, Jan 09 2017 *)
  • PARI
    T(n,k)=if(n
    				
  • Sage
    def A117401(n, k, m): return (m+2)^(k*(n-k))
    flatten([[A117401(n, k, 0) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Jun 28 2021

Formula

G.f.: A(x,y) = Sum_{n>=0} x^n/(1 - 2^n*x*y).
G.f. satisfies: A(x,y) = 1/(1 - x*y) + x*A(x,2*y).
Equals ConvOffsStoT transform of the 2^n series: (1, 2, 4, 8, ...); e.g., ConvOffs transform of (1, 2, 4, 8) = (1, 8, 16, 8, 1). - Gary W. Adamson, Apr 21 2008
T(n,k) = (1/n)*( 2^(n-k)*k*T(n-1,k-1) + 2^k*(n-k)*T(n-1,k) ), where T(i,j)=0 if j>i. - Tom Edgar, Feb 20 2014
Let E(x) = Sum_{n>=0} x^n/2^C(n,2). Then E(x)*E(y*x) = Sum_{n>=0} Sum_{k=0..n} T(n,k)*y^k*x^n/2^C(n,2). - Geoffrey Critzer, May 31 2020
T(n, k, m) = (m+2)^(k*(n-k)) with m = 0. - G. C. Greubel, Jun 28 2021

A055165 Number of invertible n X n matrices with entries equal to 0 or 1.

Original entry on oeis.org

1, 1, 6, 174, 22560, 12514320, 28836612000, 270345669985440, 10160459763342013440
Offset: 0

Views

Author

Ulrich Hermisson (uhermiss(AT)server1.rz.uni-leipzig.de), Jun 18 2000

Keywords

Comments

All eigenvalues are nonzero.

Examples

			For n=2 the 6 matrices are {{{0, 1}, {1, 0}}, {{0, 1}, {1, 1}}, {{1, 0}, {0, 1}}, {{1, 0}, {1, 1}}, {{1, 1}, {0, 1}}, {{1, 1}, {1, 0}}}.
		

Crossrefs

Cf. A056990, A056989, A046747, A055165, A002416, A003024 (positive definite matrices).
A046747(n) + a(n) = 2^(n^2) = total number of n X n (0, 1) matrices = sequence A002416.
Main diagonal of A064230.

Programs

  • PARI
    a(n)=sum(t=0,2^n^2-1,!!matdet(matrix(n,n,i,j,(t>>(i*n+j-n-1))%2))) \\ Charles R Greathouse IV, Feb 09 2016
    
  • Python
    from itertools import product
    from sympy import Matrix
    def A055165(n): return sum(1 for s in product([0,1],repeat=n**2) if Matrix(n,n,s).det() != 0) # Chai Wah Wu, Sep 24 2021

Formula

For an asymptotic estimate see A046747. A002884 is a lower bound. A002416 is an upper bound.
a(n) = n! * A088389(n). - Gerald McGarvey, Oct 20 2007

Extensions

More terms from Miodrag Zivkovic (ezivkovm(AT)matf.bg.ac.rs), Feb 28 2006
Description improved by Jeffrey Shallit, Feb 17 2016
a(0)=1 prepended by Alois P. Heinz, Jun 18 2022

A046747 Number of n X n rational {0,1}-matrices of determinant 0.

Original entry on oeis.org

1, 10, 338, 42976, 21040112, 39882864736, 292604283435872, 8286284310367538176
Offset: 1

Views

Author

Günter M. Ziegler (ziegler(AT)math.tu-berlin.de)

Keywords

Examples

			a(2)=10: the matrix of all 0's, 4 matrices with 2 zeros in the same row or column, 4 matrices with 3 zeros and the all-1 matrix.
		

Crossrefs

Programs

  • Mathematica
    Sum[KroneckerDelta[Det[Array[Mod[Floor[k/(2^(n*(#1-1)+#2-1))],2]&,{n,n}]],0],{k,0,(2^(n^2))-1}] (* John M. Campbell, Jun 24 2011 *)
    Count[Det /@ Tuples[{0, 1}, {n, n}], 0] (* David Trimas, Sep 23 2024 *)
  • PARI
    A046747(n) = m=matrix(n,n); ct=0; for(x=0,2^(n*n)-1,a=binary(x+2^(n*n)); for(i=1,n, for(j=1,n,m[i,j]=a[n*i+j+1-n])); if(matdet(m)==0,ct=ct+1,); ); ct \\ Randall L Rathbun
    
  • PARI
    a(n)=sum(i=0,2^n^2-1,matdet(matrix(n,n,x,y,(i>>(n*x+y-n-1))%2))==0) \\ Charles R Greathouse IV, Feb 21 2015

Formula

a(n) = 2^(n^2) - n! * binomial(2^n -1, n) + n! * A000410(n).
a(n) + A055165(n) = 2^(n^2) = total number of n X n (0, 1) matrices.
The probability that a random n X n {0,1}-matrix is singular is conjectured to be asymptotic to C(n+1, 2)*(1/2)^(n-1). [Corrected by N. J. A. Sloane, Jan 02 2007]

Extensions

a(8) from Vladeta Jovovic, Mar 28 2006

A326209 Number of nesting labeled digraphs with vertices {1..n}.

Original entry on oeis.org

0, 0, 4, 408, 64528
Offset: 0

Views

Author

Gus Wiseman, Jun 19 2019

Keywords

Comments

Two edges (a,b), (c,d) are nesting if a < c and b > d or a > c and b < d.
Also unsortable digraphs with vertices {1..n}, where a digraph is sortable if, when the edges are listed in lexicographic order, their targets are weakly increasing.
Also the number of semicrossing digraphs with vertices {1..n}, where two edges (a,b), (c,d) are semicrossing if a < c and b < d or a > c and b > d. For example, the a(2) = 4 semicrossing digraph edge-sets are:
{11,22}
{11,12,22}
{11,21,22}
{11,12,21,22}

Examples

			The a(2) = 4 nesting digraph edge-sets:
  {12,21}
  {11,12,21}
  {12,21,22}
  {11,12,21,22}
		

Crossrefs

Non-nesting digraphs are A326237.
Nesting set partitions are A016098.
MM-numbers of nesting multiset partitions are A326256.
MM-numbers of unsortable multiset partitions are A326258.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Tuples[Range[n],2]],!OrderedQ[Last/@#]&]],{n,4}]

Formula

A002416(n) = a(n) + A326237(n).
Showing 1-10 of 126 results. Next