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 11-20 of 28 results. Next

A337255 Irregular triangle read by rows where T(n,k) is the number of strict length-k chains of divisors starting with n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 3, 2, 1, 1, 1, 3, 3, 1, 1, 2, 1, 1, 3, 2, 1, 1, 1, 5, 7, 3, 1, 1, 1, 3, 2, 1, 3, 2, 1, 4, 6, 4, 1, 1, 1, 1, 5, 7, 3, 1, 1, 1, 5, 7, 3, 1, 3, 2, 1, 3, 2, 1, 1, 1, 7, 15, 13, 4, 1, 2, 1, 1, 3, 2, 1, 3, 3, 1, 1, 5, 7, 3, 1, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Aug 23 2020

Keywords

Examples

			Sequence of rows begins:
     1: {1}           16: {1,4,6,4,1}
     2: {1,1}         17: {1,1}
     3: {1,1}         18: {1,5,7,3}
     4: {1,2,1}       19: {1,1}
     5: {1,1}         20: {1,5,7,3}
     6: {1,3,2}       21: {1,3,2}
     7: {1,1}         22: {1,3,2}
     8: {1,3,3,1}     23: {1,1}
     9: {1,2,1}       24: {1,7,15,13,4}
    10: {1,3,2}       25: {1,2,1}
    11: {1,1}         26: {1,3,2}
    12: {1,5,7,3}     27: {1,3,3,1}
    13: {1,1}         28: {1,5,7,3}
    14: {1,3,2}       29: {1,1}
    15: {1,3,2}       30: {1,7,12,6}
Row n = 24 counts the following chains:
  24  24/1   24/2/1   24/4/2/1   24/8/4/2/1
      24/2   24/3/1   24/6/2/1   24/12/4/2/1
      24/3   24/4/1   24/6/3/1   24/12/6/2/1
      24/4   24/4/2   24/8/2/1   24/12/6/3/1
      24/6   24/6/1   24/8/4/1
      24/8   24/6/2   24/8/4/2
      24/12  24/6/3   24/12/2/1
             24/8/1   24/12/3/1
             24/8/2   24/12/4/1
             24/8/4   24/12/4/2
             24/12/1  24/12/6/1
             24/12/2  24/12/6/2
             24/12/3  24/12/6/3
             24/12/4
             24/12/6
		

Crossrefs

A008480 gives rows ends.
A067824 gives row sums.
A073093 gives row lengths.
A334996 appears to be the case of chains ending with 1.
A337071 is the sum of row n!.
A000005 counts divisors.
A001055 counts factorizations.
A001222 counts prime factors with multiplicity.
A067824 counts chains of divisors starting with n.
A074206 counts chains of divisors from n to 1.
A122651 counts chains of divisors summing to n.
A167865 counts chains of divisors > 1 summing to n.
A251683 counts chains of divisors from n to 1 by length.
A253249 counts nonempty chains of divisors.
A337070 counts chains of divisors starting with A006939(n).
A337256 counts chains of divisors.

