A015617 -id:A015617 - cp's OEIS frontend

A185953 Number of pairwise coprime triples of positive integers with largest element n (i.e., A015617(n) - A015617(n-1)).

0, 0, 1, 1, 5, 1, 11, 6, 12, 5, 31, 6, 45, 14, 21, 24, 79, 15, 101, 27, 49, 37, 149, 28, 126, 56, 109, 57, 241, 28, 277, 106, 137, 100, 183, 62, 395, 128, 193, 102, 489, 65, 541, 165, 203, 189, 649, 116, 551, 170, 347, 231, 829, 147, 506, 234, 434, 307, 1027, 119, 1101, 364, 450, 412, 727

Offset: 1

N. J. A. Sloane, Feb 07 2011

Robert Israel, Posting to Sequence Fans Mailing List, Feb 06, 2011

Alois P. Heinz, Table of n, a(n) for n = 1..10000

Cf. A015617. Row sums of triangle A186230. Column 3 of triangle A186972.

a[n_] := Sum[Boole[GCD[a, n] == 1 && GCD[b, n] == 1 && GCD[a, b] == 1], {a, 1, n-2}, {b, a+1, n-1}]; Array[a, 100] (* Jean-François Alcover, Mar 05 2019, from PARI *)

a(n)=sum(a=1,n-2,sum(b=a+1,n-1,gcd(a,n)==1&&gcd(b,n)==1&&gcd(a,b)==1)) \\ Charles R Greathouse IV, Apr 28 2015

A220377 Number of partitions of n into three distinct and mutually relatively prime parts.

Original entry on oeis.org

1, 0, 2, 1, 3, 1, 6, 1, 7, 3, 7, 3, 14, 3, 15, 6, 14, 6, 25, 6, 22, 10, 25, 9, 42, 8, 34, 15, 37, 15, 53, 13, 48, 22, 53, 17, 78, 17, 65, 30, 63, 24, 99, 24, 88, 35, 84, 30, 126, 34, 103, 45, 103, 38, 166, 35, 124, 57, 128, 51, 184, 44, 150, 67, 172, 52, 218

Offset: 6

Carl Najafi, Dec 13 2012

nonn

The Heinz numbers of these partitions are the intersection of A005117 (strict), A014612 (triples), and A302696 (coprime). - Gus Wiseman, Oct 14 2020

			For n=10 we have three such partitions: 1+2+7, 1+4+5 and 2+3+5.
From _Gus Wiseman_, Oct 14 2020: (Start)
The a(6) = 1 through a(20) = 15 triples (empty column indicated by dot, A..H = 10..17):
321  .  431  531  532  731  543  751  743  753  754  971  765  B53  875
        521       541       651       752  951  853  B51  873  B71  974
                  721       732       761  B31  871  D31  954  D51  A73
                            741       851       952       972       A91
                            831       941       B32       981       B54
                            921       A31       B41       A71       B72
                                      B21       D21       B43       B81
                                                          B52       C71
                                                          B61       D43
                                                          C51       D52
                                                          D32       D61
                                                          D41       E51
                                                          E31       F41
                                                          F21       G31
                                                                    H21
(End)

Fausto A. C. Cariboni, Table of n, a(n) for n = 6..10000 (terms 6..1000 from Seiichi Manyama)

Cf. A015617, A300815.
A023022 is the 2-part version.
A101271 is the relative prime instead of pairwise coprime version.
A220377*6 is the ordered version.
A305713 counts these partitions of any length, with Heinz numbers A302797.
A307719 is the non-strict version.
A337461 is the non-strict ordered version.
A337563 is the case with no 1's.
A337605 is the pairwise non-coprime instead of pairwise coprime version.
A001399(n-6) counts strict 3-part partitions, with Heinz numbers A007304.
A008284 counts partitions by sum and length, with strict case A008289.
A318717 counts pairwise non-coprime strict partitions.
A326675 ranks pairwise coprime sets.
A327516 counts pairwise coprime partitions.
A337601 counts 3-part partitions whose distinct parts are pairwise coprime.
Cf. A000217, A007360, A023023, A051424, A078374, A087087, A302696, A333227, A337485, A337561.

