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 10 results.

A051026 Number of primitive subsequences of {1, 2, ..., n}.

Original entry on oeis.org

1, 2, 3, 5, 7, 13, 17, 33, 45, 73, 103, 205, 253, 505, 733, 1133, 1529, 3057, 3897, 7793, 10241, 16513, 24593, 49185, 59265, 109297, 163369, 262489, 355729, 711457, 879937, 1759873, 2360641, 3908545, 5858113, 10534337, 12701537, 25403073, 38090337, 63299265, 81044097, 162088193, 205482593, 410965185, 570487233, 855676353
Offset: 0

Views

Author

Keywords

Comments

a(n) counts all subsequences of {1, ..., n} in which no term divides any other. If n is a prime a(n) = 2*a(n-1)-1 because for each subsequence s counted by a(n-1) two different subsequences are counted by a(n): s and s,n. There is only one exception: 1,n is not a primitive subsequence because 1 divides n. For all n>1: a(n) < 2*a(n-1). - Alois P. Heinz, Mar 07 2011
Maximal primitive subsets are counted by A326077. - Gus Wiseman, Jun 07 2019

Examples

			a(4) = 7, the primitive subsequences (including the empty sequence) are: (), (1), (2), (3), (4), (2,3), (3,4).
a(5) = 13 = 2*7-1, the primitive subsequences are: (), (5), (1), (2), (2,5), (3), (3,5), (4), (4,5), (2,3), (2,3,5), (3,4), (3,4,5).
From _Gus Wiseman_, Jun 07 2019: (Start)
The a(0) = 1 through a(5) = 13 primitive (pairwise indivisible) subsets:
  {}  {}   {}   {}     {}     {}
      {1}  {1}  {1}    {1}    {1}
           {2}  {2}    {2}    {2}
                {3}    {3}    {3}
                {2,3}  {4}    {4}
                       {2,3}  {5}
                       {3,4}  {2,3}
                              {2,5}
                              {3,4}
                              {3,5}
                              {4,5}
                              {2,3,5}
                              {3,4,5}
a(n) is also the number of subsets of {1..n} containing all of their pairwise products <= n as well as any quotients of divisible elements. For example, the a(0) = 1 through a(5) = 13 subsets are:
  {}  {}   {}     {}       {}         {}
      {1}  {1}    {1}      {1}        {1}
           {1,2}  {1,2}    {1,3}      {1,3}
                  {1,3}    {1,4}      {1,4}
                  {1,2,3}  {1,2,4}    {1,5}
                           {1,3,4}    {1,2,4}
                           {1,2,3,4}  {1,3,4}
                                      {1,3,5}
                                      {1,4,5}
                                      {1,2,3,4}
                                      {1,2,4,5}
                                      {1,3,4,5}
                                      {1,2,3,4,5}
Also the number of subsets of {1..n} containing all of their multiples <= n. For example, the a(0) = 1 through a(5) = 13 subsets are:
  {}  {}   {}     {}       {}         {}
      {1}  {2}    {2}      {3}        {3}
           {1,2}  {3}      {4}        {4}
                  {2,3}    {2,4}      {5}
                  {1,2,3}  {3,4}      {2,4}
                           {2,3,4}    {3,4}
                           {1,2,3,4}  {3,5}
                                      {4,5}
                                      {2,3,4}
                                      {2,4,5}
                                      {3,4,5}
                                      {2,3,4,5}
                                      {1,2,3,4,5}
(End)
From _Gus Wiseman_, Mar 12 2024: (Start)
Also the number of subsets of {1..n} containing all divisors of the elements. For example, the a(0) = 1 through a(6) = 17 subsets are:
  {}  {}   {}     {}       {}         {}
      {1}  {1}    {1}      {1}        {1}
           {1,2}  {1,2}    {1,2}      {1,2}
                  {1,3}    {1,3}      {1,3}
                  {1,2,3}  {1,2,3}    {1,5}
                           {1,2,4}    {1,2,3}
                           {1,2,3,4}  {1,2,4}
                                      {1,2,5}
                                      {1,3,5}
                                      {1,2,3,4}
                                      {1,2,3,5}
                                      {1,2,4,5}
                                      {1,2,3,4,5}
