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-9 of 9 results.

A007426 d_4(n), or tau_4(n), the number of ordered factorizations of n as n = rstu.

Original entry on oeis.org

1, 4, 4, 10, 4, 16, 4, 20, 10, 16, 4, 40, 4, 16, 16, 35, 4, 40, 4, 40, 16, 16, 4, 80, 10, 16, 20, 40, 4, 64, 4, 56, 16, 16, 16, 100, 4, 16, 16, 80, 4, 64, 4, 40, 40, 16, 4, 140, 10, 40, 16, 40, 4, 80, 16, 80, 16, 16, 4, 160, 4, 16, 40, 84, 16, 64, 4, 40, 16, 64, 4, 200, 4, 16, 40, 40, 16
Offset: 1

Views

Author

Keywords

Comments

Inverse Möbius transform applied thrice to all 1's sequence; or, Dirichlet convolution of d(n) (A000005).
Let n = Product p_i^e_i. tau (A000005) is tau_2, A007425 is tau_3, this sequence is tau_4, where tau_k(n) (also written as d_k(n)) = Product_i binomial(k-1+e_i, k-1) is the k-th Piltz function. It gives the number of ordered factorizations of n as a product of k terms.
Appears to equal the number of solid partitions of n that can be extended in exactly 4 ways to a solid partition of n + 1 by adding one element. - Wouter Meeussen, Sep 11 2004
Equals row sums of A127172. - Gary W. Adamson, Nov 05 2007

References

  • A. Ivic, The Riemann Zeta-Function, Wiley, NY, 1985, see p. xv.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A007425.
Cf. A127172, A051731, A061202 (partial sums).
Column k=4 of A077592.

