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 11-20 of 107 results. Next

A292886 Number of knapsack factorizations of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 3, 2, 2, 1, 4, 1, 2, 2, 4, 1, 4, 1, 4, 2, 2, 1, 7, 2, 2, 3, 4, 1, 5, 1, 6, 2, 2, 2, 8, 1, 2, 2, 7, 1, 5, 1, 4, 4, 2, 1, 11, 2, 4, 2, 4, 1, 7, 2, 7, 2, 2, 1, 11, 1, 2, 4, 7, 2, 5, 1, 4, 2, 5, 1, 15, 1, 2, 4, 4, 2, 5, 1, 11, 4, 2, 1, 11, 2
Offset: 1

Views

Author

Gus Wiseman, Sep 26 2017

Keywords

Comments

A knapsack factorization is a finite multiset of positive integers greater than one such that every distinct submultiset has a different product.
The sequence giving the number of factorizations of n is described as "the multiplicative partition function" (see A001055), so knapsack factorizations are a multiplicative generalization of knapsack partitions. - Gus Wiseman, Oct 24 2017

Examples

			The a(36) = 8 factorizations are 2*2*3*3, 2*2*9, 2*18, 3*3*4, 3*12, 4*9, 6*6, 36. The factorization 2*3*6 is not knapsack.
		

Crossrefs

Programs

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

A303707 Number of factorizations of n using elements of A007916 (numbers that are not perfect powers).

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, 8, 1, 2, 3, 3, 2, 5, 1, 5, 1, 2, 1, 9, 2, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Apr 29 2018

Keywords

Comments

First differs from A081707 at a(60) = 9, A081707(60) = 8.

Examples

			The a(60) = 9 factorizations are (2*2*3*5), (2*2*15), (2*3*10), (2*5*6), (2*30), (3*20), (5*12), (6*10), (60).
		

Crossrefs

Programs

  • Mathematica
    radQ[n_]:=Or[n===1,GCD@@FactorInteger[n][[All,2]]===1];
    facsr[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facsr[n/d],Min@@#>=d&]],{d,Select[Rest[Divisors[n]],radQ]}]];
    Table[Length[facsr[n]],{n,100}]

Formula

Dirichlet g.f.: Product_{n in A007916} 1/(1 - n^s).

A317791 Number of non-isomorphic multiset partitions of the multiset of prime indices of n (row n of A112798).

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Aug 07 2018

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n.
a(n) depends only on prime signature of n (cf. A025487). - Antti Karttunen, Dec 03 2018
Are any terms of the complement known? In particular, does this sequence contain 6? - Gus Wiseman, Oct 21 2022

Examples

			Non-isomorphic representatives of the a(42) = 3 multiset partitions are {{1,2,4}}, {{1},{2,4}}, {{1},{2},{4}}.
Non-isomorphic representatives of the a(60) = 9 multiset partitions:
  {1123},
  {1}{123}, {2}{113}, {11}{23}, {12}{13},
  {1}{1}{23}, {1}{2}{13}, {2}{3}{11},
  {1}{1}{2}{3}.
Missing from this list are {3}{112} and {1}{3}{12}, which are isomorphic to {2}{113} and {1}{2}{13} respectively.
For n = 180 = 2^2 * 3^2 * 5, there are A001055(180) = 26 different factorizations to one or more factors larger than 1. Of these 18 are such that by swapping 2 and 3 in each factor of that factorization the result is another, different factorization of 180, while the other 8 cases are such that 2 <-> 3 swap doesn't change the factorization. Thus a(180) = 18/2 + 8 = 17. - _Antti Karttunen_, Dec 03 2018
		

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]]]];
    sysnorm[{}] := {};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]]}]])]];
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Length[Union[sysnorm/@mps[primeMS[n]]]],{n,100}]

Formula

For all n, a(n) <= A001055(n). - Antti Karttunen, Dec 01 2018
If n is squarefree with k prime factors, or if n = p^k for p prime, we have a(n) = A000041(k).
a(n) = A318285(A181819(n)). - Andrew Howroyd, Jan 17 2023

Extensions

Terms corrected by Gus Wiseman, Dec 04 2018

A316978 Number of factorizations of n into factors > 1 with no equivalent primes.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 4, 1, 1, 1, 5, 1, 4, 1, 4, 1, 1, 1, 7, 2, 1, 3, 4, 1, 1, 1, 7, 1, 1, 1, 7, 1, 1, 1, 7, 1, 1, 1, 4, 4, 1, 1, 12, 2, 4, 1, 4, 1, 7, 1, 7, 1, 1, 1, 7, 1, 1, 4, 11, 1, 1, 1, 4, 1, 1, 1, 16, 1, 1, 4, 4, 1, 1, 1, 12, 5, 1, 1, 7, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Jul 18 2018