Programs

  • Maple
    b:= proc(n) option remember; expand(x*(1 +
          add(b(d), d=numtheory[divisors](n) minus {n})))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=1..degree(p)))(b(n)):
    seq(T(n), n=1..50);  # Alois P. Heinz, Aug 23 2020
  • Mathematica
    chss[n_]:=Prepend[Join@@Table[Prepend[#,n]&/@chss[d],{d,Most[Divisors[n]]}],{n}];
    Table[Length[Select[chss[n],Length[#]==k&]],{n,30},{k,1+PrimeOmega[n]}]

A337107 Irregular triangle read by rows where T(n,k) is the number of strict length-k chains of divisors from n! to 1.

Original entry on oeis.org

1, 0, 1, 0, 1, 2, 0, 1, 6, 9, 4, 0, 1, 14, 45, 52, 20, 0, 1, 28, 183, 496, 655, 420, 105, 0, 1, 58, 633, 2716, 5755, 6450, 3675, 840, 0, 1, 94, 1659, 11996, 46235, 106806, 155869, 145384, 84276, 27720, 3960
Offset: 1

Views

Author

Gus Wiseman, Aug 23 2020

Keywords

Comments

Row n > 1 appears to be row n! of A334996.

Examples

			Triangle begins:
    1
    0    1
    0    1    2
    0    1    6    9    4
    0    1   14   45   52   20
    0    1   28  183  496  655  420  105
    0    1   58  633 2716 5755 6450 3675  840
Row n = 4 counts the following chains:
  24/1  24/2/1   24/4/2/1   24/8/4/2/1
        24/3/1   24/6/2/1   24/12/4/2/1
        24/4/1   24/6/3/1   24/12/6/2/1
        24/6/1   24/8/2/1   24/12/6/3/1
        24/8/1   24/8/4/1
        24/12/1  24/12/2/1
                 24/12/3/1
                 24/12/4/1
                 24/12/6/1
		

Crossrefs

A097805 is the restriction to powers of 2.
A325617 is the maximal case.
A337105 gives row sums.
A337106 is column k = 3.
A000005 counts divisors.
A000142 lists factorial numbers.
A001055 counts factorizations.
A074206 counts chains of divisors from n to 1.
A027423 counts divisors of factorial numbers.
A067824 counts chains of divisors starting with n.
A076716 counts factorizations of factorial numbers.
A253249 counts chains of divisors.
A337071 counts chains starting with n!.

Programs

  • Maple
    b:= proc(n) option remember; expand(x*(`if`(n=1, 1, 0) +
          add(b(d), d=numtheory[divisors](n) minus {n})))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=1..degree(p)))(b(n!)):
    seq(T(n), n=1..10);  # Alois P. Heinz, Aug 23 2020
  • Mathematica
    nv=5;
    chnsc[n_]:=Select[Prepend[Join@@Table[Prepend[#,n]&/@chnsc[d],{d,DeleteCases[Divisors[n],n]}],{n}],MemberQ[#,1]&];
    Table[Length[Select[chnsc[n!],Length[#]==k&]],{n,nv},{k,1+PrimeOmega[n!]}]

A343662 Irregular triangle read by rows where T(n,k) is the number of strict length k chains of divisors of n, 0 <= k <= Omega(n) + 1.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, May 01 2021

Keywords

Examples

			Triangle begins:
   1:  1  1
   2:  1  2  1
   3:  1  2  1
   4:  1  3  3  1
   5:  1  2  1
   6:  1  4  5  2
   7:  1  2  1
   8:  1  4  6  4  1
   9:  1  3  3  1
  10:  1  4  5  2
  11:  1  2  1
  12:  1  6 12 10  3
  13:  1  2  1
  14:  1  4  5  2
  15:  1  4  5  2
  16:  1  5 10 10  5  1
For example, row n = 12 counts the following chains:
  ()  (1)   (2/1)   (4/2/1)   (12/4/2/1)
      (2)   (3/1)   (6/2/1)   (12/6/2/1)
      (3)   (4/1)   (6/3/1)   (12/6/3/1)
      (4)   (4/2)   (12/2/1)
      (6)   (6/1)   (12/3/1)
      (12)  (6/2)   (12/4/1)
            (6/3)   (12/4/2)
            (12/1)  (12/6/1)
            (12/2)  (12/6/2)
            (12/3)  (12/6/3)
            (12/4)
            (12/6)
		

Crossrefs

Column k = 1 is A000005.
Row ends are A008480.
Row lengths are A073093.
Column k = 2 is A238952.
The case from n to 1 is A334996 or A251683 (row sums: A074206).
A non-strict version is A334997 (transpose: A077592).
The case starting with n is A337255 (row sums: A067824).
Row sums are A337256 (nonempty: A253249).
A001055 counts factorizations.
A001221 counts distinct prime factors.
A001222 counts prime factors with multiplicity.
A097805 counts compositions by sum and length.
A122651 counts strict chains of divisors summing to n.
A146291 counts divisors of n with k prime factors (with multiplicity).
A163767 counts length n - 1 chains of divisors of n.
A167865 counts strict chains of divisors > 1 summing to n.
A337070 counts strict chains of divisors starting with superprimorials.

Programs

  • Mathematica
    Table[Length[Select[Reverse/@Subsets[Divisors[n],{k}],And@@Divisible@@@Partition[#,2,1]&]],{n,15},{k,0,PrimeOmega[n]+1}]

A343657 Sum of number of divisors of x^y for each x >= 1, y >= 0, x + y = n.

Original entry on oeis.org

1, 2, 4, 7, 12, 18, 27, 39, 56, 77, 103, 134, 174, 223, 283, 356, 445, 547, 666, 802, 959, 1139, 1344, 1574, 1835, 2128, 2454, 2815, 3213, 3648, 4126, 4653, 5239, 5888, 6608, 7407, 8298, 9288, 10385, 11597, 12936, 14408, 16025, 17799, 19746, 21882, 24221
Offset: 1

Views

Author

Gus Wiseman, Apr 29 2021

Keywords

Examples

			The a(7) = 27 divisors:
  1  32  81  64  25  6  1
     16  27  32  5   3
     8   9   16  1   2
     4   3   8       1
     2   1   4
     1       2
             1
		

Crossrefs

Antidiagonal row sums (row sums of the triangle) of A343656.
Dominated by A343661.
A000005(n) counts divisors of n.
A000312(n) = n^n.
A007318(n,k) counts k-sets of elements of {1..n}.
A009998(n,k) = n^k (as an array, offset 1).
A059481(n,k) counts k-multisets of elements of {1..n}.
A343658(n,k) counts k-multisets of divisors of n.

Programs

  • Mathematica
    Total/@Table[DivisorSigma[0,k^(n-k)],{n,30},{k,n}]

Formula

a(n) = Sum_{k=1..n} A000005(k^(n-k)).

A200221 Ordered factorizations of n with 3 parts.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 3, 0, 0, 0, 3, 0, 3, 0, 3, 0, 0, 0, 9, 0, 0, 1, 3, 0, 6, 0, 6, 0, 0, 0, 12, 0, 0, 0, 9, 0, 6, 0, 3, 3, 0, 0, 18, 0, 3, 0, 3, 0, 9, 0, 9, 0, 0, 0, 21, 0, 0, 3, 10, 0, 6, 0, 3, 0, 6, 0, 27, 0, 0, 3, 3, 0, 6, 0, 18, 3, 0, 0, 21
Offset: 1

Views

Author

Peter Luschny, Nov 14 2011

Keywords

Examples

			a(24) = 9 = card({{4,3,2}, {4,2,3}, {3,4,2}, {3,2,4}, {2,4,3}, {2,3,4}, {6,2,2},{2,6,2}, {2,2,6}}).
		

References

  • Arnold Knopfmacher and Michael Mays, Ordered and Unordered Factorizations of Integers, The Mathematica Journal, Vol 10 (1).

Crossrefs

Cf. A200214.
Column k=3 of A251683.

Programs

  • Maple
    with(numtheory):
    b:= proc(n) option remember; expand((`if`(isprime(n), 0,
          add(b(n/d), d=divisors(n) minus {1, n}))+1)*x)
        end:
    a:= n-> coeff(b(n), x, 3):
    seq(a(n), n=1..100);  # Alois P. Heinz, Dec 07 2014
  • Mathematica
    OrderedFactorizations[1] = {{}}; OrderedFactorizations[n_?PrimeQ] := {{n}}; OrderedFactorizations[n_] := OrderedFactorizations[n] = Flatten[Function[d, Prepend[#, d] & /@ OrderedFactorizations[n/d]] /@ Rest[Divisors[n]], 1]; a[n_] := With[{of3 = Sort /@ Select[OrderedFactorizations[n], Length[#] == 3 &] // Union}, Length[Permutations /@ of3 // Flatten[#, 1] &]]; Table[a[n], {n, 1, 84}] (* Jean-François Alcover, Jul 02 2013, copied and adapted from The Mathematica Journal *)
    nn = 200; f[list_, i_] := list[[i]]; a = Prepend[Table[1, {nn}], 0];
    c = Table[DirichletConvolve[f[a, n], f[a, n], n, m], {m, 1, nn}];
    Table[DirichletConvolve[f[a, n], f[c, n], n, m], {m, 1, nn}] (* Geoffrey Critzer, Apr 06 2020 *)

Formula

Dirichlet g.f.: (zeta(s)-1)^3. - Geoffrey Critzer, Apr 06 2020
Sum_{k=1..n} a(k) ~ n*(log(n)^2/2 + (3*gamma - 4)*log(n) + 3*gamma^2 - 9*gamma - 3*sg1 + 7), where gamma is the Euler-Mascheroni constant A001620 and sg1 is the first Stieltjes constant (see A082633). - Vaclav Kotesovec, Apr 07 2020

A254578 Number of ordered factorizations into distinct factors.

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 1, 3, 1, 3, 1, 5, 1, 3, 3, 3, 1, 5, 1, 5, 3, 3, 1, 13, 1, 3, 3, 5, 1, 13, 1, 5, 3, 3, 3, 13, 1, 3, 3, 13, 1, 13, 1, 5, 5, 3, 1, 21, 1, 5, 3, 5, 1, 13, 3, 13, 3, 3, 1, 29, 1, 3, 5, 11, 3, 13, 1, 5, 3, 13, 1, 29, 1, 3, 5, 5, 3, 13, 1, 21, 3, 3
Offset: 1

Views

Author

Geoffrey Critzer, Feb 01 2015

Keywords

Examples

			a(20)=5 because there are 5 ordered factorizations of 20 into distinct factors: 2*10, 4*5, 5*4, 10*2, 20.
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    b:= proc(n, i, p) option remember; `if`(n<=i, (p+1)!, 0)+add(
          b(n/d, d-1, p+1), d=select(x->x<=i, divisors(n)minus{1, n}))
        end:
    a:= n-> b(n$2, 0):
    seq(a(n), n=1..100);  # Alois P. Heinz, Feb 01 2015
  • Mathematica
    f[n_] := f[n] = Level[Table[Map[Prepend[#, d] &, f[n/d]], {d,Rest[Divisors[n]]}], {2}];
    f[1] = {{}};
    Map[Length,Map[Select[#, Apply[Unequal, #] &] &, Table[f[n], {n, 1, 60}]]]

A343935 Number of ways to choose a multiset of n divisors of n.

Original entry on oeis.org

1, 3, 4, 15, 6, 84, 8, 165, 55, 286, 12, 6188, 14, 680, 816, 4845, 18, 33649, 20, 53130, 2024, 2300, 24, 2629575, 351, 3654, 4060, 237336, 30, 10295472, 32, 435897, 7140, 7770, 8436, 177232627, 38, 10660, 11480, 62891499, 42, 85900584, 44, 1906884, 2118760
Offset: 1

Views

Author

Gus Wiseman, May 05 2021

Keywords

Examples

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

Crossrefs

Diagonal n = k of A343658.
Choosing n divisors of n - 1 gives A343936.
The version for chains of divisors is A343939.
A000005 counts divisors.
A000312 = n^n.
A007318 counts k-sets of elements of {1..n}.
A009998 = n^k (as an array, offset 1).
A059481 counts k-multisets of elements of {1..n}.
A146291 counts divisors of n with k prime factors (with multiplicity).
A253249 counts nonempty chains of divisors of n.
Strict chains of divisors:
- A067824 counts strict chains of divisors starting with n.
- A074206 counts strict chains of divisors from n to 1.
- A251683 counts strict length k + 1 chains of divisors from n to 1.
- A334996 counts strict length-k chains of divisors from n to 1.
- A337255 counts strict length-k chains of divisors starting with n.
- A337256 counts strict chains of divisors of n.
- A343662 counts strict length-k chains of divisors.

Programs

  • Mathematica
    multchoo[n_,k_]:=Binomial[n+k-1,k];
    Table[multchoo[DivisorSigma[0,n],n],{n,25}]
  • Python
    from math import comb
    from sympy import divisor_count
    def A343935(n): return comb(divisor_count(n)+n-1,n) # Chai Wah Wu, Jul 05 2024

Formula

a(n) = ((sigma(n), n)) = binomial(sigma(n) + n - 1, n) where sigma = A000005 and binomial = A007318.

A124433 Irregular array {a(n,m)} read by rows where (sum{n>=1} sum{m=1 to A001222(n)+1} a(n,m)*y^m/n^x) = 1/(zeta(x)-1+1/y) for all x and y where the double sum converges.

Original entry on oeis.org

1, 0, -1, 0, -1, 0, -1, 1, 0, -1, 0, -1, 2, 0, -1, 0, -1, 2, -1, 0, -1, 1, 0, -1, 2, 0, -1, 0, -1, 4, -3, 0, -1, 0, -1, 2, 0, -1, 2, 0, -1, 3, -3, 1, 0, -1, 0, -1, 4, -3, 0, -1, 0, -1, 4, -3, 0, -1, 2, 0, -1, 2, 0, -1, 0, -1, 6, -9, 4, 0, -1, 1, 0, -1, 2, 0, -1, 2, -1, 0, -1, 4, -3, 0, -1, 0, -1, 6, -6, 0, -1, 0, -1, 4, -6, 4, -1, 0, -1, 2, 0, -1, 2, 0, -1
Offset: 1

Views

Author

Leroy Quet, Dec 15 2006

Keywords

Comments

Row n has A001222(n)+1 terms. The polynomial P_n(y) = (sum{m=1 to A001222(n)+1} a(n,m)*y^m) is a generalization of the Mobius (Moebius) function, where P_n(1) = A008683(n).
From Gus Wiseman, Aug 24 2020: (Start)
Up to sign, also the number of strict length-k chains of divisors from n to 1, 1 <= k <= 1 + A001222(n). For example, row n = 36 counts the following chains (empty column indicated by dot):
. 36/1 36/2/1 36/4/2/1 36/12/4/2/1
36/3/1 36/6/2/1 36/12/6/2/1
36/4/1 36/6/3/1 36/12/6/3/1
36/6/1 36/9/3/1 36/18/6/2/1
36/9/1 36/12/2/1 36/18/6/3/1
36/12/1 36/12/3/1 36/18/9/3/1
36/18/1 36/12/4/1
36/12/6/1
36/18/2/1
36/18/3/1
36/18/6/1
36/18/9/1
(End)

Examples

			1/(zeta(x) - 1 + 1/y) = y - y^2/2^x - y^2/3^x + ( - y^2 + y^3)/4^x - y^2/5^x + ( - y^2 + 2y^3)/6^x - y^2/7^x + ...
From _Gus Wiseman_, Aug 24 2020: (Start)
The sequence of rows begins:
     1: 1              16: 0 -1 3 -3 1     31: 0 -1
     2: 0 -1           17: 0 -1            32: 0 -1 4 -6 4 -1
     3: 0 -1           18: 0 -1 4 -3       33: 0 -1 2
     4: 0 -1 1         19: 0 -1            34: 0 -1 2
     5: 0 -1           20: 0 -1 4 -3       35: 0 -1 2
     6: 0 -1 2         21: 0 -1 2          36: 0 -1 7 -12 6
     7: 0 -1           22: 0 -1 2          37: 0 -1
     8: 0 -1 2 -1      23: 0 -1            38: 0 -1 2
     9: 0 -1 1         24: 0 -1 6 -9 4     39: 0 -1 2
    10: 0 -1 2         25: 0 -1 1          40: 0 -1 6 -9 4
    11: 0 -1           26: 0 -1 2          41: 0 -1
    12: 0 -1 4 -3      27: 0 -1 2 -1       42: 0 -1 6 -6
    13: 0 -1           28: 0 -1 4 -3       43: 0 -1
    14: 0 -1 2         29: 0 -1            44: 0 -1 4 -3
    15: 0 -1 2         30: 0 -1 6 -6       45: 0 -1 4 -3
(End)
		

Crossrefs

A008480 gives rows ends (up to sign).
A008683 gives row sums (the Moebius function).
A073093 gives row lengths.
A074206 gives unsigned row sums.
A097805 is the restriction to powers of 2 (up to sign).
A251683 is the unsigned version with zeros removed.
A334996 is the unsigned version (except with a(1) = 0).
A334997 is an unsigned non-strict version.
A337107 is the restriction to factorial numbers.
A001055 counts factorizations.
A001222 counts prime factors with multiplicity.
A067824 counts strict chains of divisors starting with n.
A074206 counts strict chains of divisors from n to 1.
A122651 counts strict chains of divisors summing to n.
A167865 counts strict chains of divisors > 1 summing to n.
A253249 counts strict chains of divisors.
A337105 counts strict chains of divisors from n! to 1.

Programs

  • Mathematica
    f[l_List] := Block[{n = Length[l] + 1, c},c = Plus @@ Last /@ FactorInteger[n];Append[l, Prepend[ -Plus @@ Pick[PadRight[ #, c] & /@ l, Mod[n, Range[n - 1]], 0],0]]];Nest[f, {{1}}, 34] // Flatten(* Ray Chandler, Feb 13 2007 *)
    chnsc[n_]:=If[n==1,{{}},Prepend[Join@@Table[Prepend[#,n]&/@chnsc[d],{d,DeleteCases[Divisors[n],1|n]}],{n}]];
    Table[(-1)^k*Length[Select[chnsc[n],Length[#]==k&]],{n,30},{k,0,PrimeOmega[n]}] (* Gus Wiseman, Aug 24 2020 *)

Formula

a(1,1)=1. a(n,1) = 0 for n>=2. a(n,m+1) = -sum{k|n,k < n} a(k,m), where, for the purpose of this sum, a(k,m) = 0 if m > A001222(k)+1.

Extensions

Extended by Ray Chandler, Feb 13 2007

A254577 Total number of factors over all ordered factorizations of n.

Original entry on oeis.org

1, 1, 1, 3, 1, 5, 1, 8, 3, 5, 1, 18, 1, 5, 5, 20, 1, 18, 1, 18, 5, 5, 1, 56, 3, 5, 8, 18, 1, 31, 1, 48, 5, 5, 5, 75, 1, 5, 5, 56, 1, 31, 1, 18, 18, 5, 1, 160, 3, 18, 5, 18, 1, 56, 5, 56, 5, 5, 1, 132, 1, 5, 18, 112, 5, 31, 1, 18, 5, 31, 1, 264, 1, 5, 18, 18, 5
Offset: 1

Views

Author

Geoffrey Critzer, Feb 01 2015

Keywords

Comments

What is the limit log(Sum_{k=1..n} a(k)) / log(n) ?. - Vaclav Kotesovec, Feb 03 2019

Examples

			a(20)=18 because in the ordered factorizations of twenty: 20, 2*10, 4*5, 5*4, 10*2, 2*2*5, 2*5*2, 5*2*2 there are a total of 18 factors.
		

Crossrefs

Cf. A074206.

Programs

  • Maple
    with(numtheory):
    b:= proc(n) option remember; expand(x*(1+
          add(b(n/d), d=divisors(n) minus {1, n})))
        end:
    a:= n-> (p-> add(coeff(p, x, i)*i, i=1..degree(p)))(b(n)):
    seq(a(n), n=1..100);  # Alois P. Heinz, Feb 01 2015
  • Mathematica
    f[n_] := f[n] =Level[Table[Map[Prepend[#, d] &, f[n/d]], {d,Rest[Divisors[n]]}], {2}];
    f[1] = {{}};
    g[list_] := Sum[list[[i]] i, {i, 1, Length[list]}];
    Prepend[Rest[Map[g,Map[Table[Count[#, i], {i, 1, Max[#]}] &,Map[Length, Map[Sort, Table[f[n], {n, 1, 60}]], {2}]]]], 1]

Formula

Dirichlet generating function: zeta(s)/(1 - zeta(s))^2.
a(n) = Sum_{k>=1} A251683(n,k)*k.

A343939 Number of n-chains of divisors of n.

Original entry on oeis.org

1, 3, 4, 15, 6, 49, 8, 165, 55, 121, 12, 1183, 14, 225, 256, 4845, 18, 3610, 20, 4851, 484, 529, 24, 73125, 351, 729, 4060, 12615, 30, 29791, 32, 435897, 1156, 1225, 1296, 494209, 38, 1521, 1600, 505981, 42, 79507, 44, 46575, 49726, 2209, 48
Offset: 1

Views

Author

Gus Wiseman, May 05 2021

Keywords

Examples

			The a(1) = 1 through a(5) = 6 chains:
  (1)  (1/1)  (1/1/1)  (1/1/1/1)  (1/1/1/1/1)
       (2/1)  (3/1/1)  (2/1/1/1)  (5/1/1/1/1)
       (2/2)  (3/3/1)  (2/2/1/1)  (5/5/1/1/1)
              (3/3/3)  (2/2/2/1)  (5/5/5/1/1)
                       (2/2/2/2)  (5/5/5/5/1)
                       (4/1/1/1)  (5/5/5/5/5)
                       (4/2/1/1)
                       (4/2/2/1)
                       (4/2/2/2)
                       (4/4/1/1)
                       (4/4/2/1)
                       (4/4/2/2)
                       (4/4/4/1)
                       (4/4/4/2)
                       (4/4/4/4)
		

Crossrefs

Diagonal n = k - 1 of the array A077592.
Chains of length n - 1 are counted by A163767.
Diagonal n = k of the array A334997.
The version counting all multisets of divisors (not just chains) is A343935.
A000005(n) counts divisors of n.
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-1) counts strict k-chains of divisors from n to 1.
A253249(n) counts nonempty chains of divisors of n.
A334996(n,k) counts strict k-chains of divisors from n to 1.
A337255(n,k) counts strict k-chains of divisors starting with n.
A343658(n,k) counts k-multisets of divisors of n.
A343662(n,k) counts strict k-chains of divisors of n (row sums: A337256).

Programs

  • Mathematica
    Table[Length[Select[Tuples[Divisors[n],n],OrderedQ[#]&&And@@Divisible@@@Reverse/@Partition[#,2,1]&]],{n,10}]
Previous Showing 11-20 of 28 results. Next