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-10 of 18 results. Next

A038041 Number of ways to partition an n-set into subsets of equal size.

Original entry on oeis.org

1, 2, 2, 5, 2, 27, 2, 142, 282, 1073, 2, 32034, 2, 136853, 1527528, 4661087, 2, 227932993, 2, 3689854456, 36278688162, 13749663293, 2, 14084955889019, 5194672859378, 7905858780927, 2977584150505252, 13422745388226152, 2, 1349877580746537123, 2
Offset: 1

Views

Author

Keywords

Comments

a(n) = 2 iff n is prime with a(p) = card{ 1|2|3|...|p-1|p, 123...p } = 2. - Bernard Schott, May 16 2019

Examples

			a(4) = card{ 1|2|3|4, 12|34, 14|23, 13|24, 1234 } = 5.
From _Gus Wiseman_, Jul 12 2019: (Start)
The a(6) = 27 set partitions:
  {{1}{2}{3}{4}{5}{6}}  {{12}{34}{56}}  {{123}{456}}  {{123456}}
                        {{12}{35}{46}}  {{124}{356}}
                        {{12}{36}{45}}  {{125}{346}}
                        {{13}{24}{56}}  {{126}{345}}
                        {{13}{25}{46}}  {{134}{256}}
                        {{13}{26}{45}}  {{135}{246}}
                        {{14}{23}{56}}  {{136}{245}}
                        {{14}{25}{36}}  {{145}{236}}
                        {{14}{26}{35}}  {{146}{235}}
                        {{15}{23}{46}}  {{156}{234}}
                        {{15}{24}{36}}
                        {{15}{26}{34}}
                        {{16}{23}{45}}
                        {{16}{24}{35}}
                        {{16}{25}{34}}
(End)
		

Crossrefs

Cf. A061095 (same but with labeled boxes), A005225, A236696, A055225, A262280, A262320.
Column k=1 of A208437.
Row sums of A200472 and A200473.
Cf. A000110, A007837 (different lengths), A035470 (equal sums), A275780, A317583, A320324, A322794, A326512 (equal averages), A326513.

