A076078 -id:A076078 - cp's OEIS frontend

A212171 Prime signature of n (nonincreasing version): row n of table lists positive exponents in canonical prime factorization of n, in nonincreasing order.

1, 1, 2, 1, 1, 1, 1, 3, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 4, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 3, 1, 2, 1, 1, 3, 2, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 4, 1, 2, 2, 1, 1, 1, 2, 1, 1, 3, 1, 1, 1, 3, 1

Offset: 2

Matthew Vandermast, Jun 03 2012

Length of row n equals A001221(n).
The multiset of positive exponents in n's prime factorization completely determines a(n) for a host of OEIS sequences, including several "core" sequences.  Of those not cross-referenced here or in A212172, many can be found by searching the database for A025487.
(Note: Differing opinions may exist about whether the prime signature of n should be defined as this multiset itself, or as a symbol or collection of symbols that identify or "signify" this multiset. The definition of this sequence is designed to be compatible with either view, as are the original comments. When n >= 2, the customary ways to signify the multiset of exponents in n's prime factorization are to list the constituent exponents in either nonincreasing or nondecreasing order; this table gives the nonincreasing version.)
Table lists exponents in the order in which they appear in the prime factorization of a member of A025487.  This ordering is common in database comments (e.g., A008966).
Each possible multiset of an integer's positive prime factorization exponents corresponds to a unique partition that contains the same elements (cf. A000041). This includes the multiset of 1's positive exponents, { } (the empty multiset), which corresponds to the partition of 0.
Differs from A124010 from a(23) on, corresponding to the factorization of 18 = 2^1*3^2 which is here listed as row 18 = [2, 1], but as [1, 2] (in the order of the prime factors) in A124010 and also in A118914 which lists the prime signatures in nondecreasing order (so that row 12 = 2^2*3^1 is also [1, 2]). - M. F. Hasler, Apr 08 2022

			First rows of table read:
  1;
  1;
  2;
  1;
  1,1;
  1;
  3;
  2;
  1,1;
  1;
  2,1;
  ...
The multiset of positive exponents in the prime factorization of 6 = 2*3 is {1,1} (1s are often left implicit as exponents). The prime signature of 6 is therefore {1,1}.
12 = 2^2*3 has positive exponents 2 and 1 in its prime factorization, as does 18 = 2*3^2. Rows 12 and 18 of the table both read {2,1}.

Jason Kimberley, Table of i, a(i) for i = 2..24301 (n = 2..10000)

Cf. A025487, A001221 (row lengths), A001222 (row sums). A118914 gives the nondecreasing version. A124010 lists exponents in n's prime factorization in natural order, with A124010(1) = 0.
A212172 cross-references over 20 sequences that depend solely on n's prime exponents >= 2, including the "core" sequence A000688.  Other sequences determined by the exponents in the prime factorization of n include:
Multiplicative: A000005, A007425, A008683, A008836, A034444, A037445, A181819.
Additive: A001221, A001222, A056169.
Other: A001055, A008480, A010553, A038548, A050320, A051707, A071625, A074206, A076078, A085082, A088873.
A highly incomplete selection of sequences, each definable by the set of prime signatures possessed by its members: A000040, A000290, A000578, A000583, A000961, A001248, A001358, A001597, A001694, A002808, A004709, A005117, A006881, A013929, A030059, A030229, A052486.

&cat[Reverse(Sort([pe[2]:pe in Factorisation(n)])):n in[1..76]]; // Jason Kimberley, Jun 13 2012

apply( {A212171_row(n)=vecsort(factor(n)[,2]~,,4)}, [1..40])\\ M. F. Hasler, Apr 19 2022

Row n of A118914, reversed.

Row n of A124010 for n > 1, with exponents sorted in nonincreasing order. Equivalently, row A046523(n) of A124010 for n > 1.

A303838 Number of z-forests with least common multiple n > 1.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 3, 1, 2, 2, 1, 1, 3, 1, 3, 2, 2, 1, 4, 1, 2, 1, 3, 1, 8, 1, 1, 2, 2, 2, 5, 1, 2, 2, 4, 1, 8, 1, 3, 3, 2, 1, 5, 1, 3, 2, 3, 1, 4, 2, 4, 2, 2, 1, 16, 1, 2, 3, 1, 2, 8, 1, 3, 2, 8, 1, 7, 1, 2, 3, 3, 2, 8, 1, 5, 1, 2, 1, 16, 2, 2

Offset: 1

Gus Wiseman, May 19 2018

nonn

Given a finite set S of positive integers greater than 1, let G(S) be the simple labeled graph with vertex set S and edges between any two vertices that have a common divisor greater than 1. For example, G({6,14,15,35}) is a 4-cycle. A set S is said to be connected if G(S) is a connected graph. The clutter density of S is defined to be Sum_{s in S} (omega(s) - 1) - omega(LCM(S)), where omega = A001221 and LCM is least common multiple. A z-forest is a finite set of pairwise indivisible positive integers greater than 1 such that all connected components are z-trees, meaning they have clutter density -1.
This is a generalization to multiset systems of the usual definition of hyperforest (viz. hypergraph F such that two distinct hyperedges of F intersect in at most a common vertex and such that every cycle of F is contained in a hyperedge).
If n is squarefree with k prime factors, then a(n) = A134954(k).
Differs from A324837 at positions {1, 180, 210, ...}. For example, a(210) = 55, A324837(210) = 49.

			The a(60) = 16 z-forests together with the corresponding multiset systems (see A112798, A302242) are the following.
       (60): {{1,1,2,3}}
     (3,20): {{2},{1,1,3}}
     (4,15): {{1,1},{2,3}}
     (4,30): {{1,1},{1,2,3}}
     (5,12): {{3},{1,1,2}}
     (6,20): {{1,2},{1,1,3}}
    (10,12): {{1,3},{1,1,2}}
    (12,15): {{1,1,2},{2,3}}
    (12,20): {{1,1,2},{1,1,3}}
    (15,20): {{2,3},{1,1,3}}
    (3,4,5): {{2},{1,1},{3}}
   (3,4,10): {{2},{1,1},{1,3}}
    (4,5,6): {{1,1},{3},{1,2}}
   (4,6,10): {{1,1},{1,2},{1,3}}
   (4,6,15): {{1,1},{1,2},{2,3}}
  (4,10,15): {{1,1},{1,3},{2,3}}

Gus Wiseman, Table of n, a(n) for n = 1..250
Roland Bacher, On the enumeration of labelled hypertrees and of labelled bipartite trees, arXiv:1102.2708 [math.CO], 2011.

Cf. A006126, A030019, A048143, A076078, A112798, A134954, A275307, A285572, A286518, A286520, A293993, A293994, A302242, A303837, A304118.

zsm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[Less@@#,GCD@@s[[#]]]>1&]},If[c=={},s,zsm[Union[Append[Delete[s,List/@c[[1]]],LCM@@s[[c[[1]]]]]]]]];
zensity[s_]:=Total[(PrimeNu[#]-1&)/@s]-PrimeNu[LCM@@s];
Table[Length[Select[Rest[Subsets[Rest[Divisors[n]]]],Function[s,LCM@@s==n&&And@@Table[zensity[Select[s,Divisible[m,#]&]]==-1,{m,zsm[s]}]&&Select[Tuples[s,2],UnsameQ@@#&&Divisible@@#&]=={}]]],{n,100}]

A324744 Number of maximal subsets of {1...n} containing no element whose prime indices all belong to the subset.

Original entry on oeis.org

1, 1, 2, 2, 3, 4, 4, 5, 6, 8, 8, 11, 11, 22, 22, 22, 22, 28, 28, 44, 44, 52, 52, 76, 76, 88, 88, 96, 96, 184, 184, 240, 240, 264, 264, 296, 296, 592, 592, 592, 592, 728, 728, 1456, 1456, 1456, 1456, 2912, 2912, 3168, 3168, 3168, 3168, 5568, 5568, 5568, 5568

Offset: 0

Gus Wiseman, Mar 15 2019

nonn

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 a(1) = 1 through a(8) = 6 maximal subsets:
  {1}  {1}  {2}    {1,3}  {1,3}    {1,3,6}    {3,4,6}    {1,3,6,7}
       {2}  {1,3}  {2,4}  {1,5}    {1,5,6}    {1,3,6,7}  {1,5,6,7}
                   {3,4}  {3,4}    {3,4,6}    {1,5,6,7}  {3,4,6,8}
                          {2,4,5}  {2,4,5,6}  {2,4,5,6}  {3,6,7,8}
                                              {2,5,6,7}  {2,4,5,6,8}
                                                         {2,5,6,7,8}

Andrew Howroyd, Table of n, a(n) for n = 0..100

The non-maximal case is A324738. The case for subsets of {2...n} is A324762.
Cf. A000720, A001462, A007097, A076078, A084422, A085945, A112798, A276625, A290822, A304360, A306844, A320426, A324764.
Cf. A324694, A324736, A324743, A324749, A324754, A324759.

maxim[s_]:=Complement[s,Last/@Select[Tuples[s,2],UnsameQ@@#&&SubsetQ@@#&]];
Table[Length[maxim[Select[Subsets[Range[n]],!MemberQ[#,k_/;SubsetQ[#,PrimePi/@First/@FactorInteger[k]]]&]]],{n,0,10}]

pset(n)={my(b=0, f=factor(n)[, 1]); sum(i=1, #f, 1<<(primepi(f[i])))}
a(n)={my(p=vector(n, k, if(k==1, 1, pset(k))), d=0); for(i=1, #p, d=bitor(d, p[i]));
my(ismax(b)=for(k=1, #p, if(!bittest(b,k) && bitnegimply(p[k], b), my(e=bitor(b, 1<#p, ismax(b), my(f=bitnegimply(p[k], b)); if(!f || bittest(d, k), self()(k+1, b)) + if(f, self()(k+1, b+(1<Andrew Howroyd, Aug 27 2019

Terms a(16) and beyond from Andrew Howroyd, Aug 27 2019

A304717 Number of connected strict integer partitions of n with pairwise indivisible parts.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 2, 2, 1, 3, 2, 4, 3, 5, 2, 5, 4, 6, 3, 7, 6, 9, 5, 9, 8, 13, 10, 15, 9, 15, 13, 18, 14, 22, 21, 26, 19, 29, 24, 36, 31, 40, 35, 45, 38, 54, 55, 59, 55, 70, 69, 84, 74, 89, 86, 107, 103, 119, 115, 143, 143, 159

Offset: 1

Gus Wiseman, May 17 2018

nonn

Given a finite set S of positive integers greater than one, let G(S) be the simple labeled graph with vertex set S and edges between any two vertices with a common divisor. For example, G({6,14,15,35}) is a 4-cycle. A multiset S is said to be connected if G(S) is a connected graph.

			The a(34) = 13 connected strict integer partitions with pairwise indivisible parts are (34), (18,16), (20,14), (22,12), (24,10), (26,8), (28,6), (30,4), (14,12,8), (15,10,9), (20,8,6), (14,10,6,4), (15,9,6,4). Their corresponding multiset multisystems (see A112798, A302242) are the following.
         (34): {{1,7}}
       (30 4): {{1,2,3},{1,1}}
       (28 6): {{1,1,4},{1,2}}
       (26 8): {{1,6},{1,1,1}}
      (24 10): {{1,1,1,2},{1,3}}
      (22 12): {{1,5},{1,1,2}}
      (20 14): {{1,1,3},{1,4}}
     (20 8 6): {{1,1,3},{1,1,1},{1,2}}
      (18 16): {{1,2,2},{1,1,1,1}}
    (15 10 9): {{2,3},{1,3},{2,2}}
   (15 9 6 4): {{2,3},{2,2},{1,2},{1,1}}
    (14 12 8): {{1,4},{1,1,2},{1,1,1}}
  (14 10 6 4): {{1,4},{1,3},{1,2},{1,1}}

Cf. A000009, A003963, A006126, A048143, A051424, A054921, A076078, A259936, A281116, A285572, A285573, A286518, A286520, A293993, A302242, A303362, A304714, A304716.

zsm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[Less@@#,GCD@@s[[#]]]>1&]},If[c==={},s,zsm[Union[Append[Delete[s,List/@c[[1]]],LCM@@s[[c[[1]]]]]]]]];
Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&Length[zsm[#]]===1&&Select[Tuples[#,2],UnsameQ@@#&&Divisible@@#&]==={}&]],{n,30}]

A333226 Least common multiple of the n-th composition in standard order.

Original entry on oeis.org

1, 2, 1, 3, 2, 2, 1, 4, 3, 2, 2, 3, 2, 2, 1, 5, 4, 6, 3, 6, 2, 2, 2, 4, 3, 2, 2, 3, 2, 2, 1, 6, 5, 4, 4, 3, 6, 6, 3, 4, 6, 2, 2, 6, 2, 2, 2, 5, 4, 6, 3, 6, 2, 2, 2, 4, 3, 2, 2, 3, 2, 2, 1, 7, 6, 10, 5, 12, 4, 4, 4, 12, 3, 6, 6, 3, 6, 6, 3, 10, 4, 6, 6, 6, 2, 2

Offset: 1

Gus Wiseman, Mar 26 2020

nonn

The k-th composition in standard order (row k of A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again.

The version for binary indices is A271410.
The version for prime indices is A290103.
Positions of first appearances are A333225.
Let q(k) be the k-th composition in standard order:
- The terms of q(k) are row k of A066099.
- The sum of q(k) is A070939(k).
- The product of q(k) is A124758(k).
- The GCD of q(k) is A326674(k).
- The LCM of q(k) is A333226(k).
Cf. A000120, A029931, A048793, A074971, A076078, A233564, A285572, A289508, A289509, A324837, A333227, A333492.

stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
Table[LCM@@stc[n],{n,100}]

A324738 Number of subsets of {1...n} containing no element > 1 whose prime indices all belong to the subset.

Original entry on oeis.org

1, 2, 3, 5, 8, 13, 26, 42, 72, 120, 232, 376, 752, 1128, 2256, 4512, 8256, 13632, 27264, 42048, 82944, 158976, 313344, 497664, 995328, 1700352, 3350016, 5815296, 11630592, 17491968, 34983936, 56954880, 108933120, 210788352, 418258944, 804667392, 1609334784

Offset: 0

Gus Wiseman, Mar 13 2019

nonn

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 a(0) = 1 through a(6) = 26 subsets:
  {}  {}   {}   {}     {}     {}       {}
      {1}  {1}  {1}    {1}    {1}      {1}
           {2}  {2}    {2}    {2}      {2}
                {3}    {3}    {3}      {3}
                {1,3}  {4}    {4}      {4}
                       {1,3}  {5}      {5}
                       {2,4}  {1,3}    {6}
                       {3,4}  {1,5}    {1,3}
                              {2,4}    {1,5}
                              {2,5}    {1,6}
                              {3,4}    {2,4}
                              {4,5}    {2,5}
                              {2,4,5}  {2,6}
                                       {3,4}
                                       {3,6}
                                       {4,5}
                                       {4,6}
                                       {5,6}
                                       {1,3,6}
                                       {1,5,6}
                                       {2,4,5}
                                       {2,4,6}
                                       {2,5,6}
                                       {3,4,6}
                                       {4,5,6}
                                       {2,4,5,6}

Andrew Howroyd, Table of n, a(n) for n = 0..100

The maximal case is A324744. The case of subsets of {2...n} is A324739. The strict integer partition version is A324749. The integer partition version is A324754. The Heinz number version is A324759. An infinite version is A324694.
Cf. A000720, A001221, A001462, A007097, A076078, A084422, A085945, A112798, A276625, A279861, A290689, A290822, A304360, A306844.
Cf. A324695, A324736, A324741, A324750, A324755, A324760.

Table[Length[Select[Subsets[Range[n]],!MemberQ[#,k_/;SubsetQ[#,PrimePi/@First/@FactorInteger[k]]]&]],{n,0,10}]

pset(n)={my(b=0,f=factor(n)[,1]); sum(i=1, #f, 1<<(primepi(f[i])))}
a(n)={my(p=vector(n,k,if(k==1, 1, pset(k))), d=0); for(i=1, #p, d=bitor(d, p[i]));
((k,b)->if(k>#p, 1, my(t=self()(k+1,b)); if(bitnegimply(p[k], b), t+=if(bittest(d,k), self()(k+1, b+(1<Andrew Howroyd, Aug 16 2019

Terms a(21) and beyond from Andrew Howroyd, Aug 16 2019

A324742 Number of subsets of {2...n} containing no prime indices of the elements.

Original entry on oeis.org

1, 2, 3, 6, 10, 16, 24, 48, 84, 144, 228, 420, 648, 1080, 1800, 3600, 5760, 11136, 16704, 31104, 53568, 90624, 136896, 269952, 515712, 862080, 1708800, 3171840, 4832640, 9325440, 14890752, 29781504, 52245504, 88418304, 166017024, 331628544, 497645568, 829409280

Offset: 1

Gus Wiseman, Mar 15 2019

nonn

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 a(1) = 1 through a(6) = 16 subsets:
  {}  {}   {}   {}     {}       {}
      {2}  {2}  {2}    {2}      {2}
           {3}  {3}    {3}      {3}
                {4}    {4}      {4}
                {2,4}  {5}      {5}
                {3,4}  {2,4}    {6}
                       {2,5}    {2,4}
                       {3,4}    {2,5}
                       {4,5}    {3,4}
                       {2,4,5}  {3,6}
                                {4,5}
                                {4,6}
                                {5,6}
                                {2,4,5}
                                {3,4,6}
                                {4,5,6}
An example for n = 20 is {4,5,6,12,17,18,19}, with prime indices:
   4: {1,1}
   5: {3}
   6: {1,2}
  12: {1,1,2}
  17: {7}
  18: {1,2,2}
  19: {8}
None of these prime indices {1,2,3,7,8} belong to the set, as required.

Andrew Howroyd, Table of n, a(n) for n = 1..100

The maximal case is A324763. The version for subsets of {1...n} is A324741. The strict integer partition version is A324752. The integer partition version is A324757. The Heinz number version is A324761. An infinite version is A304360.
Cf. A000720, A001462, A007097, A076078, A084422, A085945, A112798, A276625, A290689, A290822, A306844, A324764.
Cf. A324695, A324737, A324743, A324751, A324756, A324758.

Table[Length[Select[Subsets[Range[2,n]],Intersection[#,PrimePi/@First/@Join@@FactorInteger/@#]=={}&]],{n,10}]

pset(n)={my(b=0,f=factor(n)[,1]); sum(i=1, #f, 1<<(primepi(f[i])))}
a(n)={my(p=vector(n-1,k,pset(k+1)>>1), d=0); for(i=1, #p, d=bitor(d, p[i]));
((k,b)->if(k>#p, 1, my(t=self()(k+1,b)); if(!bitand(p[k], b), t+=if(bittest(d,k), self()(k+1, b+(1<Andrew Howroyd, Aug 16 2019

Terms a(21) and beyond from Andrew Howroyd, Aug 16 2019

A305566 Number of finite sets of relatively prime positive integers > 1 with least common multiple n.

Original entry on oeis.org

0, 0, 0, 0, 0, 2, 0, 0, 0, 2, 0, 10, 0, 2, 2, 0, 0, 10, 0, 10, 2, 2, 0, 44, 0, 2, 0, 10, 0, 84, 0, 0, 2, 2, 2, 122, 0, 2, 2, 44, 0, 84, 0, 10, 10, 2, 0, 184, 0, 10, 2, 10, 0, 44, 2, 44, 2, 2, 0, 1590, 0, 2, 10, 0, 2, 84, 0, 10, 2, 84, 0, 1156, 0, 2, 10, 10, 2

Offset: 1

Gus Wiseman, Jun 05 2018

nonn

From Robert Israel, Jun 06 2018: (Start)
a(n) depends only on the prime signature of n.
If n is in A000961, a(n)=0.
If n is in A006881, a(n)=2. (End)
If n = p^k*q, where p and q are distinct primes and k >= 1, then a(n) = 3*4^(k-1)-2^(k-1). - Robert Israel, Jun 07 2018

			The a(12) = 10 sets:
{3,4},
{2,3,4}, {2,3,12}, {3,4,6}, {3,4,12},
{2,3,4,6}, {2,3,4,12}, {2,3,6,12}, {3,4,6,12},
{2,3,4,6,12}.

Robert Israel, Table of n, a(n) for n = 1..10000

Cf. A000961, A006881, A076078, A181819, A281116, A285572, A290103, A304818, A305563, A305564, A305565, A305567.

f:= proc(n) g(sort(map(t -> t[2],ifactors(n)[2]))) end proc:
f(1):= 0:
g:= proc(L) option remember;
  local nL, Cands, nC, Cons, i;
  nL:= nops(L);
  Cands:= [[]];
  for i from 1 to nL do
    Cands:= [seq(seq([op(s),t],t=0..L[i]),s=Cands)];
  od:
  Cands:= remove(t -> max(t)=0, Cands);
  nC:= nops(Cands);
  Cons:= [seq(select(t -> Cands[t][i]=0, {$1..nC}),i=1..nL),
          seq(select(t -> Cands[t][i]=L[i], {$1..nC}), i=1..nL)];
  h(Cons, {$1..nC})
end proc:
h:= proc(Cons, Cands)
  local t,i,Consi, Candsi;
  if Cons = [] then return 2^nops(Cands) fi;
  t:= 0;
  for i from 1 to nops(Cons[1]) do
    Consi:= map(proc(t) if member(Cons[1][i],t) then NULL else t minus Cons[1][1..i-1] fi end proc, Cons[2..-1]);
    if member({},Consi) then next fi;
    Candsi:= Cands minus Cons[1][1..i];
    t:= t + procname(Consi, Candsi)
  od;
  t
end proc:
map(f, [$1..100]); # Robert Israel, Jun 07 2018

Table[Length[Select[Subsets[Rest[Divisors[n]]],And[GCD@@#==1,LCM@@#==n]&]],{n,100}]

A225520 The number of subsets of the set of divisors of n in which elements are pairwise coprime.

Original entry on oeis.org

2, 4, 4, 6, 4, 10, 4, 8, 6, 10, 4, 16, 4, 10, 10, 10, 4, 16, 4, 16, 10, 10, 4, 22, 6, 10, 8, 16, 4, 30, 4, 12, 10, 10, 10, 26, 4, 10, 10, 22, 4, 30, 4, 16, 16, 10, 4, 28, 6, 16, 10, 16, 4, 22, 10, 22, 10, 10, 4, 50, 4, 10, 16, 14, 10, 30, 4, 16, 10, 30, 4, 36

Offset: 1

R. J. Mathar, May 09 2013

nonn

Note that this is not 1+A048691(n); n=30 is a counterexample.

The number of all subsets of the set of divisors (without the restriction) is 2^A000005(n), which therefore is an upper bound of the current sequence.

			For n=6, the set of divisors is {1,2,3,6} and the a(6)=10 subsets with pairwise coprime entries are {}, {1}, {2}, {3}, {6}, {1,2}, {1,3}, {1,6}, {2,3} and {1,2,3}.

T. D. Noe, Table of n, a(n) for n = 1..359

Cf. A076078 (subsets with lcm equal to n), A084422 (subsets of 1 through n).

paircoprime := proc(s)
    local L,i,j ;
    L := convert(s,list) ;
    for i from 1 to nops(L)-1 do
        for j from i+1 to nops(L) do
            if igcd(op(i,L),op(j,L)) <> 1 then
                return false;
            end if;
        end do:
    end do:
    return true;
end proc:
A225520 := proc(n)
    local dvs,a,p ;
    dvs := numtheory[divisors](n) ;
    a := 0 ;
    for p in combinat[powerset](dvs) do
        if paircoprime(p) then
            a := a+1 ;
        end if;
    end do:
    a ;
end proc:

Table[Length[Select[Subsets[Divisors[n]], If[Length[#] < 2, True, If[Length[#] == 2, CoprimeQ @@ #, And @@ CoprimeQ @@ #]] &]], {n, 100}] (* T. D. Noe, May 09 2013 *)

A324763 Number of maximal subsets of {2...n} containing no prime indices of the elements.

Original entry on oeis.org

1, 1, 2, 2, 2, 3, 6, 6, 6, 6, 10, 10, 16, 16, 16, 16, 24, 24, 48, 48, 48, 48, 84, 84, 84, 84, 84, 84, 144, 144, 228, 228, 228, 228, 228, 228, 420, 420, 420, 420, 648, 648, 1080, 1080, 1080, 1080, 1800, 1800, 1800, 1800, 1800, 1800, 3600, 3600, 3600, 3600, 3600

Offset: 1

Gus Wiseman, Mar 17 2019

nonn

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 a(1) = 1 through a(9) = 6 subsets:
  {}  {2}  {2}  {2,4}  {3,4}    {2,4,5}  {2,4,5}  {2,4,5,8}  {2,4,5,8}
           {3}  {3,4}  {2,4,5}  {3,4,6}  {2,5,7}  {2,5,7,8}  {2,5,7,8}
                                {4,5,6}  {3,4,6}  {3,4,6,8}  {3,4,6,8,9}
                                         {3,6,7}  {3,6,7,8}  {3,6,7,8,9}
                                         {4,5,6}  {4,5,6,8}  {4,5,6,8,9}
                                         {5,6,7}  {5,6,7,8}  {5,6,7,8,9}

Andrew Howroyd, Table of n, a(n) for n = 1..100

The non-maximal version is A324742.
The version for subsets of {1...n} is A324741.
An infinite version is A304360.
Cf. A076078, A084422, A085945, A112798, A276625, A306844, A324764.
Cf. A324695, A324743, A324751, A324756, A324758, A324762.

maxim[s_]:=Complement[s,Last/@Select[Tuples[s,2],UnsameQ@@#&&SubsetQ@@#&]];
Table[Length[maxim[Select[Subsets[Range[2,n]],Intersection[#,PrimePi/@First/@Join@@FactorInteger/@#]=={}&]]],{n,10}]

pset(n)={my(b=0, f=factor(n)[, 1]); sum(i=1, #f, 1<<(primepi(f[i])))}
a(n)={my(p=vector(n-1, k, pset(k+1)>>1), d=0); for(i=1, #p, d=bitor(d, p[i]));
my(ismax(b)=my(e=0); forstep(k=#p, 1, -1, if(bittest(b,k), e=bitor(e,p[k]), if(!bittest(e,k) && !bitand(p[k], b), return(0)) )); 1);
((k, b)->if(k>#p, ismax(b), my(f=!bitand(p[k], b)); if(!f || bittest(d, k), self()(k+1, b)) + if(f, self()(k+1, b+(1<Andrew Howroyd, Aug 26 2019

Terms a(16) and beyond from Andrew Howroyd, Aug 26 2019

cp's OEIS Frontend

A212171 Prime signature of n (nonincreasing version): row n of table lists positive exponents in canonical prime factorization of n, in nonincreasing order.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

PARI

Formula

A303838 Number of z-forests with least common multiple n > 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A324744 Number of maximal subsets of {1...n} containing no element whose prime indices all belong to the subset.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A304717 Number of connected strict integer partitions of n with pairwise indivisible parts.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A333226 Least common multiple of the n-th composition in standard order.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

A324738 Number of subsets of {1...n} containing no element > 1 whose prime indices all belong to the subset.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A324742 Number of subsets of {2...n} containing no prime indices of the elements.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A305566 Number of finite sets of relatively prime positive integers > 1 with least common multiple n.

Views

Author

Keywords