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

A067824 a(1) = 1; for n > 1, a(n) = 1 + Sum_{0 < d < n, d|n} a(d).

Original entry on oeis.org

1, 2, 2, 4, 2, 6, 2, 8, 4, 6, 2, 16, 2, 6, 6, 16, 2, 16, 2, 16, 6, 6, 2, 40, 4, 6, 8, 16, 2, 26, 2, 32, 6, 6, 6, 52, 2, 6, 6, 40, 2, 26, 2, 16, 16, 6, 2, 96, 4, 16, 6, 16, 2, 40, 6, 40, 6, 6, 2, 88, 2, 6, 16, 64, 6, 26, 2, 16, 6, 26, 2, 152, 2, 6, 16, 16, 6, 26, 2, 96, 16, 6, 2, 88, 6, 6, 6, 40, 2, 88, 6, 16, 6, 6, 6, 224, 2, 16, 16, 52
Offset: 1

Views

Author

Reinhard Zumkeller, Feb 08 2002

Keywords

Comments

By a result of Karhumaki and Lifshits, this is also the number of polynomials p(x) with coefficients in {0,1} that divide x^n-1 and such that (x^n-1)/ {(x-1)p(x)} has all coefficients in {0,1}.
The number of tiles of a discrete interval of length n (an interval of Z). - Eric H. Rivals (rivals(AT)lirmm.fr), Mar 13 2007
Bodini and Rivals proved this is the number of tiles of a discrete interval of length n and also is the number (A107067) of polynomials p(x) with coefficients in {0,1} that divide x^n-1 and such that (x^n-1)/ {(x-1)p(x)} has all coefficients in {0,1} (Bodini, Rivals, 2006). This structure of such tiles is also known as Krasner's factorization (Krasner and Ranulac, 1937). The proof also gives an algorithm to recognize if a set is a tile in optimal time and in this case, to compute the smallest interval it can tile (Bodini, Rivals, 2006). - Eric H. Rivals (rivals(AT)lirmm.fr), Mar 13 2007
Number of lone-child-avoiding rooted achiral (or generalized Bethe) trees with positive integer leaves summing to n, where a rooted tree is lone-child-avoiding if all terminal subtrees have at least two branches, and achiral if all branches directly under any given node are equal. For example, the a(6) = 6 trees are 6, (111111), (222), ((11)(11)(11)), (33), ((111)(111)). - Gus Wiseman, Jul 13 2018. Updated Aug 22 2020.
From Gus Wiseman, Aug 20 2020: (Start)
Also the number of strict chains of divisors starting with n. For example, the a(n) chains for n = 1, 2, 4, 6, 8, 12 are:
1 2 4 6 8 12
2/1 4/1 6/1 8/1 12/1
4/2 6/2 8/2 12/2
4/2/1 6/3 8/4 12/3
6/2/1 8/2/1 12/4
6/3/1 8/4/1 12/6
8/4/2 12/2/1
8/4/2/1 12/3/1
12/4/1
12/4/2
12/6/1
12/6/2
12/6/3
12/4/2/1
12/6/2/1
12/6/3/1
(End)
a(n) is the number of chains including n of the divisor lattice of divisors of n, which is to say, a(n) is the number of (d_1,d_2,...,d_k) such that d_1 < d_2 < ... < d_k = n and d_i divides d_{i+1} for 1 <= i <= k-1. Using this definition, the recurrence a(n) = 1 + Sum_{0 < d < n, d|n} a(d) is evident by enumerating the preceding element of n in the chains. If we count instead the chains whose LCM of members is n, then a(1) would be 2 because the empty chain is included, and we would obtain 2*A074206(n). - Jianing Song, Aug 21 2024

Examples

			a(12) = 1 + a(6) + a(4) + a(3) + a(2) + a(1)
= 1+(1+a(3)+a(2)+a(1))+(1+a(2)+a(1))+(1+a(1))+(1+a(1))+(1)
= 1+(1+(1+a(1))+(1+a(1))+1)+(1+(1+a(1))+1)+(1+1)+(1+1)+(1)
= 1+(1+(1+1)+(1+1)+1)+(1+(1+1)+1)+(1+1)+(1+1)+(1)
= 1 + 6 + 4 + 2 + 2 + 1 = 16.
		

