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

A074206 Kalmár's [Kalmar's] problem: number of ordered factorizations of n.

Original entry on oeis.org

0, 1, 1, 1, 2, 1, 3, 1, 4, 2, 3, 1, 8, 1, 3, 3, 8, 1, 8, 1, 8, 3, 3, 1, 20, 2, 3, 4, 8, 1, 13, 1, 16, 3, 3, 3, 26, 1, 3, 3, 20, 1, 13, 1, 8, 8, 3, 1, 48, 2, 8, 3, 8, 1, 20, 3, 20, 3, 3, 1, 44, 1, 3, 8, 32, 3, 13, 1, 8, 3, 13, 1, 76, 1, 3, 8, 8, 3, 13, 1, 48, 8, 3, 1, 44, 3, 3, 3, 20, 1, 44, 3, 8, 3, 3, 3, 112
Offset: 0

Views

Author

N. J. A. Sloane, Apr 29 2003

Keywords

Comments

a(0)=0, a(1)=1; thereafter a(n) is the number of ordered factorizations of n as a product of integers greater than 1.
Kalmár (1931) seems to be the earliest reference that mentions this sequence (as opposed to A002033). - N. J. A. Sloane, May 05 2016
a(n) is the permanent of the n-1 X n-1 matrix A with (i,j) entry = 1 if j|i+1 and = 0 otherwise. This is because ordered factorizations correspond to nonzero elementary products in the permanent. For example, with n=6, 3*2 -> 1,3,6 [partial products] -> 6,3,1 [reverse list] -> (6,3)(3,1) [partition into pairs with offset 1] -> (5,3)(2,1) [decrement first entry] -> (5,3)(2,1)(1,2)(3,4)(4,5) [append pairs (i,i+1) to get a permutation] -> elementary product A(1,2)A(2,1)A(3,4)A(4,5)A(5,3). - David Callan, Oct 19 2005
This sequence is important in describing the amount of energy in all wave structures in the Universe according to harmonics theory. - Ray Tomes (ray(AT)tomes.biz), Jul 22 2007
a(n) appears to be the number of permutation matrices contributing to the Moebius function. See A008683 for more information. Also a(n) appears to be the Moebius transform of A067824. Furthermore it appears that except for the first term a(n)=A067824(n)*(1/2). Are there other sequences such that when the Moebius transform is applied, the new sequence is also a constant factor times the starting sequence? - Mats Granvik, Jan 01 2009
Numbers divisible by n distinct primes appear to have ordered factorization values that can be found in an n-dimensional summatory Pascal triangle. For example, the ordered factorization values for numbers divisible by two distinct primes can be found in table A059576. - Mats Granvik, Sep 06 2009
Inverse Mobius transform of A174725 and also except for the first term, inverse Mobius transform of A174726. - Mats Granvik, Mar 28 2010
a(n) is a lower bound on the worst-case number of solutions to the probed partial digest problem for n fragments of DNA; see the Newberg & Naor reference, below. - Lee A. Newberg, Aug 02 2011
All integers greater than 1 are perfect numbers over this sequence (for definition of A-perfect numbers, see comment to A175522). - Vladimir Shevelev, Aug 03 2011
If n is squarefree, then a(n) = A000670(A001221(n)) = A000670(A001222(n)). - Vladimir Shevelev and Franklin T. Adams-Watters, Aug 05 2011
A034776 lists the values taken by this sequence. - Robert G. Wilson v, Jun 02 2012
From Gus Wiseman, Aug 25 2020: (Start)
Also the number of strict chains of divisors from n to 1. For example, the a(n) chains for n = 1, 2, 4, 6, 8, 12, 30 are:
1 2/1 4/1 6/1 8/1 12/1 30/1
4/2/1 6/2/1 8/2/1 12/2/1 30/2/1
6/3/1 8/4/1 12/3/1 30/3/1
8/4/2/1 12/4/1 30/5/1
12/6/1 30/6/1
12/4/2/1 30/10/1
12/6/2/1 30/15/1
12/6/3/1 30/6/2/1
30/6/3/1
30/10/2/1
30/10/5/1
30/15/3/1
30/15/5/1
(End)
a(n) is also the number of ways to tile a strip of length log(n) with tiles having lengths {log(k) : k>=2}. - David Bevan, Jan 07 2025

Examples

			G.f. = x + x^2 + x^3 + 2*x^4 + x^5 + 3*x^6 + x^7 + 4*x^8 + 2*x^9 + 3*x^10 + ...
Number of ordered factorizations of 8 is 4: 8 = 2*4 = 4*2 = 2*2*2.
		

References

  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 126, see #27.
  • R. Honsberger, Mathematical Gems III, M.A.A., 1985, p. 141.
  • Kalmár, Laszlo, A "factorisatio numerorum" problemajarol [Hungarian], Matemat. Fizik. Lapok, 38 (1931), 1-15.
  • J. Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 124.

Crossrefs

Apart from initial term, same as A002033.
a(A002110) = A000670, row sums of A251683.
A173382 (and A025523) gives partial sums.
A124433 has these as unsigned row sums.
A334996 has these as row sums.
A001055 counts factorizations.
A001222 counts prime factors with multiplicity.
A008480 counts ordered prime factorizations.
A067824 counts strict chains of divisors starting with n.
A122651 counts strict chains of divisors summing to n.
A253249 counts strict chains of divisors.

