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

A101271 Number of partitions of n into 3 distinct and relatively prime parts.

Original entry on oeis.org

1, 1, 2, 3, 4, 5, 6, 8, 9, 12, 12, 16, 15, 21, 20, 26, 25, 33, 28, 40, 36, 45, 42, 56, 44, 65, 56, 70, 64, 84, 66, 96, 81, 100, 88, 120, 90, 133, 110, 132, 121, 161, 120, 175, 140, 176, 156, 208, 153, 220, 180, 222, 196, 261, 184, 280, 225, 270, 240, 312, 230, 341, 272
Offset: 6

Views

Author

Vladeta Jovovic, Dec 19 2004

Keywords

Comments

The Heinz numbers of these partitions are the intersection of A289509 (relatively prime), A005117 (strict), and A014612 (triple). - Gus Wiseman, Oct 15 2020

Examples

			For n=10 we have 4 such partitions: 1+2+7, 1+3+6, 1+4+5 and 2+3+5.
From _Gus Wiseman_, Oct 13 2020: (Start)
The a(6) = 1 through a(18) = 15 triples (A..F = 10..15):
  321  421  431  432  532  542  543  643  653  654  754  764  765
            521  531  541  632  651  652  743  753  763  854  873
                 621  631  641  732  742  752  762  853  863  954
                      721  731  741  751  761  843  871  872  972
                           821  831  832  851  852  943  953  981
                                921  841  932  861  952  962  A53
                                     931  941  942  961  971  A71
                                     A21  A31  951  A51  A43  B43
                                          B21  A32  B32  A52  B52
                                               A41  B41  A61  B61
                                               B31  C31  B42  C51
                                               C21  D21  B51  D32
                                                         C32  D41
                                                         C41  E31
                                                         D31  F21
                                                         E21
(End)
		

Crossrefs

A000741 is the ordered non-strict version.
A001399(n-6) does not require relative primality.
A023022 counts pairs instead of triples.
A023023 is the not necessarily strict version.
A078374 counts these partitions of any length, with Heinz numbers A302796.
A101271*6 is the ordered version.
A220377 is the pairwise coprime instead of relatively prime version.
A284825 counts the case that is pairwise non-coprime also.
A337605 is the pairwise non-coprime instead of relatively prime version.
A008289 counts strict partitions by sum and length.
A007304 gives the Heinz numbers of 3-part strict partitions.
A307719 counts 3-part pairwise coprime partitions.
A337601 counts 3-part partitions whose distinct parts are pairwise coprime.