References

  • Olivier Bodini and Eric Rivals. Tiling an Interval of the Discrete Line. In M. Lewenstein and G. Valiente, editors, Proc. of the 17th Annual Symposium on Combinatorial Pattern Matching (CPM), volume 4009 of Lecture Notes in Computer Science, pages 117-128. Springer Verlag, 2006.
  • Juhani Karhumaki, Yury Lifshits and Wojciech Rytter, Tiling Periodicity, in Combinatorial Pattern Matching, Lecture Notes in Computer Science, Volume 4580/2007, Springer-Verlag.

Crossrefs

Cf. A122408 (fixed points).
Inverse Möbius transform of A074206.
A001055 counts factorizations.
A008480 counts maximal chains of divisors starting with n.
A074206 counts chains of divisors from n to 1.
A253249 counts nonempty chains of divisors.
A337070 counts chains of divisors starting with A006939(n).
A337071 counts chains of divisors starting with n!.
A337256 counts chains of divisors.
Cf. A001221, A001222, A002033, A124010, A337074, A337105, A378223, A378225 (Dirichlet inverse).

Programs

  • Haskell
    a067824 n = 1 + sum (map a067824 [d | d <- [1..n-1], mod n d == 0])
    -- Reinhard Zumkeller, Oct 13 2011
    
  • Maple
    a:= proc(n) option remember;
          1+add(a(d), d=numtheory[divisors](n) minus {n})
        end:
    seq(a(n), n=1..100);  # Alois P. Heinz, Apr 17 2021
  • Mathematica
    a[1]=1; a[n_] := a[n] = 1+Sum[If[Mod[n,d]==0, a[d], 0], {d, 1, n-1}]; Array[a,100] (* Jean-François Alcover, Apr 28 2011 *)
  • PARI
    A=vector(100);A[1]=1; for(n=2,#A,A[n]=1+sumdiv(n,d,A[d])); A \\ Charles R Greathouse IV, Nov 20 2012

Formula

a(n) = 2*A074206(n), n>1. - Vladeta Jovovic, Jul 03 2005
a(p^k) = 2^k for primes p. - Reinhard Zumkeller, Sep 03 2006
a(n) = Sum_{d|n} A002033(d-1) = Sum_{d|n} A074206(d). - Gus Wiseman, Jul 13 2018
Dirichlet g.f.: zeta(s) / (2 - zeta(s)). - Álvar Ibeas, Dec 30 2018
G.f. A(x) satisfies: A(x) = x/(1 - x) + Sum_{k>=2} A(x^k). - Ilya Gutkovskiy, May 18 2019

Extensions

Entry revised by N. J. A. Sloane, Aug 27 2006

A251683 Irregular triangular array: T(n,k) is the number of ordered factorizations of n with exactly k factors, n >= 1, 1 <= k <= A086436(n).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 4, 3, 1, 1, 2, 1, 2, 1, 3, 3, 1, 1, 1, 4, 3, 1, 1, 4, 3, 1, 2, 1, 2, 1, 1, 6, 9, 4, 1, 1, 1, 2, 1, 2, 1, 1, 4, 3, 1, 1, 6, 6, 1, 1, 4, 6, 4, 1, 1, 2, 1, 2, 1, 2, 1, 7, 12, 6, 1, 1, 2, 1, 2, 1, 6, 9, 4
Offset: 1

Views

Author

Geoffrey Critzer, Dec 06 2014

Keywords

Comments

Row sums = A074206.
Row lengths give A086436.
T(n,2) = A070824(n).
T(n,3) = A200221(n).
Sum_{k>=1} k*T(n,k) = A254577.
For all n > 1, Sum_{k=1..A086436(n)} (-1)^k*T(n,k) = A008683(n). - Geoffrey Critzer, May 25 2018
From Gus Wiseman, Aug 21 2020: (Start)
Also the number of strict length k + 1 chains of divisors from n to 1. For example, row n = 24 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
(End)

Examples

			Triangle T(n,k) begins:
  1;
  1;
  1;
  1, 1;
  1;
  1, 2;
  1;
  1, 2, 1;
  1, 1;
  1, 2;
  1;
  1, 4, 3;
  1;
  1, 2;
  1, 2;
  ...
There are 8 ordered factorizations of the integer 12: 12, 6*2, 4*3, 3*4, 2*6, 3*2*2, 2*3*2, 2*2*3.  So T(12,1)=1, T(12,2)=4, and T(12,3)=3.
		

Crossrefs

A008480 gives rows ends.
A086436 gives row lengths.
A124433 is the same except for signs and zeros.
A334996 is the same except for zeros.
A337107 is the restriction to factorial numbers (but with zeros).
A000005 counts divisors.
A001055 counts factorizations.
A001222 counts prime factors with multiplicity.
A074206 counts strict chains of divisors from n to 1.
A067824 counts strict chains of divisors starting with n.
A122651 counts strict chains of divisors summing to n.
A167865 counts strict chains of divisors > 1 summing to n.
A253249 counts strict nonempty chains of divisors of n.
A337071 counts strict chains of divisors starting with n!.
A337256 counts strict chains of divisors of n.

Programs

  • Maple
    with(numtheory):
    b:= proc(n) option remember; expand(x*(1+
          add(b(n/d), d=divisors(n) minus {1, n})))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=1..degree(p)))(b(n)):
    seq(T(n), n=1..100);  # Alois P. Heinz, Dec 07 2014
  • Mathematica
    f[1] = {{}};
    f[n_] := f[n] =
      Level[Table[
        Map[Prepend[#, d] &, f[n/d]], {d, Rest[Divisors[n]]}], {2}];
    Prepend[Map[Select[#, # > 0 &] &,
      Drop[Transpose[
        Table[Map[Count[#, k] &,
          Map[Length, Table[f[n], {n, 1, 40}], {2}]], {k, 1, 10}]],
       1]],{1}] // Grid
    (* Second program: *)
    b[n_] := b[n] = x(1+Sum[b[n/d], {d, Divisors[n]~Complement~{1, n}}]);
    T[n_] := CoefficientList[b[n]/x, x];
    Array[T, 100] // Flatten (* Jean-François Alcover, Nov 17 2020, after Alois P. Heinz *)

Formula

Dirichlet g.f.: 1/(1 - y*(zeta(x)-1)).

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

A336424 Number of factorizations of n where each factor belongs to A130091 (numbers with distinct prime multiplicities).

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 3, 1, 1, 1, 5, 1, 3, 1, 3, 1, 1, 1, 5, 2, 1, 3, 3, 1, 1, 1, 7, 1, 1, 1, 6, 1, 1, 1, 5, 1, 1, 1, 3, 3, 1, 1, 9, 2, 3, 1, 3, 1, 5, 1, 5, 1, 1, 1, 4, 1, 1, 3, 11, 1, 1, 1, 3, 1, 1, 1, 11, 1, 1, 3, 3, 1, 1, 1, 9, 5, 1, 1, 4, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Aug 03 2020

Keywords

Comments

A number's prime signature (row n of A124010) is the sequence of positive exponents in its prime factorization, so a number has distinct prime multiplicities iff all the exponents in its prime signature are distinct.

Examples

			The a(n) factorizations for n = 2, 4, 8, 60, 16, 36, 32, 48:
  2  4    8      5*12     16       4*9      32         48
     2*2  2*4    3*20     4*4      3*12     4*8        4*12
          2*2*2  3*4*5    2*8      3*3*4    2*16       3*16
                 2*2*3*5  2*2*4    2*18     2*4*4      3*4*4
                          2*2*2*2  2*2*9    2*2*8      2*24
                                   2*2*3*3  2*2*2*4    2*3*8
                                            2*2*2*2*2  2*2*12
                                                       2*2*3*4
                                                       2*2*2*2*3
		

Crossrefs

A327523 is the case when n is restricted to belong to A130091 also.
A001055 counts factorizations.
A007425 counts divisors of divisors.
A045778 counts strict factorizations.
A074206 counts ordered factorizations.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
A253249 counts nonempty chains of divisors.
A281116 counts factorizations with no common divisor.
A302696 lists numbers whose prime indices are pairwise coprime.
A305149 counts stable factorizations.
A320439 counts factorizations using A289509.
A327498 gives the maximum divisor with distinct prime multiplicities.
A336500 counts divisors of n in A130091 with quotient also in A130091.
A336568 = not a product of two numbers with distinct prime multiplicities.
A336569 counts maximal chains of elements of A130091.
A337256 counts chains of divisors.

Programs

  • Mathematica
    facsusing[s_,n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facsusing[Select[s,Divisible[n/d,#]&],n/d],Min@@#>=d&]],{d,Select[s,Divisible[n,#]&]}]];
    Table[Length[facsusing[Select[Range[2,n],UnsameQ@@Last/@FactorInteger[#]&],n]],{n,100}]

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

A336423 Number of strict chains of divisors from n to 1 using terms of A130091 (numbers with distinct prime multiplicities).

Original entry on oeis.org

1, 1, 1, 2, 1, 0, 1, 4, 2, 0, 1, 5, 1, 0, 0, 8, 1, 5, 1, 5, 0, 0, 1, 14, 2, 0, 4, 5, 1, 0, 1, 16, 0, 0, 0, 0, 1, 0, 0, 14, 1, 0, 1, 5, 5, 0, 1, 36, 2, 5, 0, 5, 1, 14, 0, 14, 0, 0, 1, 0, 1, 0, 5, 32, 0, 0, 1, 5, 0, 0, 1, 35, 1, 0, 5, 5, 0, 0, 1, 36, 8, 0, 1, 0
Offset: 1

Views

Author

Gus Wiseman, Jul 27 2020

Keywords

Comments

A number's prime signature (row n of A124010) is the sequence of positive exponents in its prime factorization, so a number has distinct prime multiplicities iff all the exponents in its prime signature are distinct.

Examples

			The a(n) chains for n = 4, 8, 12, 16, 24, 32:
  4/1    8/1      12/1      16/1        24/1         32/1
  4/2/1  8/2/1    12/2/1    16/2/1      24/2/1       32/2/1
         8/4/1    12/3/1    16/4/1      24/3/1       32/4/1
         8/4/2/1  12/4/1    16/8/1      24/4/1       32/8/1
                  12/4/2/1  16/4/2/1    24/8/1       32/16/1
                            16/8/2/1    24/12/1      32/4/2/1
                            16/8/4/1    24/4/2/1     32/8/2/1
                            16/8/4/2/1  24/8/2/1     32/8/4/1
                                        24/8/4/1     32/16/2/1
                                        24/12/2/1    32/16/4/1
                                        24/12/3/1    32/16/8/1
                                        24/12/4/1    32/8/4/2/1
                                        24/8/4/2/1   32/16/4/2/1
                                        24/12/4/2/1  32/16/8/2/1
                                                     32/16/8/4/1
                                                     32/16/8/4/2/1
		

Crossrefs

A336569 is the maximal case.
A336571 does not require n itself to have distinct prime multiplicities.
A000005 counts divisors.
A007425 counts divisors of divisors.
A074206 counts strict chains of divisors from n to 1.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
A253249 counts nonempty strict chains of divisors.
A327498 gives the maximum divisor with distinct prime multiplicities.
A336422 counts divisible pairs of divisors, both in A130091.
A336424 counts factorizations using A130091.
A336500 counts divisors of n in A130091 with quotient also in A130091.
A337256 counts strict chains of divisors.

Programs

  • Mathematica
    strchns[n_]:=If[n==1,1,If[!UnsameQ@@Last/@FactorInteger[n],0,Sum[strchns[d],{d,Select[Most[Divisors[n]],UnsameQ@@Last/@FactorInteger[#]&]}]]];
    Table[strchns[n],{n,100}]

A163767 a(n) = tau_{n}(n) = number of ordered n-factorizations of n.

Original entry on oeis.org

1, 2, 3, 10, 5, 36, 7, 120, 45, 100, 11, 936, 13, 196, 225, 3876, 17, 3078, 19, 4200, 441, 484, 23, 62400, 325, 676, 3654, 11368, 29, 27000, 31, 376992, 1089, 1156, 1225, 443556, 37, 1444, 1521, 459200, 41, 74088, 43, 43560, 46575, 2116, 47, 11995200, 1225
Offset: 1

Views

Author

Paul D. Hanna, Aug 04 2009

Keywords

Comments

Also the number of length n - 1 chains of divisors of n. - Gus Wiseman, May 07 2021

Examples

			Successive Dirichlet self-convolutions of the all 1's sequence begin:
(1),1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,... (A000012)
1,(2),2,3,2,4,2,4,3,4,2,6,2,4,4,5,... (A000005)
1,3,(3),6,3,9,3,10,6,9,3,18,3,9,9,15,... (A007425)
1,4,4,(10),4,16,4,20,10,16,4,40,4,16,16,35,... (A007426)
1,5,5,15,(5),25,5,35,15,25,5,75,5,25,25,70,... (A061200)
1,6,6,21,6,(36),6,56,21,36,6,126,6,36,36,126,... (A034695)
1,7,7,28,7,49,(7),84,28,49,7,196,7,49,49,210,... (A111217)
1,8,8,36,8,64,8,(120),36,64,8,288,8,64,64,330,... (A111218)
1,9,9,45,9,81,9,165,(45),81,9,405,9,81,81,495,... (A111219)
1,10,10,55,10,100,10,220,55,(100),10,550,10,100,... (A111220)
1,11,11,66,11,121,11,286,66,121,(11),726,11,121,... (A111221)
1,12,12,78,12,144,12,364,78,144,12,(936),12,144,... (A111306)
...
where the main diagonal forms this sequence.
From _Gus Wiseman_, May 07 2021: (Start)
The a(1) = 1 through a(5) = 5 chains of divisors:
  ()  (1)  (1/1)  (1/1/1)  (1/1/1/1)
      (2)  (3/1)  (2/1/1)  (5/1/1/1)
           (3/3)  (2/2/1)  (5/5/1/1)
                  (2/2/2)  (5/5/5/1)
                  (4/1/1)  (5/5/5/5)
                  (4/2/1)
                  (4/2/2)
                  (4/4/1)
                  (4/4/2)
                  (4/4/4)
(End)
		

Crossrefs

Main diagonal of A077592.
Diagonal n = k + 1 of the array A334997.
The version counting all multisets of divisors (not just chains) is A343935.
A000005 counts divisors.
A001055 counts factorizations (strict: A045778, ordered: A074206).
A001221 counts distinct prime factors.
A001222 counts prime factors with multiplicity.
A067824 counts strict chains of divisors starting with n.
A122651 counts strict chains of divisors summing to n.
A146291 counts divisors of n with k prime factors (with multiplicity).
A167865 counts strict chains of divisors > 1 summing to n.
A253249 counts nonempty strict chains of divisors of n.
A251683/A334996 count strict nonempty length-k divisor chains from n to 1.
A337255 counts strict length-k chains of divisors starting with n.
A339564 counts factorizations with a selected factor.
A343662 counts strict length-k chains of divisors (row sums: A337256).
Cf. A060690.

Programs

  • Mathematica
    Table[Times@@(Binomial[#+n-1,n-1]&/@FactorInteger[n][[All,2]]),{n,1,50}] (* Enrique Pérez Herrero, Dec 25 2013 *)
  • PARI
    {a(n,m=n)=if(n==1,1,if(m==1,1,sumdiv(n,d,a(d,1)*a(n/d,m-1))))}
    
  • Python
    from math import prod, comb
    from sympy import factorint
    def A163767(n): return prod(comb(n+e-1,e) for e in factorint(n).values()) # Chai Wah Wu, Jul 05 2024

Formula

a(p) = p for prime p.
a(n) = n^k when n is the product of k distinct primes (conjecture).
a(n) = n-th term of the n-th Dirichlet self-convolution of the all 1's sequence.
a(2^n) = A060690(n). - Alois P. Heinz, Jun 12 2024

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

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

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.
Showing 1-10 of 13 results. Next