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.

A007837 Number of partitions of n-set with distinct block sizes.

Original entry on oeis.org

1, 1, 1, 4, 5, 16, 82, 169, 541, 2272, 17966, 44419, 201830, 802751, 4897453, 52275409, 166257661, 840363296, 4321172134, 24358246735, 183351656650, 2762567051857, 10112898715063, 62269802986835, 343651382271526, 2352104168848091, 15649414071734847
Offset: 0

Views

Author

Keywords

Comments

Conjecture: the Gauss congruences a(n*p^k) == a(n*p^(k-1)) (mod p^k) hold for all primes p and positive integers n and k. Cf. A185895. - Peter Bala, Mar 17 2022

Examples

			From _Gus Wiseman_, Jul 13 2019: (Start)
The a(1) = 1 through a(5) = 16 set partitions with distinct block sizes:
  {{1}}  {{1,2}}  {{1,2,3}}    {{1,2,3,4}}    {{1,2,3,4,5}}
                  {{1},{2,3}}  {{1},{2,3,4}}  {{1},{2,3,4,5}}
                  {{1,2},{3}}  {{1,2,3},{4}}  {{1,2},{3,4,5}}
                  {{1,3},{2}}  {{1,2,4},{3}}  {{1,2,3},{4,5}}
                               {{1,3,4},{2}}  {{1,2,3,4},{5}}
                                              {{1,2,3,5},{4}}
                                              {{1,2,4},{3,5}}
                                              {{1,2,4,5},{3}}
                                              {{1,2,5},{3,4}}
                                              {{1,3},{2,4,5}}
                                              {{1,3,4},{2,5}}
                                              {{1,3,4,5},{2}}
                                              {{1,3,5},{2,4}}
                                              {{1,4},{2,3,5}}
                                              {{1,4,5},{2,3}}
                                              {{1,5},{2,3,4}}
(End)
		

Crossrefs

