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

A122848 Exponential Riordan array (1, x(1+x/2)).

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 0, 3, 1, 0, 0, 3, 6, 1, 0, 0, 0, 15, 10, 1, 0, 0, 0, 15, 45, 15, 1, 0, 0, 0, 0, 105, 105, 21, 1, 0, 0, 0, 0, 105, 420, 210, 28, 1, 0, 0, 0, 0, 0, 945, 1260, 378, 36, 1, 0, 0, 0, 0, 0, 945, 4725, 3150, 630, 45, 1, 0, 0, 0, 0, 0, 0, 10395, 17325, 6930, 990, 55, 1, 0, 0
Offset: 0

Views

Author

Paul Barry, Sep 14 2006

Keywords

Comments

Entries are Bessel polynomial coefficients. Row sums are A000085. Diagonal sums are A122849. Inverse is A122850. Product of A007318 and A122848 gives A100862.
T(n,k) is the number of self-inverse permutations of {1,2,...,n} having exactly k cycles. - Geoffrey Critzer, May 08 2012
Bessel numbers of the second kind. For relations to the Hermite polynomials and the Catalan (A033184 and A009766) and Fibonacci (A011973, A098925, and A092865) matrices, see Yang and Qiao. - Tom Copeland, Dec 18 2013.
Also the inverse Bell transform of the double factorial of odd numbers Product_{k= 0..n-1} (2*k+1) (A001147). For the definition of the Bell transform see A264428 and for cross-references A265604. - Peter Luschny, Dec 31 2015

Examples

			Triangle begins:
    1
    0    1
    0    1    1
    0    0    3    1
    0    0    3    6    1
    0    0    0   15   10    1
    0    0    0   15   45   15    1
    0    0    0    0  105  105   21    1
    0    0    0    0  105  420  210   28    1
    0    0    0    0    0  945 1260  378   36    1
From _Gus Wiseman_, Jan 12 2021: (Start)
As noted above, a(n) is the number of set partitions of {1..n} into k singletons or pairs. This is also the number of set partitions of subsets of {1..n} into n - k pairs. In the first case, row n = 5 counts the following set partitions:
  {{1},{2,3},{4,5}}  {{1},{2},{3},{4,5}}  {{1},{2},{3},{4},{5}}
  {{1,2},{3},{4,5}}  {{1},{2},{3,4},{5}}
  {{1,2},{3,4},{5}}  {{1},{2,3},{4},{5}}
  {{1,2},{3,5},{4}}  {{1,2},{3},{4},{5}}
  {{1},{2,4},{3,5}}  {{1},{2},{3,5},{4}}
  {{1},{2,5},{3,4}}  {{1},{2,4},{3},{5}}
  {{1,3},{2},{4,5}}  {{1},{2,5},{3},{4}}
  {{1,3},{2,4},{5}}  {{1,3},{2},{4},{5}}
  {{1,3},{2,5},{4}}  {{1,4},{2},{3},{5}}
  {{1,4},{2},{3,5}}  {{1,5},{2},{3},{4}}
  {{1,4},{2,3},{5}}
  {{1,4},{2,5},{3}}
  {{1,5},{2},{3,4}}
  {{1,5},{2,3},{4}}
  {{1,5},{2,4},{3}}
In the second case, we have:
  {{1,2},{3,4}}  {{1,2}}  {}
  {{1,2},{3,5}}  {{1,3}}
  {{1,2},{4,5}}  {{1,4}}
  {{1,3},{2,4}}  {{1,5}}
  {{1,3},{2,5}}  {{2,3}}
  {{1,3},{4,5}}  {{2,4}}
  {{1,4},{2,3}}  {{2,5}}
  {{1,4},{2,5}}  {{3,4}}
  {{1,4},{3,5}}  {{3,5}}
  {{1,5},{2,3}}  {{4,5}}
  {{1,5},{2,4}}
  {{1,5},{3,4}}
  {{2,3},{4,5}}
  {{2,4},{3,5}}
  {{2,5},{3,4}}
(End)
		

Crossrefs

Row sums are A000085.
Column sums are A001515.
Same as A049403 but with a first column k = 0.
The same set partitions counted by number of pairs are A100861.
Reversing rows gives A111924 (without column k = 0).
A047884 counts standard Young tableaux by size and greatest row length.
A238123 counts standard Young tableaux by size and least row length.
A320663/A339888 count unlabeled multiset partitions into singletons/pairs.
A322661 counts labeled covering half-loop-graphs.
A339742 counts factorizations into distinct primes or squarefree semiprimes.

