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 50 results. Next

A326083 Number of subsets of {1..n} containing all of their pairwise sums <= n.

Original entry on oeis.org

1, 2, 3, 5, 7, 12, 16, 27, 37, 58, 80, 131, 171, 277, 380, 580, 785, 1250, 1655, 2616, 3516, 5344, 7257, 11353, 14931, 23204, 31379, 47511, 63778, 98681, 130503, 201357, 270038, 407429, 548090, 840171, 1110429, 1701872, 2284325, 3440337, 4601656
Offset: 0

Views

Author

Gus Wiseman, Jun 05 2019

Keywords

Comments

The summands are allowed to be equal. The case where they must be distinct is A326080. If A007865 counts sum-free sets, this sequence counts sum-closed sets. This is different from sum-full sets (A093971).
From Gus Wiseman, Jul 08 2019: (Start)
Also the number of subsets of {1..n} containing no sum of any multiset of the elements. For example, the a(0) = 1 through a(6) = 16 subsets are:
{} {} {} {} {} {} {}
{1} {1} {1} {1} {1} {1}
{2} {2} {2} {2} {2}
{3} {3} {3} {3}
{2,3} {4} {4} {4}
{2,3} {5} {5}
{3,4} {2,3} {6}
{2,5} {2,3}
{3,4} {2,5}
{3,5} {3,4}
{4,5} {3,5}
{3,4,5} {4,5}
{4,6}
{5,6}
{3,4,5}
{4,5,6}
(End)

Examples

			The a(0) = 1 through a(6) = 16 subsets:
  {}  {}   {}     {}       {}         {}           {}
      {1}  {2}    {2}      {3}        {3}          {4}
           {1,2}  {3}      {4}        {4}          {5}
                  {2,3}    {2,4}      {5}          {6}
                  {1,2,3}  {3,4}      {2,4}        {3,6}
                           {2,3,4}    {3,4}        {4,5}
                           {1,2,3,4}  {3,5}        {4,6}
                                      {4,5}        {5,6}
                                      {2,4,5}      {2,4,6}
                                      {3,4,5}      {3,4,6}
                                      {2,3,4,5}    {3,5,6}
                                      {1,2,3,4,5}  {4,5,6}
                                                   {2,4,5,6}
                                                   {3,4,5,6}
                                                   {2,3,4,5,6}
                                                   {1,2,3,4,5,6}