Programs

  • Maple
    m:=3: with(numtheory): g:=sum(mobius(k)*x^(m*(m+1)/2*k)/Product(1-x^(i*k),i=1..m),k=1..20): gser:=series(g,x=0,80): seq(coeff(gser,x^n),n=6..77); # Emeric Deutsch, May 31 2005
  • Mathematica
    Table[Length[Select[IntegerPartitions[n,{3}],UnsameQ@@#&&GCD@@#==1&]],{n,6,50}] (* Gus Wiseman, Oct 13 2020 *)

Formula

G.f. for the number of partitions of n into m distinct and relatively prime parts is Sum(moebius(k)*x^(m*(m+1)/2*k)/Product(1-x^(i*k), i=1..m), k=1..infinity).

Extensions

More terms from Emeric Deutsch, May 31 2005

A000742 Number of compositions of n into 4 ordered relatively prime parts.

Original entry on oeis.org

1, 4, 10, 20, 34, 56, 80, 120, 154, 220, 266, 360, 420, 560, 614, 816, 884, 1120, 1210, 1540, 1572, 2020, 2080, 2544, 2638, 3276, 3200, 4060, 4040, 4840, 4896, 5960, 5710, 7140, 6954, 8216, 8136, 9880, 9244, 11480, 11010, 12824, 12650, 15180, 14024, 17276
Offset: 4

Views

Author

Keywords

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Magma
    [&+[MoebiusMu(n div d)*Binomial(d-1, 3):d in Divisors(n)]:n in[4..49]]; // Marius A. Burtea, Feb 08 2020
  • Maple
    with(numtheory):
    a:= n-> add(mobius(n/d)*binomial(d-1, 3), d=divisors(n)):
    seq(a(n), n=4..50);  # Alois P. Heinz, Feb 05 2020
  • Mathematica
    a[n_] := Sum[Boole[Divisible[n, k]] MoebiusMu[n/k] Binomial[k - 1, 3], {k, 1, n}]; Table[a[n], {n, 4, 51}] (* Jean-François Alcover, Feb 11 2016 *)

Formula

Möbius transform of C(n-1,3).
G.f.: Sum_{k>=1} mu(k) * x^(4*k) / (1 - x^k)^4. - Ilya Gutkovskiy, Feb 05 2020

Extensions

Offset changed to 4 by Ilya Gutkovskiy, Feb 05 2020

A000743 Number of compositions of n into 5 ordered relatively prime parts.

Original entry on oeis.org

1, 5, 15, 35, 70, 125, 210, 325, 495, 700, 1000, 1330, 1820, 2305, 3060, 3750, 4830, 5775, 7315, 8490, 10625, 12155, 14880, 16835, 20475, 22620, 27405, 30100, 35750, 39100, 46360, 49655, 58905, 62985, 73320, 78340, 91390, 95720, 111930, 117425
Offset: 5

Views

Author

Keywords

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

  • Magma
    [&+[MoebiusMu(n div d)*Binomial(d-1, 4):d in Divisors(n)]:n in[5..44]]; // Marius A. Burtea, Feb 08 2020
  • Maple
    with(numtheory):
    a:= n-> add(mobius(n/d)*binomial(d-1, 4), d=divisors(n)):
    seq(a(n), n=5..50);  # Alois P. Heinz, Feb 05 2020
  • Mathematica
    a[n_] := Sum[Boole[Divisible[n, k]] MoebiusMu[n/k] Binomial[k - 1, 4], {k, 1, n}]; Table[a[n], {n, 5, 52}] (* Jean-François Alcover, Feb 11 2016 *)

Formula

Möbius transform of binomial(n-1, 4).
G.f.: Sum_{k>=1} mu(k) * x^(5*k) / (1 - x^k)^5. - Ilya Gutkovskiy, Feb 05 2020

Extensions

Offset changed to 5 by Ilya Gutkovskiy, Feb 05 2020

A023032 Number of compositions of n into 7 ordered relatively prime parts.

Original entry on oeis.org

1, 7, 28, 84, 210, 462, 924, 1715, 3003, 4998, 8008, 12348, 18564, 27048, 38759, 54054, 74613, 100478, 134596, 176176, 230202, 294294, 376740, 471933, 593775, 731276, 905982, 1099560, 1344903, 1610322, 1947792, 2306220, 2759757, 3235484, 3838380
Offset: 7

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    [&+[MoebiusMu(n div d)*Binomial(d-1,6):d in Divisors(n)]:n in[7..41]]; // Marius A. Burtea, Feb 07 2020
  • Maple
    with(numtheory):
    a:= n-> add(mobius(n/d)*binomial(d-1, 6), d=divisors(n)):
    seq(a(n), n=7..50);  # Alois P. Heinz, Feb 05 2020
  • Mathematica
    a[n_]:=DivisorSum[n, Binomial[#-1, 6] MoebiusMu[n/#]&]; Array[a, 37, 7] (* or *) a[n_]:=Sum[Boole[Divisible[n, k]] MoebiusMu[n/k] Binomial[k-1, 6], {k, 1, n}]; Table[a[n], {n, 7, 45}] (* Vincenzo Librandi, Feb 07 2020 *)

Formula

Moebius transform of C(n-1,6).
G.f.: Sum_{k>=1} mu(k) * x^(7*k) / (1 - x^k)^7. - Ilya Gutkovskiy, Feb 05 2020

A023035 Number of compositions of n into 10 ordered relatively prime parts.

Original entry on oeis.org

1, 10, 55, 220, 715, 2002, 5005, 11440, 24310, 48620, 92377, 167960, 293920, 497420, 817135, 1307504, 2042755, 3124550, 4686110, 6906900, 10013002, 14307150, 20155070, 28048790, 38555660, 52451256, 70583095, 94143280, 124355000
Offset: 10

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    [&+[MoebiusMu(n div d)*Binomial(d-1, 9):d in Divisors(n)]:n in[10..42]]; // Marius A. Burtea, Feb 07 2020
  • Maple
    with(numtheory):
    a:= n-> add(mobius(n/d)*binomial(d-1, 9), d=divisors(n)):
    seq(a(n), n=10..50);  # Alois P. Heinz, Feb 05 2020
  • Mathematica
    a[n_]:=Sum[Boole[Divisible[n, k]] MoebiusMu[n/k] Binomial[k -1, 9], {k, 1, n}]; Table[a[n], {n, 10, 45}] (* or *) Table[a[n],{n, 9, 45}] a[n_]:=DivisorSum[n, Binomial[#-1,9] MoebiusMu[n/#]&]; Array[a, 37, 10] (* Vincenzo Librandi, Feb 08 2020 *)

Formula

Moebius transform of C(n-1,9).
G.f.: Sum_{k>=1} mu(k) * x^(10*k) / (1 - x^k)^10. - Ilya Gutkovskiy, Feb 05 2020

A101391 Triangle read by rows: T(n,k) is the number of compositions of n into k parts x_1, x_2, ..., x_k such that gcd(x_1,x_2,...,x_k) = 1 (1<=k<=n).

Original entry on oeis.org

1, 0, 1, 0, 2, 1, 0, 2, 3, 1, 0, 4, 6, 4, 1, 0, 2, 9, 10, 5, 1, 0, 6, 15, 20, 15, 6, 1, 0, 4, 18, 34, 35, 21, 7, 1, 0, 6, 27, 56, 70, 56, 28, 8, 1, 0, 4, 30, 80, 125, 126, 84, 36, 9, 1, 0, 10, 45, 120, 210, 252, 210, 120, 45, 10, 1, 0, 4, 42, 154, 325, 461, 462, 330, 165, 55, 11, 1, 0, 12, 66, 220, 495, 792, 924, 792, 495, 220, 66, 12, 1
Offset: 1

Views

Author

Emeric Deutsch, Jan 26 2005

Keywords

Comments

If instead we require that the individual parts (x_i,x_j) be relatively prime, we get A282748. This is the question studied by Shonhiwa (2006). - N. J. A. Sloane, Mar 05 2017.

Examples

			T(6,3)=9 because we have 411,141,114 and the six permutations of 123 (222 does not qualify).
T(8,3)=18 because binomial(0,2)*mobius(8/1)+binomial(1,2)*mobius(8/2)+binomial(3,2)*mobius(8/4)+binomial(7,2)*mobius(8/8)=0+0+(-3)+21=18.
Triangle begins:
   1;
   0,  1;
   0,  2,  1;
   0,  2,  3,   1;
   0,  4,  6,   4,   1;
   0,  2,  9,  10,   5,   1;
   0,  6, 15,  20,  15,   6,   1;
   0,  4, 18,  34,  35,  21,   7,   1;
   0,  6, 27,  56,  70,  56,  28,   8,   1;
   0,  4, 30,  80, 125, 126,  84,  36,   9,   1;
   0, 10, 45, 120, 210, 252, 210, 120,  45,  10,  1;
   0,  4, 42, 154, 325, 461, 462, 330, 165,  55, 11,  1;
   0, 12, 66, 220, 495, 792, 924, 792, 495, 220, 66, 12, 1;
  ...
From _Gus Wiseman_, Oct 19 2020: (Start)
Row n = 6 counts the following compositions:
  (15)  (114)  (1113)  (11112)  (111111)
  (51)  (123)  (1122)  (11121)
        (132)  (1131)  (11211)
        (141)  (1212)  (12111)
        (213)  (1221)  (21111)
        (231)  (1311)
        (312)  (2112)
        (321)  (2121)
        (411)  (2211)
               (3111)
Missing are: (42), (24), (33), (222).
(End)
		

Crossrefs

Mirror image of A039911.
Row sums are A000740.
A000837 counts relatively prime partitions.
A135278 counts compositions by length.
A282748 is the pairwise coprime instead of relatively prime version.
A282750 is the unordered version.
A291166 ranks these compositions (evidently).
T(2n+1,n+1) gives A000984.

Programs

  • Maple
    with(numtheory): T:=proc(n,k) local d, j, b: d:=divisors(n): for j from 1 to tau(n) do b[j]:=binomial(d[j]-1,k-1)*mobius(n/d[j]) od: sum(b[i],i=1..tau(n)) end: for n from 1 to 14 do seq(T(n,k),k=1..n) od; # yields the sequence in triangular form
    # second Maple program:
    b:= proc(n, g) option remember; `if`(n=0, `if`(g=1, 1, 0),
          expand(add(b(n-j, igcd(g, j))*x, j=1..n)))
        end:
    T:= (n, k)-> coeff(b(n,0),x,k):
    seq(seq(T(n,k), k=1..n), n=1..14);  # Alois P. Heinz, May 05 2025
  • Mathematica
    t[n_, k_] := Sum[Binomial[d-1, k-1]*MoebiusMu[n/d], {d, Divisors[n]}]; Table[t[n, k], {n, 2, 14}, {k, 2, n}] // Flatten (* Jean-François Alcover, Jan 20 2014 *)
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n,{k}],GCD@@#==1&]],{n,10},{k,2,n}] (* change {k,2,n} to {k,1,n} for the version with zeros. - Gus Wiseman, Oct 19 2020 *)
  • PARI
    T(n, k) = sumdiv(n, d, binomial(d-1, k-1)*moebius(n/d)); \\ Michel Marcus, Mar 09 2016

Formula

T(n,k) = Sum_{d|n} binomial(d-1,k-1)*mobius(n/d).
Sum_{k=1..n} k * T(n,k) = A085411(n). - Alois P. Heinz, May 05 2025

Extensions

Definition clarified by N. J. A. Sloane, Mar 05 2017
Edited by Alois P. Heinz, May 05 2025

A023033 Number of compositions of n into 8 ordered relatively prime parts.

Original entry on oeis.org

1, 8, 36, 120, 330, 792, 1716, 3432, 6434, 11440, 19440, 31824, 50352, 77520, 116160, 170544, 244826, 346104, 479908, 657792, 886314, 1184040, 1557312, 2035800, 2623140, 3365736, 4260608, 5379616, 6704742, 8347680, 10263648, 12619464
Offset: 8

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    [&+[MoebiusMu(n div d)*Binomial(d-1,7):d in Divisors(n)]:n in[8..40]]; // Marius A. Burtea, Feb 07 2020
  • Maple
    with(numtheory):
    a:= n-> add(mobius(n/d)*binomial(d-1, 7), d=divisors(n)):
    seq(a(n), n=8..50);  # Alois P. Heinz, Feb 05 2020
  • Mathematica
    a[n_]:=DivisorSum[n, Binomial[# - 1, 7] MoebiusMu[n/#]&]; Array[a, 37, 8] (* or *) a[n_]:=Sum[Boole[Divisible[n, k]] MoebiusMu[n/k] Binomial[k - 1, 7], {k, 1, n}]; Table[a[n], {n, 8, 45}] (* Vincenzo Librandi, Feb 07 2020 *)

Formula

Moebius transform of C(n-1,7).
G.f.: Sum_{k>=1} mu(k) * x^(8*k) / (1 - x^k)^8. - Ilya Gutkovskiy, Feb 05 2020

A023034 Number of compositions of n into 9 ordered relatively prime parts.

Original entry on oeis.org

1, 9, 45, 165, 495, 1287, 3003, 6435, 12870, 24309, 43758, 75573, 125970, 203445, 319770, 490149, 735471, 1081080, 1562274, 2218788, 3108105, 4289133, 5852925, 7882290, 10518255, 13871286, 18156204, 23511345, 30260340, 38564262, 48902997
Offset: 9

Views

Author

Keywords

Crossrefs

Programs

  • Magma
    [&+[MoebiusMu(n div d)*Binomial(d-1,8):d in Divisors(n)]:n in[9..39]]; // Marius A. Burtea, Feb 07 2020
  • Maple
    with(numtheory):
    a:= n-> add(mobius(n/d)*binomial(d-1, 8), d=divisors(n)):
    seq(a(n), n=9..50);  # Alois P. Heinz, Feb 05 2020
  • Mathematica
    Table[a[n],{n,9,45}]a[n_]:=DivisorSum[n, Binomial[#-1, 8] MoebiusMu[n/#]&]; Array[a, 37, 9] (* or *) a[n_]:=Sum[Boole[Divisible[n,k]] MoebiusMu[n/k] Binomial[k-1,8],{k,1,n}];Table[a[n],{n,9,45}] (* Vincenzo Librandi, Feb 08 2020 *)

Formula

Moebius transform of C(n-1,8).
G.f.: Sum_{k>=1} mu(k) * x^(9*k) / (1 - x^k)^9. - Ilya Gutkovskiy, Feb 05 2020

A341870 Number of partitions of n into 6 distinct and relatively prime parts.

Original entry on oeis.org

1, 1, 2, 3, 5, 7, 11, 14, 20, 26, 35, 44, 58, 71, 90, 110, 136, 163, 199, 235, 282, 330, 391, 453, 532, 610, 709, 808, 931, 1052, 1206, 1353, 1540, 1718, 1945, 2158, 2432, 2682, 3009, 3305, 3692, 4035, 4493, 4891, 5427, 5883, 6510, 7033, 7758, 8352, 9192, 9862, 10829, 11584, 12687, 13539
Offset: 21

Views

Author

Ilya Gutkovskiy, Feb 23 2021

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 76; CoefficientList[Series[Sum[MoebiusMu[k] x^(21 k)/Product[1 - x^(j k), {j, 1, 6}], {k, 1, nmax}], {x, 0, nmax}], x] // Drop[#, 21] &

Formula

G.f.: Sum_{k>=1} mu(k)* x^(21*k) / Product_{j=1..6} (1 - x^(j*k)).
Showing 1-9 of 9 results.