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.

Previous Showing 51-60 of 102 results. Next

A330462 Triangle read by rows where T(n,k) is the number of k-element sets of nonempty sets of positive integers with total sum n.

Original entry on oeis.org

1, 0, 1, 0, 1, 0, 0, 2, 1, 0, 0, 2, 2, 0, 0, 0, 3, 4, 0, 0, 0, 0, 4, 6, 2, 0, 0, 0, 0, 5, 11, 3, 0, 0, 0, 0, 0, 6, 16, 8, 0, 0, 0, 0, 0, 0, 8, 25, 15, 1, 0, 0, 0, 0, 0, 0, 10, 35, 28, 4, 0, 0, 0, 0, 0, 0, 0, 12, 52, 46, 9, 0, 0, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Gus Wiseman, Dec 18 2019

Keywords

Examples

			Triangle begins:
  1
  0  1
  0  1  0
  0  2  1  0
  0  2  2  0  0
  0  3  4  0  0  0
  0  4  6  2  0  0  0
  0  5 11  3  0  0  0  0
  0  6 16  8  0  0  0  0  0
  0  8 25 15  1  0  0  0  0  0
  0 10 35 28  4  0  0  0  0  0  0
  ...
Row n = 7 counts the following set-systems:
  {{7}}      {{1},{6}}      {{1},{2},{4}}
  {{1,6}}    {{2},{5}}      {{1},{2},{1,3}}
  {{2,5}}    {{3},{4}}      {{1},{3},{1,2}}
  {{3,4}}    {{1},{1,5}}
  {{1,2,4}}  {{1},{2,4}}
             {{2},{1,4}}
             {{2},{2,3}}
             {{3},{1,3}}
             {{4},{1,2}}
             {{1},{1,2,3}}
             {{1,2},{1,3}}
		

Crossrefs

