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.

Previous Showing 21-30 of 75 results. Next

A305149 Number of factorizations of n whose distinct factors are pairwise indivisible and greater than 1.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, May 26 2018

Keywords

Examples

			The a(60) = 8 factorizations are (2*2*3*5), (2*2*15), (3*4*5), (3*20), (4*15), (5*12), (6*10), (60). Missing from this list are (2*3*10), (2*5*6), (2*30).
		

Crossrefs

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[n],Select[Tuples[Union[#],2],UnsameQ@@#&&Divisible@@#&]=={}&]],{n,100}]
  • PARI
    pairwise_indivisible(v) = { for(i=1,#v,for(j=i+1,#v,if(!(v[j]%v[i]),return(0)))); (1); };
    A305149(n, m=n, facs=List([])) = if(1==n, pairwise_indivisible(Set(facs)), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s += A305149(n/d, d, newfacs))); (s)); \\ Antti Karttunen, Oct 08 2018

Extensions

More terms from Antti Karttunen, Oct 08 2018

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

Views

Author

Gus Wiseman, May 19 2018

Keywords

Comments

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.

Examples

			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}}
		

Crossrefs

Programs

  • Mathematica
    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}]

A087086 Primitive sets of integers, each subset mapped onto a unique binary integer, values here shown in decimal.

Original entry on oeis.org

0, 1, 2, 4, 6, 8, 12, 16, 18, 20, 22, 24, 28, 32, 40, 48, 56, 64, 66, 68, 70, 72, 76, 80, 82, 84, 86, 88, 92, 96, 104, 112, 120, 128, 132, 144, 148, 160, 176, 192, 196, 208, 212, 224, 240, 256, 258, 264, 272, 274, 280, 288, 296, 304, 312, 320, 322, 328, 336, 338, 344
Offset: 0

Views

Author

Alan Sutcliffe (alansut(AT)ntlworld.com), Aug 14 2003

Keywords

Comments

A primitive set of integers has no pair of elements one of which divides the other. Each element i in a subset contributes 2^(i-1) to the binary value for that subset. The integers missing from the sequence correspond to nonprimitive subsets.

Examples

			a(10)=22 since the 10th primitive set counting from 0 is {5,3,2}, which maps onto 10110 binary = 22 decimal.
From _Gus Wiseman_, Oct 31 2019: (Start)
The sequence of terms together with their binary expansions and binary indices begins:
   0:       0 ~ {}
   1:       1 ~ {1}
   2:      10 ~ {2}
   4:     100 ~ {3}
   6:     110 ~ {2,3}
   8:    1000 ~ {4}
  12:    1100 ~ {3,4}
  16:   10000 ~ {5}
  18:   10010 ~ {2,5}
  20:   10100 ~ {3,5}
  22:   10110 ~ {2,3,5}
  24:   11000 ~ {4,5}
  28:   11100 ~ {3,4,5}
(End)
		

References

  • Alan Sutcliffe, Divisors and Common Factors in Sets of Integers, awaiting publication

Crossrefs

A051026 gives the number of primitive subsets of the integers 1 to n.
The version for prime indices (rather than binary indices) is A316476.
The relatively prime case is A328671.
Partitions with no consecutive divisible parts are A328171.
Compositions without consecutive divisible parts are A328460.
A ranking of antichains is A326704.