Keywords

Comments

In a factorization, two primes are equivalent if each factor has in its prime factorization the same multiplicity of both primes.

Examples

			The a(36) = 7 factorizations are (2*2*3*3), (2*2*9), (2*3*6), (3*3*4), (2*18), (3*12), (4*9). Missing from this list are (6*6) and (36).
		

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]]}]];
    dual[eds_]:=Table[First/@Position[eds,x],{x,Union@@eds}];
    Table[Length[Select[facs[n],UnsameQ@@dual[primeMS/@#]&]],{n,100}]

Formula

a(prime^n) = A000041(n).
a(squarefree) = 1.

A114592 Sum_{n>=1} a(n)/n^s = Product_{k>=2} (1 - 1/k^s).

Original entry on oeis.org

1, -1, -1, -1, -1, 0, -1, 0, -1, 0, -1, 1, -1, 0, 0, 0, -1, 1, -1, 1, 0, 0, -1, 1, -1, 0, 0, 1, -1, 1, -1, 1, 0, 0, 0, 1, -1, 0, 0, 1, -1, 1, -1, 1, 1, 0, -1, 1, -1, 1, 0, 1, -1, 1, 0, 1, 0, 0, -1, 1, -1, 0, 1, 0, 0, 1, -1, 1, 0, 1, -1, 1, -1, 0, 1, 1, 0, 1
Offset: 1

Views

Author

Leroy Quet, Dec 11 2005

Keywords

Comments

For n >= 2, Sum_{k|n} A001055(n/k) * a(k) = 0. A114591(n) = Sum_{k|n} a(k).
First entry greater than 1 in absolute value is a(360) = -2. - Gus Wiseman, Sep 15 2018

Examples

			24 can be factored into distinct integers (each >= 2) as 24; as 4*6, 3*8 and 2*12; and as 2*3*4. (A045778(24) = 5).
So a(24) = (-1)^1 + 3*(-1)^2 + (-1)^3 = 1, where the 1 exponent is due to the 1 factor of the 24 = 24 factorization and the 2 exponent is due to the 3 cases of 2 factors each of the 24 = 4*6 = 3*8 = 2*12 factorizations and the 3 exponent is due to the 24 = 2*3*4 factorization.
		

Crossrefs

Programs

  • Mathematica
    strfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strfacs[n/d],Min@@#>d&]],{d,Rest[Divisors[n]]}]];
    Table[Sum[(-1)^Length[f],{f,strfacs[n]}],{n,100}] (* Gus Wiseman, Sep 15 2018 *)
  • PARI
    A114592aux(n, k) = if(1==n, 1, sumdiv(n, d, if(d > 1 && d <= k && d < n, (-1)*A114592aux(n/d, d-1))) - (n<=k)); \\ After code in A045778.
    A114592(n) = A114592aux(n,n); \\ Antti Karttunen, Jul 23 2017

Formula

a(1) = 1; for n>= 2, a(n) = sum, over ways to factor n into any number of distinct integers >= 2, of (-1)^(number of integers in a factorization). (See example.)

Extensions

More terms from Antti Karttunen, Jul 23 2017

A316441 a(n) = Sum (-1)^k where the sum is over all factorizations of n into factors > 1 and k is the number of factors.

Original entry on oeis.org

1, -1, -1, 0, -1, 0, -1, -1, 0, 0, -1, 0, -1, 0, 0, 1, -1, 0, -1, 0, 0, 0, -1, 1, 0, 0, -1, 0, -1, 1, -1, -1, 0, 0, 0, 1, -1, 0, 0, 1, -1, 1, -1, 0, 0, 0, -1, 0, 0, 0, 0, 0, -1, 1, 0, 1, 0, 0, -1, 1, -1, 0, 0, 1, 0, 1, -1, 0, 0, 1, -1, 0, -1, 0, 0, 0, 0, 1, -1
Offset: 1

Views

Author

Gus Wiseman, Jul 03 2018

Keywords

Comments

First term greater than 1 in absolute value is a(256) = 2.

Examples

			The factorizations of 24 are (2*2*2*3), (2*2*6), (2*3*4), (2*12), (3*8), (4*6), (24); so a(24) = 1 - 2 + 3 - 1 = 1.
		

Crossrefs

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Sum[(-1)^Length[f],{f,facs[n]}],{n,200}]
  • PARI
    A316441(n, m=n, k=0) = if(1==n, (-1)^k, my(s=0); fordiv(n, d, if((d>1)&&(d<=m), s += A316441(n/d, d, k+1))); (s)); \\ Antti Karttunen, Sep 08 2018, after Michael B. Porter's code for A001055

Formula

Dirichlet g.f.: Product_{n > 1} 1/(1 + 1/n^s).

Extensions

Secondary offset added by Antti Karttunen, Sep 08 2018

A336424 Number of factorizations of n where each factor belongs to A130091 (numbers with distinct prime multiplicities).

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Aug 03 2020

Keywords

Comments

A number's prime signature (row n of A124010) is the sequence of positive exponents in its prime factorization, so a number has distinct prime multiplicities iff all the exponents in its prime signature are distinct.

Examples

			The a(n) factorizations for n = 2, 4, 8, 60, 16, 36, 32, 48:
  2  4    8      5*12     16       4*9      32         48
     2*2  2*4    3*20     4*4      3*12     4*8        4*12
          2*2*2  3*4*5    2*8      3*3*4    2*16       3*16
                 2*2*3*5  2*2*4    2*18     2*4*4      3*4*4
                          2*2*2*2  2*2*9    2*2*8      2*24
                                   2*2*3*3  2*2*2*4    2*3*8
                                            2*2*2*2*2  2*2*12
                                                       2*2*3*4
                                                       2*2*2*2*3
		

Crossrefs

A327523 is the case when n is restricted to belong to A130091 also.
A001055 counts factorizations.
A007425 counts divisors of divisors.
A045778 counts strict factorizations.
A074206 counts ordered factorizations.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
A253249 counts nonempty chains of divisors.
A281116 counts factorizations with no common divisor.
A302696 lists numbers whose prime indices are pairwise coprime.
A305149 counts stable factorizations.
A320439 counts factorizations using A289509.
A327498 gives the maximum divisor with distinct prime multiplicities.
A336500 counts divisors of n in A130091 with quotient also in A130091.
A336568 = not a product of two numbers with distinct prime multiplicities.
A336569 counts maximal chains of elements of A130091.
A337256 counts chains of divisors.

Programs

  • Mathematica
    facsusing[s_,n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facsusing[Select[s,Divisible[n/d,#]&],n/d],Min@@#>=d&]],{d,Select[s,Divisible[n,#]&]}]];
    Table[Length[facsusing[Select[Range[2,n],UnsameQ@@Last/@FactorInteger[#]&],n]],{n,100}]

A339564 Number of ways to choose a distinct factor in a factorization of n (pointed factorizations).

Original entry on oeis.org

0, 1, 1, 2, 1, 3, 1, 4, 2, 3, 1, 7, 1, 3, 3, 7, 1, 7, 1, 7, 3, 3, 1, 14, 2, 3, 4, 7, 1, 10, 1, 12, 3, 3, 3, 17, 1, 3, 3, 14, 1, 10, 1, 7, 7, 3, 1, 26, 2, 7, 3, 7, 1, 14, 3, 14, 3, 3, 1, 25, 1, 3, 7, 19, 3, 10, 1, 7, 3, 10, 1, 36, 1, 3, 7, 7, 3, 10, 1, 26, 7, 3
Offset: 1

Views

Author

Gus Wiseman, Apr 10 2021

Keywords

Examples

			The pointed factorizations of n for n = 2, 4, 6, 8, 12, 24, 30:
  ((2))  ((4))    ((6))    ((8))      ((12))     ((24))       ((30))
         ((2)*2)  ((2)*3)  ((2)*4)    ((2)*6)    ((3)*8)      ((5)*6)
                  (2*(3))  (2*(4))    (2*(6))    (3*(8))      (5*(6))
                           ((2)*2*2)  ((3)*4)    ((4)*6)      ((2)*15)
                                      (3*(4))    (4*(6))      (2*(15))
                                      ((2)*2*3)  ((2)*12)     ((3)*10)
                                      (2*2*(3))  (2*(12))     (3*(10))
                                                 ((2)*2*6)    ((2)*3*5)
                                                 (2*2*(6))    (2*(3)*5)
                                                 ((2)*3*4)    (2*3*(5))
                                                 (2*(3)*4)
                                                 (2*3*(4))
                                                 ((2)*2*2*3)
                                                 (2*2*2*(3))
		

Crossrefs

The additive version is A000070 (strict: A015723).
The unpointed version is A001055 (strict: A045778, ordered: A074206, listed: A162247).
Allowing point (1) gives A057567.
Choosing a position instead of value gives A066637.
The ordered additive version is A336875.
A000005 counts divisors.
A001787 count normal multisets with a selected position.
A001792 counts compositions with a selected position.
A006128 counts partitions with a selected position.
A066186 count strongly normal multisets with a selected position.
A254577 counts ordered factorizations with a selected position.

Programs

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

Formula

a(n) = A057567(n) - A001055(n).
a(n) = Sum_{d|n, d>1} A001055(n/d).

A259936 Number of ways to express the integer n as a product of its unitary divisors (A034444).

Original entry on oeis.org

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, 5, 1, 1, 2, 2, 2, 2, 1, 2, 2, 2, 1, 5, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 5, 1, 2, 2, 1, 2, 5, 1, 2, 2, 5, 1, 2, 1, 2, 2, 2, 2, 5, 1, 2, 1, 2, 1, 5, 2, 2, 2, 2, 1, 5, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 1, 5, 1, 2, 5
Offset: 1

Views

Author

Geoffrey Critzer, Jul 09 2015

Keywords

Comments

Equivalently, a(n) is the number of ways to express the cyclic group Z_n as a direct sum of its Hall subgroups. A Hall subgroup of a finite group G is a subgroup whose order is coprime to its index.
a(n) is the number of ways to partition the set of distinct prime factors of n.
Also the number of singleton or pairwise coprime factorizations of n. - Gus Wiseman, Sep 24 2019

Examples

			a(60) = 5 because we have: 60 = 4*3*5 = 4*15 = 3*20 = 5*12.
For n = 36, its unitary divisors are 1, 4, 9, 36. From these we obtain 36 either as 1*36 or 4*9, thus a(36) = 2. - _Antti Karttunen_, Oct 21 2017
		

Crossrefs

Differs from A050320 for the first time at n=36.
Differs from A354870 for the first time at n=210, where a(210) = 15, while A354870(210) = 12.
Related classes of factorizations:
- No conditions: A001055
- Strict: A045778
- Constant: A089723
- Distinct multiplicities: A255231
- Singleton or coprime: A259936
- Relatively prime: A281116
- Aperiodic: A303386
- Stable (indivisible): A305149
- Connected: A305193
- Strict relatively prime: A318721
- Uniform: A319269
- Intersecting: A319786
- Constant or distinct factors coprime: A327399
- Constant or relatively prime: A327400
- Coprime: A327517
- Not relatively prime: A327658
- Distinct factors coprime: A327695

Programs

  • Maple
    map(combinat:-bell @ nops @ numtheory:-factorset, [$1..100]); # Robert Israel, Jul 09 2015
  • Mathematica
    Table[BellB[PrimeNu[n]], {n, 1, 75}]
    (* second program *)
    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],Length[#]==1||CoprimeQ@@#&]],{n,100}] (* Gus Wiseman, Sep 24 2019 *)
  • PARI
    a(n) = my(t=omega(n), x='x, m=contfracpnqn(matrix(2, t\2, y, z, if( y==1, -z*x^2, 1 - (z+1)*x)))); polcoeff(1/(1 - x + m[2, 1]/m[1, 1]) + O(x^(t+1)), t) \\ Charles R Greathouse IV, Jun 30 2017

Formula

a(n) = A000110(A001221(n)).
a(n > 1) = A327517(n) + 1. - Gus Wiseman, Sep 24 2019

Extensions

Incorrect comment removed by Antti Karttunen, Jun 11 2022

A319765 Number of non-isomorphic intersecting multiset partitions of weight n whose dual is also an intersecting multiset partition.

Original entry on oeis.org

1, 1, 3, 6, 15, 31, 74, 156, 358, 792, 1821
Offset: 0

Views

Author

Gus Wiseman, Sep 27 2018

Keywords

Comments

The dual of a multiset partition has, for each vertex, one part consisting of the indices (or positions) of the parts containing that vertex, counted with multiplicity. For example, the dual of {{1,2},{2,2}} is {{1},{1,2,2}}.
A multiset partition is intersecting iff no two parts are disjoint. The dual of a multiset partition is intersecting iff every pair of distinct vertices appear together in some part.
The weight of a multiset partition 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(4) = 15 multiset partitions:
1: {{1}}
2: {{1,1}}
   {{1,2}}
   {{1},{1}}
3: {{1,1,1}}
   {{1,2,2}}
   {{1,2,3}}
   {{1},{1,1}}
   {{2},{1,2}}
   {{1},{1},{1}}
4: {{1,1,1,1}}
   {{1,1,2,2}}
   {{1,2,2,2}}
   {{1,2,3,3}}
   {{1,2,3,4}}
   {{1},{1,1,1}}
   {{1},{1,2,2}}
   {{2},{1,2,2}}
   {{3},{1,2,3}}
   {{1,1},{1,1}}
   {{1,2},{1,2}}
   {{1,2},{2,2}}
   {{1},{1},{1,1}}
   {{2},{2},{1,2}}
   {{1},{1},{1},{1}}
		

Crossrefs

Previous Showing 11-20 of 107 results. Next