Programs

  • Mathematica
    ppl[n_,k_]:=Switch[k,0,{n},1,IntegerPartitions[n],_,Join@@Table[Union[Sort/@Tuples[ppl[#,k-1]&/@ptn]],{ptn,IntegerPartitions[n]}]];
    Table[Length[Select[ppl[n,2],And[UnsameQ@@#,And@@UnsameQ@@@#,Length[#]==k]&]],{n,0,10},{k,0,n}]
  • PARI
    L(n)={eta(x^2 + O(x*x^n))/eta(x + O(x*x^n))}
    A(n)={my(c=L(n), v=Vec(prod(k=1, n, (1 + x^k*y + O(x*x^n))^polcoef(c,k)))); vector(#v, n, Vecrev(v[n],n))}
    {my(T=A(12)); for(n=1, #T, print(T[n]))} \\ Andrew Howroyd, Dec 29 2019

Formula

G.f.: Product_{j>=1} (1 + y*x^j)^A000009(j). - Andrew Howroyd, Dec 29 2019

A356939 MM-numbers of multisets of intervals. Products of primes indexed by members of A073485.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 8, 9, 10, 11, 12, 13, 15, 16, 17, 18, 20, 22, 24, 25, 26, 27, 30, 31, 32, 33, 34, 36, 39, 40, 41, 44, 45, 47, 48, 50, 51, 52, 54, 55, 59, 60, 62, 64, 65, 66, 67, 68, 72, 75, 78, 80, 81, 82, 83, 85, 88, 90, 93, 94, 96, 99, 100, 102, 104, 108
Offset: 1

Views

Author

Gus Wiseman, Sep 12 2022

Keywords

Comments

An interval such as {3,4,5} is a set of positive integers with all differences of adjacent elements equal to 1.
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.
We define the multiset of multisets with MM-number n to be formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. The size of this multiset of multisets is A302242(n). For example, the prime indices of 78 are {1,2,6}, so the multiset of multisets with MM-number 78 is {{},{1},{1,2}}.

Examples

			The initial terms and corresponding multisets of multisets:
   1: {}
   2: {{}}
   3: {{1}}
   4: {{},{}}
   5: {{2}}
   6: {{},{1}}
   8: {{},{},{}}
   9: {{1},{1}}
  10: {{},{2}}
  11: {{3}}
  12: {{},{},{1}}
  13: {{1,2}}
  15: {{1},{2}}
  16: {{},{},{},{}}
		

Crossrefs

The initial version is A356940.
Intervals are counted by A000012, A001227, ranked by A073485.
Other types: A107742, A356936, A356937, A356938.
Other conditions: A302478, A302492, A356930, A356935, A356944, A356955.
A000041 counts integer partitions, strict A000009.
A000688 counts factorizations into prime powers.
A001055 counts factorizations.
A001221 counts prime divisors, sum A001414.
A001222 counts prime factors with multiplicity.
A056239 adds up prime indices, row sums of A112798.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    chQ[y_]:=Or[Length[y]<=1,Union[Differences[y]]=={1}];
    Select[Range[100],And@@chQ/@primeMS/@primeMS[#]&]

A381719 Numbers whose prime indices cannot be partitioned into sets with a common sum.

Original entry on oeis.org

12, 18, 20, 24, 28, 40, 44, 45, 48, 50, 52, 54, 56, 60, 63, 68, 72, 75, 76, 80, 84, 88, 90, 92, 96, 98, 99, 104, 108, 112, 116, 117, 120, 124, 126, 132, 135, 136, 140, 144, 147, 148, 152, 153, 156, 160, 162, 164, 168, 171, 172, 175, 176, 184, 188, 189, 192
Offset: 1

Views

Author

Gus Wiseman, Apr 22 2025

Keywords

Comments

Differs from A059404, A323055, A376250 in lacking 150.
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, sum A056239.
Also numbers that cannot be factored into squarefree numbers with a common sum of prime indices (A056239).

Examples

			The prime indices of 150 are {1,2,3,3}, and {{3},{3},{1,2}} is a partition into sets with a common sum, so 150 is not in the sequence.
		

Crossrefs

Twice-partitions of this type (sets with a common sum) are counted by A279788.
These multiset partitions (sets with a common sum) are ranked by A326534 /\ A302478.
For distinct block-sums we have A381806, counted by A381990 (complement A381992).
For constant blocks we have A381871 (zeros of A381995), counted by A381993.
Partitions of this type are counted by A381994.
These are the zeros of A382080.
Normal multiset partitions of this type are counted by A382429, see A326518.
The complement counted by A383308.
A000041 counts integer partitions, strict A000009.
A001055 counts factorizations, strict A045778.
A050320 counts factorizations into squarefree numbers, see A381078, A381454.
A050326 counts factorizations into distinct squarefree numbers.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A317141 counts coarsenings of prime indices, refinements A300383.
A381633 counts set systems with distinct sums, see A381634, A293243.
Set multipartitions: A089259, A116540, A270995, A296119, A318360.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    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]]]];
    Select[Range[100],Select[mps[prix[#]], SameQ@@Total/@#&&And@@UnsameQ@@@#&]=={}&]

A382428 Number of normal multiset partitions of weight n into sets with distinct sizes.

Original entry on oeis.org

1, 1, 1, 6, 8, 35, 292, 673, 2818, 16956, 219772, 636748, 3768505, 20309534, 183403268, 3227600747, 12272598308, 81353466578, 561187259734, 4416808925866, 50303004612136, 1238783066956740, 5566249468690291, 44970939483601100, 330144217684933896, 3131452652308459402
Offset: 0

Views

Author

Gus Wiseman, Mar 29 2025

Keywords

Comments

We call a multiset or multiset partition normal iff it covers an initial interval of positive integers. The weight of a multiset partition is the sum of sizes of its blocks.

Examples

			The a(1) = 1 through a(4) = 8 multiset partitions:
  {{1}}  {{1,2}}  {{1,2,3}}    {{1,2,3,4}}
                  {{1},{1,2}}  {{1},{1,2,3}}
                  {{1},{2,3}}  {{1},{2,3,4}}
                  {{2},{1,2}}  {{2},{1,2,3}}
                  {{2},{1,3}}  {{2},{1,3,4}}
                  {{3},{1,2}}  {{3},{1,2,3}}
                               {{3},{1,2,4}}
                               {{4},{1,2,3}}
		

Crossrefs

For distinct sums instead of sizes we have A116539, see A050326.
Without distinct lengths we have A116540 (normal set multipartitions).
Without strict blocks we have A326517, for sum instead of size A326519.
For equal instead of distinct sizes we have A331638.
Twice-partitions of this type are counted by A358830.
For distinct sums instead of sizes we have A381718.
For equal instead of distinct sizes we have A382429.
A000670 counts patterns, ranked by A055932 and A333217, necklace A019536.
A001055 count factorizations, strict A045778.
Normal multiset partitions: A034691, A035310, A255906.
Set multipartitions: A089259, A270995, A296119, A318360.

Programs

  • Mathematica
    allnorm[n_Integer]:=Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1];
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[mset_]:=Union[Sort[Sort/@(#/.x_Integer:>mset[[x]])]&/@sps[Range[Length[mset]]]];
    Table[Length[Join@@(Select[mps[#],UnsameQ@@Length/@#&&And@@UnsameQ@@@#&]&/@allnorm[n])],{n,0,5}]
  • PARI
    R(n, k)={Vec(prod(j=1, n, 1 + binomial(k, 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, Mar 31 2025

Extensions

a(10) onwards from Andrew Howroyd, Mar 31 2025

A330460 Triangle read by rows where T(n,k) is the number of set partitions with k blocks and total sum n.

Original entry on oeis.org

1, 0, 1, 0, 1, 0, 0, 2, 1, 0, 0, 2, 1, 0, 0, 0, 3, 2, 0, 0, 0, 0, 4, 5, 1, 0, 0, 0, 0, 5, 6, 1, 0, 0, 0, 0, 0, 6, 9, 2, 0, 0, 0, 0, 0, 0, 8, 13, 3, 0, 0, 0, 0, 0, 0, 0, 10, 23, 10, 1, 0, 0, 0, 0, 0, 0, 0, 12, 27, 11, 1, 0, 0, 0, 0, 0, 0, 0, 0, 15, 40, 19, 2, 0, 0, 0, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Gus Wiseman, Dec 18 2019

Keywords

Examples

			Triangle begins:
  1
  0  1
  0  1  0
  0  2  1  0
  0  2  1  0  0
  0  3  2  0  0  0
  0  4  5  1  0  0  0
  0  5  6  1  0  0  0  0
  0  6  9  2  0  0  0  0  0
  0  8 13  3  0  0  0  0  0  0
  0 10 23 10  1  0  0  0  0  0  0
  0 12 27 11  1  0  0  0  0  0  0  0
  0 15 40 19  2  0  0  0  0  0  0  0  0
Row n = 8 counts the following set partitions:
  {{8}}      {{1},{7}}    {{1},{2},{5}}
  {{3,5}}    {{2},{6}}    {{1},{3},{4}}
  {{2,6}}    {{3},{5}}
  {{1,7}}    {{1},{3,4}}
  {{1,3,4}}  {{1},{2,5}}
  {{1,2,5}}  {{2},{1,5}}
             {{3},{1,4}}
             {{4},{1,3}}
             {{5},{1,2}}
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i, k) option remember; `if`(i*(i+1)/2 k*
             b(n-i, t, k)+b(n-i, t, k+1))(min(n-i, i-1))))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=0..n))(b(n$2, 0)):
    seq(T(n), n=0..15);  # Alois P. Heinz, Dec 29 2019
  • Mathematica
    ppl[n_,k_]:=Switch[k,0,{n},1,IntegerPartitions[n],_,Join@@Table[Union[Sort/@Tuples[ppl[#,k-1]&/@ptn]],{ptn,IntegerPartitions[n]}]];
    Table[Length[Select[ppl[n,2],Length[#]==k&&And[UnsameQ@@#,UnsameQ@@Join@@#]&]],{n,0,10},{k,0,n}]
    (* Second program: *)
    b[n_, i_, k_] := b[n, i, k] = If[i(i+1)/2 < n, 0, If[n == 0, x^k, b[n, i-1, k] + Function[t, k*b[n-i, t, k] + b[n-i, t, k + 1]][Min[n-i, i-1]]]];
    T[n_] := PadRight[CoefficientList[b[n, n, 0], x], n + 1];
    T /@ Range[0, 15] // Flatten (* Jean-François Alcover, May 16 2021, after Alois P. Heinz *)
  • PARI
    A(n)={my(v=Vec(prod(k=1, n, 1 + x^k*y + O(x*x^n)))); vector(#v, n, my(p=v[n]); vector(n, k, sum(i=k, n, polcoef(p,i-1)*stirling(i-1, k-1, 2))))}
    {my(T=A(12)); for(n=1, #T, print(T[n]))} \\ Andrew Howroyd, Dec 29 2019

Formula

T(n,k) = Sum_{k <= i <= n} A060016(n,i) * A008277(i,k).
For n > 0, T(n,2) = Sum_{k = 1..n} (2^(k - 1) -1) * A060016(n,k).

A356944 MM-numbers of multisets of gapless multisets of positive integers. Products of primes indexed by elements of A073491.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70
Offset: 1

Views

Author

Gus Wiseman, Sep 12 2022

Keywords

Comments

A multiset is gapless if it covers an interval of positive integers. For example, {2,3,3,4} is gapless but {1,1,3,3} is not.
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.
We define the multiset of multisets with MM-number n to be formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. The size of this multiset of multisets is A302242(n). For example, the prime indices of 78 are {1,2,6}, so the multiset of multisets with MM-number 78 is {{},{1},{1,2}}.

Examples

			The initial terms and corresponding multiset partitions:
   1: {}
   2: {{}}
   3: {{1}}
   4: {{},{}}
   5: {{2}}
   6: {{},{1}}
   7: {{1,1}}
   8: {{},{},{}}
   9: {{1},{1}}
  10: {{},{2}}
  11: {{3}}
  12: {{},{},{1}}
  13: {{1,2}}
  14: {{},{1,1}}
  15: {{1},{2}}
  16: {{},{},{},{}}
		

Crossrefs

Gapless multisets are counted by A034296, ranked by A073491.
The initial version is A356955.
Other types: A356233, A356941, A356942, A356943.
Other conditions: A302478, A302492, A356930, A356935, A356939, A356940.
A000041 counts integer partitions, strict A000009.
A000688 counts factorizations into prime powers.
A001055 counts factorizations.
A001221 counts prime divisors, sum A001414.
A001222 counts prime factors with multiplicity.
A011782 counts multisets covering an initial interval.
A056239 adds up prime indices, row sums of A112798.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    nogapQ[m_]:=Or[m=={},Union[m]==Range[Min[m],Max[m]]];
    Select[Range[100],And@@nogapQ/@primeMS/@primeMS[#]&]

A387115 Number of ways to choose a sequence of distinct strict integer partitions, one of each prime index of n.

Original entry on oeis.org

1, 1, 1, 0, 2, 1, 2, 0, 0, 2, 3, 0, 4, 2, 2, 0, 5, 0, 6, 0, 2, 3, 8, 0, 2, 4, 0, 0, 10, 2, 12, 0, 3, 5, 4, 0, 15, 6, 4, 0, 18, 2, 22, 0, 0, 8, 27, 0, 2, 2, 5, 0, 32, 0, 6, 0, 6, 10, 38, 0, 46, 12, 0, 0, 8, 3, 54, 0, 8, 4, 64, 0, 76, 15, 2, 0, 6, 4, 89, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Aug 20 2025

Keywords

Comments

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 axiom of choice says that, given any sequence of nonempty sets, it is possible to choose a sequence containing an element from each. In the strict version, the elements of this sequence must be distinct, meaning none is chosen more than once.

Examples

			The prime indices of 15 are (2,3), and there are a(15) = 2 choices:
  ((2),(3))
  ((2),(2,1))
The prime indices of 121 are (5,5), and there are a(121) = 6 choices:
  ((5),(4,1))
  ((5),(3,2))
  ((4,1),(5))
  ((4,1),(3,2))
  ((3,2),(5))
  ((3,2),(4,1))
		

Crossrefs

For divisors instead of partitions we have A355739, see A355740, A355733, A355734.
Allowing repeated partitions gives A357982, see A299200, A357977, A357978.
Twice-partitions of this type are counted by A358914, strict case of A270995.
The disjoint case is A383706.
Allowing non-strict partitions gives A387110, for prime factors A387133.
For initial intervals instead of strict partitions we have A387111.
For constant instead of strict partitions we have A387120.
Positions of 0 are A387176 (non-choosable), complement A387177 (choosable).
A000041 counts integer partitions, strict A000009.
A003963 multiplies together the prime indices of n.
A112798 lists prime indices, row sums A056239 or A066328, lengths A001222.
A120383 lists numbers divisible by all of their prime indices.
A289509 lists numbers with relatively prime prime indices.

Programs

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

A336343 Number of ways to choose a strict partition of each part of a strict composition of n.

Original entry on oeis.org

1, 1, 1, 4, 6, 11, 26, 39, 78, 142, 320, 488, 913, 1558, 2798, 5865, 9482, 16742, 28474, 50814, 82800, 172540, 266093, 472432, 790824, 1361460, 2251665, 3844412, 7205416, 11370048, 19483502, 32416924, 54367066, 88708832, 149179800, 239738369, 445689392
Offset: 0

Views

Author

Gus Wiseman, Jul 19 2020

Keywords

Comments

A strict composition of n (A032020) is a finite sequence of distinct positive integers summing to n.
Is there a simple generating function?

Examples

			The a(1) = 1 through a(5) = 11 ways:
  (1)  (2)  (3)      (4)        (5)
            (2,1)    (3,1)      (3,2)
            (1),(2)  (1),(3)    (4,1)
            (2),(1)  (3),(1)    (1),(4)
                     (1),(2,1)  (2),(3)
                     (2,1),(1)  (3),(2)
                                (4),(1)
                                (1),(3,1)
                                (2,1),(2)
                                (2),(2,1)
                                (3,1),(1)
		

Crossrefs

Multiset partitions of partitions are A001970.
Strict compositions are counted by A032020, A072574, and A072575.
Splittings of strict partitions are A072706.
Set partitions of strict partitions are A294617.
Splittings of partitions with distinct sums are A336131.
Partitions:
- Partitions of each part of a partition are A063834.
- Compositions of each part of a partition are A075900.
- Strict partitions of each part of a partition are A270995.
- Strict compositions of each part of a partition are A336141.
Strict partitions:
- Partitions of each part of a strict partition are A271619.
- Compositions of each part of a strict partition are A304961.
- Strict partitions of each part of a strict partition are A279785.
- Strict compositions of each part of a strict partition are A336142.
Compositions:
- Partitions of each part of a composition are A055887.
- Compositions of each part of a composition are A133494.
- Strict partitions of each part of a composition are A304969.
- Strict compositions of each part of a composition are A307068.
Strict compositions:
- Partitions of each part of a strict composition are A336342.
- Compositions of each part of a strict composition are A336127.
- Strict partitions of each part of a strict composition are A336343.
- Strict compositions of each part of a strict composition are A336139.

Programs

  • Mathematica
    strptn[n_]:=Select[IntegerPartitions[n],UnsameQ@@#&];
    Table[Length[Join@@Table[Tuples[strptn/@ctn],{ctn,Join@@Permutations/@strptn[n]}]],{n,0,10}]
  • PARI
    \\ here Q(N) gives A000009 as a vector.
    Q(n) = {Vec(eta(x^2 + O(x*x^n))/eta(x + O(x*x^n)))}
    seq(n)={my(b=Q(n)); [subst(serlaplace(p),y,1) | p<-Vec(prod(k=1, n, 1 + y*x^k*b[1+k] + O(x*x^n)))]} \\ Andrew Howroyd, Apr 16 2021

Formula

G.f.: Sum_{k>=0} k! * [y^k](Product_{j>=1} 1 + y*x^j*A000009(j)). - Andrew Howroyd, Apr 16 2021

A356934 Number of multisets of odd-size multisets whose multiset union is a size-n multiset covering an initial interval with weakly decreasing multiplicities.

Original entry on oeis.org

1, 1, 2, 6, 17, 46, 166, 553, 2093
Offset: 0

Views

Author

Gus Wiseman, Sep 09 2022

Keywords

Examples

			The a(1) = 1 through a(4) = 17 multiset partitions:
  {{1}}  {{1},{1}}  {{1,1,1}}      {{1},{1,1,1}}
         {{1},{2}}  {{1,1,2}}      {{1},{1,1,2}}
                    {{1,2,3}}      {{1},{1,2,2}}
                    {{1},{1},{1}}  {{1},{1,2,3}}
                    {{1},{1},{2}}  {{1},{2,3,4}}
                    {{1},{2},{3}}  {{2},{1,1,1}}
                                   {{2},{1,1,2}}
                                   {{2},{1,1,3}}
                                   {{2},{1,3,4}}
                                   {{3},{1,1,2}}
                                   {{3},{1,2,4}}
                                   {{4},{1,2,3}}
                                   {{1},{1},{1},{1}}
                                   {{1},{1},{1},{2}}
                                   {{1},{1},{2},{2}}
                                   {{1},{1},{2},{3}}
                                   {{1},{2},{3},{4}}
		

Crossrefs

A000041 counts integer partitions, strict A000009.
A000670 counts patterns, ranked by A333217, necklace A019536.
A011782 counts multisets covering an initial interval.
Odd-size multisets are counted by A000302, A027193, A058695, ranked by A026424.
Other conditions: A035310, A063834, A330783, A356938, A356943, A356954.
Other types: A050330, A356932, A356933, A356935.

Programs

  • 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]]]];
    strnorm[n_]:=Flatten[MapIndexed[Table[#2,{#1}]&,#]]&/@IntegerPartitions[n];
    Table[Length[Select[Join@@mps/@strnorm[n],OddQ[Times@@Length/@#]&]],{n,0,5}]

A356937 Number of multisets of intervals whose multiset union is of size n and covers an initial interval of positive integers.

Original entry on oeis.org

1, 1, 3, 9, 29, 94, 310, 1026, 3411, 11360, 37886, 126442, 422203, 1410189, 4711039, 15740098, 52593430, 175742438, 587266782, 1962469721, 6558071499, 21915580437, 73237274083, 244744474601, 817889464220, 2733235019732, 9133973730633, 30524096110942, 102006076541264
Offset: 0

Views

Author

Gus Wiseman, Sep 08 2022

Keywords

Comments

An interval such as {3,4,5} is a set with all differences of adjacent elements equal to 1.

Examples

			The a(1) = 1 through a(3) = 9 set multipartitions (multisets of sets):
  {{1}}  {{1,2}}    {{1,2,3}}
         {{1},{1}}  {{1},{1,2}}
         {{1},{2}}  {{1},{2,3}}
                    {{2},{1,2}}
                    {{3},{1,2}}
                    {{1},{1},{1}}
                    {{1},{1},{2}}
                    {{1},{2},{2}}
                    {{1},{2},{3}}
		

Crossrefs

A000041 counts integer partitions, strict A000009.
A000670 counts patterns, ranked by A333217, necklace A019536.
A011782 counts multisets covering an initial interval.
Intervals are counted by A000012, A001227, ranked by A073485.
Other conditions: A034691, A116540, A255906, A356933, A356942.
Other types: A107742, A356936, A356938, A356939.

Programs

  • Mathematica
    allnorm[n_]:=If[n<=0,{{}},Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1]];
    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]]]];
    chQ[y_]:=Or[Length[y]<=1,Union[Differences[y]]=={1}];
    Table[Length[Select[Join@@mps/@allnorm[n],And@@chQ/@#&]],{n,0,5}]
  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
    R(n,k) = {EulerT(vector(n, j, max(0, 1+k-j)))}
    seq(n) = {my(A=1+O(y*y^n)); for(k = 1, n, A += x^k*(1 + y*Ser(R(n,k), y) - polcoef(1/(1 - x*A) + O(x^(k+2)), k+1))); Vec(subst(A,x,1))} \\ Andrew Howroyd, Jan 01 2023

Extensions

Terms a(10) and beyond from Andrew Howroyd, Jan 01 2023
Previous Showing 51-60 of 102 results. Next