Programs

  • Maple
    A007426 := proc(n) local e,j; e := ifactors(n)[2]: product(binomial(3+e[j][2],3), j=1..nops(e)); end;
  • Mathematica
    tau[n_, 1] = 1; tau[n_, k_] := tau[n, k] = Plus @@ (tau[ #, k - 1] & /@ Divisors[n]); Table[ tau[n, 4], {n, 77}] (* Robert G. Wilson v, Nov 02 2005 *)
    a[n_] := DivisorSum[n, DivisorSigma[0, n/#]*DivisorSigma[0, #]&]; Array[a, 80] (* Jean-François Alcover, Dec 01 2015 *)
    tau[1, k_] := 1; tau[n_, k_] := Times @@ (Binomial[Last[#]+k-1, k-1]& /@ FactorInteger[n]); Table[tau[n, 4], {n, 1, 100}] (* Amiram Eldar, Sep 13 2020 *)
  • PARI
    for(n=1,100,print1(sumdiv(n,k,sumdiv(k,x,numdiv(x))),","))
    
  • PARI
    a(n)=sumdiv(n,d,numdiv(n/d)*numdiv(d))
    
  • PARI
    a(n, f=factor(n))=f=f[, 2]; prod(i=1, #f, binomial(f[i]+3, 3)) \\ Charles R Greathouse IV, Oct 28 2017
    
  • PARI
    for(n=1, 100, print1(numerator(direuler(p=2, n, 1/(1-X)^4)[n]), ", ")) \\ Vaclav Kotesovec, May 06 2025
    
  • Python
    from math import prod, comb
    from sympy import factorint
    def A007426(n): return prod(comb(3+e,3) for e in factorint(n).values()) # Chai Wah Wu, Dec 22 2024

Formula

a(n) = Sum_{d dividing n} tau(d)*tau(n/d). - Benoit Cloitre, May 12 2003
Dirichlet g.f.: zeta^4(x).
G.f.: Sum_{k>=1} tau_3(k)*x^k/(1 - x^k). - Ilya Gutkovskiy, Oct 30 2018

Extensions

More terms from Robert G. Wilson v, Nov 02 2005

A061201 Partial sums of A007425: (tau<=)_3(n).

Original entry on oeis.org

1, 4, 7, 13, 16, 25, 28, 38, 44, 53, 56, 74, 77, 86, 95, 110, 113, 131, 134, 152, 161, 170, 173, 203, 209, 218, 228, 246, 249, 276, 279, 300, 309, 318, 327, 363, 366, 375, 384, 414, 417, 444, 447, 465, 483, 492, 495, 540, 546, 564, 573, 591, 594, 624, 633, 663
Offset: 1

Views

Author

Vladeta Jovovic, Apr 21 2001

Keywords

Comments

(tau<=)_k(n) = |{(x_1,x_2,...,x_k): x_1*x_2*...*x_k<=n}|, i.e., tau<=_k(n) is number of solutions to x_1*x_2*...*x_k<=n, x_i > 0.
A061201(n) is the number of 4-tuples (w,x,y,z) having all terms in {1,...,n} and w=x*y*z; see A211795 for a list of related counting sequences. - Clark Kimberling, Apr 28 2012
The formula for Sum_{k=1..n} d3(k) in the Benoit Cloitre article on page 15 is incorrect. For correct asymptotic formula see below or generate it in the Mathematica: Residue[Zeta[s]^3 * n^s/s, {s, 1}] // Expand. - Vaclav Kotesovec, Aug 19 2021

References

  • M. N. Huxley, Area, Lattice Points and Exponential Sums, Oxford, 1996; p. 239.

Crossrefs

Cf. tau_2(n): A000005, tau_3(n): A007425, tau_4(n): A007426, tau_5(n): A061200, tau_6(n): A034695, (unordered) 2-factorizations of n: A038548, (unordered) 3-factorizations of n: A034836, A001055, (tau<=)_2(n): A006218, (tau<=)_4(n): A061202, (tau<=)_5(n): A061203, (tau<=)_6(n): A061204.

Programs

  • Magma
    [&+[NumberOfDivisors(k)*Floor(n/k): k in [1..n]]: n in [1..56]];  // Bruno Berselli, Apr 13 2011
    
  • Maple
    b:= proc(k, n) option remember; uses numtheory;
         `if`(k=1, 1, add(b(k-1, d), d=divisors(n)))
        end:
    a:= proc(n) option remember; `if`(n=0, 0, b(3, n)+a(n-1)) end:
    seq(a(n), n=1..76);  # Alois P. Heinz, Oct 23 2023
  • Mathematica
    a[n_] := Sum[ DivisorSigma[0, k]*Floor[n/k], {k, 1, n}]; Table[a[n], {n, 1, 56}] (* Jean-François Alcover, Sep 20 2011, after Benoit Cloitre *)
    (* Asymptotics: *) n*(Log[n]^2/2 + (3*EulerGamma - 1)*Log[n] + 3*EulerGamma^2 - 3*EulerGamma - 3*StieltjesGamma[1] + 1) (* Vaclav Kotesovec, Sep 09 2018 *)
    Accumulate[a[n_]:=DivisorSum[n, DivisorSigma[0, #]&]; Array[a, 60]] (* Vincenzo Librandi, Jan 12 2020 *)
  • PARI
    a(n)=sum(k=1,n,numdiv(k)*floor(n/k)) \\ Benoit Cloitre, Apr 19 2007
    
  • PARI
    { for (n=1, 1000, write("b061201.txt", n, " ", sum(k=1, n, numdiv(k)*(n\k))) ) } \\ Harry J. Smith, Jul 18 2009
    
  • PARI
    my(N=60, x='x+O('x^N)); Vec(sum(k=1, N, numdiv(k)*x^k/(1-x^k))/(1-x)) \\ Seiichi Manyama, Jul 24 2022
    
  • Python
    from math import isqrt
    from sympy import integer_nthroot
    def A061201(n): return (m:=integer_nthroot(n,3)[0])**3+3*sum(-(s:=isqrt(r:=n//i))**2+(sum(r//k for k in range(1,s+1))<<1)-sum(n//(i*j) for j in range(1,m+1)) for i in range(1,m+1)) # Chai Wah Wu, Oct 23 2023

Formula

(tau<=)k(n) = Sum{i=1..n} tau_k(i).
a(n) = n * ( log(n)^2/2 + (3*g-1)*log(n) + 3*g^2-3*g-3*g1+1 ) + O(sqrt(n)), where g is the Euler-Mascheroni number ~ 0.57721... (see A001620), and g1 is the first Stieltjes constant ~ -0.072816 (see A082633). The determination of the precise size of the error term is an unsolved problem - see references. - Andrew Lelechenko, Apr 15 2011 [corrected by Vaclav Kotesovec, Sep 09 2018]
a(n) = Sum_{k=1..n} A000005(k)*floor(n/k). - Benoit Cloitre, Apr 19 2007
To compute a(n) for huge n (see A180365) in sublinear use a(n) = 3*Sum_{i=1..n3} A006218(n/i) - Sum_{j=1..n3} floor(n/(i*j)) + n3^3, where n3 = floor(n^(1/3)). - Andrew Lelechenko, Apr 15 2011
a(n) = Sum_{k=1..n} Sum_{i=1..n} floor(n/(i*k)). - Wesley Ivan Hurt, Sep 14 2017
G.f.: (1/(1-x)) * Sum_{k>=1} A000005(k) * x^k/(1 - x^k). - Seiichi Manyama, Jul 24 2022

A061200 tau_5(n) = number of ordered 5-factorizations of n.

Original entry on oeis.org

1, 5, 5, 15, 5, 25, 5, 35, 15, 25, 5, 75, 5, 25, 25, 70, 5, 75, 5, 75, 25, 25, 5, 175, 15, 25, 35, 75, 5, 125, 5, 126, 25, 25, 25, 225, 5, 25, 25, 175, 5, 125, 5, 75, 75, 25, 5, 350, 15, 75, 25, 75, 5, 175, 25, 175, 25, 25, 5, 375, 5, 25, 75, 210, 25, 125, 5, 75, 25, 125, 5
Offset: 1

Views

Author

Vladeta Jovovic, Apr 21 2001

Keywords

Crossrefs

Cf. tau_2(n): A000005, tau_3(n): A007425, tau_4(n): A007426, tau_6(n): A034695, (unordered) 2-factorization of n: A038548, (unordered) 3-factorization of n: A034836, A001055, A006218, A061201, A061202, A061203 (partial sums), A061204.
Column k=5 of A077592.

Programs

  • Mathematica
    tau[n_, 1] = 1; tau[n_, k_] := tau[n, k] = Plus @@ (tau[ #, k - 1] & /@ Divisors[n]); Table[ tau[n, 5], {n, 77}] (* Robert G. Wilson v *)
    tau[1, k_] := 1; tau[n_, k_] := Times @@ (Binomial[Last[#]+k-1, k-1]& /@ FactorInteger[n]); Table[tau[n, 5], {n, 1, 100}] (* Amiram Eldar, Sep 13 2020 *)
  • PARI
    for(n=1,100,print1(sumdiv(n,k,sumdiv(k,x,sumdiv(x,y,numdiv(y)))),","))
    
  • PARI
    a(n)=sumdivmult(n,k,sumdivmult(k,x,sumdivmult(x,y,numdiv(y)))) \\ Charles R Greathouse IV, Sep 09 2014
    
  • PARI
    a(n, f=factor(n))=f=f[, 2]; prod(i=1, #f, binomial(f[i]+4, 4)) \\ Charles R Greathouse IV, Oct 28 2017
    
  • PARI
    for(n=1, 100, print1(numerator(direuler(p=2, n, 1/(1-X)^5)[n]), ", ")) \\ Vaclav Kotesovec, May 06 2025
    
  • Python
    from math import prod, comb
    from sympy import factorint
    def A061200(n): return prod(comb(4+e,4) for e in factorint(n).values()) # Chai Wah Wu, Dec 22 2024

Formula

tau_k(n) = |{(x_1,x_2,...,x_k): x_1*x_2*...*x_k=n}|, number of ordered k-factorizations of n.
tau_k(p^m) = (-1)^(k-1)*binomial(-m-1,k-1), p prime.
limit(tau_k(n)/n^epsilon, n=infinity) = 0, for any epsilon>0.
tau_k(n) = Sum_{d|n} tau_(k-1)(d), tau_1(n)=1.
Dirichlet g.f.: (zeta(s))^k.
For explicit formula, see A007425.
G.f.: Sum_{k>=1} tau_4(k)*x^k/(1 - x^k). - Ilya Gutkovskiy, Oct 30 2018

A061203 (tau<=)_5(n).

Original entry on oeis.org

1, 6, 11, 26, 31, 56, 61, 96, 111, 136, 141, 216, 221, 246, 271, 341, 346, 421, 426, 501, 526, 551, 556, 731, 746, 771, 806, 881, 886, 1011, 1016, 1142, 1167, 1192, 1217, 1442, 1447, 1472, 1497, 1672, 1677, 1802, 1807, 1882, 1957, 1982, 1987, 2337, 2352
Offset: 1

Views

Author

Vladeta Jovovic, Apr 21 2001

Keywords

Comments

(tau<=)_k(n) = |{(x_1,x_2,...,x_k): x_1*x_2*...*x_k <= n}|, i.e., (tau<=)_k(n) is number of solutions to x_1*x_2*...*x_k <= n, x_i > 0.
Partial sums of A061200.
Equals row sums of triangle A140705. - Gary W. Adamson, May 24 2008

Crossrefs

Cf. tau_2(n): A000005, tau_3(n): A007425, tau_4(n): A007426, tau_5(n): A061200, tau_6(n): A034695, (unordered) 2-factorizations of n: A038548, (unordered) 3-factorizations of n: A034836, A001055, (tau<=)_2(n): A006218, (tau<=)_3(n): A061201, (tau<=)_4(n): A061202, (tau<=)_6(n): A061204.
Cf. A140705.

Programs

  • Maple
    b:= proc(k, n) option remember; uses numtheory;
         `if`(k=1, 1, add(b(k-1, d), d=divisors(n)))
        end:
    a:= proc(n) option remember; `if`(n=0, 0, b(5, n)+a(n-1)) end:
    seq(a(n), n=1..49);  # Alois P. Heinz, Feb 13 2022
  • Mathematica
    nmax = 50;
    tau4 = Table[DivisorSum[n, DivisorSigma[0, n/#]*DivisorSigma[0, #] &], {n, 1, nmax}];
    Accumulate[Table[Sum[tau4[[d]], {d, Divisors[n]}], {n, nmax}]] (* Vaclav Kotesovec, Sep 10 2018 *)

Formula

(tau<=)k(n) = Sum{i=1..n} tau_k(i).
a(n) = Sum_{k=1..n} tau_{4}(k) * floor(n/k), where tau_{4} is A007426. - Enrique Pérez Herrero, Jan 23 2013
a(n) ~ n*(log(n)^4/24 + (5*g/6 - 1/6)*log(n)^3 + 10*g1^2 + (5*g^2 - 5*g/2 - 5*g1/2 + 1/2)*log(n)^2 + (10*g^3 - 10*g^2 + (5 - 20*g1)*g + 5*g1 + 5*g2/2 - 1)*log(n) + 5*g^4 - 10*g^3 + (10 - 30*g1)*g^2 + (20*g1 + 10*g2 - 5)*g - 5*g1 - 5*g2/2 - 5*g3/6 + 1), where g is the Euler-Mascheroni constant A001620 and g1, g2, g3 are the Stieltjes constants, see A082633, A086279 and A086280. - Vaclav Kotesovec, Sep 10 2018

A061204 (tau<=)_6(n).

Original entry on oeis.org

1, 7, 13, 34, 40, 76, 82, 138, 159, 195, 201, 327, 333, 369, 405, 531, 537, 663, 669, 795, 831, 867, 873, 1209, 1230, 1266, 1322, 1448, 1454, 1670, 1676, 1928, 1964, 2000, 2036, 2477, 2483, 2519, 2555, 2891, 2897, 3113, 3119, 3245, 3371, 3407, 3413
Offset: 1

Views

Author

Vladeta Jovovic, Apr 21 2001

Keywords

Comments

(tau<=)_k(n) = |{(x_1,x_2,...,x_k): x_1*x_2*...*x_k<=n}|, i.e. (tau<=)_k(n) is number of solutions to x_1*x_2*...*x_k<=n, x_i>0.

Crossrefs

Cf. tau_2(n): A000005, tau_3(n): A007425, tau_4(n): A007426, tau_5(n): A061200, tau_6(n): A034695, (unordered) 2-factorizations of n: A038548, (unordered) 3-factorizations of n: A034836, A001055, (tau<=)_2(n): A006218, (tau<=)_3(n): A061201, (tau<=)_4(n): A061202, (tau<=)_5(n): A061203.

Programs

  • Mathematica
    nmax = 50; tau4 = Table[DivisorSum[n, DivisorSigma[0, n/#]*DivisorSigma[0, #] &], {n, 1, nmax}]; tau5 = Table[Sum[tau4[[d]], {d, Divisors[n]}], {n, nmax}]; Accumulate[Table[Sum[tau5[[d]], {d, Divisors[n]}], {n, nmax}]] (* Vaclav Kotesovec, Sep 10 2018 *)

Formula

(tau<=)k(n)=Sum{i=1..n} tau_k(i). a(n)=partial sums of A034695.
a(n) = Sum_{k=1..n} tau_{5}(k) * floor(n/k), where tau_{5} is A061200. - Enrique Pérez Herrero, Jan 23 2013
a(n) ~ n*(log(n)^5/120 + (g/4 - 1/24)*log(n)^4 + (5*g^2/2 - g - g1 + 1/6)*log(n)^3 + (10*g^3 - 15*g^2/2 + (3 - 15*g1)*g + 3*g1 + 3*g2/2 - 1/2)*log(n)^2 + (15*g^4 - 20*g^3 + (15 - 60*g1)*g^2 + (30*g1 + 15*g2 - 6)*g + 15*g1^2 - 6*g1 - 3*g2 - g3 + 1)*log(n) + 6*g^5 - 15*g^4 + (20 - 60*g1)*g^3 + (60*g1 + 30*g2 - 15)*g^2 + (60*g1^2 - 30*g1 - 15*g2 - 5*g3 + 6)*g - 15*g1^2 + g1*(6 - 15*g2) + 3*g2 + g3 + g4/4 - 1), where g is the Euler-Mascheroni constant A001620 and g1, g2, g3, g4 are the Stieltjes constants, see A082633, A086279, A086280 and A086281. - Vaclav Kotesovec, Sep 10 2018

A077593 Table by antidiagonals where T(n,k) = Sum_{i=1..n} T(floor(n/i),k-1) starting with T(n,0)=1 if n>0 and T(0,0)=0.

Original entry on oeis.org

0, 0, 1, 0, 1, 1, 0, 1, 2, 1, 0, 1, 3, 3, 1, 0, 1, 4, 5, 4, 1, 0, 1, 5, 7, 8, 5, 1, 0, 1, 6, 9, 13, 10, 6, 1, 0, 1, 7, 11, 19, 16, 14, 7, 1, 0, 1, 8, 13, 26, 23, 25, 16, 8, 1, 0, 1, 9, 15, 34, 31, 39, 28, 20, 9, 1, 0, 1, 10, 17, 43, 40, 56, 43, 38, 23, 10, 1, 0, 1, 11, 19, 53, 50, 76, 61, 63
Offset: 0

Views

Author

Henry Bottomley, Nov 08 2002

Keywords

Examples

			Rows start:
 0,0,0,0,0,0...;
 1,1,1,1,1,1...;
 1,2,3,4,5,6...;
 1,3,5,7,9,11...;
 1,4,8,13,19,26,...;
 ...
		

Crossrefs

Rows include (with offsets) A000004, A000012, A000027, A005408, A034856, A052905.
Cf. A077593.

Formula

T(n, k) = T(n-1, k) + A077592(n, k). Writing m as Sum_{i} p_i^e_i, T(n, k) = Sum_{m=1..n} Product_{i} C(k+e_i-1, e_i).

A140705 A000012 * A051731^4.

Original entry on oeis.org

1, 5, 1, 9, 1, 1, 19, 5, 1, 1, 23, 5, 1, 1, 1, 39, 9, 5, 1, 1, 1, 43, 9, 5, 1, 1, 1, 1, 63, 19, 5, 5, 1, 1, 1, 1, 73, 19, 9, 5, 1, 1, 1, 1, 1, 89, 23, 9, 5, 5, 1, 1, 1, 1, 1, 93, 23, 9, 5, 5, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Gary W. Adamson, May 24 2008

Keywords

Comments

Row sums = A061203: (1, 6, 11, 26, 31, 56,...).
Left column = A061202: (1, 5, 9, 19, 23, 39,...).

Examples

			First few rows of the triangle are:
1;
5, 1;
9, 1, 1;
19, 5, 1, 1;
23, 5, 1, 1, 1;
39, 9, 5, 1, 1, 1;
43, 9, 5, 1, 1, 1, 1;
63, 19, 5, 5, 1, 1, 1, 1;
...
		

Crossrefs

Formula

A000012 * A051731^4 as infinite lower triangular matrices, where A051731 = the inverse Mobius transform and A000012 = an infinite lower triangular matrix with all 1's.

Extensions

a(60) split into 5,1 by Georg Fischer, Aug 27 2023

A140703 A000012 * A051731^3.

Original entry on oeis.org

1, 4, 1, 7, 1, 1, 13, 4, 1, 1, 16, 4, 1, 1, 1, 25, 7, 4, 1, 1, 1, 28, 7, 4, 1, 1, 1, 1, 38, 13, 4, 4, 1, 1, 1, 1, 44, 13, 7, 4, 1, 1, 1, 1, 1, 53, 16, 7, 4, 4, 1, 1, 1, 1, 1, 56, 16, 7, 4, 4, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Gary W. Adamson, May 24 2008

Keywords

Comments

Row sums = A061202: (1, 5, 9, 19, 23, 39,...).
Leftmost column = A061201: (1, 4, 7, 13, 16, 25, 28,...).

Examples

			First few rows of the triangle are:
1;
4, 1;
7, 1, 1;
13, 4, 1, 1;
16, 4, 1, 1, 1;
25, 7, 4, 1, 1, 1;
28, 7, 4, 1, 1, 1, 1;
38, 13, 4, 4, 1, 1, 1, 1;
...
		

Crossrefs

Formula

A000012 * A051731^3 as infinite lower triangular matrices, where A000012 = [1; 1,1; 1,1,1;...] and A051731 = the inverse Mobius transform.

A189020 a(n) = Sum_{k=1..10^n} tau_4(k), where tau_4 is the number of ordered factorizations into 4 factors (A007426).

Original entry on oeis.org

1, 89, 3575, 93237, 1951526, 35270969, 578262093, 8840109380, 128217432396, 1784942188189, 24045237260214, 315312623543840, 4042957241191810, 50862246063060180, 629513636928477232, 7681900592647818929, 92587253467765253144, 1103781870246459696784, 13031388731053572679450, 152516435040764735691556, 1771079109308495896176156
Offset: 0

Views

Author

Andrew Lelechenko, Apr 15 2011

Keywords

Comments

Using that tau_4 = tau_2 ** tau_2, where ** means Dirichlet convolution and tau_2 is (A000005), one can calculate a(n) faster than in O(10^n) operations - namely in O(10^(3n/4)) or even in O(10^(2n/3)). See links for details.

Crossrefs

Cf. A057494 - partial sums up to 10^n of the divisors function tau_2 (A000005), A180361 - of the unitary divisors function tau_2* (A034444), A180365 - of the 3-divisors function tau_3 (A007425).
Also see A072692 for such sums of the sum of divisors function (A000203), A084237 for sums of Moebius function (A008683), A064018 for sums of Euler totient function (A000010).

Formula

a(n) = A061202(10^n) = Sum_{k=1..10^n} A007426(n).

Extensions

a(16)-a(20) from Henri Lifchitz, Feb 05 2025
Showing 1-9 of 9 results.