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

A171873 Column sums of triangle A171871.

Original entry on oeis.org

1, 1, 2, 10, 280, 1173468
Offset: 0

Views

Author

Robert Munafo, Jan 21 2010

Keywords

Comments

Next term is known to be greater than 220146725295227, based on link between A171871 and A039754.

A171875 Subdiagonal of triangle A171871: Classifications of N elements containing exactly N-2 binary partitions.

Original entry on oeis.org

0, 0, 1, 3, 17, 74, 358, 1631, 7563, 34751, 160807
Offset: 2

Views

Author

Robert Munafo, Jan 21 2010

Keywords

Crossrefs

Cf. A171871.

A000055 Number of trees with n unlabeled nodes.

Original entry on oeis.org

1, 1, 1, 1, 2, 3, 6, 11, 23, 47, 106, 235, 551, 1301, 3159, 7741, 19320, 48629, 123867, 317955, 823065, 2144505, 5623756, 14828074, 39299897, 104636890, 279793450, 751065460, 2023443032, 5469566585, 14830871802, 40330829030, 109972410221, 300628862480, 823779631721, 2262366343746, 6226306037178
Offset: 0

Views

Author

Keywords

Comments

Also, number of unlabeled 2-gonal 2-trees with n-1 2-gons, for n>0. [Corrected by Andrei Zabolotskii, Jul 29 2025]
Main diagonal of A054924.
Left border of A157905. - Gary W. Adamson, Mar 08 2009
From Robert Munafo, Jan 24 2010: (Start)
Also counts classifications of n items that require exactly n-1 binary partitions; see Munafo link at A005646, also A171871 and A171872.
The 11 trees for n = 7 are illustrated at the Munafo web link.
Link to A171871/A171872 conjectured by Robert Munafo, then proved by Andrew Weimholt and Franklin T. Adams-Watters on Dec 29 2009. (End)
This is also "Number of tree perfect graphs on n nodes" [see Hougardy]. - N. J. A. Sloane, Dec 04 2015
For n > 0, a(n) is the number of ways to arrange n-1 unlabeled non-intersecting circles on a sphere. - Vladimir Reshetnikov, Aug 25 2016
All trees for n=1 through n=12 are depicted in Chapter 1 of the Steinbach reference. On p. 6 appear encircled two trees (with n=10) which seem inequivalent only when considered as ordered (planar) trees. Earlier instances of such possibly (in)equivalent trees could appear from n=6 on (and from n=9 on without equivalence modulo plane symmetry) but are not drawn separately there. - M. F. Hasler, Aug 29 2017

Examples

			a(1) = 1 [o]; a(2) = 1 [o-o]; a(3) = 1 [o-o-o];
a(4) = 2 [o-o-o and o-o-o-o]
            |
            o
G.f. = 1 + x + x^2 + x^3 + 2*x^4 + 3*x^5 + 6*x^6 + 11*x^7 + 23*x^8 + ...
		

References

  • F. Bergeron, G. Labelle and P. Leroux, Combinatorial Species and Tree-Like Structures, Camb. 1998, p. 279.
  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 55.
  • N. L. Biggs et al., Graph Theory 1736-1936, Oxford, 1976, p. 49.
  • A. Cayley, On the analytical forms called trees, with application to the theory of chemical combinations, Reports British Assoc. Advance. Sci. 45 (1875), 257-305 = Math. Papers, Vol. 9, 427-460 (see p. 459).
  • S. R. Finch, Mathematical Constants, Cambridge, 2003, pp. 295-316.
  • J. L. Gross and J. Yellen, eds., Handbook of Graph Theory, CRC Press, 2004; p. 526.
  • F. Harary, Graph Theory. Addison-Wesley, Reading, MA, 1969, p. 232.
  • F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973, p. 58 and 244.
  • D. E. Knuth, Fundamental Algorithms, 3d Ed. 1997, pp. 386-88.
  • R. C. Read and R. J. Wilson, An Atlas of Graphs, Oxford, 1998.
  • J. Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 138.
  • 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. A000676 (centered trees), A000677 (bicentered trees), A027416 (trees with a centroid), A102011 (trees with a bicentroid), A034853 (refined by diameter), A238414 (refined by maximum vertex degree).