Programs

  • Haskell
    a074206 n | n <= 1 = n
    | otherwise = 1 + (sum $ map (a074206 . (div n)) $
    tail $ a027751_row n)
    -- Reinhard Zumkeller, Oct 01 2012
    
  • Maple
    a := array(1..150): for k from 1 to 150 do a[k] := 0 od: a[1] := 1: for j from 2 to 150 do for m from 1 to j-1 do if j mod m = 0 then a[j] := a[j]+a[m] fi: od: od: for k from 1 to 150 do printf(`%d,`,a[k]) od: # James Sellers, Dec 07 2000
    A074206:= proc(n) option remember; if n > 1 then `+`(op(apply(A074206, numtheory[divisors](n)[1..-2]))) else n fi end: # M. F. Hasler, Oct 12 2018
  • Mathematica
    a[0] = 0; a[1] = 1; a[n_] := a[n] = a /@ Most[Divisors[n]] // Total; a /@ Range[20000] (* N. J. A. Sloane, May 04 2016, based on program in A002033 *)
    ccc[n_]:=Switch[n,0,{},1,{{1}},,Join@@Table[Prepend[#,n]&/@ccc[d],{d,Most[Divisors[n]]}]]; Table[Length[ccc[n]],{n,0,100}] (* _Gus Wiseman, Aug 25 2020 *)
  • PARI
    A=vector(100);A[1]=1; for(n=2,#A,A[n]=1+sumdiv(n,d,A[d])); A/=2; A[1]=1; concat(0,A) \\ Charles R Greathouse IV, Nov 20 2012
    
  • PARI
    {a(n) = if( n<2, n>0, my(A = divisors(n)); sum(k=1, #A-1, a(A[k])))}; /* Michael Somos, Dec 26 2016 */
    
  • PARI
    A074206(n)=if(n>1, sumdiv(n, i, if(iA074206(i))),n) \\ M. F. Hasler, Oct 12 2018
    
  • PARI
    A74206=[1]; A074206(n)={if(#A74206A074206(i)))} \\ Use memoization for computing many values. - M. F. Hasler, Oct 12 2018
    
  • PARI
    first(n) = {my(res = vector(n, i, 1)); for(i = 2, n, for(j = 2, n \ i, res[i*j] += res[i])); concat(0, res)} \\ David A. Corneth, Oct 13 2018
    
  • PARI
    first(n) = {my(res = vector(n, i, 1)); for(i = 2, n, d = divisors(i); res[i] += sum(j = 1, #d-1, res[d[j]])); concat(0, res)} \\ somewhat faster than progs above for finding first terms of n. \\ David A. Corneth, Oct 12 2018
    
  • PARI
    a(n)={if(!n, 0, my(sig=factor(n)[,2], m=vecsum(sig)); sum(k=0, m, prod(i=1, #sig, binomial(sig[i]+k-1, k-1))*sum(r=k, m, binomial(r,k)*(-1)^(r-k))))} \\ Andrew Howroyd, Aug 30 2020
    
  • Python
    from math import prod
    from functools import lru_cache
    from sympy import divisors, factorint, prime
    @lru_cache(maxsize=None)
    def A074206(n): return sum(A074206(d) for d in divisors(prod(prime(i+1)**e for i,e in enumerate(sorted(factorint(n).values(),reverse=True))),generator=True,proper=True)) if n > 1 else n # Chai Wah Wu, Sep 16 2022
  • SageMath
    @cached_function
    def minus_mu(n):
        if n < 2: return n
        return sum(minus_mu(d) for d in divisors(n)[:-1])
    # Note that changing the sign of the sum gives the Möbius function A008683.
    print([minus_mu(n) for n in (0..96)]) # Peter Luschny, Dec 26 2016
    

Formula

With different offset: a(n) = sum of all a(i) such that i divides n and i < n. - Clark Kimberling
a(p^k) = 2^(k-1) if k>0 and p is a prime.
Dirichlet g.f.: 1/(2-zeta(s)). - Herbert S. Wilf, Apr 29 2003
a(n) = A067824(n)/2 for n>1; a(A122408(n)) = A122408(n)/2. - Reinhard Zumkeller, Sep 03 2006
If p,q,r,... are distinct primes, then a(p*q)=3, a(p^2*q)=8, a(p*q*r)=13, a(p^3*q)=20, etc. - Vladimir Shevelev, Aug 03 2011 [corrected by Charles R Greathouse IV, Jun 02 2012]
a(0) = 0, a(1) = 1; a(n) = [x^n] Sum_{k=1..n-1} a(k)*x^k/(1 - x^k). - Ilya Gutkovskiy, Dec 11 2017
a(n) = a(A046523(n)); a(A025487(n)) = A050324(n): a(n) depends only on the nonzero exponents in the prime factorization of n, more precisely prime signature of n, cf. A124010 and A320390. - M. F. Hasler, Oct 12 2018
a(n) = A000670(A001221(n)) for squarefree n. In particular a(A002110(n)) = A000670(n). - Amiram Eldar, May 13 2019
a(n) = A050369(n)/n, for n>=1. - Ridouane Oudra, Aug 31 2019
a(n) = A361665(A181819(n)). - Pontus von Brömssen, Mar 25 2023
From Ridouane Oudra, Nov 02 2023: (Start)
If p,q are distinct primes, and n,m>0 then we have:
a(p^n*q^m) = Sum_{k=0..min(n,m)} 2^(n+m-k-1)*binomial(n,k)*binomial(m,k);
More generally: let tau[k](n) denote the number of ordered factorizations of n as a product of k terms, also named the k-th Piltz function (see A007425), then we have for n>1:
a(n) = Sum_{j=1..bigomega(n)} Sum_{k=1..j} (-1)^(j-k)*binomial(j,k)*tau[k](n), or
a(n) = Sum_{j=1..bigomega(n)} Sum_{k=0..j-1} (-1)^k*binomial(j,k)*tau[j-k](n). (End)

Extensions

Originally this sequence was merged with A002033, the number of perfect partitions. Herbert S. Wilf suggested that it warrants an entry of its own.

A167865 Number of partitions of n into distinct parts greater than 1, with each part divisible by the next.

Original entry on oeis.org

1, 0, 1, 1, 1, 1, 2, 1, 2, 2, 2, 1, 4, 1, 3, 3, 3, 1, 5, 1, 5, 4, 3, 1, 6, 2, 5, 4, 5, 1, 9, 1, 6, 4, 4, 4, 8, 1, 6, 6, 7, 1, 11, 1, 8, 8, 4, 1, 10, 3, 10, 5, 8, 1, 11, 4, 10, 7, 6, 1, 13, 1, 10, 11, 7, 6, 15, 1, 9, 5, 11, 1, 14, 1, 9, 12, 8, 5, 15, 1, 16, 9, 8, 1, 18, 5, 12, 7, 10, 1, 21, 7, 13, 11, 5
Offset: 0

Views

Author

Max Alekseyev, Nov 13 2009

Keywords

Comments

Number of lone-child-avoiding achiral rooted trees with n + 1 vertices, 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 vertex are equal. The Matula-Goebel numbers of these trees are given by A331967. - Gus Wiseman, Feb 07 2020

Examples

			a(12) = 4: [12], [10,2], [9,3], [8,4].
a(14) = 3: [14], [12,2], [8,4,2].
a(18) = 5: [18], [16,2], [15,3], [12,6], [12,4,2].
From _Gus Wiseman_, Jul 13 2018: (Start)
The a(36) = 8 lone-child-avoiding achiral rooted trees with 37 vertices:
  (oooooooooooooooooooooooooooooooooooo)
  ((oo)(oo)(oo)(oo)(oo)(oo)(oo)(oo)(oo)(oo)(oo)(oo))
  ((ooo)(ooo)(ooo)(ooo)(ooo)(ooo)(ooo)(ooo)(ooo))
  ((ooooo)(ooooo)(ooooo)(ooooo)(ooooo)(ooooo))
  ((oooooooo)(oooooooo)(oooooooo)(oooooooo))
  (((ooo)(ooo))((ooo)(ooo))((ooo)(ooo))((ooo)(ooo)))
  ((ooooooooooo)(ooooooooooo)(ooooooooooo))
  ((ooooooooooooooooo)(ooooooooooooooooo))
(End)
		

Crossrefs

The semi-achiral version is A320268.
Matula-Goebel numbers of these trees are A331967.
The semi-lone-child-avoiding version is A331991.
Achiral rooted trees are counted by A003238.

Programs

  • Maple
    with(numtheory):
    a:= proc(n) option remember;
          `if`(n=0, 1, add(a((n-d)/d), d=divisors(n) minus{1}))
        end:
    seq(a(n), n=0..200);  # Alois P. Heinz, Mar 28 2011
  • Mathematica
    a[0] = 1; a[n_] := a[n] = DivisorSum[n, a[(n-#)/#]&, #>1&]; Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Oct 07 2015 *)
  • PARI
    { A167865(n) = if(n==0,return(1)); sumdiv(n,d, if(d>1, A167865((n-d)\d) ) ) }

Formula

a(0) = 1 and for n>=1, a(n) = Sum_{d|n, d>1} a((n-d)/d).
G.f. A(x) satisfies: A(x) = 1 + x^2*A(x^2) + x^3*A(x^3) + x^4*A(x^4) + ... - Ilya Gutkovskiy, May 09 2019

A069916 Number of log-concave compositions (ordered partitions) of n.

Original entry on oeis.org

1, 1, 2, 4, 6, 9, 14, 20, 26, 36, 47, 60, 80, 102, 127, 159, 194, 236, 291, 355, 425, 514, 611, 718, 856, 1009, 1182, 1381, 1605, 1861, 2156, 2496, 2873, 3299, 3778, 4301, 4902, 5574, 6325, 7176, 8116, 9152, 10317, 11610, 13028, 14611, 16354, 18259, 20365
Offset: 0

Views

Author

Pontus von Brömssen, Apr 24 2002

Keywords

Comments

These are compositions with weakly decreasing first quotients, where the first quotients of a sequence are defined as if the sequence were an increasing divisor chain, so for example the first quotients of (6,3,1) are (1/2,1/3). - Gus Wiseman, Mar 16 2021

Examples

			Out of the 8 compositions of 4, only 2+1+1 and 1+1+2 are not log-concave, so a(4)=6.
From _Gus Wiseman_, Mar 15 2021: (Start)
The a(1) = 1 through a(6) = 14 compositions:
  (1)  (2)   (3)    (4)     (5)      (6)
       (11)  (12)   (13)    (14)     (15)
             (21)   (22)    (23)     (24)
             (111)  (31)    (32)     (33)
                    (121)   (41)     (42)
                    (1111)  (122)    (51)
                            (131)    (123)
                            (221)    (132)
                            (11111)  (141)
                                     (222)
                                     (231)
                                     (321)
                                     (1221)
                                     (111111)
(End)
		

Crossrefs

The version for differences instead of quotients is A070211.
A000005 counts constant compositions.
A000009 counts strictly increasing (or strictly decreasing) compositions.
A000041 counts weakly increasing (or weakly decreasing) compositions.
A001055 counts factorizations.
A002843 counts compositions with adjacent parts x <= 2y.
A003238 counts chains of divisors summing to n-1, with strict case A122651.
A003242 counts anti-run compositions.
A074206 counts ordered factorizations.
A167865 counts strict chains of divisors summing to n.

Programs

  • Mathematica
    (* This program is not suitable for computing a large number of terms *)
    compos[n_] := Permutations /@ IntegerPartitions[n] // Flatten[#, 1]&;
    logConcaveQ[p_] := And @@ Table[p[[i]]^2 >= p[[i-1]]*p[[i+1]], {i, 2, Length[p]-1}]; a[n_] := Count[compos[n], p_?logConcaveQ]; Table[an = a[n]; Print["a(", n, ") = ", an]; a[n], {n, 0, 25}] (* Jean-François Alcover, Feb 29 2016 *)
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],GreaterEqual@@Divide@@@Reverse/@Partition[#,2,1]&]],{n,0,15}] (* Gus Wiseman, Mar 15 2021 *)
  • Sage
    def A069916(n) : return sum(all(p[i]^2 >= p[i-1] * p[i+1] for i in range(1, len(p)-1)) for p in Compositions(n)) # Eric M. Schmidt, Sep 29 2013

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

A334996 Irregular triangle read by rows: T(n, m) is the number of ways to distribute Omega(n) objects into precisely m distinct boxes, with no box empty (Omega(n) >= m).

Original entry on oeis.org

0, 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
Offset: 1

Views

Author

Stefano Spezia, May 19 2020

Keywords

Comments

n is the specification number for a set of Omega(n) objects (see Theorem 3 in Beekman's article).
The specification number of a multiset is also called its Heinz number. - Gus Wiseman, Aug 25 2020
From Gus Wiseman, Aug 25 2020: (Start)
For n > 1, T(n,k) is also the number of ordered factorizations of n into k factors > 1. For example, row n = 24 counts the following ordered factorizations (the first column is empty):
24 3*8 2*2*6 2*2*2*3
4*6 2*3*4 2*2*3*2
6*4 2*4*3 2*3*2*2
8*3 2*6*2 3*2*2*2
12*2 3*2*4
2*12 3*4*2
4*2*3
4*3*2
6*2*2
For n > 1, T(n,k) is also the number of strict length-k chains of divisors from n to 1. For example, row n = 36 counts the following chains (the first column is empty):
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

			The triangle T(n, m) begins
  n\m| 0     1     2     3     4
  ---+--------------------------
   1 | 0
   2 | 0     1
   3 | 0     1
   4 | 0     1     1
   5 | 0     1
   6 | 0     1     2
   7 | 0     1
   8 | 0     1     2     1
   9 | 0     1     1
  10 | 0     1     2
  11 | 0     1
  12 | 0     1     4     3
  13 | 0     1
  14 | 0     1     2
  15 | 0     1     2
  16 | 0     1     3     3     1
  ...
From _Gus Wiseman_, Aug 25 2020: (Start)
Row n = 36 counts the following distributions of {1,1,2,2} (the first column is empty):
  {1122}  {1}{122}  {1}{1}{22}  {1}{1}{2}{2}
          {11}{22}  {1}{12}{2}  {1}{2}{1}{2}
          {112}{2}  {11}{2}{2}  {1}{2}{2}{1}
          {12}{12}  {1}{2}{12}  {2}{1}{1}{2}
          {122}{1}  {12}{1}{2}  {2}{1}{2}{1}
          {2}{112}  {1}{22}{1}  {2}{2}{1}{1}
          {22}{11}  {12}{2}{1}
                    {2}{1}{12}
                    {2}{11}{2}
                    {2}{12}{1}
                    {2}{2}{11}
                    {22}{1}{1}
(End)
		

References

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

Crossrefs

Cf. A000007 (1st column), A000012 (2nd column), A001222 (Omega function), A002033 (row sums shifted left), A007318.
A008480 gives rows ends.
A073093 gives row lengths.
A074206 gives row sums.
A112798 constructs the multiset with each specification number.
A124433 is a signed version.
A251683 is the version with zeros removed.
A334997 is the non-strict version.
A337107 is the restriction to factorial numbers.
A001055 counts factorizations.
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 chains of divisors.
A337105 counts strict chains of divisors from n! to 1.

Programs

  • Mathematica
    tau[n_,k_]:=If[n==1,1,Product[Binomial[Extract[Extract[FactorInteger[n],i],2]+k,k],{i,1,Length[FactorInteger[n]]}]]; (* A334997 *)
    T[n_,m_]:=Sum[(-1)^k*Binomial[m,k]*tau[n,m-k-1],{k,0,m-1}]; Table[T[n,m],{n,1,30},{m,0,PrimeOmega[n]}]//Flatten
    (* second program *)
    chc[n_]:=If[n==1,{{}},Prepend[Join@@Table[Prepend[#,n]&/@chc[d],{d,DeleteCases[Divisors[n],1|n]}],{n}]]; (* change {{}} to {} if a(1) = 0 *)
    Table[Length[Select[chc[n],Length[#]==k&]],{n,30},{k,0,PrimeOmega[n]}] (* Gus Wiseman, Aug 25 2020 *)
  • PARI
    TT(n, k) = if (k==0, 1, sumdiv(n, d, TT(d, k-1))); \\ A334997
    T(n, m) = sum(k=0, m-1, (-1)^k*binomial(m, k)*TT(n, m-k-1));
    tabf(nn) = {for (n=1, nn, print(vector(bigomega(n)+1, k, T(n, k-1))););} \\ Michel Marcus, May 20 2020

Formula

T(n, m) = Sum_{k=0..m-1} (-1)^k*binomial(m,k)*tau_{m-k-1}(n), where tau_s(r) = A334997(r, s) (see Theorem 3, Lemma 1 and Lemma 2 in Beekman's article).
Conjecture: Sum_{m=0..Omega(n)} T(n, m) = A002033(n-1) for n > 1.
The above conjecture is true since T(n, m) is also the number of ordered factorizations of n into m factors (see Comments) and A002033(n-1) is the number of ordered factorizations of n. - Stefano Spezia, Aug 21 2025

A342083 Number of chains of strictly inferior divisors from n to 1.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Feb 28 2021

Keywords

Comments

We define a divisor d|n to be strictly inferior if d < n/d. Strictly inferior divisors are counted by A056924 and listed by A341674.
These chains have first-quotients (in analogy with first-differences) that are term-wise > their decapitation (maximum element removed). Equivalently, x > y^2 for all adjacent x, y. For example, the divisor chain q = 60/6/2/1 has first-quotients (10,3,2), which are > (6,2,1), so q is counted under a(60).
Also the number of factorizations of n where each factor is greater than the product of all previous factors.

Examples

			The a(n) chains for n = 2, 6, 12, 24, 42, 48, 60, 72:
  2/1  6/1    12/1    24/1    42/1      48/1      60/1      72/1
       6/2/1  12/2/1  24/2/1  42/2/1    48/2/1    60/2/1    72/2/1
              12/3/1  24/3/1  42/3/1    48/3/1    60/3/1    72/3/1
                      24/4/1  42/6/1    48/4/1    60/4/1    72/4/1
                              42/6/2/1  48/6/1    60/5/1    72/6/1
                                        48/6/2/1  60/6/1    72/8/1
                                                  60/6/2/1  72/6/2/1
                                                            72/8/2/1
The a(n) factorizations for n = 2, 6, 12, 24, 42, 48, 60, 72:
  2  6    12   24    42     48     60      72
     2*3  2*6  3*8   6*7    6*8    2*30    8*9
          3*4  4*6   2*21   2*24   3*20    2*36
               2*12  3*14   3*16   4*15    3*24
                     2*3*7  4*12   5*12    4*18
                            2*3*8  6*10    6*12
                                   2*3*10  2*4*9
                                           2*3*12
		

Crossrefs

The restriction to powers of 2 is A040039.
Not requiring strict inferiority gives A074206 (ordered factorizations).
The weakly inferior version is A337135.
The strictly superior version is A342084.
The weakly superior version is A342085.
The additive version is A342098, or A000929 allowing equality.
A000005 counts divisors.
A001055 counts factorizations.
A003238 counts chains of divisors summing to n-1, with strict case A122651.
A038548 counts inferior (or superior) divisors.
A056924 counts strictly inferior (or strictly superior) divisors.
A067824 counts strict chains of divisors starting with n.
A167865 counts strict chains of divisors > 1 summing to n.
A207375 lists central divisors.
A253249 counts strict chains of divisors.
A334996 counts ordered factorizations by product and length.
A334997 counts chains of divisors of n by length.
A342086 counts chains of divisors with strictly increasing quotients > 1.
- Inferior: A033676, A066839, A072499, A161906.
- Superior: A033677, A070038, A161908.
- Strictly Inferior: A060775, A070039, A333806, A341674.
- Strictly Superior: A048098, A064052, A140271, A238535, A341673.

Programs

  • Mathematica
    cen[n_]:=If[n==1,{{1}},Prepend[#,n]&/@Join@@cen/@Select[Divisors[n],#
    				

Formula

G.f.: x + Sum_{k>=1} a(k) * x^(k*(k + 1)) / (1 - x^k). - Ilya Gutkovskiy, Nov 03 2021

A342084 Number of chains of distinct strictly superior divisors starting with n.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 4, 1, 2, 2, 3, 1, 4, 1, 4, 2, 2, 1, 9, 1, 2, 2, 4, 1, 7, 1, 6, 2, 2, 2, 10, 1, 2, 2, 9, 1, 6, 1, 4, 4, 2, 1, 19, 1, 4, 2, 4, 1, 8, 2, 9, 2, 2, 1, 20, 1, 2, 4, 10, 2, 6, 1, 4, 2, 7, 1, 29, 1, 2, 4, 4, 2, 6, 1, 19, 3, 2, 1, 19, 2
Offset: 1

Views

Author

Gus Wiseman, Feb 28 2021

Keywords

Comments

We define a divisor d|n to be strictly superior if d > n/d. Strictly superior divisors are counted by A056924 and listed by A341673.
These chains have first-quotients (in analogy with first-differences) that are term-wise < their decapitation (maximum element removed). Equivalently, x < y^2 for all adjacent x, y. For example, the divisor chain q = 30/6/3 has first-quotients (5,2), which are < (6,3), so q is counted under a(30).
Also the number of ordered factorizations of n where each factor is less than the product of all previous factors.

Examples

			The a(n) chains for n = 2, 6, 12, 16, 24, 30, 32, 36:
  2  6    12      16      24         30       32         36
     6/3  12/4    16/8    24/6       30/6     32/8       36/9
          12/6    16/8/4  24/8       30/10    32/16      36/12
          12/6/3          24/12      30/15    32/8/4     36/18
                          24/6/3     30/6/3   32/16/8    36/12/4
                          24/8/4     30/10/5  32/16/8/4  36/12/6
                          24/12/4    30/15/5             36/18/6
                          24/12/6                        36/18/9
                          24/12/6/3                      36/12/6/3
                                                         36/18/6/3
The a(n) ordered factorizations for n = 2, 6, 12, 16, 24, 30, 32, 36:
  2  6    12     16     24       30     32       36
     3*2  4*3    8*2    6*4      6*5    8*4      9*4
          6*2    4*2*2  8*3      10*3   16*2     12*3
          3*2*2         12*2     15*2   4*2*4    18*2
                        3*2*4    3*2*5  8*2*2    4*3*3
                        4*2*3    5*2*3  4*2*2*2  6*2*3
                        4*3*2    5*3*2           6*3*2
                        6*2*2                    9*2*2
                        3*2*2*2                  3*2*2*3
                                                 3*2*3*2
		

Crossrefs

The restriction to powers of 2 is A045690, with reciprocal version A040039.
The inferior version is A337135.
The strictly inferior version is A342083.
The superior version is A342085.
The additive version allowing equality is A342094 or A342095.
The additive version is A342096 or A342097.
A000005 counts divisors.
A001055 counts factorizations.
A003238 counts divisibility chains summing to n-1, with strict case A122651.
A038548 counts inferior (or superior) divisors.
A056924 counts strictly inferior (or strictly superior) divisors.
A067824 counts strict chains of divisors starting with n.
A074206 counts strict chains of divisors from n to 1 (also ordered factorizations).
A167865 counts strict chains of divisors > 1 summing to n.
A207375 lists central divisors.
A253249 counts strict chains of divisors.
A334996 counts ordered factorizations by product and length.
A334997 counts chains of divisors of n by length.
- Superior: A033677, A070038, A161908, A341591.
- Strictly Inferior: A060775, A070039, A333806, A341674.
- Strictly Superior: A064052/A048098, A140271, A238535, A341642, A341673.

Programs

  • Mathematica
    ceo[n_]:=Prepend[Prepend[#,n]&/@Join@@ceo/@Select[Most[Divisors[n]],#>n/#&],{n}];
    Table[Length[ceo[n]],{n,100}]

Formula

a(2^n) = A045690(n).

A169594 Number of divisors of n, counting divisor multiplicity in n.

Original entry on oeis.org

1, 2, 2, 4, 2, 4, 2, 6, 4, 4, 2, 7, 2, 4, 4, 9, 2, 7, 2, 7, 4, 4, 2, 10, 4, 4, 6, 7, 2, 8, 2, 11, 4, 4, 4, 12, 2, 4, 4, 10, 2, 8, 2, 7, 7, 4, 2, 14, 4, 7, 4, 7, 2, 10, 4, 10, 4, 4, 2, 13, 2, 4, 7, 15, 4, 8, 2, 7, 4, 8, 2, 16, 2, 4, 7, 7, 4, 8, 2, 14, 9, 4, 2, 13, 4, 4, 4, 10, 2, 13, 4, 7, 4, 4, 4, 17, 2, 7
Offset: 1

Views

Author

Joseph L. Pe, Dec 02 2009

Keywords

Comments

The multiplicity of a divisor d > 1 in n is defined as the largest power i for which d^i divides n; and for d = 1 it is defined as 1.
a(n) is also the sum of the multiplicities of the divisors of n.
In other words, a(n) = 1 + sum of the highest exponents e_i for which each number k_i in range 2 .. n divide n, as {k_i}^{e_i} | n. For nondivisors of n this exponent e_i is 0, for n itself it is 1. - Antti Karttunen, May 20 2017
From Gus Wiseman, Mar 25 2021: (Start)
Also the number of strict chains of divisors ending with n and having constant (equal) first quotients. The case starting with 1 is A089723. For example, the a(1) = 1 through a(12) = 7 chains are:
1 2 3 4 5 6 7 8 9 10 11 12
1|2 1|3 1|4 1|5 1|6 1|7 1|8 1|9 1|10 1|11 1|12
2|4 2|6 2|8 3|9 2|10 2|12
1|2|4 3|6 4|8 1|3|9 5|10 3|12
2|4|8 4|12
1|2|4|8 6|12
3|6|12
(End)
a(n) depends only on the prime signature of n. - David A. Corneth, Mar 28 2021

Examples

			The divisors of 8 are 1, 2, 4, 8 of multiplicity 1, 3, 1, 1, respectively. So a(8) = 1 + 3 + 1 + 1 = 6.
		

Crossrefs

Cf. A168512.
Row sums of A286561, A286563 and A286564.
A001055 counts factorizations (strict: A045778, ordered: A074206).
A057567 counts chains of divisors with weakly increasing first quotients.
A067824 counts strict chains of divisors ending with n.
A253249 counts strict chains of divisors.
A334997 counts chains of divisors of n by length.
A342086 counts chains of divisors with strictly increasing first quotients.
A342496 counts partitions with equal first quotients (strict: A342515, ranking: A342522, ordered: A342495).
A342530 counts chains of divisors with distinct first quotients.
First differences of A078651.

Programs

  • Maple
    a := n -> ifelse(n < 2, 1, 1 + add(padic:-ordp(n, k), k = 2..n)):
    seq(a(n), n = 1..98);  # Peter Luschny, Apr 10 2025
  • Mathematica
    divmult[d_, n_] := Module[{output, i}, If[d == 1, output = 1, If[d == n, output = 1, i = 0; While[Mod[n, d^(i + 1)] == 0, i = i + 1]; output = i]]; output]; dmt0[n_] := Module[{divs, l}, divs = Divisors[n]; l = Length[divs]; Sum[divmult[divs[[i]], n], {i, 1, l}]]; Table[dmt0[i], {i, 1, 40}]
    Table[1 + DivisorSum[n, IntegerExponent[n, #] &, # > 1 &], {n, 98}] (* Michael De Vlieger, May 20 2017 *)
  • PARI
    A286561(n,k) = { my(i=1); if(1==k, 1, while(!(n%(k^i)), i = i+1); (i-1)); };
    A169594(n) = sumdiv(n,d,A286561(n,d)); \\ Antti Karttunen, May 20 2017
    
  • PARI
    a(n) = { if(n == 1, return(1)); my(f = factor(n), u = vecmax(f[, 2]), cf = f, res = numdiv(f) - u + 1); for(i = 2, u, cf[, 2] = f[, 2]\i; res+=numdiv(factorback(cf)) ); res } \\ David A. Corneth, Mar 29 2021
    
  • PARI
    A169594(n) = {my(s=0, k=2); while(k<=n, s+=valuation(n, k); k=k+1); s + 1} \\ Zhuorui He, Aug 28 2025
    
  • Python
    def a286561(n, k):
        i=1
        if k==1: return 1
        while n%(k**i)==0:
            i+=1
        return i-1
    def a(n): return sum([a286561(n, d) for d in divisors(n)]) # Indranil Ghosh, May 20 2017
  • Scheme
    (define (A169594 n) (add (lambda (k) (A286561bi n k)) 1 n))
    ;; Implements sum_{i=lowlim..uplim} intfun(i)
    (define (add intfun lowlim uplim) (let sumloop ((i lowlim) (res 0)) (cond ((> i uplim) res) (else (sumloop (1+ i) (+ res (intfun i)))))))
    ;; For A286561bi see A286561. - Antti Karttunen, May 20 2017
    

Formula

From Friedjof Tellkamp, Feb 29 2024: (Start)
a(n) = A309891(n) + 1.
G.f.: x/(1-x) + Sum_{k>=2, j>=1} x^(k^j)/(1-x^(k^j)).
Dirichlet g.f.: zeta(s) * (1 + Sum_{k>=1} (zeta(k*s) - 1)).
Sum_{n>=1} a(n)/n^2 = (7/24) * Pi^2. (End)

Extensions

Extended by Ray Chandler, Dec 08 2009

A337135 a(1) = 1; for n > 1, a(n) = Sum_{d|n, d <= sqrt(n)} a(d).

Original entry on oeis.org

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

Views

Author

Ilya Gutkovskiy, Nov 21 2020

Keywords

Comments

From Gus Wiseman, Mar 05 2021: (Start)
This sequence counts all of the following essentially equivalent things:
1. Chains of distinct inferior divisors from n to 1, where a divisor d|n is inferior if d <= n/d. Inferior divisors are counted by A038548 and listed by A161906.
2. Chains of divisors from n to 1 whose first-quotients (in analogy with first-differences) are term-wise greater than or equal to their decapitation (maximum element removed). For example, the divisor chain q = 60/4/2/1 has first-quotients (15,2,2), which are >= (4,2,1), so q is counted under a(60).
3. Chains of divisors from n to 1 such that x >= y^2 for all adjacent x, y.
4. Factorizations of n where each factor is greater than or equal to the product of all previous factors.
(End)

Examples

			From _Gus Wiseman_, Mar 05 2021: (Start)
The a(n) chains for n = 1, 2, 4, 12, 16, 24, 36, 60:
  1  2/1  4/1    12/1    16/1      24/1      36/1      60/1
          4/2/1  12/2/1  16/2/1    24/2/1    36/2/1    60/2/1
                 12/3/1  16/4/1    24/3/1    36/3/1    60/3/1
                         16/4/2/1  24/4/1    36/4/1    60/4/1
                                   24/4/2/1  36/6/1    60/5/1
                                             36/4/2/1  60/6/1
                                             36/6/2/1  60/4/2/1
                                                       60/6/2/1
The a(n) factorizations for n = 2, 4, 12, 16, 24, 36, 60:
    2  4    12   16     24     36     60
       2*2  2*6  2*8    3*8    4*9    2*30
            3*4  4*4    4*6    6*6    3*20
                 2*2*4  2*12   2*18   4*15
                        2*2*6  3*12   5*12
                               2*2*9  6*10
                               2*3*6  2*2*15
                                      2*3*10
(End)
		

Crossrefs

Cf. A002033, A008578 (positions of 1's), A068108.
The restriction to powers of 2 is A018819.
Not requiring inferiority gives A074206 (ordered factorizations).
The strictly inferior version is A342083.
The strictly superior version is A342084.
The weakly superior version is A342085.
The additive version is A000929, or A342098 forbidding equality.
A000005 counts divisors, with sum A000203.
A001055 counts factorizations.
A003238 counts chains of divisors summing to n-1, with strict case A122651.
A038548 counts inferior (or superior) divisors.
A056924 counts strictly inferior (or strictly superior) divisors.
A067824 counts strict chains of divisors starting with n.
A167865 counts strict chains of divisors > 1 summing to n.
A207375 lists central divisors.
A253249 counts strict chains of divisors.
A334996 counts ordered factorizations by product and length.
A334997 counts chains of divisors of n by length.
A342086 counts strict factorizations of divisors.
- Inferior: A033676, A066839, A072499, A161906.
- Superior: A033677, A070038, A161908.
- Strictly Inferior: A060775, A070039, A333806, A341674.
- Strictly Superior: A048098, A064052, A140271, A238535, A341673.

Programs

  • Maple
    a:= proc(n) option remember; `if`(n=1, 1, add(
          `if`(d<=n/d, a(d), 0), d=numtheory[divisors](n)))
        end:
    seq(a(n), n=1..128);  # Alois P. Heinz, Jun 24 2021
  • Mathematica
    a[1] = 1; a[n_] := a[n] = DivisorSum[n, a[#] &, # <= Sqrt[n] &]; Table[a[n], {n, 95}]
    (* second program *)
    asc[n_]:=Prepend[#,n]&/@Prepend[Join@@Table[asc[d],{d,Select[Divisors[n],#Gus Wiseman, Mar 05 2021 *)

Formula

G.f.: Sum_{k>=1} a(k) * x^(k^2) / (1 - x^k).
a(2^n) = A018819(n). - Gus Wiseman, Mar 08 2021

A342087 Number of chains of divisors starting with n and having no adjacent parts x <= y^2.

Original entry on oeis.org

1, 2, 2, 2, 2, 4, 2, 4, 2, 4, 2, 6, 2, 4, 4, 4, 2, 6, 2, 6, 4, 4, 2, 8, 2, 4, 4, 6, 2, 8, 2, 6, 4, 4, 4, 8, 2, 4, 4, 8, 2, 10, 2, 6, 6, 4, 2, 12, 2, 6, 4, 6, 2, 10, 4, 8, 4, 4, 2, 14, 2, 4, 6, 6, 4, 10, 2, 6, 4, 8, 2, 16, 2, 4, 6, 6, 4, 10, 2, 12, 4, 4, 2, 14
Offset: 1

Views

Author

Gus Wiseman, Mar 05 2021

Keywords

Comments

An alternative wording: Number of chains of divisors starting with n and having all adjacent parts x > y^2.

Examples

			The chains for n = 1, 2, 6, 12, 24, 42, 48:
   1    2      6        12        24        42          48
        2/1    6/1      12/1      24/1      42/1        48/1
               6/2      12/2      24/2      42/2        48/2
               6/2/1    12/3      24/3      42/3        48/3
                        12/2/1    24/4      42/6        48/4
                        12/3/1    24/2/1    42/2/1      48/6
                                  24/3/1    42/3/1      48/2/1
                                  24/4/1    42/6/1      48/3/1
                                            42/6/2      48/4/1
                                            42/6/2/1    48/6/1
                                                        48/6/2
                                                        48/6/2/1
		

Crossrefs

The restriction to powers of 2 is A018819.
Not requiring strict inferiority gives A067824.
The weakly inferior version is twice A337135.
The case ending with 1 is counted by A342083.
The strictly superior version is A342084.
The weakly superior version is A342085.
The additive version is A342098, or A000929 allowing equality.
A000005 counts divisors, with sum A000203.
A001055 counts factorizations.
A003238 counts chains of divisors summing to n-1, with strict case A122651.
A038548 counts inferior (or superior) divisors.
A056924 counts strictly inferior (or strictly superior) divisors.
A067824 counts strict chains of divisors starting with n.
A074206 counts ordered factorizations.
A167865 counts strict chains of divisors > 1 summing to n.
A253249 counts strict chains of divisors.
A334997 counts chains of divisors of n by length.

Programs

  • Mathematica
    cem[n_]:=Prepend[Prepend[#,n]&/@Join@@cem/@Most[Divisors[n]],{n}];
    Table[Length[Select[cem[n],And@@Thread[Divide@@@Partition[#,2,1]>Rest[#]]&]],{n,30}]

Formula

For n > 1, a(n) = 2*A342083(n).
Showing 1-10 of 43 results. Next