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

A320887 Number of multiset partitions of factorizations of n into factors > 1 such that all the parts have the same product.

Original entry on oeis.org

1, 1, 1, 3, 1, 2, 1, 4, 3, 2, 1, 4, 1, 2, 2, 9, 1, 4, 1, 4, 2, 2, 1, 7, 3, 2, 4, 4, 1, 5, 1, 8, 2, 2, 2, 12, 1, 2, 2, 7, 1, 5, 1, 4, 4, 2, 1, 12, 3, 4, 2, 4, 1, 7, 2, 7, 2, 2, 1, 11, 1, 2, 4, 22, 2, 5, 1, 4, 2, 5, 1, 16, 1, 2, 4, 4, 2, 5, 1, 12, 9, 2, 1, 11, 2, 2, 2, 7, 1, 11, 2, 4, 2, 2, 2, 19, 1, 4, 4, 12, 1, 5, 1, 7, 5
Offset: 1

Views

Author

Gus Wiseman, Oct 23 2018

Keywords

Examples

			The a(36) = 12 multiset partitions:
  (2*2*3*3)    (6)*(2*3)  (6)*(6)  (36)
  (2*3)*(2*3)  (2*2*9)    (2*18)
               (2*3*6)    (3*12)
               (3*3*4)    (4*9)
                          (6*6)
		

Crossrefs

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[With[{g=GCD@@FactorInteger[n][[All,2]]},Sum[Binomial[Length[facs[n^(1/d)]]+d-1,d],{d,Divisors[g]}]],{n,100}]
  • PARI
    A001055(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1)&&(d<=m), s += A001055(n/d, d))); (s));
    A052409(n) = { my(k=ispower(n)); if(k, k, n>1); }; \\ From A052409
    A320887(n) = if(1==n,n,my(r); sumdiv(A052409(n), d, binomial(A001055(sqrtnint(n,d)) + d - 1, d))); \\ Antti Karttunen, Nov 17 2019

Formula

a(n) = Sum_{d|A052409(n)} binomial(A001055(n^(1/d)) + d - 1, d).
a(n) = a(A046523(n)). - Antti Karttunen, Nov 17 2019

Extensions

Data section extended up to term a(105) by Antti Karttunen, Nov 17 2019

A320888 Number of set multipartitions (multisets of sets) of factorizations of n into factors > 1 such that all the parts have the same product.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 3, 2, 2, 1, 3, 1, 2, 2, 4, 1, 3, 1, 3, 2, 2, 1, 5, 2, 2, 3, 3, 1, 5, 1, 4, 2, 2, 2, 8, 1, 2, 2, 5, 1, 5, 1, 3, 3, 2, 1, 7, 2, 3, 2, 3, 1, 5, 2, 5, 2, 2, 1, 9, 1, 2, 3, 9, 2, 5, 1, 3, 2, 5, 1, 9, 1, 2, 3, 3, 2, 5, 1, 7, 4, 2, 1, 9, 2, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Oct 23 2018

Keywords

Examples

			The a(144) = 20 set multipartitions:
  (2*3*4*6)    (2*8*9)     (2*72)     (144)
  (2*6)*(2*6)  (3*6*8)     (3*48)
  (2*6)*(3*4)  (2*3*24)    (4*36)
  (3*4)*(3*4)  (2*4*18)    (6*24)
               (2*6*12)    (8*18)
               (3*4*12)    (9*16)
               (12)*(2*6)  (12)*(12)
               (12)*(3*4)
		

Crossrefs

Programs

  • Mathematica
    strfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strfacs[n/d],Min@@#>d&]],{d,Rest[Divisors[n]]}]];
    Table[With[{g=GCD@@FactorInteger[n][[All,2]]},Sum[Binomial[Length[strfacs[n^(1/d)]]+d-1,d],{d,Divisors[g]}]],{n,100}]

Formula

a(n) = Sum_{d|A052409(n)} binomial(A045778(n^(1/d)) + d - 1, d).

A320889 Number of set partitions of strict factorizations of n into factors > 1 such that all the blocks have the same product.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Oct 23 2018

Keywords

Examples

			The a(144) = 17 set partitions:
  (2*3*4*6)    (2*8*9)     (2*72)  (144)
  (2*6)*(3*4)  (3*6*8)     (3*48)
               (2*3*24)    (4*36)
               (2*4*18)    (6*24)
               (2*6*12)    (8*18)
               (3*4*12)    (9*16)
               (2*6)*(12)
               (3*4)*(12)
		

Crossrefs

Programs

  • Mathematica
    strfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strfacs[n/d],Min@@#>d&]],{d,Rest[Divisors[n]]}]];
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    Table[Length[Join@@Table[Select[sps[fac],SameQ@@Times@@@#&],{fac,strfacs[n]}]],{n,100}]
Showing 1-3 of 3 results.