A000005 -id:A000005 - cp's OEIS frontend

A101509 Binomial transform of tau(n) (see A000005).

1, 3, 7, 16, 35, 75, 159, 334, 696, 1442, 2976, 6123, 12562, 25706, 52492, 107014, 217877, 443061, 899957, 1826078, 3701783, 7498261, 15178255, 30706320, 62085915, 125465715, 253415981, 511608490, 1032427637, 2082680887, 4199956101, 8467124805, 17064784905, 34382825363, 69256687719, 139465867773

Offset: 0

Paul Barry, Dec 05 2004

Row sums of A101508.

Also: Number of matrices with positive integer coefficients such that the sum of all entries equals n+1, cf. link "Partitions and A101509". - M. F. Hasler, Jan 14 2009

			From _Gus Wiseman_, Jan 16 2019: (Start)
The a(3) = 16 ways to arrange the parts of an integer partition of 4 into a matrix:
  [4] [1 3] [3 1] [2 2] [1 1 2] [1 2 1] [2 1 1] [1 1 1 1]
.
  [1] [3] [2] [1 1]
  [3] [1] [2] [1 1]
.
  [1] [1] [2]
  [1] [2] [1]
  [2] [1] [1]
.
  [1]
  [1]
  [1]
  [1]
(End)

M. F. Hasler, Table of n, a(n) for n = 0..500
L. Manor, M. F. Hasler, Partitions and A101509. SeqFan list, Jan 14 2009
N. J. A. Sloane, Transforms

Cf. A000005 (tau), A101508, A160399.

Cf. A000219, A047966, A053529, A319066, A323300, A323301, A323307, A323429.

bintr:= proc(p) proc(n) add(p(k) *binomial(n, k), k=0..n) end end:
a:= bintr(n-> numtheory[tau](n+1)):
seq(a(n), n=0..40);  # Alois P. Heinz, Jan 30 2011

a[n_] := Sum[DivisorSigma[0, k+1]*Binomial[n, k], {k, 0, n}]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Feb 18 2017 *)

A101509(n) = sum( k=0,n, numdiv(k+1)*binomial(n,k)) \\ M. F. Hasler, Jan 14 2009

a(n) = Sum_{k=0..n, Sum_{i=0..n, if(mod(i+1, k+1)=0, binomial(n, i), 0)}}.
G.f.: 1/x * Sum_{n>=1} z^n/(1-z^n) (Lambert series) where z=x/(1-x). - Joerg Arndt, Jan 30 2011
a(n) ~ 2^n * (log(n/2) + 2*gamma), where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, Mar 07 2020

A184389 a(n) = Sum_{k=1..tau(n)} k, where tau is the number of divisors of n (A000005).

Original entry on oeis.org

1, 3, 3, 6, 3, 10, 3, 10, 6, 10, 3, 21, 3, 10, 10, 15, 3, 21, 3, 21, 10, 10, 3, 36, 6, 10, 10, 21, 3, 36, 3, 21, 10, 10, 10, 45, 3, 10, 10, 36, 3, 36, 3, 21, 21, 10, 3, 55, 6, 21, 10, 21, 3, 36, 10, 36, 10, 10, 3, 78, 3, 10, 21, 28, 10, 36, 3, 21, 10, 36, 3, 78

Offset: 1

Jaroslav Krizek, Jan 12 2011

nonn

Length of row n in triangle A187207. - Omar E. Pol, Aug 07 2011

Number of pairs of even divisors of 2n, (d1,d2), such that d1<=d2. - Wesley Ivan Hurt, Aug 24 2020

			For n = 4; tau(4) = 3; a(4) = 1+2+3 = 6.

Enrique Pérez Herrero, Table of n, a(n) for n = 1..5000

Cf. A000005 (tau), A000217 (triangular numbers).
Cf. A184387, A184388, A184390, A184391, A130674.
Cf. A035116, A066446.
Cf. A187207, A135539, A337362, A129308.

a184389 = a000217 . a000005'  -- Reinhard Zumkeller, Sep 08 2015

A184389 := proc(n)
    A000217(numtheory[tau](n)) ;