Programs

  • Mathematica
    stableQ[u_,Q_]:=!Apply[Or,Outer[#1=!=#2&&Q[#1,#2]&,u,u,1],{0,1}];
    Select[Range[0,100],stableQ[Join@@Position[Reverse[IntegerDigits[#,2]],1],Divisible]&] (* Gus Wiseman, Oct 31 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

Views

Author

Gus Wiseman, May 17 2018

Keywords

Comments

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.

Examples

			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}}
		

Crossrefs

Programs

  • Mathematica
    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}]

A320632 Numbers k such that there exists a pair of factorizations of k into factors > 1 where no factor of one divides any factor of the other.

Original entry on oeis.org

36, 60, 72, 84, 90, 100, 108, 120, 126, 132, 140, 144, 150, 156, 168, 180, 196, 198, 200, 204, 210, 216, 220, 225, 228, 234, 240, 252, 260, 264, 270, 276, 280, 288, 294, 300, 306, 308, 312, 315, 324, 330, 336, 340, 342, 348, 350, 360, 364, 372, 378, 380, 390
Offset: 1

Views

Author

Gus Wiseman, Dec 09 2018

Keywords

Comments

Positions of nonzero terms in A322437 or A322438.
Mats Granvik has conjectured that these are all the positive integers k such that sigma_0(k) - 2 > (bigomega(k) - 1) * omega(k), where sigma_0 = A000005, omega = A001221, and bigomega = A001222. - Gus Wiseman, Nov 12 2019
Numbers with more semiprime divisors than prime divisors. - Wesley Ivan Hurt, Jun 10 2021

Examples

			An example of such a pair for 36 is (4*9)|(6*6).
		

Crossrefs

The following are additional cross-references relating to Granvik's conjecture.
bigomega(n) * omega(n) is A113901(n).
(bigomega(n) - 1) * omega(n) is A307409(n).
sigma_0(n) - bigomega(n) * omega(n) is A328958(n).
sigma_0(n) - 2 - (omega(n) - 1) * nu(n) is A328959(n).

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Select[Range[100],Select[Subsets[facs[#],{2}],And[!Or@@Divisible@@@Tuples[#],!Or@@Divisible@@@Reverse/@Tuples[#]]&]!={}&]
  • PARI
    factorizations(n, m=n, f=List([]), z=List([])) = if(1==n, listput(z,Vec(f)); z, my(newf); fordiv(n, d, if((d>1)&&(d<=m), newf = List(f); listput(newf,d); z = factorizations(n/d, d, newf, z))); (z));
    is_ndf_pair(fac1,fac2) = { for(i=1,#fac1,for(j=1,#fac2,if(!(fac1[i]%fac2[j])||!(fac2[j]%fac1[i]),return(0)))); (1); };
    has_at_least_one_ndfpair(z) = { for(i=1,#z,for(j=i+1,#z,if(is_ndf_pair(z[i],z[j]),return(1)))); (0); };
    isA320632(n) = has_at_least_one_ndfpair(Vec(factorizations(n))); \\ Antti Karttunen, Dec 10 2020

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

Views

Author

Gus Wiseman, Mar 26 2020

Keywords

Comments

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.

Crossrefs

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).

Programs

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

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).

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

Views

Author

Gus Wiseman, Jun 05 2018

Keywords

Comments

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

Examples

			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}.
		

Crossrefs

Programs

  • Maple
    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
  • Mathematica
    Table[Length[Select[Subsets[Rest[Divisors[n]]],And[GCD@@#==1,LCM@@#==n]&]],{n,100}]

A074971 Number of partitions of n into distinct parts of order n.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 3, 1, 2, 1, 1, 1, 4, 1, 1, 1, 2, 1, 4, 1, 1, 1, 1, 1, 6, 1, 1, 1, 3, 1, 3, 1, 1, 1, 1, 1, 6, 1, 1, 1, 1, 1, 4, 1, 2, 1, 1, 1, 32, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 25, 1, 1, 1, 1, 1, 2, 1, 4, 1, 1, 1, 24, 1, 1, 1, 2, 1, 24, 1, 1, 1, 1, 1, 12, 1, 1, 1, 3, 1, 2
Offset: 1

Views

Author

Vladeta Jovovic, Oct 05 2002

Keywords

Comments

Order of partition is lcm of its parts.

Examples

			The a(36) = 6 partitions are (36), (18,12,6), (18,12,4,2), (18,12,3,2,1), (18,9,4,3,2), (12,9,6,4,3,2). - _Gus Wiseman_, Aug 01 2018
		

Crossrefs

Programs

  • PARI
    A074971(n) = { my(q=0); fordiv(n,i,my(p=1); fordiv(i,j,p *= (1 + 'x^j)); q += moebius(n/i)*p); polcoeff(q,n); }; \\ Antti Karttunen, Dec 19 2018

Formula

Coefficient of x^n in expansion of Sum_{i divides n} mu(n/i)*Product_{j divides i} (1+x^j).
Previous Showing 21-30 of 75 results. Next