Programs

  • Maple
    A038041 := proc(n) local d;
    add(n!/(d!*(n/d)!^d), d = numtheory[divisors](n)) end:
    seq(A038041(n),n = 1..29); # Peter Luschny, Apr 16 2011
  • Mathematica
    a[n_] := Block[{d = Divisors@ n}, Plus @@ (n!/(#! (n/#)!^#) & /@ d)]; Array[a, 29] (* Robert G. Wilson v, Apr 16 2011 *)
    Table[Sum[n!/((n/d)!*(d!)^(n/d)), {d, Divisors[n]}], {n, 1, 31}] (* Emanuele Munarini, Jan 30 2014 *)
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    Table[Length[Select[sps[Range[n]],SameQ@@Length/@#&]],{n,0,8}] (* Gus Wiseman, Jul 12 2019 *)
  • Maxima
    a(n):= lsum(n!/((n/d)!*(d!)^(n/d)),d,listify(divisors(n)));
    makelist(a(n),n,1,40); /* Emanuele Munarini, Feb 03 2014 */
    
  • PARI
    /* compare to A061095 */
    mnom(v)=
    /* Multinomial coefficient s! / prod(j=1, n, v[j]!) where
      s= sum(j=1, n, v[j]) and n is the number of elements in v[]. */
    sum(j=1, #v, v[j])! / prod(j=1, #v, v[j]!)
    A038041(n)={local(r=0);fordiv(n,d,r+=mnom(vector(d,j,n/d))/d!);return(r);}
    vector(33,n,A038041(n)) /* Joerg Arndt, Apr 16 2011 */
    
  • Python
    import math
    def a(n):
        count = 0
        for k in range(1, n + 1):
            if n % k == 0:
                count += math.factorial(n) // (math.factorial(k) ** (n // k) * math.factorial(n // k))
        return count # Paul Muljadi, Sep 25 2024

Formula

a(n) = Sum_{d divides n} (n!/(d!*((n/d)!)^d)).
E.g.f.: Sum_{k >= 1} (exp(x^k/k!)-1).

Extensions

More terms from Erich Friedman

A321455 Number of ways to factor n into factors > 1 all having the same sum of prime indices.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 2, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 2, 1, 2, 1, 1, 1, 3, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 3, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 4, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 3
Offset: 1

Views

Author

Gus Wiseman, Nov 10 2018

Keywords

Comments

Also the number of multiset partitions of the multiset of prime indices of n with equal block-sums.
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 sum of prime indices of n is A056239(n).

Examples

			The a(1440) = 6 factorizations into factors all having the same sum of prime indices:
  (10*12*12)
  (5*6*6*8)
  (9*10*16)
  (30*48)
  (36*40)
  (1440)
The a(900) = 5 multiset partitions with equal block-sums:
  {{1,1,2,2,3,3}}
  {{3,3},{1,1,2,2}}
  {{1,2,3},{1,2,3}}
  {{1,3},{1,3},{2,2}}
  {{3},{3},{1,2},{1,2}}
		

Crossrefs

Positions of 1's are A321453. Positions of terms > 1 are A321454.

Programs

  • Mathematica
    hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
    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],SameQ@@hwt/@#&]],{n,100}]
  • PARI
    A056239(n) = if(1==n, 0, my(f=factor(n)); sum(i=1, #f~, f[i, 2] * primepi(f[i, 1])));
    all_have_same_sum_of_pis(facs) = if(!#facs, 1, (#Set(apply(A056239,facs)) == 1));
    A321455(n, m=n, facs=List([])) = if(1==n, all_have_same_sum_of_pis(facs), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s += A321455(n/d, d, newfacs))); (s)); \\ Antti Karttunen, Jan 20 2025

Extensions

Data section extended to a(108) by Antti Karttunen, Jan 20 2025

A321469 Number of factorizations of n into factors > 1 with different sums of prime indices. Number of multiset partitions of the multiset of prime indices of n with distinct block-sums.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Nov 11 2018

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 sum of prime indices of n is A056239(n).

Examples

			The a(72) = 8 multiset partitions with distinct block-sums:
    {{1,1,1,2,2}}
   {{1},{1,1,2,2}}
   {{2},{1,1,1,2}}
   {{1,1},{1,2,2}}
   {{1,2},{1,1,2}}
   {{2,2},{1,1,1}}
  {{1},{2},{1,1,2}}
  {{1},{1,1},{2,2}}
Missing from this list are:
    {{1},{1},{1,2,2}}
    {{1},{1,2},{1,2}}
    {{2},{2},{1,1,1}}
    {{2},{1,1},{1,2}}
   {{1},{1},{1},{2,2}}
   {{1},{1},{2},{1,2}}
   {{1},{2},{2},{1,1}}
  {{1},{1},{1},{2},{2}}
		

Crossrefs

Programs

  • Mathematica
    primeMS[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]]]];
    Table[Length[Select[mps[primeMS[n]],UnsameQ@@Sort[Total/@#]&]],{n,100}]
  • PARI
    A056239(n) = if(1==n, 0, my(f=factor(n)); sum(i=1, #f~, f[i, 2] * primepi(f[i, 1])));
    all_have_different_sum_of_pis(facs) = if(!#facs, 1, (#Set(apply(A056239,facs)) == #facs));
    A321469(n, m=n, facs=List([])) = if(1==n, all_have_different_sum_of_pis(facs), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s += A321469(n/d, d, newfacs))); (s)); \\ Antti Karttunen, Jan 20 2025

Extensions

Data section extended to a(105) by Antti Karttunen, Jan 20 2025

A320324 Numbers of which each prime index has the same number of prime factors, counted with multiplicity.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 13, 15, 16, 17, 19, 23, 25, 27, 29, 31, 32, 33, 37, 41, 43, 45, 47, 49, 51, 53, 55, 59, 61, 64, 67, 71, 73, 75, 79, 81, 83, 85, 89, 91, 93, 97, 99, 101, 103, 107, 109, 113, 121, 123, 125, 127, 128, 131, 135, 137, 139, 149, 151, 153
Offset: 1

Views

Author

Gus Wiseman, Oct 10 2018

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n.

Examples

			The terms together with their corresponding multiset multisystems (A302242):
   1: {}
   2: {{}}
   3: {{1}}
   4: {{},{}}
   5: {{2}}
   7: {{1,1}}
   8: {{},{},{}}
   9: {{1},{1}}
  11: {{3}}
  13: {{1,2}}
  15: {{1},{2}}
  16: {{},{},{},{}}
  17: {{4}}
  19: {{1,1,1}}
  23: {{2,2}}
  25: {{2},{2}}
  27: {{1},{1},{1}}
  29: {{1,3}}
  31: {{5}}
  32: {{},{},{},{},{}}
  33: {{1},{3}}
  37: {{1,1,2}}
  41: {{6}}
  43: {{1,4}}
  45: {{1},{1},{2}}
  47: {{2,3}}
  49: {{1,1},{1,1}}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[100],SameQ@@PrimeOmega/@PrimePi/@First/@FactorInteger[#]&]
  • PARI
    is(n) = #Set(apply(p -> bigomega(primepi(p)), factor(n)[,1]~))<=1 \\ Rémy Sigrist, Oct 11 2018

A326028 Number of factorizations of n into factors > 1 with integer geometric mean.

Original entry on oeis.org

0, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2
Offset: 1

Views

Author

Gus Wiseman, Jul 15 2019

Keywords

Comments

First differs from A294336 and A316782 at a(36) = 5.

Examples

			The a(4) = 2 through a(36) = 5 factorizations (showing only the cases where n is a perfect power).
  (4)    (8)      (9)    (16)       (25)   (27)     (32)         (36)
  (2*2)  (2*2*2)  (3*3)  (2*8)      (5*5)  (3*3*3)  (2*2*2*2*2)  (4*9)
                         (4*4)                                   (6*6)
                         (2*2*2*2)                               (2*18)
                                                                 (3*12)
		

Crossrefs

Positions of terms > 1 are the perfect powers A001597.
Partitions with integer geometric mean are A067539.
Subsets with integer geometric mean are A326027.
Factorizations with integer average and geometric mean are A326647.

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],IntegerQ[GeometricMean[#]]&]],{n,2,100}]
  • PARI
    A326028(n, m=n, facmul=1, facnum=0) = if(1==n,facnum>0 && ispower(facmul,facnum), my(s=0); fordiv(n, d, if((d>1)&&(d<=m), s += A326028(n/d, d, facmul*d, facnum+1))); (s)); \\ Antti Karttunen, Nov 10 2024

Formula

a(2^n) = A067538(n).

Extensions

a(89) onwards from Antti Karttunen, Nov 10 2024

A317583 Number of multiset partitions of normal multisets of size n such that all blocks have the same size.

Original entry on oeis.org

1, 4, 8, 30, 32, 342, 128, 3754, 11360, 56138, 2048, 3834670, 8192, 27528494, 577439424, 2681075210, 131072, 238060300946, 524288, 11045144602614, 115488471132032, 49840258213638, 8388608, 152185891301461434, 140102945910265344, 124260001149229146, 85092642310351607968
Offset: 1

Views

Author

Gus Wiseman, Aug 01 2018

Keywords

Comments

A multiset is normal if it spans an initial interval of positive integers.
a(n) is the number of nonnegative integer matrices with total sum n, nonzero rows and each column with the same sum with columns in nonincreasing lexicographic order. - Andrew Howroyd, Jan 15 2020

Examples

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

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]]]];
    allnorm[n_]:=Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1];
    Table[Length[Select[Join@@mps/@allnorm[n],SameQ@@Length/@#&]],{n,8}]
  • PARI
    \\ here U(n,m) gives number for m blocks of size n.
    U(n,m)={sum(k=1, n*m, binomial(binomial(k+n-1, n)+m-1, m)*sum(r=k, n*m, binomial(r, k)*(-1)^(r-k)) )}
    a(n)={sumdiv(n, d, U(d, n/d))} \\ Andrew Howroyd, Sep 15 2018

Formula

a(p) = 2^p for prime p. - Andrew Howroyd, Sep 15 2018
a(n) = Sum_{d|n} A331315(n/d, d). - Andrew Howroyd, Jan 15 2020

Extensions

Terms a(9) and beyond from Andrew Howroyd, Sep 15 2018

A326515 Number of factorizations of n into factors > 1 where every factor has the same average of prime indices.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jul 12 2019

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.

Examples

			The a(900) = 9 factorizations:
  (3*3*10*10),
  (3*3*100), (3*10*30), (9*10*10),
  (3*300), (9*100), (10*90), (30*30),
  (900).
		

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]]}]];
    Table[Length[Select[facs[n],SameQ@@Mean/@primeMS/@#&]],{n,100}]
  • PARI
    avgpis(n) = { my(f=factor(n)); f[,1] = apply(primepi,f[,1]); (1/bigomega(n))*sum(i=1,#f~,f[i,2]*f[i,1]); };
    has_same_average_of_pis(facs) = if(!#facs, 1, my(avg=0); for(i=1,#facs,if(!avg, avg=avgpis(facs[i]), if(avg!=avgpis(facs[i]), return(0)))); (1));
    A326515(n, m=n, facs=List([])) = if(1==n, has_same_average_of_pis(facs), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s += A326515(n/d, d, newfacs))); (s)); \\ Antti Karttunen, Jan 20 2025

Extensions

Data section extended to a(105) by Antti Karttunen, Jan 20 2025

A326516 Number of factorizations of n into factors > 1 where each factor has a different average of prime indices.

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

Views

Author

Gus Wiseman, Jul 12 2019

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.

Examples

			The a(60) = 8 factorizations: (2*5*6), (3*4*5), (2*30), (3*20), (4*15), (5*12), (6*10), (60).
		

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]]}]];
    Table[Length[Select[facs[n],UnsameQ@@Mean/@primeMS/@#&]],{n,100}]
  • PARI
    avgpis(n) = { my(f=factor(n)); f[,1] = apply(primepi,f[,1]); (1/bigomega(n))*sum(i=1,#f~,f[i,2]*f[i,1]); };
    all_have_different_average_of_pis(facs) = if(!#facs, 1, (#Set(apply(avgpis,facs)) == #facs));
    A326516(n, m=n, facs=List([])) = if(1==n, all_have_different_average_of_pis(facs), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s += A326516(n/d, d, newfacs))); (s)); \\ Antti Karttunen, Jan 20 2025

Extensions

Data section extended to a(105) by Antti Karttunen, Jan 20 2025

A326647 Number of factorizations of n into factors > 1 with integer average and integer geometric mean.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Jul 16 2019

Keywords

Examples

			The a(216) = 5 factorizations:
  (2*4*27)
  (3*3*24)
  (3*6*12)
  (6*6*6)
  (216)
The a(729) = 8 factorizations:
  (3*3*3*3*3*3)
  (3*3*81)
  (3*9*27)
  (3*243)
  (9*9*9)
  (9*81)
  (27*27)
  (729)
		

Crossrefs

Positions of terms > 1 are the perfect powers A001597.
Factorizations with integer average are A326622.
Factorizations with integer geometric mean are A326028.
Partitions with integer average and geometric mean are A326641.

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],IntegerQ[Mean[#]]&&IntegerQ[GeometricMean[#]]&]],{n,2,100}]

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}]
Showing 1-10 of 18 results. Next