end proc: # R. J. Mathar, Oct 04 2014

((#+1)#)/2&/@DivisorSigma[0,Range[80]] (* Harvey P. Dale, Feb 27 2013 *)

a(n) = my(nd=numdiv(n)); nd*(nd+1)/2; \\ Michel Marcus, Jun 25 2016

a(n) = A000217(A000005(n)) = (1/2)*A000005(n)*(A000005(n)+1).
a(n) = A066446(n) + A000005(n) = A035116(n) - A066446(n). - Reinhard Zumkeller, Sep 08 2015
Dirichlet g.f.: zeta(s)^2*(zeta(s)^2 + zeta(2*s))/(2*zeta(2*s)). - Ilya Gutkovskiy, Jun 25 2016
a(n) = Sum_{d1|(2*n), d2|(2*n), d1 and d2 even, d1<=d2} 1. - Wesley Ivan Hurt, Aug 24 2020
a(n) = Sum_{d|n} A018892(d). - Daniel Suteu, Jan 08 2021
a(n) = Sum_{d|n} A135539(n,d). - Ridouane Oudra, May 29 2025
a(n) = A337362(n) + A129308(n). - Ridouane Oudra, May 30 2025

A204270 a(n) = tau(n)*Pell(n), where tau(n) = A000005(n), the number of divisors of n.

Original entry on oeis.org

1, 4, 10, 36, 58, 280, 338, 1632, 2955, 9512, 11482, 83160, 66922, 323128, 780100, 2354160, 2273378, 16465260, 13250218, 95966568, 154455860, 372889432, 450117362, 4346717760, 3935214363, 12667263848, 30581480180, 110745336312, 89120964298

Offset: 1

Paul D. Hanna, Jan 14 2012

nonn

Compare g.f. to the Lambert series identity: Sum_{n>=1} x^n/(1-x^n) = Sum_{n>=1} tau(n)*x^n.
Related identities:
(1) Sum_{n>=1} n^k*Pell(n)*x^n/(1 - A002203(n)*x^n + (-1)^n*x^(2*n)) = Sum_{n>=1} sigma_{k}(n)*Pell(n)*x^n for k>=0.
(2) Sum_{n>=1} phi(n)*Pell(n)*x^n/(1 - A002203(n)*x^n + (-1)^n*x^(2*n)) = Sum_{n>=1} n*Pell(n)*x^n.
(3) Sum_{n>=1} moebius(n)*Pell(n)*x^n/(1 - A002203(n)*x^n + (-1)^n*x^(2*n)) = x.
(4) Sum_{n>=1} lambda(n)*Pell(n)*x^n/(1 - A002203(n)*x^n + (-1)^n*x^(2*n)) = Sum_{n>=1} Pell(n^2)*x^(n^2).

			G.f.: A(x) = 1 + 4*x + 10*x^2 + 36*x^3 + 58*x^4 + 280*x^5 + 338*x^6 +...
where A(x) = x/(1-2*x-x^2) + 2*x^2/(1-6*x^2+x^4) + 5*x^3/(1-14*x^3-x^6) + 12*x^4/(1-34*x^4+x^8) + 29*x^5/(1-82*x^5-x^10) + 70*x^6/(1-198*x^6+x^12) +...+ Pell(n)*x^n/(1 - A002203(n)*x^n + (-1)^n*x^(2*n)) +...

G. C. Greubel, Table of n, a(n) for n = 1..1000

Cf. A203847, A204271, A204272, A204273, A204274, A204275, A000005 (tau), A002203, A000045.

Table[DivisorSigma[0, n] Fibonacci[n, 2], {n, 1, 50}] (* G. C. Greubel, Jan 05 2018 *)

/* Subroutines used in PARI programs below: */
{Pell(n)=polcoeff(x/(1-2*x-x^2+x*O(x^n)), n)}
{A002203(n)=polcoeff(2*(1-x)/(1-2*x-x^2+x*O(x^n)), n)}

PARI
```
{a(n)=sigma(n,0)*Pell(n)}
```

{a(n)=polcoeff(sum(m=1,n,Pell(m)*x^m/(1-A002203(m)*x^m+(-1)^m*x^(2*m)+x*O(x^n))),n)}

G.f.: Sum_{n>=1} Pell(n)*x^n/(1 - A002203(n)*x^n + (-1)^n*x^(2*n)) = Sum_{n>=1} tau(n)*Pell(n)*x^n, where Pell(n) = A000129(n) and A002203 is the companion Pell numbers.

A064491 a(1) = 1, a(n+1) = a(n) + tau(a(n)), where tau(n) (A000005) is the number of divisors of n.

Original entry on oeis.org

1, 2, 4, 7, 9, 12, 18, 24, 32, 38, 42, 50, 56, 64, 71, 73, 75, 81, 86, 90, 102, 110, 118, 122, 126, 138, 146, 150, 162, 172, 178, 182, 190, 198, 210, 226, 230, 238, 246, 254, 258, 266, 274, 278, 282, 290, 298, 302, 306, 318, 326, 330, 346, 350, 362, 366, 374

Offset: 1

Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Oct 04 2001

a(n) = partial sums of A165930(n). [Jaroslav Krizek, Sep 30 2009]

Claudia Spiro, Problem proposed at West Coast Number Theory Meeting, 1977. [If you change the starting term, does the resulting sequence always join this one? Does the parity of terms change infinitely often?] - From N. J. A. Sloane, Jan 11 2013

T. D. Noe, Table of n, a(n) for n=1..1000

Cf. A000005, A165930, A165494, A165495, A165496, A165497, A165498, A165499.

a064491 n = a064491_list !! (n-1)
a064491_list = iterate a062249 1  -- Reinhard Zumkeller, Mar 29 2014

a[n_] := a[n] = a[n - 1] + DivisorSigma[0, a[n - 1]]; a[1] = 1; Table[a[n], {n, 1, 57}] (* Jean-François Alcover, Oct 11 2012 *)
NestList[#+DivisorSigma[0,#]&,1,60] (* Harvey P. Dale, Feb 05 2017 *)

{ for (n=1, 1000, if (n>1, a+=numdiv(a), a=1); write("b064491.txt", n, " ", a) ) } \\ Harry J. Smith, Sep 16 2009

from itertools import islice
from sympy import divisor_count
def A064491gen(): # generator of terms
    n = 1
    yield n
    while True:
        n += divisor_count(n)
        yield n
A064491_list = list(islice(A064491gen(),20)) # Chai Wah Wu, Dec 13 2021

It seems likely that there exist constants c_1 and c_2 such that c_1*n*log(n) < a(n) < c_2*n*log(n) for all sufficiently large n. - Franklin T. Adams-Watters, Jun 25 2008

a(n+1) = A062249(a(n)). - Reinhard Zumkeller, Mar 29 2014

Beginning of sequence corrected by T. D. Noe, Sep 13 2007

A082284 a(n) = smallest number k such that k - tau(k) = n, or 0 if no such number exists, where tau(n) = the number of divisors of n (A000005).

Original entry on oeis.org

1, 3, 6, 5, 8, 7, 9, 0, 0, 11, 14, 13, 18, 0, 20, 17, 24, 19, 22, 0, 0, 23, 25, 27, 0, 0, 32, 29, 0, 31, 34, 35, 40, 0, 38, 37, 0, 0, 44, 41, 0, 43, 46, 0, 50, 47, 49, 51, 56, 0, 0, 53, 0, 57, 58, 0, 0, 59, 62, 61, 72, 65, 68, 0, 0, 67, 0, 0, 0, 71, 74, 73, 84, 77, 0, 0, 81, 79, 82, 0, 88

Offset: 0

Amarnath Murthy, Apr 14 2003

nonn

a(p-2) = p for odd primes p.

Antti Karttunen, Table of n, a(n) for n = 0..124340

Column 1 of A265751.
Cf. A000005, A002182, A002183, A049820, A060990, A261100.
Cf. A262686 (the largest such number), A262511 (positions where these are equal and nonzero).
Cf. A266114 (same sequence sorted into ascending order, with zeros removed).
Cf. A266115 (positive numbers missing from this sequence).
Cf. A266110 (number of iterations before zero is reached), A266116 (final nonzero value reached).
Cf. also tree A263267 and its illustration.

N:= 1000: # to get a(0) .. a(N)
V:= Array(0..N):
for k from 1 to 2*(N+1) do
  v:= k - numtheory:-tau(k);
  if v <= N and V[v] = 0 then V[v]:= k fi
od:
seq(V[n],n=0..N); # Robert Israel, Dec 21 2015

Table[k = 1; While[k - DivisorSigma[0, k] != n && k <= 2 (n + 1), k++]; If[k > 2 (n + 1), 0, k], {n, 0, 80}] (* Michael De Vlieger, Dec 22 2015 *)

allocatemem(123456789);
uplim1 = 2162160 + 320; \\ = A002182(41) + A002183(41).
uplim2 = 2162160;
v082284 = vector(uplim1);
A082284 = n -> if(!n,1,v082284[n]);
for(n=1, uplim1, k = n-numdiv(n); if((0 == A082284(k)), v082284[k] = n));
for(n=0, 124340, write("b082284.txt", n, " ", A082284(n)));
\\ Antti Karttunen, Dec 21 2015

(define (A082284 n) (if (zero? n) 1 (let ((u (+ n (A002183 (+ 2 (A261100 n)))))) (let loop ((k n)) (cond ((= (A049820 k) n) k) ((> k u) 0) (else (loop (+ 1 k))))))))
;; Antti Karttunen, Dec 21 2015

Other identities and observations. For all n >= 0:

a(n) <= A262686(n).

More terms from David Wasserman, Aug 31 2004

A203847 a(n) = tau(n)*Fibonacci(n), where tau(n) = A000005(n), the number of divisors of n.

Original entry on oeis.org

1, 2, 4, 9, 10, 32, 26, 84, 102, 220, 178, 864, 466, 1508, 2440, 4935, 3194, 15504, 8362, 40590, 43784, 70844, 57314, 370944, 225075, 485572, 785672, 1906866, 1028458, 6656320, 2692538, 13069854, 14098312, 22811548, 36909860, 134373168, 48315634, 156352676, 252983944

Offset: 1

Paul D. Hanna, Jan 11 2012

nonn

Compare g.f. to the Lambert series identity: Sum_{n>=1} x^n/(1-x^n) = Sum_{n>=1} tau(n)*x^n.
Related identities:
(1) Sum_{n>=1} n^k*Fibonacci(n)*x^n/(1 - Lucas(n)*x^n + (-1)^n*x^(2*n)) = Sum_{n>=1} sigma_{k}(n)*Fibonacci(n)*x^n for k>=0.
(2) Sum_{n>=1} phi(n)*Fibonacci(n)*x^n/(1 - Lucas(n)*x^n + (-1)^n*x^(2*n)) = Sum_{n>=1} n*Fibonacci(n)*x^n.
(3) Sum_{n>=1} moebius(n)*Fibonacci(n)*x^n/(1 - Lucas(n)*x^n + (-1)^n*x^(2*n)) = x.
(4) Sum_{n>=1} lambda(n)*Fibonacci(n)*x^n/(1 - Lucas(n)*x^n + (-1)^n*x^(2*n)) = Sum_{n>=1} Fibonacci(n^2)*x^(n^2).

			G.f.: A(x) = x + 2*x^2 + 4*x^3 + 9*x^4 + 10*x^5 + 32*x^6 + 26*x^7 +...
where A(x) = x/(1-x-x^2) + x^2/(1-3*x^2+x^4) + 2*x^3/(1-4*x^3-x^6) + 3*x^4/(1-7*x^4+x^8) + 5*x^5/(1-11*x^5-x^10) + 8*x^6/(1-18*x^6+x^12) +...+ Fibonacci(n)*x^n/(1 - Lucas(n)*x^n + (-1)^n*x^(2*n)) +...

G. C. Greubel, Table of n, a(n) for n = 1..2500

Cf. A203848, A203849, A203838, A204060, A000005 (tau), A000204 (Lucas), A000045.
Cf. A205507, A205882, A205963, A205964, A205965, A205966.
Cf. A205967, A205968, A205969, A205970, A205971.
Cf. A205972, A205973, A205974, A205975, A205976.
Cf. A204291, A203801, A203850, A203860, A203861.

Table[DivisorSigma[0, n]*Fibonacci[n], {n, 50}] (* G. C. Greubel, Jul 17 2018 *)

PARI
```
{a(n)=sigma(n,0)*fibonacci(n)}
```

{Lucas(n)=fibonacci(n-1)+fibonacci(n+1)}
{a(n)=polcoeff(sum(m=1,n,fibonacci(m)*x^m/(1-Lucas(m)*x^m+(-1)^m*x^(2*m)+x*O(x^n))),n)}

a(n) = numdiv(n)*fibonacci(n); \\ Michel Marcus, Jul 18 2018

G.f.: Sum_{n>=1} Fibonacci(n)*x^n/(1 - Lucas(n)*x^n + (-1)^n*x^(2*n)) = Sum_{n>=1} tau(n)*Fibonacci(n)*x^n, where Lucas(n) = A000204(n).

A284005 a(0) = 1, and for n > 1, a(n) = (1 + A000120(n))*a(floor(n/2)); also a(n) = A000005(A283477(n)).

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 18, 24, 16, 24, 36, 48, 54, 72, 96, 120, 32, 48, 72, 96, 108, 144, 192, 240, 162, 216, 288, 360, 384, 480, 600, 720, 64, 96, 144, 192, 216, 288, 384, 480, 324, 432, 576, 720, 768, 960, 1200, 1440, 486, 648, 864, 1080, 1152, 1440, 1800, 2160, 1536, 1920, 2400, 2880, 3000

Offset: 0

Antti Karttunen, Mar 18 2017

Antti Karttunen, Table of n, a(n) for n = 0..8191
Index entries for sequences related to binary expansion of n

Cf. A000005, A283477, A284001.

Similar recurrences: A124758, A243499, A329369, A341392.

Table[DivisorSigma[0, #] &@ Apply[Times, Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e == 1 :> {Times @@ Prime@ Range@ PrimePi@ p, e}]] &[Times @@ Prime@ Flatten@ Position[#, 1] &@ Reverse@ IntegerDigits[n, 2]], {n, 0, 71}] (* Michael De Vlieger, Mar 18 2017 *)

A284005(n) = numdiv(A283477(n)); \\ edited by Michel Marcus, May 01 2019, M. F. Hasler, Nov 10 2019

a(n) = my(k=if(n,logint(n,2)),s=1); prod(i=0,k, s+=bittest(n,k-i)); \\ Kevin Ryde, Jan 20 2021

(define (A284005 n) (A000005 (A283477 n)))

a(n) = A000005(A283477(n)).
Conjecture: a(n) = 2*a(f(n)) + Sum_{k=0..floor(log_2(n))-1} a(f(n) + 2^k*(1 - T(n,k))) for n > 1 with a(0) = 1, a(1) = 2, f(n) = A053645(n), T(n,k) = floor(n/2^k) mod 2. - Mikhail Kurkov, Nov 10 2019
From Mikhail Kurkov, Aug 23 2021: (Start)
a(2n+1) = a(n) + a(2n) for n >= 0.
a(2n) = a(n) + a(2n - 2^A007814(n)) for n > 0 with a(0) = 1. (End)
Conjecture: a(n) = Sum_{k=0..n} (binomial(n, k) mod 2)*A329369(k). In other words, this sequence is modulo 2 binomial transform of A329369. - Mikhail Kurkov, Mar 10 2023
Conjecture: a(2^m*(2n+1)) = Sum_{k=0..m+1} binomial(m+1, k)*a(2^k*n) for m >= 0, n >= 0 with a(0) = 1. - Mikhail Kurkov, Apr 24 2023

Made Mikhail Kurkov's Nov 10 2019 formula the new primary name of this sequence - Antti Karttunen, Dec 30 2020

A062011 a(n) = 2tau(n) = 2A000005(n).

Original entry on oeis.org

2, 4, 4, 6, 4, 8, 4, 8, 6, 8, 4, 12, 4, 8, 8, 10, 4, 12, 4, 12, 8, 8, 4, 16, 6, 8, 8, 12, 4, 16, 4, 12, 8, 8, 8, 18, 4, 8, 8, 16, 4, 16, 4, 12, 12, 8, 4, 20, 6, 12, 8, 12, 4, 16, 8, 16, 8, 8, 4, 24, 4, 8, 12, 14, 8, 16, 4, 12, 8, 16, 4, 24, 4, 8, 12, 12, 8, 16, 4, 20, 10, 8, 4, 24, 8, 8, 8, 16

Offset: 1

Ahmed Fares (ahmedfares(AT)my-deja.com), Jul 12 2001

Old definition was "Number of cyclic subgroups of the group C_n X C_2 (where C_n is the cyclic group with n elements)."
More generally, the number of cyclic subgroups of the group C_n X C_m is Sum_{i|n, j|m} phi(i)*phi(j)/phi(lcm(i,j)), where phi=Euler totient function, cf. A000010. - Vladeta Jovovic, Jul 15 2001
Number of divisors of p*n, where p is any prime not dividing n. - Reinhard Zumkeller, May 17 2006
From Enrique Pérez Herrero, Jul 21 2011: (Start)
If p(x) is a polynomial with integer coefficients, and if r is an integer zero of p(x), then r is a divisor of the constant term c_0 of p(x). Under this theorem, p(x) can have a(c_0) possible integer roots.
a(n) is the number of integer divisors of n, while A000005(n) is the number of positive divisors. (End)
Number of solutions to the Diophantine equation i*j = n*i + j. - Robert G. Wilson v, Apr 10 2019
a(n) is also the number of times n appears in the triangle A333119, or equivalently, the number of positive integer solutions of the equation A333119(x, y) = n for y < x. - Stefano Spezia, Oct 05 2022

Harry J. Smith, Table of n, a(n) for n = 1..1000
Omar E. Pol, Illustration of initial terms
Index entries for sequences related to groups

Cf. A000005, A060648, A060710.
Cf. A087560, A119416.
Cf. A333119.

A062011[n_] := 2*DivisorSigma[0,n]; Array[A062011,50] (* Enrique Pérez Herrero, Jul 15 2011 *)

for (n=1, 1000, write("b062011.txt", n, " ", 2*numdiv(n))) \\ Harry J. Smith, Jul 29 2009

a(n) = A000005(A087560(n)) = A000005(A119416(n)). - Reinhard Zumkeller, May 17 2006

L.g.f.: -log(Product_{k>=1} (1 - x^k)^(2/k)) = Sum_{n>=1} a(n)*x^n/n. - Ilya Gutkovskiy, Mar 18 2018

More terms from Vladeta Jovovic, Jul 14 2001

Edited by N. J. A. Sloane, Sep 20 2018, replacing old definition (which was of course correct) with a simple formula.

A066843 a(n) = Product_{k=1..n} d(k); d(k) = A000005(k) is the number of positive divisors of k.

Original entry on oeis.org

1, 1, 2, 4, 12, 24, 96, 192, 768, 2304, 9216, 18432, 110592, 221184, 884736, 3538944, 17694720, 35389440, 212336640, 424673280, 2548039680, 10192158720, 40768634880, 81537269760, 652298158080, 1956894474240, 7827577896960, 31310311587840, 187861869527040

Offset: 0

Leroy Quet, Jan 20 2002

nonn

a(n) is also the determinant of the symmetric n X n matrix M defined by M(i,j) = d_3(gcd(i,j)) for 1 <= i,j <= n, where d_3(n) is A007425. - Enrique Pérez Herrero, Aug 12 2011

a(n) is the number of integer sequences of length n where a(m) divides m for every term. - Franklin T. Adams-Watters, Oct 29 2017

Alois P. Heinz, Table of n, a(n) for n = 0..1310 (terms n = 1..200 from Harry J. Smith)
Antal Bege, Hadamard product of GCD matrices, Acta Univ. Sapientiae, Mathematica, 1, 1 (2009) 43-49
Mathoverflow, Product of tau(k), 2015.
Ramanujan's Papers, Some formulas in the analytic theory of numbers, Messenger of Mathematics, XLV, 1916, 81-84, Formula (10).

Cf. A000005, A001088, A066780.

with(numtheory):seq(mul(tau(k),k=1..n), n=0..26); # Zerinvary Lajos, Jan 11 2009
with(numtheory):a[0]:=1: for n from 2 to 26 do a[n]:=a[n-1]*tau(n) od: seq(a[n], n=0..26); # Zerinvary Lajos, Mar 21 2009

A066843[n_] := Product[DivisorSigma[0,i], {i,1,n}]; Array[A066843,20] (* Enrique Pérez Herrero, Aug 12 2011 *)
FoldList[Times, Array[DivisorSigma[0, #] &, 27]] (* Michael De Vlieger, Nov 01 2017 *)

{ p=1; for (n=1, 200, p*=length(divisors(n)); write("b066843.txt", n, " ", p) ) } \\ Harry J. Smith, Apr 01 2010

a(n) = Product_{p=primes<=n} Product_{1<=k<=log(n)/log(p)} (1 +1/k)^floor(n/p^k). - Leroy Quet, Mar 20 2007

a(n) = Product_{k=1..n} Product_{p prime<=n} (v_p(k) + 1), where v_p(k) is the exponent of highest power of p dividing k. - Ridouane Oudra, Apr 15 2024

a(0)=1 prepended by Alois P. Heinz, Jul 19 2023

A262686 a(n) = largest number k such that k - d(k) = n, or 0 if no such number exists, where d(n) = the number of divisors of n (A000005).

Original entry on oeis.org

2, 4, 6, 5, 8, 7, 12, 0, 0, 11, 14, 16, 18, 0, 20, 17, 24, 21, 22, 0, 0, 23, 30, 27, 0, 0, 32, 36, 0, 33, 34, 35, 40, 0, 42, 39, 0, 0, 48, 45, 0, 43, 46, 0, 50, 47, 54, 51, 60, 0, 0, 55, 0, 57, 58, 0, 0, 64, 66, 61, 72, 65, 70, 0, 0, 69, 0, 0, 0, 75, 80, 73, 84, 77, 0, 0, 81, 79, 90, 0, 88, 85, 86, 87, 96, 0, 92, 91, 0, 93, 94, 100, 98, 99, 102, 97, 108, 105, 0, 101

Offset: 0

Antti Karttunen, Sep 28 2015

nonn

Antti Karttunen, Table of n, a(n) for n = 0..10080

Cf. A000005, A002183, A049820, A060990, A261100, A262512.

Cf. also A082284 (the smallest such number), A262511 (positions where these are equal and nonzero).

Table[k = 2 n + 3; While[Nor[k - DivisorSigma[0, k] == n, k == 0], k--]; k, {n, 0, 99}] (* Michael De Vlieger, Sep 29 2015 *)

(definec (A262686 n) (if (zero? n) 2 (let ((u (+ n (A002183 (+ 2 (A261100 n)))))) (let loop ((k u)) (cond ((= (A049820 k) n) k) ((< k n) 0) (else (loop (- k 1))))))))

cp's OEIS Frontend

A101509 Binomial transform of tau(n) (see A000005).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A184389 a(n) = Sum_{k=1..tau(n)} k, where tau is the number of divisors of n (A000005).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

A204270 a(n) = tau(n)*Pell(n), where tau(n) = A000005(n), the number of divisors of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

PARI

Formula

A064491 a(1) = 1, a(n+1) = a(n) + tau(a(n)), where tau(n) (A000005) is the number of divisors of n.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Python

Formula

Extensions

A082284 a(n) = smallest number k such that k - tau(k) = n, or 0 if no such number exists, where tau(n) = the number of divisors of n (A000005).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Scheme

Formula

Extensions

A203847 a(n) = tau(n)*Fibonacci(n), where tau(n) = A000005(n), the number of divisors of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

A062011 a(n) = 2tau(n) = 2A000005(n).