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 30 results. Next

A064576 First differences of A064573, where A064573(n) is the number of partitions of n into parts which are all powers of the same prime.

Original entry on oeis.org

1, 1, 2, 1, 3, 1, 4, 2, 5, 1, 8, 1, 7, 3, 10, 1, 13, 1, 15, 4, 15, 1, 23, 2, 21, 5, 27, 1, 33, 1, 36, 6, 37, 3, 53, 1, 47, 8, 62, 1, 68, 1, 75, 11, 75, 1, 103, 2, 97, 10, 115, 1, 126, 4, 142, 13, 141, 1, 181, 1, 167, 17, 202, 4, 218, 1, 239, 16, 243, 1, 302, 1, 285, 22, 331, 3, 349
Offset: 1

Views

Author

Marc LeBrun, Sep 20 2001

Keywords

Comments

Apparently a(n)=1 when n+1 is prime.

Crossrefs

Programs

A306017 Number of non-isomorphic multiset partitions of weight n in which all parts have the same size.

Original entry on oeis.org

1, 1, 4, 6, 17, 14, 66, 30, 189, 222, 550, 112, 4696, 202, 5612, 30914, 63219, 594, 453125, 980, 3602695, 5914580, 1169348, 2510, 299083307, 232988061, 23248212, 2669116433, 14829762423, 9130, 170677509317, 13684, 1724710753084, 2199418340875, 14184712185, 38316098104262
Offset: 0

Views

Author

Gus Wiseman, Jun 17 2018

Keywords

Comments

A multiset partition of weight n is a finite multiset of finite nonempty multisets whose sizes sum to n.
Number of distinct nonnegative integer matrices with all row sums equal and total sum n up to row and column permutations. - Andrew Howroyd, Sep 05 2018
From Gus Wiseman, Oct 11 2018: (Start)
Also the number of non-isomorphic multiset partitions of weight n in which each vertex appears the same number of times. For n = 4, non-isomorphic representatives of these 17 multiset partitions are:
{{1,1,1,1}}
{{1,1,2,2}}
{{1,2,3,4}}
{{1},{1,1,1}}
{{1},{1,2,2}}
{{1},{2,3,4}}
{{1,1},{1,1}}
{{1,1},{2,2}}
{{1,2},{1,2}}
{{1,2},{3,4}}
{{1},{1},{1,1}}
{{1},{1},{2,2}}
{{1},{2},{1,2}}
{{1},{2},{3,4}}
{{1},{1},{1},{1}}
{{1},{1},{2},{2}}
{{1},{2},{3},{4}}
(End)

Examples

			Non-isomorphic representatives of the a(4) = 17 multiset partitions:
  {{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}}
  {{1,2},{1,2}}
  {{1,2},{2,2}}
  {{1,2},{3,3}}
  {{1,2},{3,4}}
  {{1,3},{2,3}}
  {{1},{1},{1},{1}}
  {{1},{1},{2},{2}}
  {{1},{2},{2},{2}}
  {{1},{2},{3},{3}}
  {{1},{2},{3},{4}}
		

Crossrefs

Programs

  • Mathematica
    permcount[v_List] := Module[{m = 1, s = 0, k = 0, t}, For[i = 1, i <= Length[v], i++, t = v[[i]]; k = If[i > 1 && t == v[[i - 1]], k + 1, 1]; m *= t*k; s += t]; s!/m];
    K[q_List, t_, k_] := SeriesCoefficient[1/Product[g = GCD[t, q[[j]]]; (1 - x^(q[[j]]/g))^g, {j, 1, Length[q]}], {x, 0, k}];
    RowSumMats[n_, m_, k_] := Module[{s = 0}, Do[s += permcount[q]* SeriesCoefficient[Exp[Sum[K[q, t, k]/t*x^t, {t, 1, n}]], {x, 0, n}], {q, IntegerPartitions[m]}]; s/m!];
    a[n_] := a[n] = If[n==0, 1, If[PrimeQ[n], 2 PartitionsP[n], Sum[ RowSumMats[ n/d, n, d], {d, Divisors[n]}]]];
    Table[Print[n, " ", a[n]]; a[n], {n, 0, 35}] (* Jean-François Alcover, Nov 07 2019, after Andrew Howroyd *)
  • PARI
    \\ See A318951 for RowSumMats.
    a(n)={sumdiv(n,d,RowSumMats(n/d,n,d))} \\ Andrew Howroyd, Sep 05 2018

Formula

For p prime, a(p) = 2*A000041(p).
a(n) = Sum_{d|n} A331485(n/d, d). - Andrew Howroyd, Feb 09 2020

Extensions

Terms a(11) and beyond from Andrew Howroyd, Sep 05 2018

A319056 Number of non-isomorphic multiset partitions of weight n in which (1) all parts have the same size and (2) each vertex appears the same number of times.

Original entry on oeis.org