Table[Length@Select[ IntegerPartitions[ n, {3}], #[[1]] != #[[2]] != #[[3]] && GCD[#[[1]], #[[2]]] == 1 && GCD[#[[1]], #[[3]]] == 1 && GCD[#[[2]], #[[3]]] == 1 &], {n, 6, 100}]
Table[Count[IntegerPartitions[n,{3}],?(CoprimeQ@@#&&Length[ Union[#]] == 3&)],{n,6,100}] (* _Harvey P. Dale, May 22 2020 *)

a(n)=my(P=partitions(n));sum(i=1,#P,#P[i]==3&&P[i][1]Charles R Greathouse IV, Dec 14 2012

a(n > 2) = A307719(n) - 1. - Gus Wiseman, Oct 15 2020

A186974 Irregular triangle T(n,k), n>=1, 1<=k<=A036234(n), read by rows: T(n,k) is the number of k-element subsets of {1, 2, ..., n} having pairwise coprime elements.

Original entry on oeis.org

1, 2, 1, 3, 3, 1, 4, 5, 2, 5, 9, 7, 2, 6, 11, 8, 2, 7, 17, 19, 10, 2, 8, 21, 25, 14, 3, 9, 27, 37, 24, 6, 10, 31, 42, 26, 6, 11, 41, 73, 68, 32, 6, 12, 45, 79, 72, 33, 6, 13, 57, 124, 151, 105, 39, 6, 14, 63, 138, 167, 114, 41, 6, 15, 71, 159, 192, 128, 44, 6

Offset: 1

Alois P. Heinz, Mar 02 2011

T(n,k) = 0 for k > A036234(n). The triangle contains all positive values of T.

			T(5,3) = 7 because there are 7 3-element subsets of {1,2,3,4,5} having pairwise coprime elements: {1,2,3}, {1,2,5}, {1,3,4}, {1,3,5}, {1,4,5}, {2,3,5}, {3,4,5}.
Irregular Triangle T(n,k) begins:
  1;
  2,  1;
  3,  3,  1;
  4,  5,  2;
  5,  9,  7,  2;
  6, 11,  8,  2;
  7, 17, 19, 10, 2;

Alois P. Heinz, Rows n = 1..220, flattened

Columns k=1-10 give: A000027, A015614, A015617, A015623, A015698, A186982, A186983, A186984, A186985, A186986.
Row sums give A187106.
Rightmost terms of rows give A319187.
Cf. A036234, A186972.

with(numtheory):
s:= proc(m, r) option remember; mul(`if`(i pi(n) +1:
b:= proc(t, n, k) option remember; local c, d, h;
      if k=0 or k>n then 0
    elif k=1 then 1
    elif k=2 and t=n then `if`(n<2, 0, phi(n))
    else c:= 0;
         d:= 2-irem(t, 2);
         for h from 1 to n-1 by d do
           if igcd(t, h)=1 then c:= c +b(s(t*h, h), h, k-1) fi
         od; c
      fi
    end:
T:= proc(n, k) option remember;
       b(s(n, n), n, k) +`if`(n<2, 0, T(n-1, k))
    end:
seq(seq(T(n, k), k=1..a(n)), n=1..20);

s[m_, r_] := s[m, r] = Product[If[i < r, i, 1], {i, FactorInteger[m][[All, 1]]}]; a[n_] := PrimePi[n]+1; b[t_, n_, k_] := b[t, n, k] = Module[{c, d, h}, Which[k == 0 || k > n, 0, k == 1, 1, k == 2 && t == n, If[n < 2, 0, EulerPhi[n]], True, c = 0; d = 2-Mod[t, 2]; For[h = 1, h <= n-1, h = h+d, If[ GCD[t, h] == 1, c = c + b[s[t*h, h], h, k-1]]]; c]]; t[n_, k_] := t[n, k] = b[s[n, n], n, k] + If[n < 2, 0, t[n-1, k]]; Table[Table[t[n, k], { k, 1, a[n]}], {n, 1, 20}] // Flatten (* Jean-François Alcover, Dec 17 2013, translated from Maple *)

T(n,k) = Sum_{i=1..n} A186972(i,k).

A100565 a(n) = Card{(x,y,z) : x <= y <= z, x|n, y|n, z|n, gcd(x,y)=1, gcd(x,z)=1, gcd(y,z)=1}.

Original entry on oeis.org

1, 2, 2, 3, 2, 5, 2, 4, 3, 5, 2, 8, 2, 5, 5, 5, 2, 8, 2, 8, 5, 5, 2, 11, 3, 5, 4, 8, 2, 15, 2, 6, 5, 5, 5, 13, 2, 5, 5, 11, 2, 15, 2, 8, 8, 5, 2, 14, 3, 8, 5, 8, 2, 11, 5, 11, 5, 5, 2, 25, 2, 5, 8, 7, 5, 15, 2, 8, 5, 15, 2, 18, 2, 5, 8, 8, 5, 15, 2, 14, 5, 5, 2, 25, 5, 5, 5, 11, 2, 25, 5, 8, 5, 5, 5, 17

Offset: 1

Vladeta Jovovic, Nov 28 2004

First differs from A018892 at a(30) = 15, A018892(30) = 14.
First differs from A343654 at a(210) = 51, A343654(210) = 52.
Also a(n) = Card{(x,y,z) : x <= y <= z and lcm(x,y)=n, lcm(x,z)=n, lcm(y,z)=n}.
In words, a(n) is the number of pairwise coprime unordered triples of divisors of n. - Gus Wiseman, May 01 2021

			From _Gus Wiseman_, May 01 2021: (Start)
The a(n) triples for n = 1, 2, 4, 6, 8, 12, 24:
  (1,1,1)  (1,1,1)  (1,1,1)  (1,1,1)  (1,1,1)  (1,1,1)   (1,1,1)
           (1,1,2)  (1,1,2)  (1,1,2)  (1,1,2)  (1,1,2)   (1,1,2)
                    (1,1,4)  (1,1,3)  (1,1,4)  (1,1,3)   (1,1,3)
                             (1,1,6)  (1,1,8)  (1,1,4)   (1,1,4)
                             (1,2,3)           (1,1,6)   (1,1,6)
                                               (1,2,3)   (1,1,8)
                                               (1,3,4)   (1,2,3)
                                               (1,1,12)  (1,3,4)
                                                         (1,3,8)
                                                         (1,1,12)
                                                         (1,1,24)
(End)

Antti Karttunen, Table of n, a(n) for n = 1..10000

Cf. A000005, A070919, A086222, A086165, A048691, A063647.
Positions of 2's through 5's are A000040, A001248, A030078, A068993.
The version for subsets of {1..n} instead of divisors is A015617.
The version for pairs of divisors is A018892.
The ordered version is A048785.
The strict case is A066620.
The version for strict partitions is A220377.
A version for sets of divisors of any size is A225520.
The version for partitions is A307719 (no 1's: A337563).
The case of distinct parts coprime is A337600 (ordered: A337602).
A001399(n-3) = A069905(n) = A211540(n+2) counts 3-part partitions.
A007304 ranks 3-part strict partitions.
A014311 ranks 3-part compositions.
A014612 ranks 3-part partitions.
A051026 counts pairwise indivisible subsets of {1..n}.
A302696 lists Heinz numbers of pairwise coprime partitions.
A337461 counts 3-part pairwise coprime compositions.
Cf. A000961, A000977, A007360, A023022, A087087, A276187, A282935, A337601, A337603, A338331, A343652, A343654.

pwcop[y_]:=And@@(GCD@@#==1&/@Subsets[y,{2}]);
Table[Length[Select[Tuples[Divisors[n],3],LessEqual@@#&&pwcop[#]&]],{n,30}] (* Gus Wiseman, May 01 2021 *)

A100565(n) = (numdiv(n^3)+3*numdiv(n)+2)/6; \\ Antti Karttunen, May 19 2017

a(n) = (tau(n^3) + 3*tau(n) + 2)/6.

A066620 Number of unordered triples of distinct pairwise coprime divisors of n.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 2, 0, 1, 1, 0, 0, 2, 0, 2, 1, 1, 0, 3, 0, 1, 0, 2, 0, 7, 0, 0, 1, 1, 1, 4, 0, 1, 1, 3, 0, 7, 0, 2, 2, 1, 0, 4, 0, 2, 1, 2, 0, 3, 1, 3, 1, 1, 0, 13, 0, 1, 2, 0, 1, 7, 0, 2, 1, 7, 0, 6, 0, 1, 2, 2, 1, 7, 0, 4, 0, 1, 0, 13, 1, 1, 1, 3, 0, 13, 1, 2, 1, 1, 1, 5, 0, 2, 2, 4, 0, 7, 0

Offset: 1

K. B. Subramaniam (kb_subramaniambalu(AT)yahoo.com) and Amarnath Murthy, Dec 24 2001

nonn

a(m) = a(n) if m and n have same factorization structure.

			a(24) = 3: the divisors of 24 are 1, 2, 3, 4, 6, 8, 12 and 24. The triples are (1, 2, 3), (1, 2, 9), (1, 3, 4).
a(30) = 7: the triples are (1, 2, 3), (1, 2, 5), (1, 3, 5), (2, 3, 5), (1, 3, 10), (1, 5, 6), (1, 2, 15).

Amarnath Murthy, Decomposition of the divisors of a natural number into pairwise coprime sets, Smarandache Notions Journal, vol. 12, No. 1-2-3, Spring 2001.pp 303-306.

Antti Karttunen, Table of n, a(n) for n = 1..10000

Positions of zeros are A000961.
Positions of ones are A006881.
The version for subsets of {1..n} instead of divisors is A015617.
The non-strict ordered version is A048785.
The version for pairs of divisors is A063647.
The non-strict version (3-multisets) is A100565.
The version for partitions is A220377 (non-strict: A307719).
A version for sets of divisors of any size is A225520.
A000005 counts divisors.
A001399(n-3) = A069905(n) = A211540(n+2) counts 3-part partitions.
A007304 ranks 3-part strict partitions.
A014311 ranks 3-part compositions.
A014612 ranks 3-part partitions.
A018892 counts unordered pairs of coprime divisors (ordered: A048691).
A051026 counts pairwise indivisible subsets of {1..n}.
A337461 counts 3-part pairwise coprime compositions.
A338331 lists Heinz numbers of pairwise coprime partitions.
Cf. A007360, A023022, A084422, A276187, A282935, A305713, A337563, A337605, A343652, A343655.

Table[Length[Select[Subsets[Divisors[n],{3}],CoprimeQ@@#&]],{n,100}] (* Gus Wiseman, Apr 28 2021 *)

A066620(n) = (numdiv(n^3)-3*numdiv(n)+2)/6; \\ After Jovovic's formula. - Antti Karttunen, May 27 2017

from sympy import divisor_count as d
def a(n): return (d(n**3) - 3*d(n) + 2)/6 # Indranil Ghosh, May 27 2017

In the reference it is shown that if k is a squarefree number with r prime factors and m with (r+1) prime factors then a(m) = 4*a(k) + 2^k - 1.

a(n) = (tau(n^3)-3*tau(n)+2)/6. - Vladeta Jovovic, Nov 27 2004

More terms from Vladeta Jovovic, Apr 03 2003
Name corrected by Andrey Zabolotskiy, Dec 09 2020
Name corrected by Gus Wiseman, Apr 28 2021 (ordered version is 6*a(n))

A343659 Number of maximal pairwise coprime subsets of {1..n}.

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 3, 4, 7, 9, 9, 10, 10, 12, 16, 19, 19, 20, 20, 22, 28, 32, 32, 33, 54, 61, 77, 84, 84, 85, 85, 94, 112, 123, 158, 161, 161, 176, 206, 212, 212, 214, 214, 229, 241, 260, 260, 263, 417, 428, 490, 521, 521, 526, 655, 674, 764, 818, 818, 820, 820, 874, 918, 975, 1182, 1189, 1189

Offset: 1

Gus Wiseman, Apr 26 2021

nonn

For this sequence, it does not matter whether singletons are considered pairwise coprime.
For n > 2, also the number of maximal pairwise coprime subsets of {2..n}.
For each prime p <= n, p divides exactly one element of each maximal subset. - Bert Dobbelaere, May 04 2021

			The a(1) = 1 through a(9) = 7 subsets:
  {1}  {12}  {123}  {123}  {1235}  {156}   {1567}   {1567}   {1567}
                    {134}  {1345}  {1235}  {12357}  {12357}  {12357}
                                   {1345}  {13457}  {13457}  {12579}
                                                    {13578}  {13457}
                                                             {13578}
                                                             {14579}
                                                             {15789}

Bert Dobbelaere, Table of n, a(n) for n = 1..500
Bert Dobbelaere, Python program

The case of pairs is A015614.
The case of triples is A015617.
The non-maximal version counting empty sets and singletons is A084422.
The non-maximal version counting singletons is A187106.
The non-maximal version is A320426(n) = A276187(n) + 1.
The version for indivisibility instead of coprimality is A326077.
The version for sets of divisors is A343652.
The version for sets of divisors > 1 is A343660.
A018892 counts coprime unordered pairs of divisors.
A051026 counts pairwise indivisible subsets of {1..n}.
A100565 counts pairwise coprime unordered triples of divisors.
Cf. A007360, A067824, A087087, A225520, A324837, A325683, A325859, A326358, A326496, A326675, A333227, A343653, A343655.

fasmax[y_]:=Complement[y,Union@@Most@*Subsets/@y];
Table[Length[fasmax[Select[Subsets[Range[n]],CoprimeQ@@#&]]],{n,15}]

More terms from Bert Dobbelaere, May 04 2021

cp's OEIS Frontend

A185953 Number of pairwise coprime triples of positive integers with largest element n (i.e., A015617(n) - A015617(n-1)).

Views

Author

Keywords

References

Links

Crossrefs

Programs

Mathematica

PARI

A220377 Number of partitions of n into three distinct and mutually relatively prime parts.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A186974 Irregular triangle T(n,k), n>=1, 1<=k<=A036234(n), read by rows: T(n,k) is the number of k-element subsets of {1, 2, ..., n} having pairwise coprime elements.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A100565 a(n) = Card{(x,y,z) : x <= y <= z, x|n, y|n, z|n, gcd(x,y)=1, gcd(x,z)=1, gcd(y,z)=1}.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A066620 Number of unordered triples of distinct pairwise coprime divisors of n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

Extensions

A343659 Number of maximal pairwise coprime subsets of {1..n}.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Extensions