Programs

  • Maple
    a:= proc(n) option remember; `if`(n=0, 1, add(add((-d)*(-d!)^(-k/d),
           d=numtheory[divisors](k))*(n-1)!/(n-k)!*a(n-k), k=1..n))
        end:
    seq(a(n), n=0..30);  # Alois P. Heinz, Sep 06 2008
    # second Maple program:
    A007837 := proc(n) option remember; local k; `if`(n = 0, 1,
    add(binomial(n-1, k-1) * A182927(k) * A007837(n-k), k = 1..n)) end:
    seq(A007837(i),i=0..24); # Peter Luschny, Apr 25 2011
  • Mathematica
    nn=20;p=Product[1+x^i/i!,{i,1,nn}];Drop[Range[0,nn]!CoefficientList[ Series[p,{x,0,nn}],x],1]  (* Geoffrey Critzer, Sep 22 2012 *)
    a[0]=1; a[n_] := a[n] = Sum[(n-1)!/(n-k)!*DivisorSum[k, -#*(-#!)^(-k/#)&]* a[n-k], {k, 1, n}]; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Nov 23 2015, after Vladeta Jovovic *)
  • PARI
    {my(n=20); Vec(serlaplace(prod(k=1, n, (1+x^k/k!) + O(x*x^n))))} \\ Andrew Howroyd, Dec 21 2017

Formula

E.g.f.: Product_{m >= 1} (1+x^m/m!).
a(n) = Sum_{k=1..n} (n-1)!/(n-k)!*b(k)*a(n-k), where b(k) = Sum_{d divides k} (-d)*(-d!)^(-k/d) and a(0) = 1. - Vladeta Jovovic, Oct 13 2002
E.g.f.: exp(Sum_{k>=1} Sum_{j>=1} (-1)^(k+1)*x^(j*k)/(k*(j!)^k)). - Ilya Gutkovskiy, Jun 18 2018

Extensions

More terms from Christian G. Bower
a(0)=1 prepended by Alois P. Heinz, Aug 29 2015

A326517 Number of normal multiset partitions of weight n where each part has a different size.

Original entry on oeis.org

1, 1, 2, 12, 28, 140, 956, 3520, 17792, 111600, 1144400, 4884064, 34907936, 214869920, 1881044032, 25687617152, 139175009920, 1098825972608, 8770328141888, 74286112885504, 784394159958848, 15114871659653952, 92392468773724544, 889380453354852416, 7652770202041529856
Offset: 0

Views

Author

Gus Wiseman, Jul 12 2019

Keywords

Comments

A multiset partition is normal if it covers an initial interval of positive integers.

Examples

			The a(0) = 1 through a(3) = 12 normal multiset partitions:
  {}  {{1}}  {{1,1}}  {{1,1,1}}
             {{1,2}}  {{1,1,2}}
                      {{1,2,2}}
                      {{1,2,3}}
                      {{1},{1,1}}
                      {{1},{1,2}}
                      {{1},{2,2}}
                      {{1},{2,3}}
                      {{2},{1,1}}
                      {{2},{1,2}}
                      {{2},{1,3}}
                      {{3},{1,2}}
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i, k) option remember; `if`(n=0, 1, `if`(i<1, 0,
          add(b(n-i*j, i-1, k)*binomial(i+k-1, k-1)^j, j=0..min(1, n/i))))
        end:
    a:= n->add(add(b(n$2, k-i)*(-1)^i*binomial(k, i), i=0..n), k=0..n):
    seq(a(n), n=0..25);  # Alois P. Heinz, Sep 23 2023
  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    allnorm[n_]:=If[n<=0,{{}},Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1]];
    Table[Length[Select[Join@@mps/@allnorm[n],UnsameQ@@Length/@#&]],{n,0,6}]
  • PARI
    R(n, k)={Vec(prod(j=1, n, 1 + binomial(k+j-1, j)*x^j + O(x*x^n)))}
    seq(n)={sum(k=0, n, R(n, k)*sum(r=k, n, binomial(r, k)*(-1)^(r-k)))} \\ Andrew Howroyd, Feb 07 2020

Extensions

Terms a(8) and beyond from Andrew Howroyd, Feb 07 2020

A326514 Number of factorizations of n into factors > 1 where each factor has a different number of prime factors counted with multiplicity.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jul 11 2019

Keywords

Examples

			The a(96) = 8 factorizations: (2*4*12), (2*6*8), (2*48), (3*4*8), (3*32), (4*24), (6*16), (96).
		

Crossrefs

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[n],UnsameQ@@PrimeOmega/@#&]],{n,100}]

A326533 MM-numbers of multiset partitions where each part has a different length.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 13, 14, 17, 19, 21, 22, 23, 26, 29, 31, 34, 35, 37, 38, 39, 41, 42, 43, 46, 47, 53, 57, 58, 59, 61, 62, 65, 67, 69, 70, 71, 73, 74, 77, 78, 79, 82, 83, 86, 87, 89, 94, 95, 97, 101, 103, 106, 107, 109, 111, 113, 114, 115, 118, 119, 122
Offset: 1

Views

Author

Gus Wiseman, Jul 12 2019

Keywords

Comments

These are numbers where each prime index has a different Omega (number of prime factors counted with multiplicity). 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. The multiset multisystem with MM-number n is obtained by taking the multiset of prime indices of each prime index of n. For example, the prime indices of 78 are {1,2,6}, so the multiset multisystem with MM-number 78 is {{},{1},{1,2}}.

Examples

			The sequence of multiset partitions where each part has a different average preceded by their MM-numbers begins:
   1: {}
   2: {{}}
   3: {{1}}
   5: {{2}}
   6: {{},{1}}
   7: {{1,1}}
  10: {{},{2}}
  11: {{3}}
  13: {{1,2}}
  14: {{},{1,1}}
  17: {{4}}
  19: {{1,1,1}}
  21: {{1},{1,1}}
  22: {{},{3}}
  23: {{2,2}}
  26: {{},{1,2}}
  29: {{1,3}}
  31: {{5}}
  34: {{},{4}}
  35: {{2},{1,1}}
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],UnsameQ@@PrimeOmega/@primeMS[#]&]

A326570 Number of covering antichains of subsets of {1..n} with different edge-sizes.

Original entry on oeis.org

2, 1, 1, 4, 17, 186, 3292, 139161, 14224121
Offset: 0

Views

Author

Gus Wiseman, Jul 18 2019

Keywords

Comments

An antichain is a finite set of finite sets, none of which is a subset of any other. It is covering if its union is {1..n}. The edge-sizes are the numbers of vertices in each edge, so for example the edge-sizes of {{1,3},{2,5},{3,4,5}} are {2,2,3}.

Examples

			The a(0) = 2 through a(4) = 17 antichains:
  {}    {{1}}  {{1,2}}  {{1,2,3}}    {{1,2,3,4}}
  {{}}                  {{1},{2,3}}  {{1},{2,3,4}}
                        {{2},{1,3}}  {{2},{1,3,4}}
                        {{3},{1,2}}  {{3},{1,2,4}}
                                     {{4},{1,2,3}}
                                     {{1,2},{1,3,4}}
                                     {{1,2},{2,3,4}}
                                     {{1,3},{1,2,4}}
                                     {{1,3},{2,3,4}}
                                     {{1,4},{1,2,3}}
                                     {{1,4},{2,3,4}}
                                     {{2,3},{1,2,4}}
                                     {{2,3},{1,3,4}}
                                     {{2,4},{1,2,3}}
                                     {{2,4},{1,3,4}}
                                     {{3,4},{1,2,3}}
                                     {{3,4},{1,2,4}}
		

Crossrefs

Antichain covers are A006126.
Set partitions with different block sizes are A007837.
The case without singletons is A326569.
(Antichain) covers with equal edge-sizes are A306021.

Programs

  • Mathematica
    stableSets[u_,Q_]:=If[Length[u]==0,{{}},With[{w=First[u]},Join[stableSets[DeleteCases[u,w],Q],Prepend[#,w]&/@stableSets[DeleteCases[u,r_/;r==w||Q[r,w]||Q[w,r]],Q]]]];
    cleq[n_]:=Select[stableSets[Subsets[Range[n]],SubsetQ[#1,#2]||Length[#1]==Length[#2]&],Union@@#==Range[n]&];
    Table[Length[cleq[n]],{n,0,6}]

Extensions

a(8) from Andrew Howroyd, Aug 13 2019

A332260 Triangle read by rows: T(n,k) is the number of non-isomorphic multiset partitions of weight n whose union is a k-set where each part has a different size.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 2, 3, 2, 0, 2, 5, 3, 2, 0, 3, 11, 12, 6, 3, 0, 4, 26, 39, 27, 11, 4, 0, 5, 40, 79, 67, 37, 14, 5, 0, 6, 68, 170, 184, 116, 55, 19, 6, 0, 8, 122, 407, 543, 417, 219, 91, 28, 8, 0, 10, 232, 1082, 1911, 1760, 1052, 459, 159, 42, 10
Offset: 0

Views

Author

Andrew Howroyd, Feb 08 2020

Keywords

Comments

T(n,k) is the number of nonequivalent nonnegative integer matrices with total sum n and k nonzero rows with distinct column sums up to permutation of rows and columns.

Examples

			Triangle begins:
  1;
  0, 1;
  0, 1,   1;
  0, 2,   3,   2;
  0, 2,   5,   3,   2;
  0, 3,  11,  12,   6,   3;
  0, 4,  26,  39,  27,  11,   4;
  0, 5,  40,  79,  67,  37,  14,  5;
  0, 6,  68, 170, 184, 116,  55, 19,  6;
  0, 8, 122, 407, 543, 417, 219, 91, 28, 8;
  ...
The T(4,2) = 5 multiset partitions are:
  {{1,1,2,2}}, {{1,2,2,2}}, {{1},{1,2,2}}, {{1},{2,2,2}}, {{1},{1,1,2}}.
These correspond with the following matrices:
   [2]  [1]  [1 1]  [1 0]  [1 2]
   [2]  [3]  [0 2]  [0 3]  [0 1]
		

Crossrefs

Column k=1 is A000009.
Main diagonal is A000009.
Row sums are A326026.
Cf. A332253.

Programs

  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
    D(p,n)={my(v=vector(n)); for(i=1, #p, v[p[i]]++); my(u=EulerT(v)); prod(j=1, #u, 1 + u[j]*x^j + O(x*x^n))/if(!#p, 1, prod(i=1, p[#p], i^v[i]*v[i]!))}
    M(n)={my(v=vector(n+1)); for(i=0, n, my(s=0); forpart(p=i, s+=D(p,n)); v[1+i]=Col(s)); Mat(vector(#v, i, v[i]-if(i>1, v[i-1])))}
    {my(T=M(10)); for(n=1, #T~, print(T[n, ][1..n]))}

A336138 Number of set partitions of the binary indices of n with distinct block-sums.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 2, 4, 1, 2, 2, 5, 2, 4, 5, 12, 1, 2, 2, 5, 2, 5, 4, 13, 2, 4, 5, 13, 5, 13, 13, 43, 1, 2, 2, 5, 2, 5, 5, 13, 2, 5, 4, 14, 5, 13, 14, 42, 2, 4, 5, 13, 5, 14, 13, 43, 5, 13, 14, 45, 14, 44, 44, 160, 1, 2, 2, 5, 2, 5, 5, 14, 2, 5, 5, 14, 4, 13
Offset: 0

Views

Author

Gus Wiseman, Jul 12 2020

Keywords

Comments

A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.

Examples

			The a(n) set partitions for n = 3, 7, 11, 15, 23:
  {12}    {123}      {124}      {1234}        {1235}
  {1}{2}  {1}{23}    {1}{24}    {1}{234}      {1}{235}
          {13}{2}    {12}{4}    {12}{34}      {12}{35}
          {1}{2}{3}  {14}{2}    {123}{4}      {123}{5}
                     {1}{2}{4}  {124}{3}      {125}{3}
                                {13}{24}      {13}{25}
                                {134}{2}      {135}{2}
                                {1}{2}{34}    {15}{23}
                                {1}{23}{4}    {1}{2}{35}
                                {1}{24}{3}    {1}{25}{3}
                                {14}{2}{3}    {13}{2}{5}
                                {1}{2}{3}{4}  {15}{2}{3}
                                              {1}{2}{3}{5}
		

Crossrefs

The version for twice-partitions is A271619.
The version for partitions of partitions is (also) A271619.
These set partitions are counted by A275780.
The version for factorizations is A321469.
The version for normal multiset partitions is A326519.
The version for equal block-sums is A336137.
Set partitions with distinct block-lengths are A007837.
Set partitions of binary indices are A050315.
Twice-partitions with equal sums are A279787.
Partitions of partitions with equal sums are A305551.
Normal multiset partitions with equal block-lengths are A317583.
Multiset partitions with distinct block-sums are ranked by A326535.

Programs

  • Mathematica
    bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    Table[Length[Select[sps[bpe[n]],UnsameQ@@Total/@#&]],{n,0,100}]
Showing 1-7 of 7 results.