1, 1, 4, 4, 10, 4, 21, 4, 26, 13, 28, 4, 128, 4, 39, 84, 150, 4, 358, 4, 956, 513, 86, 4, 12549, 1864, 134, 9582, 52366, 4, 301086, 4, 1042038, 407140, 336, 4690369, 61738312, 4, 532, 28011397, 2674943885, 4, 819150246, 4, 54904825372, 65666759973, 1303, 4, 4319823776760
Offset: 0

Views

Author

Gus Wiseman, Oct 10 2018

Keywords

Comments

a(p) = 4 for p prime. - Charlie Neder, Oct 15 2018

Examples

			Non-isomorphic representatives of the a(1) = 1 through a(6) = 21 multiset partitions:
  (1)  (11)    (111)      (1111)        (11111)          (111111)
       (12)    (123)      (1122)        (12345)          (111222)
       (1)(1)  (1)(1)(1)  (1234)        (1)(1)(1)(1)(1)  (112233)
       (1)(2)  (1)(2)(3)  (11)(11)      (1)(2)(3)(4)(5)  (123456)
                          (11)(22)                       (111)(111)
                          (12)(12)                       (111)(222)
                          (12)(34)                       (112)(122)
                          (1)(1)(1)(1)                   (112)(233)
                          (1)(1)(2)(2)                   (123)(123)
                          (1)(2)(3)(4)                   (123)(456)
                                                         (11)(11)(11)
                                                         (11)(12)(22)
                                                         (11)(22)(33)
                                                         (11)(23)(23)
                                                         (12)(12)(12)
                                                         (12)(13)(23)
                                                         (12)(34)(56)
                                                         (1)(1)(1)(1)(1)(1)
                                                         (1)(1)(1)(2)(2)(2)
                                                         (1)(1)(2)(2)(3)(3)
                                                         (1)(2)(3)(4)(5)(6)
		

Crossrefs

Extensions

Terms a(12) and beyond from Andrew Howroyd, Feb 03 2022

A320322 Number of integer partitions of n whose product is a perfect power.

Original entry on oeis.org

1, 0, 0, 0, 2, 2, 5, 5, 9, 11, 18, 19, 28, 30, 42, 50, 68, 76, 102, 113, 146, 170, 212, 241, 312, 356, 441, 514, 628, 720, 887, 1008, 1215, 1403, 1660, 1903, 2291, 2609, 3107, 3594, 4254, 4864, 5739, 6546, 7672, 8811, 10237, 11651, 13583, 15420, 17867, 20382
Offset: 0

Views

Author

Gus Wiseman, Oct 10 2018

Keywords

Examples

			The a(4) = 2 through a(11) = 19 integer partitions:
  4   41   33    331    8       9        55        551
  22  221  42    421    44      81       82        632
           222   2221   422     333      91        821
           411   4111   2222    441      433       911
           2211  22111  3311    4221     442       4331
                        4211    22221    811       4421
                        22211   33111    3322      8111
                        41111   42111    3331      33221
                        221111  222111   4222      33311
                                411111   4411      42221
                                2211111  22222     44111
                                         42211     222221
                                         222211    422111
                                         331111    2222111
                                         421111    3311111
                                         2221111   4211111
                                         4111111   22211111
                                         22111111  41111111
                                                   221111111
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],GCD@@FactorInteger[Times@@#][[All,2]]>1&]],{n,30}]

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

A325039 Number of integer partitions of n with the same product of parts as their conjugate.

Original entry on oeis.org

1, 1, 0, 1, 1, 1, 1, 1, 6, 2, 2, 4, 3, 5, 7, 6, 5, 7, 9, 10, 11, 18, 16, 19, 19, 16, 20, 20, 28, 39, 28, 40, 53, 45, 52, 59, 71, 61, 73, 97, 102, 95, 112, 131, 137, 148, 140, 166, 199, 181, 238, 251, 255, 289, 339, 344, 381, 398, 422, 464, 541, 555, 628, 677, 732
Offset: 0

Views

Author

Gus Wiseman, Mar 25 2019

Keywords

Comments

For example, the partition (6,4,1) with product 24 has conjugate (3,2,2,2,1,1) with product also 24.
The Heinz numbers of these partitions are given by A325040.

Examples

			The a(8) = 6 through a(15) = 6 integer partitions:
  (44)    (333)    (4321)   (641)     (4422)    (4432)     (6431)
  (332)   (51111)  (52111)  (4331)    (53211)   (6421)     (8411)
  (431)                     (322211)  (621111)  (53311)    (54221)
  (2222)                    (611111)            (432211)   (433211)
  (3221)                                        (7111111)  (632111)
  (4211)                                                   (7211111)
                                                           (42221111)
		

Crossrefs

Programs

  • Mathematica
    conj[y_]:=If[Length[y]==0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]];
    Table[Length[Select[IntegerPartitions[n],Times@@#==Times@@conj[#]&]],{n,0,30}]

Extensions

More terms from Jinyuan Wang, Jun 27 2020

A322452 Number of factorizations of n into factors > 1 not including any prime powers.

Original entry on oeis.org

1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 1, 1, 1, 2, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 2, 0, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 2, 0, 1, 1, 1, 1, 1, 0, 1, 0, 1, 0, 2, 1, 1, 1, 1, 0, 2, 1, 1, 1, 1, 1, 1, 0, 1, 1, 2, 0, 1, 0, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Dec 09 2018

Keywords

Comments

Also the number of multiset partitions of the multiset of prime indices of n with no constant parts.

Examples

			The a(840) = 11 factorizations are (6*10*14), (6*140), (10*84), (12*70), (14*60), (15*56), (20*42), (21*40), (24*35), (28*30), (840).
		

Crossrefs

Positions of 0's are the prime powers A000961.

Programs

  • Mathematica
    acfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[acfacs[n/d],Min@@#>=d&]],{d,Select[Rest[Divisors[n]],!PrimePowerQ[#]&]}]];
    Table[Length[acfacs[n]],{n,100}]
  • PARI
    A322452(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1)&&(d<=m)&&(1A322452(n/d, d))); (s)); \\ Antti Karttunen, Jan 03 2019
    
  • PARI
    first(n) = my(res=vector(n)); for(i=1, n, f=factor(i); v=vecsort(f[,2] , , 4); f[, 2] = v; fb = factorback(f); if(fb==i, res[i] = A322452(i), res[i] = res[fb])); res \\ A322452 the function above \\ David A. Corneth, Jan 03 2019

Extensions

More terms from Antti Karttunen, Jan 03 2019

A319899 Numbers whose number of prime factors with multiplicity (A001222) is the number of distinct prime factors (A001221) in the product of the prime indices (A003963).

Original entry on oeis.org

1, 3, 5, 7, 11, 15, 17, 19, 23, 26, 31, 33, 35, 39, 41, 51, 53, 55, 58, 59, 65, 67, 69, 74, 77, 83, 85, 86, 87, 91, 93, 94, 95, 97, 103, 109, 111, 119, 122, 123, 127, 129, 131, 142, 146, 155, 157, 158, 161, 165, 169, 177, 178, 179, 183, 185, 187, 191, 201, 202
Offset: 1

Views

Author

Gus Wiseman, Dec 17 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 multiset multisystem with MM-number n is formed by taking the multiset of prime indices of each part of the multiset of prime indices 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}}. This sequence lists all MM-numbers of square multiset multisystems, meaning the number of edges is equal to the number of distinct vertices.

