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 32 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

A326534 MM-numbers of multiset partitions where every part has the same sum.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 13, 16, 17, 19, 23, 25, 27, 29, 31, 32, 35, 37, 41, 43, 47, 49, 53, 59, 61, 64, 67, 71, 73, 79, 81, 83, 89, 97, 101, 103, 107, 109, 113, 121, 125, 127, 128, 131, 137, 139, 143, 149, 151, 157, 163, 167, 169, 173, 175, 179, 181, 191
Offset: 1

Views

Author

Gus Wiseman, Jul 12 2019

Keywords

Comments

First differs from A298538 in lacking 187.
These are numbers where each prime index has the same sum of prime indices. 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 multiset multisystem with MM-number n is obtained by taking the multiset of prime indices of each prime index of n. For example, the prime indices of 78 are {1,2,6}, so the multiset multisystem with MM-number 78 is {{},{1},{1,2}}.

Examples

			The sequence of multiset partitions where every part has the same sum, preceded by their MM-numbers, begins:
   1: {}
   2: {{}}
   3: {{1}}
   4: {{},{}}
   5: {{2}}
   7: {{1,1}}
   8: {{},{},{}}
   9: {{1},{1}}
  11: {{3}}
  13: {{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: {{},{},{},{},{}}
  35: {{2},{1,1}}
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],SameQ@@Total/@primeMS/@primeMS[#]&]

A326535 MM-numbers of multiset partitions where each part has a different sum.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 26, 29, 30, 31, 33, 34, 37, 38, 39, 41, 42, 43, 46, 47, 51, 53, 55, 57, 58, 59, 61, 62, 65, 66, 67, 69, 71, 73, 74, 77, 78, 79, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97, 101, 102, 103, 106, 107, 109
Offset: 1

Views

Author

Gus Wiseman, Jul 12 2019

Keywords

Comments

First differs from A298540 in having 187.
These are numbers where each prime index has a different sum of prime indices. 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 multiset multisystem with MM-number n is obtained by taking the multiset of prime indices of each prime index of n. For example, the prime indices of 78 are {1,2,6}, so the multiset multisystem with MM-number 78 is {{},{1},{1,2}}.

Examples

			The sequence of multiset partitions where each part has a different sum, preceded by their MM-numbers, begins:
   1: {}
   2: {{}}
   3: {{1}}
   5: {{2}}
   6: {{},{1}}
   7: {{1,1}}
  10: {{},{2}}
  11: {{3}}
  13: {{1,2}}
  14: {{},{1,1}}
  15: {{1},{2}}
  17: {{4}}
  19: {{1,1,1}}
  21: {{1},{1,1}}
  22: {{},{3}}
  23: {{2,2}}
  26: {{},{1,2}}
  29: {{1,3}}
  30: {{},{1},{2}}
  31: {{5}}
		

Crossrefs

Programs

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

A320325 Numbers whose product of prime indices is a perfect power.

Original entry on oeis.org

7, 9, 14, 18, 19, 21, 23, 25, 27, 28, 36, 38, 42, 46, 49, 50, 53, 54, 56, 57, 63, 72, 76, 81, 84, 92, 97, 98, 100, 103, 106, 108, 112, 114, 115, 121, 125, 126, 131, 133, 144, 147, 151, 152, 159, 161, 162, 168, 169, 171, 175, 183, 184, 185, 189, 194, 195, 196
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. The multiset of prime indices of n is row n of A112798.

Examples

			The terms together with their corresponding multiset multisystems (A302242):
   7: {{1,1}}
   9: {{1},{1}}
  14: {{},{1,1}}
  18: {{},{1},{1}}
  19: {{1,1,1}}
  21: {{1},{1,1}}
  23: {{2,2}}
  25: {{2},{2}}
  27: {{1},{1},{1}}
  28: {{},{},{1,1}}
  36: {{},{},{1},{1}}
  38: {{},{1,1,1}}
  42: {{},{1},{1,1}}
  46: {{},{2,2}}
  49: {{1,1},{1,1}}
  50: {{},{2},{2}}
  53: {{1,1,1,1}}
  54: {{},{1},{1},{1}}
  56: {{},{},{},{1,1}}
  57: {{1},{1,1,1}}
  63: {{1},{1},{1,1}}
  72: {{},{},{},{1},{1}}
  76: {{},{},{1,1,1}}
  81: {{1},{1},{1},{1}}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[100],GCD@@FactorInteger[Times@@Cases[FactorInteger[#],{p_,k_}:>PrimePi[p]^k]][[All,2]]>1&]

A319169 Number of integer partitions of n whose parts all have the same number of prime factors, counted with multiplicity.

Original entry on oeis.org

1, 1, 2, 2, 3, 3, 4, 4, 6, 6, 8, 7, 11, 11, 14, 15, 20, 19, 26, 27, 34, 35, 43, 45, 59, 60, 72, 77, 94, 98, 118, 125, 148, 158, 184, 198, 233, 245, 282, 308, 353, 374, 428, 464, 525, 566, 635, 686, 779, 832, 930, 1005, 1123, 1208, 1345, 1451, 1609, 1732, 1912
Offset: 0

Views

Author

Gus Wiseman, Oct 10 2018

Keywords

Examples

			The a(1) = 1 through a(9) = 6 integer partitions:
  1  2   3    4     5      6       7        8         9
     11  111  22    32     33      52       44        72
              1111  11111  222     322      53        333
                           111111  1111111  332       522
                                            2222      3222
                                            11111111  111111111
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i, f) option remember; `if`(n=0, 1, `if`(i<1, 0,
          b(n, i-1, f)+(o-> `if`(f in {0, o}, b(n-i, min(i, n-i),
         `if`(f=0, o, f)), 0))(numtheory[bigomega](i))))
        end:
    a:= n-> b(n$2, 0):
    seq(a(n), n=0..75);  # Alois P. Heinz, Dec 15 2018
  • Mathematica
    Table[Length[Select[IntegerPartitions[n],SameQ@@PrimeOmega/@#&]],{n,30}]
    (* Second program: *)
    b[n_, i_, f_] := b[n, i, f] = If[n == 0, 1, If[i < 1, 0,
         b[n, i-1, f] + Function[o, If[f == 0 || f == o, b[n-i, Min[i, n-i],
         If[f == 0, o, f]], 0]][PrimeOmega[i]]]];
    a[n_] := b[n, n, 0];
    a /@ Range[0, 75] (* Jean-François Alcover, May 10 2021, after Alois P. Heinz *)

Extensions

a(51)-a(58) from Chai Wah Wu, Nov 12 2018

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

A375128 Irregular triangle read by rows where row n lists the minima of maximal strictly increasing runs in the weakly increasing prime indices of n.

Original entry on oeis.org

1, 2, 1, 1, 3, 1, 4, 1, 1, 1, 2, 2, 1, 5, 1, 1, 6, 1, 2, 1, 1, 1, 1, 7, 1, 2, 8, 1, 1, 2, 1, 9, 1, 1, 1, 3, 3, 1, 2, 2, 2, 1, 1, 10, 1, 11, 1, 1, 1, 1, 1, 2, 1, 3, 1, 1, 2, 12, 1, 2, 1, 1, 1, 13, 1, 14, 1, 1, 2, 2, 1, 15, 1, 1, 1, 1, 4, 4, 1, 3, 2, 1, 1, 16
Offset: 1

Views

Author

Gus Wiseman, Aug 04 2024

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 minima of strictly increasing runs in a sequence are obtained by splitting it into maximal strictly increasing subsequences and taking the first term of each.

Examples

			The prime indices of 540 are {1,1,2,2,2,3}, with strictly increasing runs ({1},{1,2},{2},{2,3}), with minima (1,1,2,2), which is row 540.
Triangle begins:
   1:
   2:  1
   3:  2
   4:  1  1
   5:  3
   6:  1
   7:  4
   8:  1  1  1
   9:  2  2
  10:  1
  11:  5
  12:  1  1
  13:  6
  14:  1
  15:  2
  16:  1  1  1  1
		

Crossrefs

Row-minima are A055396.
Row-sums are A374706.
Row-lengths are A375136.
For leaders of constant runs we have A304038, row-sums A066328.
For compositions we have A374683, row-sums of A374684 (length A124768).
A112798 lists prime indices:
- length A001222, distinct A001221
- leader A055396
- sum A056239
- reverse A296150

Programs

  • Mathematica
    Table[If[n==1,{},First/@Split[Flatten[Cases[FactorInteger[n], {p_,k_}:>Table[PrimePi[p],{k}]]],Less]],{n,100}]

A322794 Number of factorizations of n into factors > 1 where all factors have the same number of prime factors counted with multiplicity.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Dec 26 2018

Keywords

Comments

Also the number of uniform multiset partitions of the multiset of prime indices of n, where a multiset partition is uniform if all parts have the same size.

Examples

			The a(1260) = 13 factorizations:
  (1260)  (18*70)   (4*9*35)   (2*2*3*3*5*7)
          (20*63)   (6*6*35)
          (28*45)   (4*15*21)
          (30*42)   (6*10*21)
          (12*105)  (6*14*15)
                    (9*10*14)
The a(1260) = 13 multiset partitions:
  {{1},{1},{2},{2},{3},{4}}
     {{1,1},{2,2},{3,4}}
     {{1,1},{2,3},{2,4}}
     {{1,2},{1,2},{3,4}}
     {{1,2},{1,3},{2,4}}
     {{1,2},{1,4},{2,3}}
     {{2,2},{1,3},{1,4}}
      {{1,1,2},{2,3,4}}
      {{1,2,2},{1,3,4}}
      {{1,1,3},{2,2,4}}
      {{1,1,4},{2,2,3}}
      {{1,2,3},{1,2,4}}
       {{1,1,2,2,3,4}}
		

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],SameQ@@PrimeOmega/@#&]],{n,100}]

A375136 Number of maximal strictly increasing runs in the weakly increasing prime factors of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Aug 04 2024

Keywords

Comments

For n > 1, this is one more than the number of adjacent equal terms in the multiset of prime factors of n.

Examples

			The prime factors of 540 are {2,2,3,3,3,5}, with maximal strictly increasing runs ({2},{2,3},{3},{3,5}), so a(540) = 4.
		

Crossrefs

For compositions we have A124768, row-lengths of A374683, sum A374684.
For sum of prime indices we have A374706.
Row-lengths of A375128.
A112798 lists prime indices:
- distinct A001221
- length A001222
- leader A055396
- sum A056239
- reverse A296150

Programs

  • Mathematica
    Table[Length[Split[Flatten[ConstantArray@@@FactorInteger[n]],Less]],{n,100}]

Formula

For n > 1, a(n) = A046660(n) + 1 = A001222(n) - A001221(n) + 1.

A326533 MM-numbers of multiset partitions where each part has a different length.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 13, 14, 17, 19, 21, 22, 23, 26, 29, 31, 34, 35, 37, 38, 39, 41, 42, 43, 46, 47, 53, 57, 58, 59, 61, 62, 65, 67, 69, 70, 71, 73, 74, 77, 78, 79, 82, 83, 86, 87, 89, 94, 95, 97, 101, 103, 106, 107, 109, 111, 113, 114, 115, 118, 119, 122
Offset: 1

Views

Author

Gus Wiseman, Jul 12 2019

Keywords

Comments

These are numbers where each prime index has a different Omega (number of prime factors counted with multiplicity). 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 multiset multisystem with MM-number n is obtained by taking the multiset of prime indices of each prime index of n. For example, the prime indices of 78 are {1,2,6}, so the multiset multisystem with MM-number 78 is {{},{1},{1,2}}.

Examples

			The sequence of multiset partitions where each part has a different average preceded by their MM-numbers begins:
   1: {}
   2: {{}}
   3: {{1}}
   5: {{2}}
   6: {{},{1}}
   7: {{1,1}}
  10: {{},{2}}
  11: {{3}}
  13: {{1,2}}
  14: {{},{1,1}}
  17: {{4}}
  19: {{1,1,1}}
  21: {{1},{1,1}}
  22: {{},{3}}
  23: {{2,2}}
  26: {{},{1,2}}
  29: {{1,3}}
  31: {{5}}
  34: {{},{4}}
  35: {{2},{1,1}}
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],UnsameQ@@PrimeOmega/@primeMS[#]&]
Showing 1-10 of 32 results. Next