The a(7) = 27 subsets:
  {}  {4}  {36}  {246}  {2467}  {24567}  {234567}  {1234567}
      {5}  {45}  {356}  {3467}  {34567}
      {6}  {46}  {367}  {3567}
      {7}  {47}  {456}  {4567}
           {56}  {457}
           {57}  {467}
           {67}  {567}
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[Subsets[Range[n]],SubsetQ[#,Select[Plus@@@Tuples[#,2],#<=n&]]&]],{n,0,10}]

Formula

For n > 0, a(n) = A103580(n) + 1.

A093971 Number of sum-full subsets of {1,...,n}; subsets A such that there is a solution to x+y=z for x,y,z in A.

Original entry on oeis.org

0, 1, 2, 7, 16, 40, 86, 195, 404, 873, 1795, 3727, 7585, 15537, 31368, 63582, 127933, 257746, 517312, 1038993, 2081696, 4173322, 8355792, 16731799, 33484323, 67014365, 134069494, 268234688, 536562699, 1073326281, 2146849378, 4294117419, 8588623348, 17178130162
Offset: 1

Views

Author

T. D. Noe, Apr 20 2004

Keywords

Comments

In sumset notation, number of subsets A of {1,...,n} such that the intersection of A and 2A is nonempty.
A variation of binary sum-full sets where parts can be re-used, this sequence counts subsets of {1..n} containing a part equal to the sum of two other (possibly equal) parts. The complement is counted by A007865. The non-binary version is A364914. For non-re-usable parts we have A088809. - Gus Wiseman, Aug 14 2023

Examples

			The a(1) = 0 through a(5) = 16 subsets:
  .  {1,2}  {1,2}    {1,2}      {1,2}
            {1,2,3}  {2,4}      {2,4}
                     {1,2,3}    {1,2,3}
                     {1,2,4}    {1,2,4}
                     {1,3,4}    {1,2,5}
                     {2,3,4}    {1,3,4}
                     {1,2,3,4}  {1,4,5}
                                {2,3,4}
                                {2,3,5}
                                {2,4,5}
                                {1,2,3,4}
                                {1,2,3,5}
                                {1,2,4,5}
                                {1,3,4,5}
                                {2,3,4,5}
                                {1,2,3,4,5}
		

Crossrefs

The complement is counted by A007865.
The version without re-usable parts is A088809 (differences A364756), complement A085489 (differences A364755).
The non-binary version is A364914, complement A326083.
The non-binary version w/o re-usable parts is A364534, complement A151897.
The version for partitions is A363225:
- ranks A364348,
- strict A363226,
- non-binary A364839,
- without re-usable parts A237113,
- non-binary without re-usable parts A237668.
The complement for partitions is A364345:
- ranks A364347,
- strict A364346,
- non-binary A364350,
- without re-usable parts A236912,
- non-binary without re-usable parts A237667.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Range[n]],Intersection[#,Total/@Tuples[#,2]]!={}&]],{n,0,10}] (* Gus Wiseman, Aug 14 2023 *)

Formula

a(n) = 2^n - A007865(n).

Extensions

Terms a(31) and beyond from Fausto A. C. Cariboni, Oct 01 2020

A088809 Number of subsets of {1, ..., n} that are not sum-free.

Original entry on oeis.org

0, 0, 0, 1, 3, 10, 27, 67, 154, 350, 763, 1638, 3450, 7191, 14831, 30398, 61891, 125557, 253841, 511818, 1029863, 2069341, 4153060, 8327646, 16687483, 33422562, 66916342, 133936603, 268026776, 536277032, 1072886163, 2146245056, 4293187682, 8587371116
Offset: 0

Views

Author

Reinhard Zumkeller, Oct 19 2003

Keywords

Comments

a(n) = 2^n - A085489(n); a non-sum-free subset contains at least one subset {u,v, w} with w=u+v.
A variation of binary sum-full sets where parts cannot be re-used, this sequence counts subsets of {1..n} with an element equal to the sum of two distinct others. The complement is counted by A085489. The non-binary version is A364534. For re-usable parts we have A093971. - Gus Wiseman, Aug 10 2023

Examples

			From _Gus Wiseman_, Aug 10 2023: (Start)
The set S = {1,3,6,8} has pair-sums {4,7,9,11,14}, which are all missing from S, so it is not counted under a(8).
The set {1,4,6,7} has pair-sum 1 + 6 = 7, so is counted under a(7).
The a(1) = 0 through a(5) = 10 sets:
  .  .  {1,2,3}  {1,2,3}    {1,2,3}
                 {1,3,4}    {1,3,4}
                 {1,2,3,4}  {1,4,5}
                            {2,3,5}
                            {1,2,3,4}
                            {1,2,3,5}
                            {1,2,4,5}
                            {1,3,4,5}
                            {2,3,4,5}
                            {1,2,3,4,5}
(End)
		

Crossrefs

The complement is counted by A085489, differences A364755.
With re-usable parts we have A093971, for partitions A363225.
The complement for partitions is A236912:
non-binary A237667,
ranks A364461,
strict A364533.
The version for partitions is A237113:
non-binary A237668,
ranks A364462,
strict A364670.
The non-binary version is A364534, complement A151897.
First differences are A364756.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Range[n]],Intersection[#,Total/@Subsets[#,{2}]]!={}&]],{n,0,10}] (* Gus Wiseman, Aug 10 2023 *)

Extensions

Terms a(32) and beyond from Fausto A. C. Cariboni, Sep 28 2020

A085489 a(n) is the number of subsets of {1,...,n} containing no solutions to x+y=z with x and y distinct (one version of "sum-free subsets").

Original entry on oeis.org

1, 2, 4, 7, 13, 22, 37, 61, 102, 162, 261, 410, 646, 1001, 1553, 2370, 3645, 5515, 8303, 12470, 18713, 27811, 41244, 60962, 89733, 131870, 192522, 281125, 408680, 593880, 855661, 1238592, 1779614, 2563476, 3660084, 5255913, 7473380, 10696444, 15137517
Offset: 0

Views

Author

Eric W. Weisstein, Jul 02 2003

Keywords

Comments

First differs from A151897 at a(7) = 61, A151897(7) = 60. The one subset counted under a(7) but not under A151897(7) is {1,2,4,7}. - Gus Wiseman, Jun 07 2019

Examples

			From _Gus Wiseman_, Jun 07 2019: (Start)
The a(0) = 1 through a(4) = 13 subsets:
  {}  {}   {}     {}     {}
      {1}  {1}    {1}    {1}
           {2}    {2}    {2}
           {1,2}  {3}    {3}
                  {1,2}  {4}
                  {1,3}  {1,2}
                  {2,3}  {1,3}
                         {1,4}
                         {2,3}
                         {2,4}
                         {3,4}
                         {1,2,4}
                         {2,3,4}
The a(5) = 22 subsets:
  {}  {1}  {1,2}  {1,2,4}
      {2}  {1,3}  {1,2,5}
      {3}  {1,4}  {1,3,5}
      {4}  {1,5}  {2,3,4}
      {5}  {2,3}  {2,4,5}
           {2,4}  {3,4,5}
           {2,5}
           {3,4}
           {3,5}
           {4,5}
(End)
		

Crossrefs

See A007865 for another version.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Range[n]],Intersection[ #,Select[ Plus@@@ Subsets[ #,{2}],#<=n&]]=={}&]],{n,0,10}] (* Gus Wiseman, Jun 07 2019 *)

Formula

a(n) = 2^n - A088809(n). - Reinhard Zumkeller, Oct 19 2003

Extensions

More terms from Reinhard Zumkeller, Jul 13 2003
Edited by David Wasserman, Apr 16 2008
a(0) = 1 prepended by Gus Wiseman, Jun 07 2019

A326080 Number of subsets of {1..n} containing the sum of every subset whose sum is <= n.

Original entry on oeis.org

1, 2, 4, 7, 12, 19, 31, 47, 73, 110, 168, 247, 375, 546, 817, 1193, 1769, 2552, 3791, 5445, 8012, 11517, 16899, 24144, 35391, 50427, 73614, 104984, 152656, 216802, 315689, 447473, 648813, 920163, 1332991, 1884735, 2728020, 3853437, 5568644, 7868096, 11347437
Offset: 0

Views

Author

Gus Wiseman, Jun 05 2019

Keywords

Comments

Equivalently, a(n) is the number of subsets of {1..n} containing the sum of any two distinct elements whose sum is <= n.
The summands must be distinct. The case where they are allowed to be equal is A326083.
If A151897 counts sum-free sets, this sequence counts sum-closed sets. This is different from sum-full sets (A093971).

Examples

			The a(0) = 1 through a(5) = 19 subsets:
  {}  {}   {}     {}       {}         {}
      {1}  {1}    {1}      {1}        {1}
           {2}    {2}      {2}        {2}
           {1,2}  {3}      {3}        {3}
                  {1,3}    {4}        {4}
                  {2,3}    {1,4}      {5}
                  {1,2,3}  {2,3}      {1,5}
                           {2,4}      {2,4}
                           {3,4}      {2,5}
                           {1,3,4}    {3,4}
                           {2,3,4}    {3,5}
                           {1,2,3,4}  {4,5}
                                      {1,4,5}
                                      {2,3,5}
                                      {2,4,5}
                                      {3,4,5}
                                      {1,3,4,5}
                                      {2,3,4,5}
                                      {1,2,3,4,5}
The a(6) = 31 subsets:
  {}  {1}  {1,6}  {1,5,6}  {1,4,5,6}  {1,3,4,5,6}  {1,2,3,4,5,6}
      {2}  {2,5}  {2,3,5}  {2,3,5,6}  {2,3,4,5,6}
      {3}  {2,6}  {2,4,6}  {2,4,5,6}
      {4}  {3,4}  {2,5,6}  {3,4,5,6}
      {5}  {3,5}  {3,4,5}
      {6}  {3,6}  {3,4,6}
           {4,5}  {3,5,6}
           {4,6}  {4,5,6}
           {5,6}
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[Subsets[Range[n]],SubsetQ[#,Select[Plus@@@Subsets[#,{2}],#<=n&]]&]],{n,0,10}]
  • PARI
    a(n)={
       my(recurse(k, b)=
          if( k > n, 1,
              my(t=self()(k + 1, b + (1<Andrew Howroyd, Aug 30 2019

Extensions

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

A103580 Number of nonempty subsets S of {1,2,3,...,n} that have the property that no element x of S is a nonnegative integer linear combination of elements of S-{x}.

Original entry on oeis.org

1, 2, 4, 6, 11, 15, 26, 36, 57, 79, 130, 170, 276, 379, 579, 784, 1249, 1654, 2615, 3515, 5343, 7256, 11352, 14930, 23203, 31378, 47510, 63777, 98680, 130502, 201356, 270037, 407428, 548089, 840170, 1110428, 1701871, 2284324, 3440336, 4601655
Offset: 1

Views

Author

Jeffrey Shallit, Mar 23 2005

Keywords

Examples

			a(4) = 6 because the only permissible subsets are {1}, {2}, {3}, {4}, {2,3}, {3,4}.
From _Gus Wiseman_, Jun 07 2019: (Start)
The a(1) = 1 through a(6) = 15 nonempty subsets of {1..n} containing none of their own non-singleton nonzero nonnegative linear combinations are:
  {1}  {1}  {1}    {1}    {1}      {1}
       {2}  {2}    {2}    {2}      {2}
            {3}    {3}    {3}      {3}
            {2,3}  {4}    {4}      {4}
                   {2,3}  {5}      {5}
                   {3,4}  {2,3}    {6}
                          {2,5}    {2,3}
                          {3,4}    {2,5}
                          {3,5}    {3,4}
                          {4,5}    {3,5}
                          {3,4,5}  {4,5}
                                   {4,6}
                                   {5,6}
                                   {3,4,5}
                                   {4,5,6}
a(n) is also the number of nonempty subsets of {1..n} containing all of their own nonzero nonnegative linear combinations <= n. For example the a(1) = 1 through a(6) = 15 subsets are:
  {1}  {2}    {2}      {3}        {3}          {4}
       {1,2}  {3}      {4}        {4}          {5}
              {2,3}    {2,4}      {5}          {6}
              {1,2,3}  {3,4}      {2,4}        {3,6}
                       {2,3,4}    {3,4}        {4,5}
                       {1,2,3,4}  {3,5}        {4,6}
                                  {4,5}        {5,6}
                                  {2,4,5}      {2,4,6}
                                  {3,4,5}      {3,4,6}
                                  {2,3,4,5}    {3,5,6}
                                  {1,2,3,4,5}  {4,5,6}
                                               {2,4,5,6}
                                               {3,4,5,6}
                                               {2,3,4,5,6}
                                               {1,2,3,4,5,6}
(End)
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[Subsets[Range[n],{1,n}],SubsetQ[#,Select[Plus@@@Tuples[#,2],#<=n&]]&]],{n,10}] (* Gus Wiseman, Jun 07 2019 *)

Formula

a(n) = A326083(n) - 1. - Gus Wiseman, Jun 07 2019

Extensions

More terms from David Wasserman, Apr 16 2008

A084422 Number of subsets of integers 1 through n (including the empty set) containing no pair of integers that share a common factor.

Original entry on oeis.org

1, 2, 4, 8, 12, 24, 28, 56, 72, 104, 116, 232, 248, 496, 544, 616, 728, 1456, 1520, 3040, 3232, 3616, 3872, 7744, 8000, 11168, 11904, 14656, 15488, 30976, 31232, 62464, 69888, 76160, 80256, 89856, 91648, 183296, 192640, 208640, 214272, 428544
Offset: 0

Views

Author

Matthew Vandermast, Jun 26 2003

Keywords

Comments

Also the number of subsets of {1,...,n} whose product of elements is equal to the least common multiple of elements. - Michel Marcus, Mar 27 2016

Examples

			Exactly 4 of the 2^4=16 subsets of the integers from 1 through 4 contain a pair of integers that share a common factor; these are {2,4}, {1,2,4}, {2,3,4} and {1,2,3,4}. The other 12 subsets do not; hence a(4)=12.
		

References

  • Alan Sutcliffe, Divisors and Common Factors in Sets of Integers, awaiting publication. [Apparently unpublished as of 2016]

Crossrefs

Cf. A051026 gives the number of primitive subsets. A087080 gives the number of elements in coprime subsets. A087081 gives the sum of the elements in coprime subsets.

Programs

  • Mathematica
    Prepend[Table[Length@ Select[Rest@ Subsets@ Range@ n, Times @@ # == LCM @@ # &], {n, 22}] + 1, 1] (* Michael De Vlieger, Mar 27 2016 *)
  • PARI
    a(n)=nb = 0; S = vector(n, k, k); for (i = 0, 2^n - 1, ss = vecextract(S, i); if (prod(k=1, #ss, ss[k]) == lcm(ss), nb++);); nb; \\ Michel Marcus, Mar 27 2016
    
  • PARI
    a(n,k=1)=if(n<2, return(n+1)); if(gcd(k,n)==1, a(n-1,n*k)) + a(n-1,k) \\ Charles R Greathouse IV, Aug 24 2016

Formula

a(n) = 1 + Sum_{k=1..A036234(n)} A186974(n,k) if n>0; a(0) = 1.

Extensions

More terms from Alan Sutcliffe (alansut(AT)ntlworld.com), Aug 12 2003

A334997 Array T read by ascending antidiagonals: T(n, k) = Sum_{d divides n} T(d, k-1) with T(n, 0) = 1.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 3, 3, 4, 1, 1, 2, 6, 4, 5, 1, 1, 4, 3, 10, 5, 6, 1, 1, 2, 9, 4, 15, 6, 7, 1, 1, 4, 3, 16, 5, 21, 7, 8, 1, 1, 3, 10, 4, 25, 6, 28, 8, 9, 1, 1, 4, 6, 20, 5, 36, 7, 36, 9, 10, 1, 1, 2, 9, 10, 35, 6, 49, 8, 45, 10, 11, 1, 1, 6, 3, 16, 15, 56, 7, 64, 9, 55, 11, 12, 1
Offset: 1

Views

Author

Stefano Spezia, May 19 2020

Keywords

Comments

T(n, k) is called the generalized divisor function (see Beekman).
As an array with offset n=1, k=0, T(n,k) is the number of length-k chains of divisors of n. For example, the T(4,3) = 10 chains are: 111, 211, 221, 222, 411, 421, 422, 441, 442, 444. - Gus Wiseman, Aug 04 2022

Examples

			From _Gus Wiseman_, Aug 04 2022: (Start)
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   9  16  25  36  49  64  81
  n=7:  1   2   3   4   5   6   7   8   9
  n=8:  1   4  10  20  35  56  84 120 165
The T(4,5) = 21 chains:
  (1,1,1,1,1)  (4,2,1,1,1)  (4,4,2,2,2)
  (2,1,1,1,1)  (4,2,2,1,1)  (4,4,4,1,1)
  (2,2,1,1,1)  (4,2,2,2,1)  (4,4,4,2,1)
  (2,2,2,1,1)  (4,2,2,2,2)  (4,4,4,2,2)
  (2,2,2,2,1)  (4,4,1,1,1)  (4,4,4,4,1)
  (2,2,2,2,2)  (4,4,2,1,1)  (4,4,4,4,2)
  (4,1,1,1,1)  (4,4,2,2,1)  (4,4,4,4,4)
The T(6,3) = 16 chains:
  (1,1,1)  (3,1,1)  (6,2,1)  (6,6,1)
  (2,1,1)  (3,3,1)  (6,2,2)  (6,6,2)
  (2,2,1)  (3,3,3)  (6,3,1)  (6,6,3)
  (2,2,2)  (6,1,1)  (6,3,3)  (6,6,6)
The triangular form T(n-k,k) gives the number of length k chains of divisors of n - k. It begins:
  1
  1  1
  1  2  1
  1  2  3  1
  1  3  3  4  1
  1  2  6  4  5  1
  1  4  3 10  5  6  1
  1  2  9  4 15  6  7  1
  1  4  3 16  5 21  7  8  1
  1  3 10  4 25  6 28  8  9  1
  1  4  6 20  5 36  7 36  9 10  1
  1  2  9 10 35  6 49  8 45 10 11  1
(End)
		

References

  • Richard Beekman, An Introduction to Number-Theoretic Combinatorics, Lulu Press 2017.

Crossrefs

Cf. A000217 (4th row), A000290 (6th row), A000292 (8th row), A000332 (16th row), A000389 (32nd row), A000537 (36th row), A000578 (30th row), A002411 (12th row), A002417 (24th row), A007318, A027800 (48th row), A335078, A335079.
Column k = 2 of the array is A007425.
Column k = 3 of the array is A007426.
Column k = 4 of the array is A061200.
The transpose of the array is A077592.
The subdiagonal n = k + 1 of the array is A163767.
The version counting all multisets of divisors (not just chains) is A343658.
The strict case is A343662 (row sums: A337256).
Diagonal n = k of the array is A343939.
Antidiagonal sums of the array (or row sums of the triangle) are A343940.
A067824(n) counts strict chains of divisors starting with n.
A074206(n) counts strict chains of divisors from n to 1.
A146291 counts divisors by Omega.
A251683(n,k) counts strict length k + 1 chains of divisors from n to 1.
A253249(n) counts nonempty chains of divisors of n.
A334996(n,k) counts strict length k chains of divisors from n to 1.
A337255(n,k) counts strict length k chains of divisors starting with n.

Programs

  • Mathematica
    T[n_,k_]:=If[n==1,1,Product[Binomial[Extract[Extract[FactorInteger[n],i],2]+k,k],{i,1,Length[FactorInteger[n]]}]]; Table[T[n-k,k],{n,1,13},{k,0,n-1}]//Flatten
  • PARI
    T(n, k) = if (k==0, 1, sumdiv(n, d, T(d, k-1)));
    matrix(10, 10, n, k, T(n, k-1)) \\ to see the array for n>=1, k >=0; \\ Michel Marcus, May 20 2020

Formula

T(n, k) = Sum_{d divides n} T(d, k-1) with T(n, 0) = 1 (see Theorem 3 in Beekman's article).
T(i*j, k) = T(i, k)*T(j, k) if i and j are coprime positive integers (see Lemma 1 in Beekman's article).
T(p^m, k) = binomial(m+k, k) for every prime p (see Lemma 2 in Beekman's article).

Extensions

Duplicate term removed by Stefano Spezia, Jun 03 2020

A077592 Table by antidiagonals of tau_k(n), the k-th Piltz function (see A007425), or n-th term of the sequence resulting from applying the inverse Möbius transform (k-1) times to the all-ones sequence.

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, 9, 2, 1, 1, 8, 7, 21, 5, 16, 3, 4, 1, 1, 9, 8, 28, 6, 25, 4, 10, 3, 1, 1, 10, 9, 36, 7, 36, 5, 20, 6, 4, 1, 1, 11, 10, 45, 8, 49, 6, 35, 10, 9, 2, 1, 1, 12, 11, 55, 9, 64, 7, 56, 15, 16, 3, 6, 1
Offset: 1

Views

Author

Henry Bottomley, Nov 08 2002

Keywords

Comments

As an array with offset n=0, k=1, also the number of length n chains of divisors of k. - Gus Wiseman, Aug 04 2022

Examples

			T(6,3) = 9 because we have: 1*1*6, 1*2*3, 1*3*2, 1*6*1, 2*1*3, 2*3*1, 3*1*2, 3*2*1, 6*1*1. - _Geoffrey Critzer_, Feb 16 2015
From _Gus Wiseman_, May 03 2021: (Start)
Array begins:
       k=1 k=2 k=3 k=4 k=5 k=6 k=7 k=8
  n=0:  1   1   1   1   1   1   1   1
  n=1:  1   2   2   3   2   4   2   4
  n=2:  1   3   3   6   3   9   3  10
  n=3:  1   4   4  10   4  16   4  20
  n=4:  1   5   5  15   5  25   5  35
  n=5:  1   6   6  21   6  36   6  56
  n=6:  1   7   7  28   7  49   7  84
  n=7:  1   8   8  36   8  64   8 120
  n=8:  1   9   9  45   9  81   9 165
The triangular form T(n,k) = A(n-k,k) gives the number of length n - k chains of divisors of k. It 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  9  2  1
  1  8  7 21  5 16  3  4  1
  1  9  8 28  6 25  4 10  3  1
  1 10  9 36  7 36  5 20  6  4  1
  1 11 10 45  8 49  6 35 10  9  2  1
(End)
		

Crossrefs

Columns include (with multiplicity and some offsets) A000012, A000027, A000027, A000217, A000027, A000290, A000027, A000292, A000217, A000290, A000027, A002411, A000027, A000290, A000290, A000332 etc.
Cf. A077593.
Row n = 2 of the array is A007425.
Row n = 3 of the array is A007426.
Row n = 4 of the array is A061200.
The diagonal n = k of the array (central column of the triangle) is A163767.
The transpose of the array is A334997.
Diagonal n = k of the array is A343939.
Antidiagonal sums of the array (or row sums of the triangle) are A343940.
A067824(n) counts strict chains of divisors starting with n.
A074206(n) counts strict chains of divisors from n to 1.
A146291(n,k) counts divisors of n with k prime factors (with multiplicity).
A251683(n,k) counts strict length k + 1 chains of divisors from n to 1.
A253249(n) counts nonempty chains of divisors of n.
A334996(n,k) counts strict length k chains of divisors from n to 1.
A337255(n,k) counts strict length k chains of divisors starting with n.

Programs

  • Maple
    with(numtheory):
    A:= proc(n,k) option remember; `if`(k=1, 1,
          add(A(d, k-1), d=divisors(n)))
        end:
    seq(seq(A(n, 1+d-n), n=1..d), d=1..14);  # Alois P. Heinz, Feb 25 2015
  • Mathematica
    tau[n_, 1] = 1; tau[n_, k_] := tau[n, k] = Plus @@ (tau[ #, k - 1] & /@ Divisors[n]); Table[tau[n - k + 1, k], {n, 14}, {k, n, 1, -1}] // Flatten (* Robert G. Wilson v *)
    tau[1, k_] := 1; tau[n_, k_] := Times @@ (Binomial[Last[#] + k - 1, k - 1] & /@ FactorInteger[n]); Table[tau[k, n - k + 1], {n, 1, 13}, {k, 1, n}] // Flatten (* Amiram Eldar, Sep 13 2020 *)
    Table[Length[Select[Tuples[Divisors[k],n-k],And@@Divisible@@@Partition[#,2,1]&]],{n,12},{k,1,n}] (* TRIANGLE, Gus Wiseman, May 03 2021 *)
    Table[Length[Select[Tuples[Divisors[k],n-1],And@@Divisible@@@Partition[#,2,1]&]],{n,6},{k,6}] (* ARRAY, Gus Wiseman, May 03 2021 *)

Formula

If n = Product_i p_i^e_i, then T(n,k) = Product_i C(k+e_i-1, e_i). T(n,k) = Sum_d{d|n} T(n-1,d) = A077593(n,k) - A077593(n-1,k).
Columns are multiplicative.
Dirichlet g.f. for column k: Zeta(s)^k. - Geoffrey Critzer, Feb 16 2015
A(n,k) = A334997(k,n). - Gus Wiseman, Aug 04 2022

Extensions

Typo in formula fixed by Geoffrey Critzer, Feb 16 2015

A370584 Number of subsets of {1..n} such that only one set can be obtained by choosing a different prime factor of each element.

Original entry on oeis.org

1, 1, 2, 4, 6, 12, 18, 36, 48, 68, 104, 208, 284, 568, 888, 1296, 1548, 3096, 3968, 7936, 10736, 15440, 24008, 48016, 58848, 73680, 114368, 132608, 176240, 352480, 449824, 899648, 994976, 1399968, 2160720, 2859584, 3296048, 6592096, 10156672, 14214576, 16892352
Offset: 0

Views

Author

Gus Wiseman, Feb 26 2024

Keywords

Comments

For example, the only choice of a different prime factor of each element of (4,5,6) is (2,5,3).

Examples

			The a(0) = 1 through a(6) = 18 subsets:
  {}  {}  {}   {}     {}     {}       {}
          {2}  {2}    {2}    {2}      {2}
               {3}    {3}    {3}      {3}
               {2,3}  {4}    {4}      {4}
                      {2,3}  {5}      {5}
                      {3,4}  {2,3}    {2,3}
                             {2,5}    {2,5}
                             {3,4}    {2,6}
                             {3,5}    {3,4}
                             {4,5}    {3,5}
                             {2,3,5}  {3,6}
                             {3,4,5}  {4,5}
                                      {4,6}
                                      {2,3,5}
                                      {2,5,6}
                                      {3,4,5}
                                      {3,5,6}
                                      {4,5,6}
		

Crossrefs

For divisors instead of factors we have A051026, cf. A368110, A355740.
The version for set-systems is A367904, ranks A367908.
Multisets of this type are ranked by A368101, cf. A368100, A355529.
For existence we have A370582, differences A370586.
For nonexistence we have A370583, differences A370587.
Maximal sets of this type are counted by A370585.
The version for partitions is A370594, cf. A370592, A370593.
For binary indices instead of factors we have A370638, cf. A370636, A370637.
The version for factorizations is A370645, cf. A368414, A368413.
For unlabeled multiset partitions we have A370646, cf. A368098, A368097.
A006530 gives greatest prime factor, least A020639.
A027746 lists prime factors, A112798 indices, length A001222.
A355741 counts ways to choose a prime factor of each prime index.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Range[n]], Length[Union[Sort/@Select[Tuples[If[#==1, {},First/@FactorInteger[#]]&/@#], UnsameQ@@#&]]]==1&]],{n,0,10}]

Extensions

More terms from Jinyuan Wang, Mar 28 2025
Showing 1-10 of 50 results. Next