Examples

			The sequence of multiset multisystems whose MM-numbers belong to the sequence begins:
   1: {}
   3: {{1}}
   5: {{2}}
   7: {{1,1}}
  11: {{3}}
  15: {{1},{2}}
  17: {{4}}
  19: {{1,1,1}}
  23: {{2,2}}
  26: {{},{1,2}}
  31: {{5}}
  33: {{1},{3}}
  35: {{2},{1,1}}
  39: {{1},{1,2}}
  41: {{6}}
  51: {{1},{4}}
  53: {{1,1,1,1}}
  55: {{2},{3}}
  58: {{},{1,3}}
  59: {{7}}
  65: {{2},{1,2}}
  67: {{8}}
  69: {{1},{2,2}}
  74: {{},{1,1,2}}
		

Crossrefs

Programs

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

A358911 Number of integer compositions 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, 4, 4, 7, 9, 12, 20, 21, 39, 49, 79, 109, 161, 236, 345, 512, 752, 1092, 1628, 2376, 3537, 5171, 7650, 11266, 16634, 24537, 36173, 53377, 78791, 116224, 171598, 253109, 373715, 551434, 814066, 1201466, 1773425, 2617744, 3864050, 5703840, 8419699
Offset: 0

Views

Author

Gus Wiseman, Dec 11 2022

Keywords

Examples

			The a(1) = 1 through a(8) = 9 compositions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)        (8)
       (11)  (111)  (22)    (23)     (33)      (25)       (35)
                    (1111)  (32)     (222)     (52)       (44)
                            (11111)  (111111)  (223)      (53)
                                               (232)      (233)
                                               (322)      (323)
                                               (1111111)  (332)
                                                          (2222)
                                                          (11111111)
		

Crossrefs

The case of partitions is A319169, ranked by A320324.
The weakly decreasing version is A358335, strictly A358901.
For sequences of partitions see A358905.
A001222 counts prime factors, distinct A001221.
A011782 counts compositions.
A358902 = compositions with weakly decreasing A001221, strictly A358903.
A358909 = partitions with weakly decreasing A001222, complement A358910.

Programs

  • Maple
    b:= proc(n, i) option remember; uses numtheory; `if`(n=0, 1, add(
         (t-> `if`(i<0 or i=t, b(n-j, t), 0))(bigomega(j)), j=1..n))
        end:
    a:= n-> b(n, -1):
    seq(a(n), n=0..44);  # Alois P. Heinz, Feb 12 2024
  • Mathematica
    Table[Length[Select[Join @@ Permutations/@IntegerPartitions[n],SameQ@@PrimeOmega/@#&]],{n,0,10}]

Extensions

a(21) and beyond from Lucas A. Brown, Dec 15 2022
Showing 1-10 of 30 results. Next