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 41-50 of 466 results. Next

A306021 Number of set-systems spanning {1,...,n} in which all sets have the same size.

Original entry on oeis.org

1, 1, 2, 6, 54, 1754, 1102746, 68715913086, 1180735735356265746734, 170141183460507906731293351306656207090, 7237005577335553223087828975127304177495735363998991435497132232365910414322
Offset: 0

Views

Author

Gus Wiseman, Jun 17 2018

Keywords

Comments

a(n) is the number of labeled uniform hypergraphs spanning n vertices. - Andrew Howroyd, Jan 16 2024

Examples

			The a(3) = 6 set-systems in which all sets have the same size:
  {{1,2,3}}
  {{1}, {2}, {3}}
  {{1,2}, {1,3}}
  {{1,2}, {2,3}}
  {{1,3}, {2,3}}
  {{1,2}, {1,3}, {2,3}}
		

Crossrefs

Row sums of A299471.
The unlabeled version is A301481.
The connected version is A299353.

Programs

  • Mathematica
    Table[Sum[(-1)^(n-k)*Binomial[n,k]*(1+Sum[2^Binomial[k,d]-1,{d,k}]),{k,0,n}],{n,12}]
  • PARI
    a(n) = if(n==0, 1, sum(k=0, n, sum(d=0, n, (-1)^(n-d)*binomial(n,d)*2^binomial(d,k)))) \\ Andrew Howroyd, Jan 16 2024

Formula

a(n) = Sum_{k = 0..n} (-1)^(n-k)*binomial(n,k)*(1 - k + Sum_{d = 1..k} 2^binomial(k, d)).
Inverse binomial transform of A306020. - Andrew Howroyd, Jan 16 2024

A319559 Number of non-isomorphic T_0 set systems of weight n.

Original entry on oeis.org

1, 1, 1, 2, 4, 7, 16, 35, 82, 200, 517, 1373, 3867, 11216, 33910, 105950
Offset: 0

Views

Author

Gus Wiseman, Sep 23 2018

Keywords

Comments

In a set system, two vertices are equivalent if in every block the presence of the first is equivalent to the presence of the second. The T_0 condition means that there are no equivalent vertices.
The weight of a set system is the sum of sizes of its parts. Weight is generally not the same as number of vertices.

Examples

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

Crossrefs

Extensions

a(11)-a(15) from Bert Dobbelaere, May 04 2025

A368413 Number of factorizations of n into positive integers > 1 such that it is not possible to choose a different prime factor of each factor.

Original entry on oeis.org

0, 0, 0, 1, 0, 0, 0, 2, 1, 0, 0, 1, 0, 0, 0, 4, 0, 1, 0, 1, 0, 0, 0, 3, 1, 0, 2, 1, 0, 0, 0, 6, 0, 0, 0, 4, 0, 0, 0, 3, 0, 0, 0, 1, 1, 0, 0, 7, 1, 1, 0, 1, 0, 3, 0, 3, 0, 0, 0, 2, 0, 0, 1, 10, 0, 0, 0, 1, 0, 0, 0, 10, 0, 0, 1, 1, 0, 0, 0, 7, 4, 0, 0, 2, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Dec 27 2023

Keywords

Comments

For example, the factorization f = 2*3*6 has two ways to choose a prime factor of each factor, namely (2,3,2) and (2,3,3), but neither of these has all different elements, so f is counted under a(36).

Examples

			The a(1) = 0 through a(24) = 3 factorizations:
 ... 2*2 ... 2*4   3*3 .. 2*2*3 ... 2*8     . 2*3*3 . 2*2*5 ... 2*2*6
             2*2*2                  4*4                         2*3*4
                                    2*2*4                       2*2*2*3
                                    2*2*2*2
		

Crossrefs