Programs

  • Maple
    # The function BellMatrix is defined in A264428.
    BellMatrix(n -> `if`(n<2,1,0), 9); # Peter Luschny, Jan 27 2016
  • Mathematica
    t[n_, k_] := k!*Binomial[n, k]/((2 k - n)!*2^(n - k)); Table[ t[n, k], {n, 0, 11}, {k, 0, n}] // Flatten
    (* Second program: *)
    rows = 12;
    t = Join[{1, 1}, Table[0, rows]];
    T[n_, k_] := BellY[n, k, t];
    Table[T[n, k], {n, 0, rows}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jun 23 2018,after Peter Luschny *)
    sbs[{}]:={{}};sbs[set:{i_,_}]:=Join@@Function[s,(Prepend[#1,s]&)/@sbs[Complement[set,s]]]/@Cases[Subsets[set],{i}|{i,_}];
    Table[Length[Select[sbs[Range[n]],Length[#]==k&]],{n,0,6},{k,0,n}] (* Gus Wiseman, Jan 12 2021 *)
  • PARI
    {T(n,k)=if(2*kn, 0, n!/(2*k-n)!/(n-k)!*2^(k-n))} /* Michael Somos, Oct 03 2006 */
    
  • Sage
    # uses[inverse_bell_transform from A265605]
    multifact_2_1 = lambda n: prod(2*k + 1 for k in (0..n-1))
    inverse_bell_matrix(multifact_2_1, 9) # Peter Luschny, Dec 31 2015

Formula

Number triangle T(n,k) = k!*C(n,k)/((2k-n)!*2^(n-k)).
T(n,k) = A001498(k,n-k). - Michael Somos, Oct 03 2006
E.g.f.: exp(y(x+x^2/2)). - Geoffrey Critzer, May 08 2012
Triangle equals the matrix product A008275*A039755. Equivalently, the n-th row polynomial R(n,x) is given by the Type B Dobinski formula R(n,x) = exp(-x/2)*Sum_{k>=0} P(n,2*k+1)*(x/2)^k/k!, where P(n,x) = x*(x-1)*...*(x-n+1) denotes the falling factorial polynomial. Cf. A113278. - Peter Bala, Jun 23 2014
From Daniel Checa, Aug 28 2022: (Start)
E.g.f. for the m-th column: (x^2/2+x)^m/m!.
T(n,k) = T(n-1,k-1) + (n-1)*T(n-2,k-1) for n>1 and k=1..n, T(0,0) = 1. (End)

A339888 Number of non-isomorphic multiset partitions of weight n into singletons or strict pairs.

Original entry on oeis.org

1, 1, 3, 5, 13, 23, 55, 104, 236, 470, 1039, 2140, 4712, 9962, 21961, 47484, 105464, 232324, 521338, 1167825, 2651453, 6031136, 13863054, 31987058, 74448415, 174109134, 410265423, 971839195, 2317827540, 5558092098, 13412360692, 32542049038, 79424450486
Offset: 0

Views

Author

Gus Wiseman, Jan 09 2021

Keywords

Examples

			Non-isomorphic representatives of the a(1) = 1 through a(4) = 13 multiset partitions:
  {{1}}  {{1,2}}    {{1},{2,3}}    {{1,2},{1,2}}
         {{1},{1}}  {{2},{1,2}}    {{1,2},{3,4}}
         {{1},{2}}  {{1},{1},{1}}  {{1,3},{2,3}}
                    {{1},{2},{2}}  {{1},{1},{2,3}}
                    {{1},{2},{3}}  {{1},{2},{1,2}}
                                   {{1},{2},{3,4}}
                                   {{1},{3},{2,3}}
                                   {{2},{2},{1,2}}
                                   {{1},{1},{1},{1}}
                                   {{1},{1},{2},{2}}
                                   {{1},{2},{2},{2}}
                                   {{1},{2},{3},{3}}
                                   {{1},{2},{3},{4}}
		

Crossrefs

The version for set partitions is A000085, with ordered version A080599.
The case of integer partitions is 1 + A004526(n), ranked by A003586.
Non-isomorphic multiset partitions are counted by A007716.
The case without singletons is A007717.
The version allowing non-strict pairs (x,x) is A320663.
A001190 counts rooted trees with out-degrees <= 2, ranked by A292050.
A339742 counts factorizations into distinct primes or squarefree semiprimes.
A339887 counts factorizations into primes or squarefree semiprimes.

Programs

  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
    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}
    gs(v) = {sum(i=2, #v, sum(j=1, i-1, my(g=gcd(v[i], v[j])); g*x^(2*v[i]*v[j]/g))) + sum(i=1, #v, my(r=v[i]); (1 + (1+r)%2)*x^r + ((r-1)\2)*x^(2*r))}
    a(n)={if(n==0, 1, my(s=0); forpart(p=n, s+=permcount(p)*EulerT(Vec(gs(p) + O(x*x^n), -n))[n]); s/n!)} \\ Andrew Howroyd, Apr 16 2021

Extensions

Terms a(11) and beyond from Andrew Howroyd, Apr 16 2021

A368598 Number of non-isomorphic n-element sets of singletons or pairs of elements of {1..n}, or unlabeled loop-graphs with n edges and up to n vertices.

Original entry on oeis.org

1, 1, 2, 6, 17, 52, 173, 585, 2064, 7520, 28265, 109501, 437394, 1799843, 7629463, 33302834, 149633151, 691702799, 3287804961, 16058229900, 80533510224, 414384339438, 2185878202630, 11811050484851, 65318772618624, 369428031895444, 2135166786135671, 12601624505404858
Offset: 0

Views

Author

Gus Wiseman, Jan 05 2024

Keywords

Comments

It doesn't matter for this sequence whether we use loops such as {x,x} or half-loops such as {x}.

Examples

			Non-isomorphic representatives of the a(0) = 1 through a(4) = 17 set-systems:
  {}  {{1}}  {{1},{2}}    {{1},{2},{3}}        {{1},{2},{3},{4}}
             {{1},{1,2}}  {{1},{2},{1,2}}      {{1},{2},{3},{1,2}}
                          {{1},{2},{1,3}}      {{1},{2},{3},{1,4}}
                          {{1},{1,2},{1,3}}    {{1},{2},{1,2},{1,3}}
                          {{1},{1,2},{2,3}}    {{1},{2},{1,2},{3,4}}
                          {{1,2},{1,3},{2,3}}  {{1},{2},{1,3},{1,4}}
                                               {{1},{2},{1,3},{2,3}}
                                               {{1},{2},{1,3},{2,4}}
                                               {{1},{3},{1,2},{2,4}}
                                               {{1},{1,2},{1,3},{1,4}}
                                               {{1},{1,2},{1,3},{2,3}}
                                               {{1},{1,2},{1,3},{2,4}}
                                               {{1},{1,2},{2,3},{3,4}}
                                               {{2},{1,2},{1,3},{1,4}}
                                               {{4},{1,2},{1,3},{2,3}}
                                               {{1,2},{1,3},{1,4},{2,3}}
                                               {{1,2},{1,3},{2,4},{3,4}}
		

Crossrefs

For any number of edges of any size we have A000612, covering A055621.
For any number of edges we have A000666, A054921, A322700.
The labeled version is A014068.
Counting by weight gives A320663, or A339888 with loops {x,x}.
The covering case is A368599.
For edges of any size we have A368731, covering A368186.
Row sums of A368836.
A000085 counts set partitions into singletons or pairs.
A001515 counts length-n set partitions into singletons or pairs.
A100861 counts set partitions into singletons or pairs by number of pairs.
A111924 counts set partitions into singletons or pairs by length.

Programs

  • Mathematica
    brute[m_]:=First[Sort[Table[Sort[Sort /@ (m/.Rule@@@Table[{(Union@@m)[[i]],p[[i]]}, {i,Length[p]}])], {p,Permutations[Range[Length[Union@@m]]]}]]];
    Table[Length[Union[brute /@ Subsets[Subsets[Range[n],{1,2}],{n}]]],{n,0,5}]
  • PARI
    a(n) = polcoef(G(n, O(x*x^n)), n) \\ G defined in A070166. - Andrew Howroyd, Jan 09 2024

Formula

a(n) = A070166(n, n). - Andrew Howroyd, Jan 09 2024

Extensions

Terms a(7) and beyond from Andrew Howroyd, Jan 09 2024

A339742 Number of factorizations of n into distinct primes or squarefree semiprimes.

Original entry on oeis.org

1, 1, 1, 0, 1, 2, 1, 0, 0, 2, 1, 1, 1, 2, 2, 0, 1, 1, 1, 1, 2, 2, 1, 0, 0, 2, 0, 1, 1, 4, 1, 0, 2, 2, 2, 1, 1, 2, 2, 0, 1, 4, 1, 1, 1, 2, 1, 0, 0, 1, 2, 1, 1, 0, 2, 0, 2, 2, 1, 3, 1, 2, 1, 0, 2, 4, 1, 1, 2, 4, 1, 0, 1, 2, 1, 1, 2, 4, 1, 0, 0, 2, 1, 3, 2, 2, 2, 0, 1, 3, 2, 1, 2, 2, 2, 0, 1, 1, 1, 1, 1, 4, 1, 0, 4
Offset: 1

Views

Author

Gus Wiseman, Dec 20 2020

Keywords

Comments

A squarefree semiprime (A006881) is a product of any two distinct prime numbers.
The following are equivalent characteristics for any positive integer n:
(1) the prime factors of n can be partitioned into distinct singletons or strict pairs, i.e., into a set of half-loops and edges;
(2) n can be factored into distinct primes or squarefree semiprimes.

Examples

			The a(n) factorizations for n = 6, 30, 60, 210, 420 are respectively 2, 4, 3, 10, 9:
  (6)    (5*6)    (6*10)    (6*35)     (2*6*35)
  (2*3)  (2*15)   (2*5*6)   (10*21)    (5*6*14)
         (3*10)   (2*3*10)  (14*15)    (6*7*10)
         (2*3*5)            (5*6*7)    (2*10*21)
                            (2*3*35)   (2*14*15)
                            (2*5*21)   (2*5*6*7)
                            (2*7*15)   (3*10*14)
                            (3*5*14)   (2*3*5*14)
                            (3*7*10)   (2*3*7*10)
                            (2*3*5*7)
		

Crossrefs

Dirichlet convolution of A008966 with A339661.
A008966 allows only primes.
A339661 does not allow primes, only squarefree semiprimes.
A339740 lists the positions of zeros.
A339741 lists the positions of positive terms.
A339839 allows nonsquarefree semiprimes.
A339887 is the non-strict version.
A001358 lists semiprimes, with squarefree case A006881.
A002100 counts partitions into squarefree semiprimes.
A013929 cannot be factored into distinct primes.
A293511 are a product of distinct squarefree numbers in exactly one way.
A320663 counts non-isomorphic multiset partitions into singletons or pairs.
A339840 cannot be factored into distinct primes or semiprimes.
A339841 have exactly one factorization into primes or semiprimes.
The following count factorizations:
- A001055 into all positive integers > 1.
- A050320 into squarefree numbers.
- A050326 into distinct squarefree numbers.
- A320655 into semiprimes.
- A320656 into squarefree semiprimes.
- A320732 into primes or semiprimes.
- A322353 into distinct semiprimes.
- A339742 [this sequence] into distinct primes or squarefree semiprimes.
- A339839 into distinct primes or semiprimes.
The following count vertex-degree partitions and give their Heinz numbers:
- A000569 counts graphical partitions (A320922).
- A058696 counts all partitions of 2n (A300061).
- A209816 counts multigraphical partitions (A320924).
- A339656 counts loop-graphical partitions (A339658).
-
The following count partitions/factorizations of even length and give their Heinz numbers:
- A027187/A339846 has no additional conditions (A028260).
- A338914/A339562 can be partitioned into edges (A320911).
- A338916/A339563 can be partitioned into distinct pairs (A320912).
- A339559/A339564 cannot be partitioned into distinct edges (A320894).
- A339560/A339619 can be partitioned into distinct edges (A339561).

Programs

  • Mathematica
    sqps[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sqps[n/d],Min@@#>d&]],{d,Select[Divisors[n],PrimeQ[#]||SquareFreeQ[#]&&PrimeOmega[#]==2&]}]];
    Table[Length[sqps[n]],{n,100}]
  • PARI
    A353471(n) = (numdiv(n)==2*omega(n));
    A339742(n, u=(1+n)) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1) && (dA353471(d), s += A339742(n/d, d))); (s)); \\ Antti Karttunen, May 02 2022

Formula

a(n) = Sum_{d|n squarefree} A339661(n/d).

Extensions

More terms from Antti Karttunen, May 02 2022

A368599 Number of non-isomorphic n-element sets of singletons or pairs of elements of {1..n} with union {1..n}, or unlabeled loop-graphs with n edges covering n vertices.

Original entry on oeis.org

1, 1, 2, 5, 13, 34, 97, 277, 825, 2486, 7643, 23772, 74989, 238933, 769488, 2500758, 8199828, 27106647, 90316944, 303182461, 1025139840, 3490606305, 11967066094, 41302863014, 143493606215, 501772078429, 1765928732426, 6254738346969, 22294413256484, 79968425399831
Offset: 0

Views

Author

Gus Wiseman, Jan 06 2024

Keywords

Comments

It doesn't matter for this sequence whether we use loops such as {x,x} or half-loops such as {x}.

Examples

			The a(0) = 1 through a(4) = 13 set-systems:
  {}  {{1}}  {{1},{2}}    {{1},{2},{3}}        {{1},{2},{3},{4}}
             {{1},{1,2}}  {{1},{2},{1,3}}      {{1},{2},{3},{1,4}}
                          {{1},{1,2},{1,3}}    {{1},{2},{1,2},{3,4}}
                          {{1},{1,2},{2,3}}    {{1},{2},{1,3},{1,4}}
                          {{1,2},{1,3},{2,3}}  {{1},{2},{1,3},{2,4}}
                                               {{1},{2},{1,3},{3,4}}
                                               {{1},{1,2},{1,3},{1,4}}
                                               {{1},{1,2},{1,3},{2,4}}
                                               {{1},{1,2},{2,3},{2,4}}
                                               {{1},{1,2},{2,3},{3,4}}
                                               {{1},{2,3},{2,4},{3,4}}
                                               {{1,2},{1,3},{1,4},{2,3}}
                                               {{1,2},{1,3},{2,4},{3,4}}
		

Crossrefs

For any number of edges we have A000666, A054921, A322700.
For any number of edges of any size we have A055621, non-covering A000612.
For edges of any size we have A368186, covering case of A368731.
The labeled version is A368597, covering case of A014068.
This is the covering case of A368598.
A000085 counts set partitions into singletons or pairs.
A001515 counts length-n set partitions into singletons or pairs.
A100861 counts set partitions into singletons or pairs by number of pairs.
A111924 counts set partitions into singletons or pairs by length.

Programs

  • Mathematica
    brute[m_]:=First[Sort[Table[Sort[Sort /@ (m/.Rule@@@Table[{(Union@@m)[[i]],p[[i]]}, {i,Length[p]}])], {p,Permutations[Range[Length[Union@@m]]]}]]];
    Table[Length[Union[brute /@ Select[Subsets[Subsets[Range[n],{1,2}],{n}], Union@@#==Range[n]&]]],{n,0,5}]
  • PARI
    a(n) = polcoef(G(n, O(x*x^n)) - if(n, G(n-1, O(x*x^n))), n) \\ G defined in A070166. - Andrew Howroyd, Jan 09 2024

Formula

a(n) = A070166(n,n) - A070166(n-1,n) for n > 0. - Andrew Howroyd, Jan 09 2024

Extensions

Terms a(7) and beyond from Andrew Howroyd, Jan 09 2024

A323437 Number of semistandard Young tableaux whose entries are the prime indices of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 1, 4, 1, 1, 2, 2, 2, 3, 1, 2, 2, 2, 1, 4, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 5, 1, 2, 2, 1, 2, 4, 1, 2, 2, 4, 1, 3, 1, 2, 2, 2, 2, 4, 1, 2, 1, 2, 1, 5, 2, 2, 2
Offset: 0

Views

Author

Gus Wiseman, Jan 15 2019

Keywords

Comments

Number of ways to fill a Young diagram with the prime indices of n such that all rows are weakly increasing and all columns are strictly increasing.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
Is this a duplicate of A339887? - R. J. Mathar, Feb 03 2021

Examples

			The a(60) = 5 tableaux:
  1123
.
  11   112   113
  23   3     2
.
  11
  2
  3
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    ptnplane[n_]:=Union[Map[primeMS,Join@@Permutations/@facs[n],{2}]];
    Table[Length[Select[ptnplane[y],And[And@@Less@@@#,And@@(LessEqual@@@Transpose[PadRight[#]/.(0->Infinity)])]&]],{y,100}]

Formula

Sum_{A056239(n) = k} a(k) = A003293(n).
Showing 1-6 of 6 results.