(End)
		

References

  • Blanchet-Sadri, Francine. Algorithmic combinatorics on partial words. Chapman & Hall/CRC, Boca Raton, FL, 2008. ii+385 pp. ISBN: 978-1-4200-6092-8; 1-4200-6092-9 MR2384993 (2009f:68142). See p. 320. - N. J. A. Sloane, Apr 06 2012

Crossrefs

Programs

  • Maple
    with(numtheory):
    b:= proc(s) option remember; local n;
          n:= max(s[]);
          `if`(n<0, 1, b(s minus {n}) + b(s minus divisors(n)))
        end:
    bb:= n-> b({$2..n} minus divisors(n)):
    sb:= proc(n) option remember; `if`(n<2, 0, bb(n) + sb(n-1)) end:
    a:= n-> `if`(n=0, 1, `if`(isprime(n), 2*a(n-1)-1, 2+sb(n))):
    seq(a(n), n=0..40);  # Alois P. Heinz, Mar 07 2011
  • Mathematica
    b[s_] := b[s] = With[{n=Max[s]}, If[n < 0, 1, b[Complement[s, {n}]] + b[Complement[s, Divisors[n]]]]];
    bb[n_] := b[Complement[Range[2, n], Divisors[n]]];
    sb[n_] := sb[n] = If[n < 2, 0, bb[n] + sb[n-1]];
    a[n_] := If[n == 0, 1, If[PrimeQ[n], 2a[n-1] - 1, 2 + sb[n]]]; Table[a[n], {n, 0, 37}]
    (* Jean-François Alcover, Jul 27 2011, converted from Maple *)
    Table[Length[Select[Subsets[Range[n]], SubsetQ[#,Select[Union@@Table[#*i,{i,n}],#<=n&]]&]],{n,10}] (* Gus Wiseman, Jun 07 2019 *)
    Table[Length[Select[Subsets[Range[n]], #==Union@@Divisors/@#&]],{n,0,10}] (* Gus Wiseman, Mar 12 2024 *)

Extensions

More terms from David Wasserman, May 02 2002
a(32)-a(37) from Donovan Johnson, Aug 11 2010

A343658 Array read by antidiagonals where A(n,k) is the number of ways to choose a multiset of k divisors of n.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 1, 4, 3, 3, 1, 1, 5, 4, 6, 2, 1, 1, 6, 5, 10, 3, 4, 1, 1, 7, 6, 15, 4, 10, 2, 1, 1, 8, 7, 21, 5, 20, 3, 4, 1, 1, 9, 8, 28, 6, 35, 4, 10, 3, 1, 1, 10, 9, 36, 7, 56, 5, 20, 6, 4, 1, 1, 11, 10, 45, 8, 84, 6, 35, 10, 10, 2, 1
Offset: 1

Views

Author

Gus Wiseman, Apr 29 2021

Keywords

Comments

First differs from A343656 at A(4,2) = 6, A343656(4,2) = 5.
As a triangle, T(n,k) = number of ways to choose a multiset of n - k divisors of k.

Examples

			Array begins:
       k=0 k=1 k=2 k=3 k=4 k=5 k=6 k=7 k=8
  n=1:  1   1   1   1   1   1   1   1   1
  n=2:  1   2   3   4   5   6   7   8   9
  n=3:  1   2   3   4   5   6   7   8   9
  n=4:  1   3   6  10  15  21  28  36  45
  n=5:  1   2   3   4   5   6   7   8   9
  n=6:  1   4  10  20  35  56  84 120 165
  n=7:  1   2   3   4   5   6   7   8   9
  n=8:  1   4  10  20  35  56  84 120 165
  n=9:  1   3   6  10  15  21  28  36  45
Triangle begins:
   1
   1   1
   1   2   1
   1   3   2   1
   1   4   3   3   1
   1   5   4   6   2   1
   1   6   5  10   3   4   1
   1   7   6  15   4  10   2   1
   1   8   7  21   5  20   3   4   1
   1   9   8  28   6  35   4  10   3   1
   1  10   9  36   7  56   5  20   6   4   1
   1  11  10  45   8  84   6  35  10  10   2   1
For example, row n = 6 counts the following multisets:
  {1,1,1,1,1}  {1,1,1,1}  {1,1,1}  {1,1}  {1}  {}
               {1,1,1,2}  {1,1,3}  {1,2}  {5}
               {1,1,2,2}  {1,3,3}  {1,4}
               {1,2,2,2}  {3,3,3}  {2,2}
               {2,2,2,2}           {2,4}
                                   {4,4}
Note that for n = 6, k = 4 in the triangle, the two multisets {1,4} and {2,2} represent the same divisor 4, so they are only counted once under A343656(4,2) = 5.
		

Crossrefs

Row k = 1 of the array is A000005.
Column n = 4 of the array is A000217.
Column n = 6 of the array is A000292.
Row k = 2 of the array is A184389.
The distinct products of these multisets are counted by A343656.
Antidiagonal sums of the array (or row sums of the triangle) are A343661.
A000312 = n^n.
A009998(n,k) = n^k (as an array, offset 1).
A007318 counts k-sets of elements of {1..n}.
A059481 counts k-multisets of elements of {1..n}.

Programs

  • Mathematica
    multchoo[n_,k_]:=Binomial[n+k-1,k];
    Table[multchoo[DivisorSigma[0,k],n-k],{n,10},{k,n}]
  • PARI
    A(n,k) = binomial(numdiv(n) + k - 1, k)
    { for(n=1, 9, for(k=0, 8, print1(A(n,k), ", ")); print ) } \\ Andrew Howroyd, Jan 11 2024

Formula

A(n,k) = ((A000005(n), k)) = A007318(A000005(n) + k - 1, k).
T(n,k) = ((A000005(k), n - k)) = A007318(A000005(k) + n - k - 1, n - k).

A343652 Number of maximal pairwise coprime sets of divisors of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 3, 2, 2, 1, 4, 1, 2, 2, 4, 1, 4, 1, 4, 2, 2, 1, 6, 2, 2, 3, 4, 1, 5, 1, 5, 2, 2, 2, 8, 1, 2, 2, 6, 1, 5, 1, 4, 4, 2, 1, 8, 2, 4, 2, 4, 1, 6, 2, 6, 2, 2, 1, 10, 1, 2, 4, 6, 2, 5, 1, 4, 2, 5, 1, 12, 1, 2, 4, 4, 2, 5, 1, 8, 4, 2, 1, 10, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Apr 25 2021

Keywords

Comments

Also the number of maximal pairwise coprime sets of divisors > 1 of n. For example, the a(n) sets for n = 12, 30, 36, 60, 120 are:
{6} {30} {6} {30} {30}
{12} {2,15} {12} {60} {60}
{2,3} {3,10} {18} {2,15} {120}
{3,4} {5,6} {36} {3,10} {2,15}
{2,3,5} {2,3} {3,20} {3,10}
{2,9} {4,15} {3,20}
{3,4} {5,6} {3,40}
{4,9} {5,12} {4,15}
{2,3,5} {5,6}
{3,4,5} {5,12}
{5,24}
{8,15}
{2,3,5}
{3,4,5}
{3,5,8}

Examples

			The a(n) sets for n = 12, 30, 36, 60, 120:
  {1,6}    {1,30}     {1,6}    {1,30}     {1,30}
  {1,12}   {1,2,15}   {1,12}   {1,60}     {1,60}
  {1,2,3}  {1,3,10}   {1,18}   {1,2,15}   {1,120}
  {1,3,4}  {1,5,6}    {1,36}   {1,3,10}   {1,2,15}
           {1,2,3,5}  {1,2,3}  {1,3,20}   {1,3,10}
                      {1,2,9}  {1,4,15}   {1,3,20}
                      {1,3,4}  {1,5,6}    {1,3,40}
                      {1,4,9}  {1,5,12}   {1,4,15}
                               {1,2,3,5}  {1,5,6}
                               {1,3,4,5}  {1,5,12}
                                          {1,5,24}
                                          {1,8,15}
                                          {1,2,3,5}
                                          {1,3,4,5}
                                          {1,3,5,8}
		

Crossrefs

The case of pairs is A063647.
The case of triples is A066620.
The non-maximal version counting empty sets and singletons is A225520.
The non-maximal version with no 1's is A343653.
The non-maximal version is A343655.
The version for subsets of {1..n} is A343659.
The case without 1's or singletons is A343660.
A018892 counts pairwise coprime unordered pairs of divisors.
A048691 counts pairwise coprime ordered pairs of divisors.
A048785 counts pairwise coprime ordered triples of divisors.
A084422, A187106, A276187, and A320426 count pairwise coprime sets.
A100565 counts pairwise coprime unordered triples of divisors.
A305713 counts pairwise coprime non-singleton strict partitions.
A324837 counts minimal subsets of {1...n} with least common multiple n.
A325683 counts maximal Golomb rulers.
A326077 counts maximal pairwise indivisible sets.

Programs

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

Formula

a(n) = A343660(n) + A005361(n).

A276187 Number of subsets of {1,..,n} of cardinality >= 2 such that the elements of each counted subset are pairwise coprime.

Original entry on oeis.org

0, 1, 4, 7, 18, 21, 48, 63, 94, 105, 220, 235, 482, 529, 600, 711, 1438, 1501, 3020, 3211, 3594, 3849, 7720, 7975, 11142, 11877, 14628, 15459, 30946, 31201, 62432, 69855, 76126, 80221, 89820, 91611, 183258, 192601, 208600, 214231, 428502, 431573, 863188, 900563
Offset: 1

Views

Author

Robert C. Lyons, Aug 23 2016

Keywords

Comments

n is prime if and only if a(n) = 2*a(n-1)+n-1. - Robert Israel, Aug 24 2016

Examples

			From _Gus Wiseman_, May 08 2021: (Start)
The a(2) = 1 through a(6) = 21 sets:
  {1,2}   {1,2}    {1,2}     {1,2}      {1,2}
          {1,3}    {1,3}     {1,3}      {1,3}
          {2,3}    {1,4}     {1,4}      {1,4}
         {1,2,3}   {2,3}     {1,5}      {1,5}
                   {3,4}     {2,3}      {1,6}
                  {1,2,3}    {2,5}      {2,3}
                  {1,3,4}    {3,4}      {2,5}
                             {3,5}      {3,4}
                             {4,5}      {3,5}
                            {1,2,3}     {4,5}
                            {1,2,5}     {5,6}
                            {1,3,4}    {1,2,3}
                            {1,3,5}    {1,2,5}
                            {1,4,5}    {1,3,4}
                            {2,3,5}    {1,3,5}
                            {3,4,5}    {1,4,5}
                           {1,2,3,5}   {1,5,6}
                           {1,3,4,5}   {2,3,5}
                                       {3,4,5}
                                      {1,2,3,5}
                                      {1,3,4,5}
(End)
		

Crossrefs

The case of pairs is A015614.
The indivisible instead of coprime version is A051026(n) - n.
Allowing empty sets and singletons gives A084422.
The relatively prime instead of pairwise coprime version is A085945(n) - 1.
Allowing all singletons gives A187106.
Allowing only the singleton {1} gives A320426.
Row sums of A320436, each minus one.
The maximal case is counted by A343659.
The version for sets of divisors is A343655(n) - 1.
A000005 counts divisors.
A186972 counts pairwise coprime k-sets containing n.
A186974 counts pairwise coprime k-sets.
A326675 ranks pairwise coprime non-singleton sets.

Programs

  • Maple
    f:= proc(S) option remember;
        local s, Sp;
        if S = {} then return 1 fi;
        s:= S[-1];
        Sp:= S[1..-2];
        procname(Sp) + procname(select(t -> igcd(t,s)=1, Sp))
    end proc:
    seq(f({$1..n}) - n - 1, n=1..50); # Robert Israel, Aug 24 2016
  • Mathematica
    f[S_] := f[S] = Module[{s, Sp}, If[S == {}, Return[1]]; s = S[[-1]]; Sp = S[[1;;-2]]; f[Sp] + f[Select[Sp, GCD[#, s] == 1&]]];
    Table[f[Range[n]] - n - 1, {n, 1, 50}] (* Jean-François Alcover, Sep 15 2022, after Robert Israel *)
  • PARI
    f(n,k=1)=if(n==1, return(2)); if(gcd(k,n)==1, f(n-1,n*k)) + f(n-1,k)
    a(n)=f(n)-n-1 \\ Charles R Greathouse IV, Aug 24 2016
  • Sage
    from sage.combinat.subsets_pairwise import PairwiseCompatibleSubsets
    def is_coprime(x, y): return gcd(x, y) == 1
    max_n = 40
    seq = []
    for n in range(1, max_n+1):
        P = PairwiseCompatibleSubsets(range(1,n+1), is_coprime)
        a_n = len([1 for s in P.list() if len(s) > 1])
        seq.append(a_n)
    print(seq)
    

Formula

a(n) = A320426(n) - 1. - Gus Wiseman, May 08 2021

Extensions

Name and example edited by Robert Israel, Aug 24 2016

A343655 Number of pairwise coprime sets of divisors of n, where a singleton is not considered pairwise coprime unless it is {1}.

Original entry on oeis.org

1, 2, 2, 3, 2, 6, 2, 4, 3, 6, 2, 10, 2, 6, 6, 5, 2, 10, 2, 10, 6, 6, 2, 14, 3, 6, 4, 10, 2, 22, 2, 6, 6, 6, 6, 17, 2, 6, 6, 14, 2, 22, 2, 10, 10, 6, 2, 18, 3, 10, 6, 10, 2, 14, 6, 14, 6, 6, 2, 38, 2, 6, 10, 7, 6, 22, 2, 10, 6, 22, 2, 24, 2, 6, 10, 10, 6, 22, 2
Offset: 1

Views

Author

Gus Wiseman, Apr 26 2021

Keywords

Comments

First differs from A015995 at a(210) = 88, A015995(210) = 86.

Examples

			For example, the a(n) subsets for n = 1, 2, 4, 6, 8, 12, 16, 24 are:
  {1}  {1}    {1}    {1}      {1}    {1}      {1}     {1}
       {1,2}  {1,2}  {1,2}    {1,2}  {1,2}    {1,2}   {1,2}
              {1,4}  {1,3}    {1,4}  {1,3}    {1,4}   {1,3}
                     {1,6}    {1,8}  {1,4}    {1,8}   {1,4}
                     {2,3}           {1,6}    {1,16}  {1,6}
                     {1,2,3}         {2,3}            {1,8}
                                     {3,4}            {2,3}
                                     {1,12}           {3,4}
                                     {1,2,3}          {3,8}
                                     {1,3,4}          {1,12}
                                                      {1,24}
                                                      {1,2,3}
                                                      {1,3,4}
                                                      {1,3,8}
		

Crossrefs

The case of pairs is A063647.
The case of triples is A066620.
The version with empty sets and singletons is A225520.
A version for prime indices is A304711.
The version for strict integer partitions is A305713.
The version for subsets of {1..n} is A320426 = A276187 + 1.
The version for binary indices is A326675.
The version for integer partitions is A327516.
The version for standard compositions is A333227.
The maximal case is A343652.
The case without 1's is A343653.
The case without 1's with singletons is A343654.
The maximal case without 1's 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.
A325683 counts maximal Golomb rulers.
A326077 counts maximal pairwise indivisible sets.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Divisors[n]],CoprimeQ@@#&]],{n,100}]

A343654 Number of pairwise coprime sets of divisors > 1 of n.

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
Offset: 1

Views

Author

Gus Wiseman, Apr 26 2021

Keywords

Comments

First differs from A100565 at a(210) = 52, A100565(210) = 51.

Examples

			The a(n) sets for n = 1, 2, 4, 6, 8, 12, 24, 30, 32, 36, 48:
  {}  {}   {}   {}     {}   {}     {}     {}       {}    {}     {}
      {2}  {2}  {2}    {2}  {2}    {2}    {2}      {2}   {2}    {2}
           {4}  {3}    {4}  {3}    {3}    {3}      {4}   {3}    {3}
                {6}    {8}  {4}    {4}    {5}      {8}   {4}    {4}
                {2,3}       {6}    {6}    {6}      {16}  {6}    {6}
                            {12}   {8}    {10}     {32}  {9}    {8}
                            {2,3}  {12}   {15}           {12}   {12}
                            {3,4}  {24}   {30}           {18}   {16}
                                   {2,3}  {2,3}          {36}   {24}
                                   {3,4}  {2,5}          {2,3}  {48}
                                   {3,8}  {3,5}          {2,9}  {2,3}
                                          {5,6}          {3,4}  {3,4}
                                          {2,15}         {4,9}  {3,8}
                                          {3,10}                {3,16}
                                          {2,3,5}
		

Crossrefs

The version for partitions is A007359.
The version for subsets of {1..n} is A084422.
The case of pairs is A089233.
The version with 1's is A225520.
The maximal case is A343652.
The case without empty sets or singletons is A343653.
The maximal case without singletons is A343660.
A018892 counts pairwise coprime unordered pairs of divisors.
A051026 counts pairwise indivisible subsets of {1..n}.
A100565 counts pairwise coprime unordered triples of divisors.
A187106, A276187, and A320426 count other types of pairwise coprime sets.
A326077 counts maximal pairwise indivisible sets.

Programs

  • Mathematica
    pwcop[y_]:=And@@(GCD@@#1==1&)/@Subsets[y,{2}];
    Table[Length[Select[Subsets[Rest[Divisors[n]]],pwcop]],{n,100}]

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

Views

Author

Gus Wiseman, Apr 26 2021

Keywords

Comments

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

Examples

			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}
		

Crossrefs

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.

Programs

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

Extensions

More terms from Bert Dobbelaere, May 04 2021

A343660 Number of maximal pairwise coprime sets of at least two divisors > 1 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, 4, 0, 0, 1, 1, 1, 4, 0, 1, 1, 3, 0, 4, 0, 2, 2, 1, 0, 4, 0, 2, 1, 2, 0, 3, 1, 3, 1, 1, 0, 8, 0, 1, 2, 0, 1, 4, 0, 2, 1, 4, 0, 6, 0, 1, 2, 2, 1, 4, 0, 4, 0, 1, 0, 8, 1, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Apr 26 2021

Keywords

Examples

			The a(n) sets for n = 6, 12, 24, 30, 36, 60, 72, 96:
  {2,3}  {2,3}  {2,3}  {5,6}    {2,3}  {5,6}    {2,3}  {2,3}
         {3,4}  {3,4}  {2,15}   {2,9}  {2,15}   {2,9}  {3,4}
                {3,8}  {3,10}   {3,4}  {3,10}   {3,4}  {3,8}
                       {2,3,5}  {4,9}  {3,20}   {3,8}  {3,16}
                                       {4,15}   {4,9}  {3,32}
                                       {5,12}   {8,9}
                                       {2,3,5}
                                       {3,4,5}
		

Crossrefs

The case of pairs is A089233.
The case with 1's is A343652.
The case with singletons is (also) A343652.
The non-maximal version is A343653.
The non-maximal version with 1's is A343655.
The version for subsets of {2..n} is A343659 (for n > 2).
A018892 counts coprime unordered pairs of divisors.
A051026 counts pairwise indivisible subsets of {1..n}.
A066620 counts pairwise coprime 3-sets of divisors.
A100565 counts pairwise coprime unordered triples of divisors.

Programs

  • Mathematica
    fasmax[y_]:=Complement[y,Union@@Most@*Subsets/@y];
    Table[Length[fasmax[Select[Subsets[Rest[Divisors[n]]],CoprimeQ@@#&]]],{n,100}]

Formula

a(n) = A343652(n) - A005361(n).

A343661 Sum of numbers of y-multisets of divisors of x for each x >= 1, y >= 0, x + y = n.

Original entry on oeis.org

1, 2, 4, 7, 12, 19, 30, 46, 70, 105, 155, 223, 316, 443, 619, 865, 1210, 1690, 2354, 3263, 4497, 6157, 8368, 11280, 15078, 19989, 26296, 34356, 44626, 57693, 74321, 95503, 122535, 157101, 201377, 258155, 330994, 424398, 544035, 696995, 892104, 1140298, 1455080
Offset: 1

Views

Author

Gus Wiseman, Apr 30 2021

Keywords

Examples

			The a(5) = 12 multisets of divisors:
  {1,1,1,1}  {1,1,1}  {1,1}  {1}  {}
             {1,1,2}  {1,3}  {2}
             {1,2,2}  {3,3}  {4}
             {2,2,2}
		

Crossrefs

Antidiagonal sums of the array A343658 (or row sums of the triangle).
Dominates A343657.
A000005 counts divisors.
A007318 counts k-sets of elements of {1..n}.
A059481 counts k-multisets of elements of {1..n}.
A343656 counts divisors of powers.

Programs

  • Mathematica
    multchoo[n_,k_]:=Binomial[n+k-1,k];
    Table[Sum[multchoo[DivisorSigma[0,k],n-k],{k,n}],{n,10}]

Formula

a(n) = Sum_{k=1..n} binomial(sigma(k) + n - k - 1, n - k).

A326082 Number of maximal sets of pairwise indivisible divisors of n.

Original entry on oeis.org

1, 2, 2, 3, 2, 3, 2, 4, 3, 3, 2, 5, 2, 3, 3, 5, 2, 5, 2, 5, 3, 3, 2, 8, 3, 3, 4, 5, 2, 7, 2, 6, 3, 3, 3, 9, 2, 3, 3, 8, 2, 7, 2, 5, 5, 3, 2, 12, 3, 5, 3, 5, 2, 8, 3, 8, 3, 3, 2, 15, 2, 3, 5, 7, 3, 7, 2, 5, 3, 7, 2, 15, 2, 3, 5, 5, 3, 7, 2, 12, 5, 3, 2, 15, 3
Offset: 1

Views

Author

Gus Wiseman, Jun 05 2019

Keywords

Comments

Depends only on prime signature.
The non-maximal case is A096827.

Examples

			The maximal sets of pairwise indivisible divisors of n = 1, 2, 4, 8, 12, 24, 30, 32, 36, 48, 60 are:
   1   1   1   1   1     1      1         1    1       1       1
       2   2   2   12    24     30        2    36      48      60
           4   4   2,3   2,3    5,6       4    2,3     2,3     2,15
               8   3,4   3,4    2,15      8    2,9     3,4     3,20
                   4,6   3,8    3,10      16   3,4     3,8     4,30
                         4,6    2,3,5     32   4,18    4,6     5,12
                         6,8    6,10,15        9,12    6,8     2,3,5
                         8,12                  12,18   3,16    3,4,5
                                               4,6,9   6,16    4,5,6
                                                       8,12    3,4,10
                                                       12,16   6,15,20
                                                       16,24   10,12,15
                                                               12,15,20
                                                               12,20,30
                                                               4,6,10,15
		

Crossrefs

Programs

  • Mathematica
    stableQ[u_,Q_]:=!Apply[Or,Outer[#1=!=#2&&Q[#1,#2]&,u,u,1],{0,1}];
    fasmax[y_]:=Complement[y,Union@@(Most[Subsets[#]]&/@y)];
    Table[Length[fasmax[Select[Rest[Subsets[Divisors[n]]],stableQ[#,Divisible]&]]],{n,100}]
Showing 1-10 of 10 results.