Cf. A000081 (rooted trees), A000272 (labeled trees), A000169 (labeled rooted trees), A212809 (radius of convergence).
Cf. A036361 (labeled 2-trees), A036362 (labeled 3-trees), A036506 (labeled 4-trees), A054581 (unlabeled 2-trees).
Cf. A157904, A157905, A005195 (Euler transform = forests), A095133 (multisets).
Column 0 of A335362 and A034799.
Related to A005646; see A171871 and A171872.

Programs

  • Haskell
    import Data.List (generic_index)
    import Math.OEIS (getSequenceByID)
    triangle x = (x * x + x) `div` 2
    a000055 n = let {r = genericIndex (fromJust (getSequenceByID "A000081")); (m, nEO) = divMod n 2}
                in  r n - sum (zipWith (*) (map r [0..m]) (map r [n, n-1..]))
                    + (1-nEO) * (triangle (r m + 1))
    -- Walt Rorie-Baety, Jun 12 2021
    
  • Magma
    N := 30; P := PowerSeriesRing(Rationals(),N+1); f := func< A | x*&*[Exp(Evaluate(A,x^k)/k) : k in [1..N]]>; G := x; for i in [1..N] do G := f(G); end for; G000081 := G; G000055 := 1 + G - G^2/2 + Evaluate(G,x^2)/2; A000055 := Eltseq(G000055); // Geoff Baileu (geoff(AT)maths.usyd.edu.au), Nov 30 2009
    
  • Maple
    G000055 := series(1+G000081-G000081^2/2+subs(x=x^2,G000081)/2,x,31); A000055 := n->coeff(G000055,x,n); # where G000081 is g.f. for A000081 starting with n=1 term
    with(numtheory): b:= proc(n) option remember; `if`(n<=1, n, (add(add(d*b(d), d=divisors(j)) *b(n-j), j=1..n-1))/ (n-1)) end: a:= n-> `if`(n=0, 1, b(n) -(add(b(k) *b(n-k), k=0..n) -`if`(irem(n, 2)=0, b(n/2), 0))/2):
    seq(a(n), n=0..50);
    # Alois P. Heinz, Aug 21 2008
    # Program to create b-file b000055.txt:
    A000081 := proc(n) option remember; local d, j;
    if n <= 1 then n else
        add(add(d*procname(d),d=numtheory[divisors](j))*procname(n-j),j=1..n-1)/(n-1);
    fi end:
    A000055 := proc(nmax) local a81, n, t, a, j, i ;
    a81 := [seq(A000081(i), i=0..nmax)] ; a := [] ;
    for n from 0 to nmax do
        if n = 0 then
            t := 1+op(n+1, a81) ;
        else
            t := op(n+1, a81) ;
        fi;
        if type(n, even) then
            t := t-op(1+n/2, a81)^2/2 ;
            t := t+op(1+n/2, a81)/2 ;
        fi;
        for j from 0 to (n-1)/2 do
            t := t-op(j+1, a81)*op(n-j+1, a81) ;
        od:
        a := [op(a), t] ;
    od:
    a end:
    L := A000055(1000) ;
    #  R. J. Mathar, Mar 06 2009
  • Mathematica
    s[n_, k_] := s[n, k] = a[n + 1 - k] + If[n < 2k, 0, s[n - k, k]]; a[1] = 1; a[n_] := a[n] = Sum[a[i] s[n-1, i] i, {i, 1, n-1}] / (n-1); Table[a[i] - Sum[a[j] a[i-j], {j, 1, i/2}] + If[OddQ[i], 0, a[i/2] (a[i/2] + 1)/2], {i, 1, 50}] (* Robert A. Russell *)
    b[0] = 0; b[1] = 1; b[n_] := b[n] = Sum[d*b[d]*b[n-j], {j, 1, n-1}, {d, Divisors[j]}]/(n-1); a[0] = 1; a[n_] := b[n] - (Sum[b[k]*b[n-k], {k, 0, n}] - If[Mod[n, 2] == 0, b[n/2], 0])/2; Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Apr 09 2014, after Alois P. Heinz *)
  • PARI
    {a(n) = local(A, A1, an, i, t); if( n<2, n>=0, an = Vec(A = A1 = 1 + O('x^n)); for(m=2, n, i=m\2; an[m] = sum(k=1, i, an[k] * an[m-k]) + (t = polcoeff( if( m%2, A *= (A1 - 'x^i)^-an[i], A), m-1))); t + if( n%2==0, binomial( -polcoeff(A, i-1), 2)))}; /* Michael Somos */
    
  • PARI
    N=66;  A=vector(N+1, j, 1);
    for (n=1, N, A[n+1] = 1/n * sum(k=1, n, sumdiv(k, d, d * A[d]) * A[n-k+1] ) );
    A000081=concat([0], A);
    H(t)=subst(Ser(A000081, 't), 't, t);
    x='x+O('x^N);
    Vec( 1 + H(x) - 1/2*( H(x)^2 - H(x^2) ) )
    \\ Joerg Arndt, Jul 10 2014
    
  • Python
    # uses function from A000081
    def A000055(n): return 1 if n == 0 else A000081(n)-sum(A000081(i)*A000081(n-i) for i in range(1,n//2+1)) + (0 if n % 2 else (A000081(n//2)+1)*A000081(n//2)//2) # Chai Wah Wu, Feb 03 2022
  • SageMath
    [len(list(graphs.trees(n))) for n in range(16)] # Peter Luschny, Mar 01 2020
    

Formula

G.f.: A(x) = 1 + T(x) - T^2(x)/2 + T(x^2)/2, where T(x) = x + x^2 + 2*x^3 + ... is the g.f. for A000081.
a(n) ~ A086308 * A051491^n * n^(-5/2). - Vaclav Kotesovec, Jan 04 2013
a(n) = A000081(n) - A217420(n+1), n > 0. - R. J. Mathar, Sep 19 2016
a(n) = A000676(n) + A000677(n). - R. J. Mathar, Aug 13 2018
a(n) = A000081(n) - (Sum_{1<=i<=j, i+j=n} A000081(i)*A000081(j)) + (1-(-1)^(n-1)) * binomial(A000081(n/2)+1,2) / 2 [Li, equation 4.2]. - Walt Rorie-Baety, Jul 05 2021

A039754 Irregular triangle read by rows: T(n,k) = number of binary codes of length n with k words (n >= 0, 0 <= k <= 2^n); also number of 0/1-polytopes with vertices from the unit n-cube; also number of inequivalent Boolean functions of n variables with exactly k nonzero values under action of Jevons group.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 3, 3, 6, 3, 3, 1, 1, 1, 1, 4, 6, 19, 27, 50, 56, 74, 56, 50, 27, 19, 6, 4, 1, 1, 1, 1, 5, 10, 47, 131, 472, 1326, 3779, 9013, 19963, 38073, 65664, 98804, 133576, 158658, 169112, 158658, 133576, 98804, 65664, 38073, 19963, 9013, 3779, 1326, 472, 131, 47, 10, 5, 1, 1
Offset: 0

Views

Author

Keywords

Comments

For N=1 through N=5, the first 2^(N-1) terms of row N are also found in triangle A171871, which is related to A005646. This was shown for all N by Andrew Weimholt, Dec 30 2009. [Robert Munafo, Jan 25 2010]

Examples

			Triangle begins:
  k  0  1  2  3   4   5   6   7   8   9  10  11  12 13 14 15 16   sums
n
0    1  1                                                            2
1    1  1  1                                                         3
2    1  1  2  1   1                                                  6
3    1  1  3  3   6   3   3   1   1                                 22
4    1  1  4  6  19  27  50  56  74  56  50  27  19  6  4  1  1    402
		

References

  • F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973, p. 112.
  • M. A. Harrison, Introduction to Switching and Automata Theory. McGraw Hill, NY, 1965, p. 150.

Crossrefs

Row sums give A000616. Cf. A052265.
Rows give A034188, A034189, A034190, etc.
Columns are A034198, A034199, A034200, etc.
Diagonal is A276412.
For other versions of this triangle see A171876, A039754, A276777.
Cf. A171871.

Programs

  • Mathematica
    P = IntegerPartitions;
    AC[d_Integer] := Module[{C, M, p},(* from W. Y. C. Chen algorithm *) M[p_List] := Plus @@ p!/(Times @@ p * Times @@ (Length /@ Split[p]!)); C[p_List, q_List] := Module[{r, m, k, x}, r = If[0 == Length[q], 1, 2*2^IntegerExponent[LCM @@ q, 2]]; m = LCM @@ Join[p/GCD[r, p], q/GCD[r, q]]; CoefficientList[Expand[Product[(1 + x^(k *r))^((Plus @@ Map[MoebiusMu[k/#]*2^Plus @@ GCD[#*r, Join[p, q]]&, Divisors[k]])/(k*r)), {k, 1, m}]], x]]; Sum[Binomial[d, p]*Plus @@ Plus @@ Outer[M[#1] M[#2] C[#1, #2]*2^(d - Length[#1] - Length[#2]) &, P[p], P[d - p], 1], {p, 0, d}]/(d! 2^d)]; AC[0]  = {1, 1};
    AC /@ Range[0, 5] // Flatten (* Jean-François Alcover, Dec 15 2019, after Robert A. Russell in A034189 *)
    Table[ CoefficientList[ CycleIndexPolynomial[ GraphData[ {"Hypercube", n}, "AutomorphismGroup"], Array[Subscript[x, ##] &, 2^n]] /. Table[ Subscript[x, i] -> 1 + x^i, {i, 1, 2^n}], x], {n, 1,8}] // Grid (* Geoffrey Critzer, Jan 10 2020 *)

Formula

Reference gives g.f.
Fripertinger gives g.f. for the number of classes of (n, m) nonlinear codes over an alphabet of size A.

Extensions

Corrected and extended by Vladeta Jovovic, Apr 20 2000
Entry revised by N. J. A. Sloane, Sep 19 2016
T(0, 1) = 1 inserted. (There are two 0-ary functions.) - Tilman Piesk, Jan 10 2023

A005646 Number of classifications of n elements.

Original entry on oeis.org

1, 1, 1, 3, 6, 26, 122, 1015, 11847, 208914, 5236991, 184321511
Offset: 1

Views

Author

Keywords

Comments

"A 'classification' is a set of n type-specimens each one of which is corralled on its own by the union of a set of binary partitions, none of which could be omitted without leaving 2 types unseparated".
From Robert Munafo, Jan 24 2010: (Start)
Extensive explanation with illustrations on Munafo web page.
This sequence gives the row sums of triangle A171871. (End)

Examples

			Illustrated examples on Munafo web page. - _Robert Munafo_, Jan 24 2010
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • P. J. Wexler, On the number of taxonomies; or the odds on 'structuralism', American Anthropologist, 73 (1971), 1258.

Crossrefs

Cf. A000055, A171872, A171873. - Robert Munafo, Jan 24 2010

Extensions

1015 term first calculated by Andrew Weimholt, Dec 15 2009
11847 term first calculated by Andrew Weimholt, Dec 19 2009
208914 term first calculated by Robert Munafo, Dec 29 2009
5236990 term (erroneous) from Robert Munafo, Dec 30 2009
Erroneous "5236990" corrected to 5236991 by Robert Munafo, Jan 01 2010
184321511 term first calculated by Robert Munafo, Jan 10 2010

A171872 Triangle read by columns: Distinct classifications of N elements containing exactly R binary partitions.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 2, 3, 3, 1, 1, 0, 3, 17, 36, 60, 56, 50, 27, 19, 6, 4, 1, 1, 0, 6, 74, 573, 2802, 10087, 26512, 55088, 91984, 130267, 157662, 168890, 158658, 133576, 98804, 65664, 38073, 19963, 9013, 3779
Offset: 0

Views

Author

Robert Munafo, Jan 21 2010

Keywords

Comments

The bottom of each column is marked by a single 0 in this sequence. Value is 0 for all (N,R) for which N is greater than 2^R.
See note on efficient computation in A171871.

Crossrefs

First term in each column is A000055(R+1).
Column 4 shares terms with A034189, column 5 with A034190.

A171876 Mutual solutions to two classification counting problems: binary block codes of wordlength J with N used words; and classifications of N elements by J partitions.

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 3, 1, 1, 4, 6, 19, 27, 50, 56, 1, 1, 5, 10, 47, 131, 472, 1326, 3779, 9013, 19963, 38073, 65664, 98804, 133576, 158658, 1, 1, 6, 16, 103, 497, 3253, 19735, 120843, 681474, 3561696
Offset: 0

Views

Author

Robert Munafo, Jan 21 2010

Keywords

Comments

This connection was conjectured by Robert Munafo, then proved by Andrew Weimholt.
A(n) counts 2-colorings of a J-dimensional hypercube with N red vertices and 2^J-N black, each edge has at most one red vertex. - Andrew Weimholt, Dec 30 2009
This sequence contains terms of A039754 that are found in A171871/A171872. They occur in blocks of length 2^(J-1) as shown here:
1
1,1
1,1,3,3
1,1,4,6,19,27,50,56
1,1,5,10,47,131,472,1326,3779,9013,19963,38073,65664,98804,133576,158658

Crossrefs

Showing 1-7 of 7 results.