For unlabeled graphs: A140637, complement A134964.
For labeled graphs: A367867, A367868, A140638, complement A133686.
For set-systems: A367903, ranks A367907, complement A367902, ranks A367906.
For non-isomorphic set-systems: A368094, A368409, complement A368095.
For non-isomorphic multiset partitions: A368097, A355529, A368411.
Complement for non-isomorphic multiset partitions: A368098, A368100.
The complement is counted by A368414.
For non-isomorphic set multipartitions: A368421, complement A368422.
For divisors instead of prime factors: A370813, complement A370814.
A001055 counts factorizations, strict A045778.
A007716 counts non-isomorphic multiset partitions, connected A007718.
A058891 counts set-systems, unlabeled A000612, connected A323818.
A283877 counts non-isomorphic set-systems, connected A300913.

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], Select[Tuples[First/@FactorInteger[#]&/@#], UnsameQ@@#&]=={}&]],{n,100}]

Formula

a(n) + A368414(n) = A001055(n).

A318361 Number of strict set multipartitions (sets of sets) of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 0, 2, 0, 1, 0, 5, 1, 0, 0, 4, 0, 0, 0, 15, 0, 5, 0, 1, 0, 0, 0, 16, 0, 0, 8, 0, 0, 2, 0, 52, 0, 0, 0, 23, 0, 0, 0, 7, 0, 0, 0, 0, 5, 0, 0, 68, 0, 1, 0, 0, 0, 40, 0, 1, 0, 0, 0, 14, 0, 0, 1, 203, 0, 0, 0, 0, 0, 0, 0, 111, 0, 0, 4, 0, 0, 0, 0, 41, 80, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Aug 24 2018

Keywords

Examples

			The a(24) = 16 sets of sets with multiset union {1,1,2,3,4}:
  {{1},{1,2,3,4}}
  {{1,2},{1,3,4}}
  {{1,3},{1,2,4}}
  {{1,4},{1,2,3}}
  {{1},{2},{1,3,4}}
  {{1},{3},{1,2,4}}
  {{1},{4},{1,2,3}}
  {{1},{1,2},{3,4}}
  {{1},{1,3},{2,4}}
  {{1},{1,4},{2,3}}
  {{2},{1,3},{1,4}}
  {{3},{1,2},{1,4}}
  {{4},{1,2},{1,3}}
  {{1},{2},{3},{1,4}}
  {{1},{2},{4},{1,3}}
  {{1},{3},{4},{1,2}}
		

Crossrefs

Programs

  • Mathematica
    nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
    sqfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sqfacs[n/d],Min@@#>d&]],{d,Select[Rest[Divisors[n]],SquareFreeQ]}]];
    Table[Length[sqfacs[Times@@Prime/@nrmptn[n]]],{n,90}]
  • 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}
    sig(n)={my(f=factor(n)); concat(vector(#f~, i, vector(f[i, 2], j, primepi(f[i, 1]))))}
    count(sig)={my(r=0, A=O(x*x^vecmax(sig))); for(n=1, vecsum(sig)+1, my(s=0); forpart(p=n, my(q=prod(i=1, #p, 1 + x^p[i] + A)); s+=prod(i=1, #sig, polcoef(q, sig[i]))*(-1)^#p*permcount(p)); r+=(-1)^n*s/n!); r/2}
    a(n)={if(n==1, 1, my(s=sig(n)); if(#s==1, s[1]==1, count(sig(n))))} \\ Andrew Howroyd, Dec 18 2018

Formula

a(n) = A050326(A181821(n)).
a(prime(n)^k) = A188445(n, k). - Andrew Howroyd, Dec 17 2018

A368097 Number of non-isomorphic multiset partitions of weight n contradicting a strict version of the axiom of choice.

Original entry on oeis.org

0, 0, 1, 3, 12, 37, 133, 433, 1516, 5209, 18555
Offset: 0

Views

Author

Gus Wiseman, Dec 25 2023

Keywords

Comments

A multiset partition is a finite multiset of finite nonempty multisets. The weight of a multiset partition is the sum of cardinalities of its elements. Weight is generally not the same as number of vertices.
The axiom of choice says that, given any set of nonempty sets Y, it is possible to choose a set containing an element from each. The strict version requires this set to have the same cardinality as Y, meaning no element is chosen more than once.

Examples

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

Crossrefs

The case of unlabeled graphs appears to be A140637, complement A134964.
These multiset partitions have ranks A355529.
The case of labeled graphs is A367867, complement A133686.
Set-systems not of this type are A367902, ranks A367906.
Set-systems of this type are A367903, ranks A367907.
For set-systems we have A368094, complement A368095.
The complement is A368098, ranks A368100, connected case A368412.
Minimal multiset partitions of this type are ranked by A368187.
The connected case is A368411.
Factorizations of this type are counted by A368413, complement A368414.
For set multipartitions we have A368421, complement A368422.
A000110 counts set partitions, non-isomorphic A000041.
A003465 counts covering set-systems, unlabeled A055621.
A007716 counts non-isomorphic multiset partitions, connected A007718.
A058891 counts set-systems, unlabeled A000612, connected A323818.
A283877 counts non-isomorphic set-systems, connected A300913.

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]& /@ sps[Complement[set,s]]] /@ Cases[Subsets[set],{i,_}];
    mpm[n_]:=Join@@Table[Union[Sort[Sort/@(#/.x_Integer:>s[[x]])]& /@ sps[Range[n]]], {s,Flatten[MapIndexed[Table[#2,{#1}]&,#]]& /@ IntegerPartitions[n]}];
    brute[m_]:=First[Sort[Table[Sort[Sort /@ (m/.Rule@@@Table[{i,p[[i]]},{i,Length[p]}])], {p,Permutations[Union@@m]}]]];
    Table[Length[Union[brute/@Select[mpm[n], Select[Tuples[#],UnsameQ@@#&]=={}&]]], {n,0,6}]

A317757 Number of non-isomorphic multiset partitions of size n such that the blocks have empty intersection.

Original entry on oeis.org

1, 0, 1, 4, 17, 56, 205, 690, 2446, 8506, 30429, 109449, 402486, 1501424, 5714194, 22132604, 87383864, 351373406, 1439320606, 6003166059, 25488902820, 110125079184, 483987225922, 2162799298162, 9823464989574, 45332196378784, 212459227340403, 1010898241558627, 4881398739414159
Offset: 0

Views

Author

Gus Wiseman, Aug 06 2018

Keywords

Examples

			Non-isomorphic representatives of the a(4) = 17 multiset partitions:
  {1}{234},{2}{111},{2}{113},{11}{22},{11}{23},{12}{34},
  {1}{1}{22},{1}{1}{23},{1}{2}{11},{1}{2}{12},{1}{2}{13},{1}{2}{34},{2}{3}{11},
  {1}{1}{1}{2},{1}{1}{2}{2},{1}{1}{2}{3},{1}{2}{3}{4}.
		

Crossrefs

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];
    sysnorm[m_]:=If[Union@@m!=Range[Max@@Flatten[m]],sysnorm[m/.Rule@@@Table[{(Union@@m)[[i]],i},{i,Length[Union@@m]}]],First[Sort[sysnorm[m,1]]]];sysnorm[m_,aft_]:=If[Length[Union@@m]<=aft,{m},With[{mx=Table[Count[m,i,{2}],{i,Select[Union@@m,#>=aft&]}]},Union@@(sysnorm[#,aft+1]&/@Union[Table[Map[Sort,m/.{par+aft-1->aft,aft->par+aft-1},{0,1}],{par,First/@Position[mx,Max[mx]]}]])]];
    Table[Length[Union[sysnorm/@Join@@Table[Select[mps[m],Intersection@@#=={}&],{m,strnorm[n]}]]],{n,6}]
  • 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}
    K(q, t, k)={EulerT(Vec(sum(j=1, #q, gcd(t, q[j])*x^lcm(t, q[j])) + O(x*x^k), -k))}
    R(q, n)={vector(n, t, x*Ser(K(q, t, n)/t))}
    a(n)={my(s=0); forpart(q=n, my(f=prod(i=1, #q, 1 - x^q[i]), u=R(q,n)); s+=permcount(q)*sum(k=0, n, my(c=polcoef(f,k)); if(c, c*polcoef(exp(sum(t=1, n\(k+1), x^(t*k)*u[t], O(x*x^n) ))/if(k,1-x^k,1), n))) ); s/n!} \\ Andrew Howroyd, May 30 2023

Extensions

a(8)-a(10) from Gus Wiseman, Sep 27 2018
a(0)=1 prepended and terms a(11) and beyond from Andrew Howroyd, May 30 2023

A330972 Sorted list containing the least number with each possible nonzero number of factorizations into factors > 1.

Original entry on oeis.org

1, 4, 8, 12, 16, 24, 36, 48, 60, 72, 96, 120, 128, 144, 180, 192, 216, 240, 256, 288, 360, 384, 420, 432, 480, 576, 720, 768, 840, 864, 900, 960, 1024, 1080, 1152, 1260, 1440, 1680, 1728, 1800, 1920, 2048, 2160, 2304, 2520, 2592, 2880, 3072, 3360, 3456, 3600
Offset: 1

Views

Author

Gus Wiseman, Jan 06 2020

Keywords

Comments

This is the sorted list of positions of first appearances in A001055 of each element of the range (A045782).

Examples

			Factorizations of n for n = 4, 8, 12, 16, 24, 36, 48, 60:
  4    8      12     16       24       36       48         60
  2*2  2*4    2*6    2*8      3*8      4*9      6*8        2*30
       2*2*2  3*4    4*4      4*6      6*6      2*24       3*20
              2*2*3  2*2*4    2*12     2*18     3*16       4*15
                     2*2*2*2  2*2*6    3*12     4*12       5*12
                              2*3*4    2*2*9    2*3*8      6*10
                              2*2*2*3  2*3*6    2*4*6      2*5*6
                                       3*3*4    3*4*4      3*4*5
                                       2*2*3*3  2*2*12     2*2*15
                                                2*2*2*6    2*3*10
                                                2*2*3*4    2*2*3*5
                                                2*2*2*2*3
		

Crossrefs

All terms belong to A025487
Includes all highly factorable numbers A033833.
Factorizations are A001055, with image A045782.
The least number with A045782(n) factorizations is A045783(n).
The least number with n factorizations is A330973(n).
The strict version is A330997.

Programs

  • Mathematica
    nn=1000;
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    nds=Length/@Array[facs,nn];
    Table[Position[nds,i][[1,1]],{i,First/@Gather[nds]}]

A304867 Number of non-isomorphic hypertrees of weight n.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 5, 6, 13, 20, 41, 70, 144, 266, 545, 1072, 2210, 4491, 9388, 19529, 41286, 87361, 186657, 399927, 862584, 1866461, 4058367, 8852686, 19384258, 42570435, 93783472, 207157172, 458805044, 1018564642, 2266475432, 5053991582, 11292781891, 25280844844
Offset: 0

Views

Author

Gus Wiseman, May 20 2018

Keywords

Comments

A hypertree E is a connected antichain of finite sets satisfying Sum_{e in E} (|e| - 1) = |U(E)| - 1. The weight of a hypertree is the sum of cardinalities of its elements. Weight is generally not the same as number of vertices (see A035053).
From Kevin Ryde, Feb 25 2020: (Start)
a(n), except at n=1, is the number of free trees of n edges (so n+1 vertices) where any two leaves are an even distance apart. All trees are bipartite graphs and this condition is equivalent to all leaves being in the same bipartite half. The diameter of a tree is always between two leaves so these trees have even diameter (A000676).
The correspondence between hypertrees and these free trees is described for instance by Bacher (start of section 1.2). In such a free tree, call a vertex "even" if it is an even distance from a leaf. The hypertree vertices are these even vertices. Each hyperedge is the set of vertices surrounding an odd vertex, so hypertree weight is the total number of edges in the free tree.
(End)

Examples

			Non-isomorphic representatives of the a(6) = 5 hypertrees are the following:
  {{1,2,3,4,5,6}}
  {{1,2},{1,3,4,5}}
  {{1,2,3},{1,4,5}}
  {{1,2},{1,3},{1,4}}
  {{1,2},{1,3},{2,4}}
Non-isomorphic representatives of the a(7) = 6 hypertrees are the following:
  {{1,2,3,4,5,6,7}}
  {{1,2},{1,3,4,5,6}}
  {{1,2,3},{1,4,5,6}}
  {{1,2},{1,3},{1,4,5}}
  {{1,2},{1,3},{2,4,5}}
  {{1,3},{2,4},{1,2,5}}
From _Kevin Ryde_, Feb 25 2020: (Start)
a(6) = 5 hypertrees of weight 6 and their corresponding free trees of 6 edges (7 vertices).  Each * is an "odd" vertex (odd distance to a leaf).  Each hyperedge is the set of "even" vertices surrounding an odd.
  {1,2,3,4,5,6}       3   2
                       \ /
                      4-*-1      (star 7)
                       / \
                      5   6
  .
  {1,2},{1,3,4,5}               /-3
                      2--*--1--*--4
                                \-5
  .
  {1,2,3},{1,4,5}     2-\       /-4
                         *--1--*
                      3-/       \-5
  .
  {1,2},{1,3},{1,4}    /-*--2
                      1--*--3
                       \-*--4
  .
  {1,2},{2,4},{1,3}   3--*--1--*--2--*--4   (path 7)
(End)
		

Crossrefs

Programs

  • Mathematica
    etr[p_] := Module[{b}, b[n_] := b[n] = If[n == 0, 1, Sum[Sum[d*p[d], {d, Divisors[j]}]*b[n - j], {j, 1, n}]/n]; b];
    EulerT[v_List] := With[{q = etr[v[[#]]&]}, q /@ Range[Length[v]]];
    ser[v_] := Sum[v[[i]] x^(i-1), {i, 1, Length[v]}] + O[x]^Length[v];
    c[n_] := Module[{v = {1}}, For[i = 1, i <= Ceiling[n/2], i++, v = Join[{1}, EulerT[Join[{0}, EulerT[v]]]]]; v];
    seq[n_] := Module[{u = c[n]}, x*ser[EulerT[u]]*(1 - x*ser[u]) + (1 - x)* ser[u] + x + O[x]^n // CoefficientList[#, x]&];
    seq[40] (* Jean-François Alcover, Feb 08 2020, after Andrew Howroyd *)
  • PARI
    EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
    c(n)={my(v=[1]); for(i=1, ceil(n/2), v=concat([1], EulerT(concat([0], EulerT(v))))); v}
    seq(n)={my(u=c(n)); Vec(x*Ser(EulerT(u))*(1-x*Ser(u)) + (1 - x)*Ser(u) + x + O(x*x^n))} \\ Andrew Howroyd, Aug 29 2018

Formula

a(n) = Sum_{k=1..floor(n/2)} A318601(n+1-k, k). - Andrew Howroyd, Aug 29 2018

Extensions

Terms a(10) and beyond from Andrew Howroyd, Aug 29 2018

A319646 Number of non-isomorphic weight-n chains of distinct multisets whose dual is also a chain of distinct multisets.

Original entry on oeis.org

1, 1, 1, 4, 4, 9, 17, 28, 41, 75, 122, 192, 314, 484, 771, 1216, 1861, 2848, 4395, 6610, 10037
Offset: 0

Views

Author

Gus Wiseman, Sep 25 2018

Keywords

Comments

The dual of a multiset partition has, for each vertex, one block consisting of the indices (or positions) of the blocks containing that vertex, counted with multiplicity. For example, the dual of {{1,2},{2,2}} is {{1},{1,2,2}}.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.
From Gus Wiseman, Jan 17 2019: (Start)
Also the number of plane partitions of n with no repeated rows or columns. For example, the a(6) = 17 plane partitions are:
6 51 42 321
.
5 4 41 31 32 31 22 221 211
1 2 1 2 1 11 2 1 11
.
3 21 21 111
2 2 11 11
1 1 1 1
(End)

Examples

			Non-isomorphic representatives of the a(1) = 1 through a(5) = 9 chains:
1: {{1}}
2: {{1,1}}
3: {{1,1,1}}
   {{1,2,2}}
   {{1},{1,1}}
   {{2},{1,2}}
4: {{1,1,1,1}}
   {{1,2,2,2}}
   {{1},{1,1,1}}
   {{2},{1,2,2}}
5: {{1,1,1,1,1}}
   {{1,1,2,2,2}}
   {{1,2,2,2,2}}
   {{1},{1,1,1,1}}
   {{2},{1,1,2,2}}
   {{2},{1,2,2,2}}
   {{1,1},{1,1,1}}
   {{1,2},{1,2,2}}
   {{2,2},{1,2,2}}
		

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[Reverse@*primeMS,Join@@Permutations/@facs[n],{2}]];
    Table[Sum[Length[Select[ptnplane[Times@@Prime/@y],And[UnsameQ@@#,UnsameQ@@Transpose[PadRight[#]],And@@GreaterEqual@@@#,And@@(GreaterEqual@@@Transpose[PadRight[#]])]&]],{y,IntegerPartitions[n]}],{n,10}] (* Gus Wiseman, Jan 18 2019 *)

Extensions

a(11)-a(17) from Gus Wiseman, Jan 18 2019
a(18)-a(21) from Robert Price, Jun 21 2021

A368414 Number of factorizations of n into positive integers > 1 such that it is possible to choose a different prime factor of each factor.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Dec 29 2023

Keywords

Comments

For example, the factorization f = 2*3*6 has two ways to choose a prime factor of each factor, namely (2,3,2) and (2,3,3), but neither of these has all different elements, so f is not counted under a(36).

Examples

			The a(n) factorizations for selected n:
  1    6      12     24      30       60        72      120
       2*3    2*6    2*12    2*15     2*30      2*36    2*60
              3*4    3*8     3*10     3*20      3*24    3*40
                     4*6     5*6      4*15      4*18    4*30
                             2*3*5    5*12      6*12    5*24
                                      6*10      8*9     6*20
                                      2*3*10            8*15
                                      2*5*6             10*12
                                      3*4*5             2*3*20
                                                        2*5*12
                                                        2*6*10
                                                        3*4*10
                                                        3*5*8
                                                        4*5*6
		

Crossrefs

For labeled graphs: A133686, complement A367867, A367868, A140638.
For unlabeled graphs: A134964, complement A140637.
For set-systems: A367902, ranks A367906, complement A367903, ranks A367907.
For non-isomorphic set-systems: A368095, complement A368094, A368409.
Complementary non-isomorphic multiset partitions: A368097, A355529, A368411.
For non-isomorphic multiset partitions: A368098, A368100.
The complement is counted by A368413.
For non-isomorphic set multipartitions: A368422, complement A368421.
For divisors instead of prime factors: A370813, complement A370814.
A001055 counts factorizations, strict A045778.
A007716 counts non-isomorphic multiset partitions, connected A007718.
A058891 counts set-systems, unlabeled A000612, connected A323818.
A283877 counts non-isomorphic set-systems, connected A300913.

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], Select[Tuples[First/@FactorInteger[#]&/@#], UnsameQ@@#&]!={}&]],{n,100}]

Formula

a(n) = A001055(n) - A368413(n).
Previous Showing 41-50 of 466 results. Next