A064987 -id:A064987 - cp's OEIS frontend

A069097 Moebius transform of A064987, n*sigma(n).

1, 5, 11, 22, 29, 55, 55, 92, 105, 145, 131, 242, 181, 275, 319, 376, 305, 525, 379, 638, 605, 655, 551, 1012, 745, 905, 963, 1210, 869, 1595, 991, 1520, 1441, 1525, 1595, 2310, 1405, 1895, 1991, 2668, 1721, 3025, 1891, 2882, 3045, 2755, 2255, 4136, 2737

Offset: 1

Benoit Cloitre, Apr 05 2002

Equals A127569 * [1, 2, 3, ...]. - Gary W. Adamson, Jan 19 2007
Equals row sums of triangle A143309 and of triangle A143312. - Gary W. Adamson, Aug 06 2008
Dirichlet convolution of A000290 and A000010 (see Jovovic formula). - R. J. Mathar, Feb 03 2011

Amiram Eldar, Table of n, a(n) for n = 1..10000
Wenchang Chu, Jordan's totient function and trigonometric sums, Studia Scientiarum Mathematicarum Hungarica 57 (1), 40-53 (2020).
Nikolai Osipov, Problem 12003, Amer. Math. Monthly, 124(8) (2017), p. 754.

Column 2 of A343510.
Cf. A018804, A033457, A068963, A064987, A127569, A143309, A143312, A360428.
For Sum_{k = 1..n} gcd(k,n)^m see A018804 (m = 1), A343497 (m = 3), A343498 (m = 4) and A343499 (m = 5).

A069097[n_]:=n^2*Plus @@((EulerPhi[#]/#^2)&/@ Divisors[n]); Array[A069097, 100] (* Enrique Pérez Herrero, Feb 25 2012 *)
f[p_, e_] := p^(e-1)*(p^e*(p+1)-1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 18 2020 *)

for(n=1,100,print1((sumdiv(n,k,k*sigma(k)*moebius(n/k))),","))

a(n) = Sum_{d|n} d^2*phi(n/d). - Vladeta Jovovic, Jul 31 2002
a(n) = Sum_{k=1..n} gcd(n, k)^2. - Vladeta Jovovic, Aug 27 2003
Dirichlet g.f.: zeta(s-2)*zeta(s-1)/zeta(s). - R. J. Mathar, Feb 03 2011
a(n) = n*Sum_{d|n} J_2(d)/d, where J_2 is A007434. - Enrique Pérez Herrero, Feb 25 2012.
G.f.: Sum_{n >= 1} phi(n)*(x^n + x^(2*n))/(1 - x^n)^3 = x + 5*x^2 + 11*x^3 + 22*x^4  + .... - Peter Bala, Dec 30 2013
Multiplicative with a(p^e) = p^(e-1)*(p^e*(p+1)-1). - R. J. Mathar, Jun 23 2018
Sum_{k=1..n} a(k) ~ Pi^2 * n^3 / (18*zeta(3)). - Vaclav Kotesovec, Sep 18 2020
a(n) = Sum_{k=1..n} (n/gcd(n,k))^2*phi(gcd(n,k))/phi(n/gcd(n,k)). - Richard L. Ollerton, May 07 2021
From Peter Bala, Dec 26 2023: (Start)
For n odd, a(n) = Sum_{k = 1..n} gcd(k,n)/cos(k*Pi/n)^2 (see Osipov and also Chu, p. 51).
It appears that for n odd, Sum_{k = 1..n} (-1)^(k+1)*gcd(k,n)/cos(k*Pi/n)^2 = n. (End)
a(n) = Sum_{1 <= i, j <= n} gcd(i, j, n). Cf. A360428. - Peter Bala, Jan 16 2024
Sum_{k=1..n} a(k)/k ~ Pi^2 * n^2 / (12*zeta(3)). - Vaclav Kotesovec, May 11 2024

A379499 Square array A(n, k) = A064987(A246278(n, k)), read by falling antidiagonals; A064987(n) = n*sigma(n), applied to the prime shift array.

Original entry on oeis.org

6, 28, 12, 72, 117, 30, 120, 360, 775, 56, 180, 1080, 1680, 2793, 132, 336, 672, 19500, 7392, 16093, 182, 336, 3510, 3960, 137200, 24024, 30927, 306, 496, 1584, 43400, 10192, 1948584, 55692, 88723, 380, 702, 9801, 5460, 368676, 40392, 5228860, 116280, 137541, 552, 840, 9300, 488125, 17136, 2928926, 69160, 25645860, 209760, 292537, 870

Offset: 1

Antti Karttunen, Jan 02 2025

Each column is strictly monotonic.

			The top left corner of the array:
k=|   1      2      3        4      5        6      7          8        9       10
2k|   2      4      6        8     10       12     14         16       18       20
--+---------------------------------------------------------------------------------
1 |   6,    28,    72,     120,   180,     336,   336,       496,     702,     840,
2 |  12,   117,   360,    1080,   672,    3510,  1584,      9801,    9300,    6552,
3 |  30,   775,  1680,   19500,  3960,   43400,  5460,    488125,   83790,  102300,
4 |  56,  2793,  7392,  137200, 10192,  368676, 17136,   6725201,  901208,  508326,
5 | 132, 16093, 24024, 1948584, 40392, 2928926, 50160, 235793305, 4082364, 4924458,

Elementwise product of arrays A246278 and A355927.

Cf. A000203, A003961, A064987, A379500.

up_to = 55;
A064987(n) = (n*sigma(n));
A246278sq(row,col) = if(1==row,2*col, my(f = factor(2*col)); for(i=1, #f~, f[i,1] = prime(primepi(f[i,1])+(row-1))); factorback(f));
A379499sq(row,col) = A064987(A246278sq(row,col));
A379499list(up_to) = { my(v = vector(up_to), i=0); for(a=1,oo, for(col=1,a, i++; if(i > up_to, return(v)); v[i] = A379499sq(col,(a-(col-1))))); (v); };
v379499 = A379499list(up_to);
A379499(n) = v379499[n];

A(n, k) = A246278(n, k) * A355927(n, k).

A038040 a(n) = n*d(n), where d(n) = number of divisors of n (A000005).

Original entry on oeis.org

1, 4, 6, 12, 10, 24, 14, 32, 27, 40, 22, 72, 26, 56, 60, 80, 34, 108, 38, 120, 84, 88, 46, 192, 75, 104, 108, 168, 58, 240, 62, 192, 132, 136, 140, 324, 74, 152, 156, 320, 82, 336, 86, 264, 270, 184, 94, 480, 147, 300, 204, 312, 106, 432, 220, 448, 228, 232, 118

Offset: 1

Christian G. Bower

Dirichlet convolution of sigma(n) (A000203) with phi(n) (A000010). - Michael Somos, Jun 08 2000
Dirichlet convolution of f(n)=n with itself. See the Apostol reference for Dirichlet convolutions. - Wolfdieter Lang, Sep 09 2008
Sum of all parts of all partitions of n into equal parts. - Omar E. Pol, Jan 18 2013

			For n = 6 the partitions of 6 into equal parts are [6], [3, 3], [2, 2, 2], [1, 1, 1, 1, 1, 1]. The sum of all parts is 6 + 3 + 3 + 2 + 2 + 2 + 1 + 1 + 1 + 1 + 1 + 1 = 24 equalling 6 times the number of divisors of 6, so a(6) = 24. - _Omar E. Pol_, May 08 2021

Tom M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, pp. 29 ff.
James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 162.

T. D. Noe, Table of n, a(n) for n = 1..1000
J. Bourgain, S. V. Konyagin and I. E. Shparlinski, Product sets of rationals, multiplicative translates of subgroups in residue rings and fixed points of the discrete logarithms, Int. Math. Res. Notices, 2008 (2008), Art. ID rnn 090, 1-29.
Jean Bourgain, Sergei Konyagin and Igor Shparlinski. Distribution on elements of cosets of small subgroups and applications, arXiv:1103.0567 [math.NT], Mar 2 2011.
Mikhail R. Gabdullin and Vitalii V. Iudelevich, Numbers of the form kf(k), arXiv:2201.09287 [math.NT] (2022).
Passawan Noppakaew and Prapanpong Pongsriiam, Product of Some Polynomials and Arithmetic Functions, J. Int. Seq. (2023) Vol. 26, Art. 23.9.1.
Paul Pollack, Analytic and Combinatorial Number Theory Course Notes, p. 147. [Broken link?]
Paul Pollack, Analytic and Combinatorial Number Theory Course Notes, p. 147.

Cf. A000005, A000010, A000203, A001620, A018804, A029935, A064987, A062952.
Cf. A127648, A127093, A127528.
Cf. A038044, A143127 (partial sums), A328722 (Dirichlet inverse).
Column 1 of A329323.

a038040 n = a000005 n * n  -- Reinhard Zumkeller, Jan 21 2014

with(numtheory): A038040 := n->tau(n)*n;

a[n_] := DivisorSigma[0, n]*n; Table[a[n], {n, 1, 60}] (* Jean-François Alcover, Sep 03 2012 *)

n*numlib::tau (n)$ n=1..90 // Zerinvary Lajos, May 13 2008

a(n)=if(n<1,0,direuler(p=2,n,1/(1-p*X)^2)[n])

a(n)=if(n<1,0,polcoeff(sum(k=1,n,k*x^k/(x^k-1)^2,x*O(x^n)),n)) /* Michael Somos, Jan 29 2005 */

a(n) = n*numdiv(n); \\ Michel Marcus, Oct 24 2020

from sympy import divisor_count as d
def a(n): return n*d(n)
print([a(n) for n in range(1, 60)]) # Michael S. Branicky, Mar 15 2022

[n*sigma(n,0) for n in range(1, 60)] # Stefano Spezia, Jul 20 2025

Dirichlet g.f.: zeta(s-1)^2.
G.f.: Sum_{n>=1} n*x^n/(1-x^n)^2. - Vladeta Jovovic, Dec 30 2001
Sum_{k=1..n} sigma(gcd(n, k)). Multiplicative with a(p^e) = (e+1)*p^e. - Vladeta Jovovic, Oct 30 2001
Equals A127648 * A127093 * the harmonic series, [1/1, 1/2, 1/3, ...]. - Gary W. Adamson, May 10 2007
Equals row sums of triangle A127528. - Gary W. Adamson, May 21 2007
a(n) = n*A000005(n) = A066186(n) - n*(A000041(n) - A000005(n)) = A066186(n) - n*A144300(n). - Omar E. Pol, Jan 18 2013
a(n) = A000203(n) * A240471(n) + A106315(n). - Reinhard Zumkeller, Apr 06 2014
L.g.f.: Sum_{k>=1} x^k/(1 - x^k) = Sum_{n>=1} a(n)*x^n/n. - Ilya Gutkovskiy, May 13 2017
a(n) = Sum_{d|n} A018804(d). - Amiram Eldar, Jun 23 2020
a(n) = Sum_{d|n} phi(d)*sigma(n/d). - Ridouane Oudra, Jan 21 2021
G.f.: Sum_{n >= 1} q^(n^2)*(n^2 + 2*n*q^n - n^2*q^(2*n))/(1 - q^n)^2. - Peter Bala, Jan 22 2021
a(n) = Sum_{k=1..n} sigma(n/gcd(n,k))*phi(gcd(n,k))/phi(n/gcd(n,k)). - Richard L. Ollerton, May 07 2021
Define f(x) = #{n <= x: a(n) <= x}. Gabdullin & Iudelevich show that f(x) ~ x/sqrt(log x). That is, there are 0 < A < B such that Ax/sqrt(log x) < f(x) < Bx/sqrt(log x). - Charles R Greathouse IV, Mar 15 2022
Sum_{k=1..n} a(k) ~ n^2*log(n)/2 + (gamma - 1/4)*n^2, where gamma is Euler's constant (A001620). - Amiram Eldar, Oct 25 2022
Mobius transform of A060640. - R. J. Mathar, Feb 07 2023

A001001 Number of sublattices of index n in generic 3-dimensional lattice.

Original entry on oeis.org

1, 7, 13, 35, 31, 91, 57, 155, 130, 217, 133, 455, 183, 399, 403, 651, 307, 910, 381, 1085, 741, 931, 553, 2015, 806, 1281, 1210, 1995, 871, 2821, 993, 2667, 1729, 2149, 1767, 4550, 1407, 2667, 2379, 4805, 1723, 5187, 1893, 4655, 4030, 3871, 2257, 8463, 2850, 5642, 3991, 6405, 2863

Offset: 1

N. J. A. Sloane

These sublattices are in 1-1 correspondence with matrices
[a b d]
[0 c e]
[0 0 f]
with acf = n, b = 0..c-1, d = 0..f-1, e = 0..f-1. The sublattice is primitive if gcd(a,b,c,d,e,f) = 1.
Total area of all distinct rectangles whose side lengths are divisors of n, and whose length is an integer multiple of the width. - Wesley Ivan Hurt, Aug 23 2020

Richard P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Problem 5.13(d), pp. 76 and 113.

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from T. D. Noe)
Michael Baake, Solution of the coincidence problem in dimensions d <= 4, in R. V. Moody, ed., Math. of Long-Range Aperiodic Order, Kluwer 1997, pp. 9-44.; arXiv:math/0605222 [math.MG], 2006.
M. Baake and N. Neumarker, A Note on the Relation Between Fixed Point and Orbit Count Sequences, JIS 12 (2009) 09.4.4, Section 3.
J. Liouville, Théorème concernant les sommes de diviseurs des nombres, Journal de mathématiques pures et appliquées 2e série, tome 2 (1857), p. 56.
V. A. Liskovets and A. Mednykh, Enumeration of subgroups in the fundamental groups of orientable circle bundles over surfaces, Commun. in Algebra, 28, No. 4 (2000), 1717-1738.
Y. Puri and T. Ward, Arithmetic and growth of periodic orbits, J. Integer Seqs., Vol. 4 (2001), #01.2.1.
J. S. Rutherford, The enumeration and symmetry-significant properties of derivative lattices, Act. Cryst. (1992) A48, 500-508
J. S. Rutherford, The enumeration and symmetry-significant properties of derivative lattices II, Acta Cryst. A49 (1993), 293-300. [_N. J. A. Sloane_, Mar 14 2009]
Tad White, Counting Free Abelian Actions, arXiv:1304.2830 [math.CO], 2013.
Index entries for sequences related to sublattices.

Column 3 of A160870.
Cf. A060983, A064987 (Mobius transform).
Cf. A061256, A226313, A301777.
Primes in this sequence are in A053183.
Cf. A038991, A038992, A038993, A038994, A038995, A038996, A038997, A038998, A038999.

a001001 n = sum [sum [k * (if k `mod` l == 0 then l else 0) | k <- [1..n], n `mod` k == 0 ] | l <- [1..n]]
a = [ a001001 n | n <- [1..53]]
putStrLn $ concat $ map (++ ", ") (map show a) -- Miles Wilson, Apr 04 2025

nmax := 100:
L12 := [seq(1,i=1..nmax) ];
L27 := [seq(i,i=1..nmax) ];
L290 := [seq(i^2,i=1..nmax) ];
DIRICHLET(L12,L27) ;
DIRICHLET(%,L290) ; # R. J. Mathar, Sep 25 2017

a[n_] := Sum[ d*DivisorSigma[1, d], {d, Divisors[n]}]; Table[ a[n], {n, 1, 42}] (* Jean-François Alcover, Jan 20 2012, after Vladeta Jovovic *)
f[p_, e_] := Product[(p^(e + k) - 1)/(p^k - 1), {k, 1, 2}]; a[1] = 1; a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* Amiram Eldar, Aug 29 2019 *)

N=17; default(seriesprecision,N); x=z+O(z^(N+1))
c=sum(j=1,N,j*x^j);
t=1/prod(j=1,N, eta(x^(j))^j)
t=log(t)
t=serconvol(t,c)
Vec(t)
/* Joerg Arndt, May 03 2008 */

a(n)=sumdiv(n,d, d * sumdiv(d,t, t ) );  /* Joerg Arndt, Oct 07 2012 */

a(n)=sumdivmult(n,d, sigma(d)*d) \\ Charles R Greathouse IV, Sep 09 2014

If n = Product p^m, a(n) = Product (p^(m + 1) - 1)(p^(m + 2) - 1)/(p - 1)(p^2 - 1). Or, a(n) = Sum_{d|n} sigma(n/d)*d^2, Dirichlet convolution of A000290 and A000203.
a(n) = Sum_{d|n} d*sigma(d). - Vladeta Jovovic, Apr 06 2001
Multiplicative with a(p^e) = ((p^(e+1)-1)(p^(e+2)-1))/((p-1)(p^2-1)). - David W. Wilson, Sep 01 2001
Dirichlet g.f.: zeta(s)*zeta(s-1)*zeta(s-2).
L.g.f.: -log(Product_{k>=1} (1 - x^k)^sigma(k)) = Sum_{n>=1} a(n)*x^n/n. - Ilya Gutkovskiy, May 23 2018
a(n) = Sum_{d1|n, d2|n, d1|d2} d1*d2. - Wesley Ivan Hurt, Aug 23 2020
Sum_{k=1..n} a(k) ~ c * n^3, where c = Pi^2*zeta(3)/18 = 0.659101... . - Amiram Eldar, Oct 19 2022
G.f.: Sum_{k>=1} Sum {l>=1} k*l^2*x^(k*l - 1)/(1 - x^(k*l)). - Miles Wilson, Apr 04 2025

A349169 Numbers k such that k * gcd(sigma(k), A003961(k)) is equal to the odd part of {sigma(k) * gcd(k, A003961(k))}, where A003961 shifts the prime factorization one step towards larger primes, and sigma is the sum of divisors function.

Original entry on oeis.org

1, 15, 105, 3003, 3465, 13923, 45045, 264537, 459459, 745875, 1541475, 5221125, 8729721, 10790325, 14171625, 29288025, 34563375, 57034575, 71430975, 99201375, 109643625, 144729585, 205016175, 255835125, 295708875, 356080725, 399242025, 419159475, 449323875, 928602675, 939495375, 1083656925, 1941623775, 1962350685, 2083228875

Offset: 1

Antti Karttunen, Nov 10 2021

Numbers k such that A348990(k) [= k/gcd(k, A003961(k))] is equal to A348992(k), which is the odd part of A349162(k), thus all terms must be odd, as A348990 preserves the parity of its argument.
Equally, numbers k for which gcd(A064987(k), A191002(k)) is equal to A000265(gcd(A064987(k), A341529(k))).
Also odd numbers k for which A348993(k) = A319627(k).
Odd terms of A336702 are given by the intersection of this sequence and A349174.
Conjectures:
  (1) After 1, all terms are multiples of 3. (Why?)
  (2) After 1, all terms are in A104210, in other words, for all n > 1, gcd(a(n), A003961(a(n))) > 1. Note that if we encountered a term k with gcd(k, A003961(k)) = 1, then we would have discovered an odd multiperfect number.
  (3) Apart from 1, 15, 105, 3003, 13923, 264537, all other terms are abundant. [These apparently are also the only terms that are not Zumkeller, A083207. Note added Dec 05 2024]
  (4) After 1, all terms are in A248150. (Cf. also A386430).
  (5) After 1, all terms are in A348748.
  (6) Apart from 1, there are no common terms with A349753.
Note: If any of the last four conjectures could be proved, it would refute the existence of odd perfect numbers at once. Note that it seems that gcd(sigma(k), A003961(k)) < k, for all k except these four: 1, 2, 20, 160.
Questions:
  (1) For any term x here, can 2*x be in A349745? (Partial answer: at least x should be in A191218 and should not be a multiple of 3). Would this then imply that x is an odd perfect number? (Which could explain the points (1) and (4) in above, assuming the nonexistence of opn's).

Cf. A000203, A003961, A007949, A064989, A083207, A104210, A161942, A191002, A191218, A228058, A319627, A322361, A326134, A329963, A342671, A348748, A348990, A348992, A348993, A349162, A349164, A349174.
Cf. also A349745, A349746, A349753.
Subsequence of A349749, and of A386430.

Select[Range[10^6], #1/GCD[#1, #3] == #2/(2^IntegerExponent[#2, 2]*GCD[#2, #3]) & @@ {#, DivisorSigma[1, #], Times @@ Map[NextPrime[#1]^#2 & @@ # &, FactorInteger[#]]} &] (* Michael De Vlieger, Nov 11 2021 *)

A000265(n) = (n >> valuation(n, 2));
A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
isA349169(n) = { my(s=sigma(n),u=A003961(n)); (n*gcd(s,u) == A000265(s)*gcd(n,u)); }; \\ (Program simplified Nov 30 2021)

For all n >= 1, A007949(A000203(a(n))) = A007949(a(n)). [sigma preserves the 3-adic valuation of the terms of this sequence] - Antti Karttunen, Nov 29 2021

Name changed and comment section rewritten by Antti Karttunen, Nov 29 2021

A143128 a(n) = Sum_{k=1..n} k*sigma(k).

Original entry on oeis.org

1, 7, 19, 47, 77, 149, 205, 325, 442, 622, 754, 1090, 1272, 1608, 1968, 2464, 2770, 3472, 3852, 4692, 5364, 6156, 6708, 8148, 8923, 10015, 11095, 12663, 13533, 15693, 16685, 18701, 20285, 22121, 23801, 27077, 28483, 30763, 32947, 36547

Offset: 1

Gary W. Adamson, Jul 26 2008

nonn

Partial sums of A064987. - Omar E. Pol, Jul 04 2014
a(n) is also the volume after n-th step of the symmetric staircase described in A244580 (see also A237593). - Omar E. Pol, Jul 31 2018
In general, for j >= 1 and m >= 0, Sum_{k=1..n} k^m * sigma_j(k) ~ n^(j+m+1) * zeta(j+1) / (j+m+1). - Daniel Suteu, Nov 21 2018

			a(4) = 47 = (1 + 6 + 12 + 28) where A064987 = (1, 6, 12, 28, 30, ...).
a(4) = 47 = sum of row 4 terms of triangle A110662 = (15 + 14 + 11 + 7).

Indranil Ghosh, Table of n, a(n) for n = 1..4267

Cf. A000203, A064987, A110662, A175254, A237593, A244580, A256533.

[(&+[k*DivisorSigma(1,k): k in [1..n]]): n in [1..50]]; // G. C. Greubel, Nov 21 2018

with(numtheory): a:=proc(n) options operator, arrow: sum(k*sigma(k), k=1..n) end proc: seq(a(n),n=1..40); # Emeric Deutsch, Aug 12 2008

Table[Sum[i*DivisorSigma[1, i], {i, n}], {n, 50}] (* Wesley Ivan Hurt, Jul 06 2014 *)

a(n)=sum(k=1,n,k*sigma(k)) \\ Charles R Greathouse IV, Apr 27 2015

f(n) = n*(n+1)*(2*n+1)/6; \\ A000330
g(n) = n*(n+1)/2; \\ A000217
a(n) = sum(k=1, sqrtint(n), k * f(n\k) + k^2 * g(n\k)) - f(sqrtint(n)) * g(sqrtint(n)); \\ Daniel Suteu, Nov 26 2020

def A143128(n): return sum(k**2*(m:=n//k)*(m+1)>>1 for k in range(1,n+1)) # Chai Wah Wu, Oct 20 2023

from math import isqrt
def A143128(n): return ((-((s:=isqrt(n))*(s+1))**2*(2*s+1)>>1) + sum((q:=n//k)*(q+1)*k*(3*k+2*q+1) for k in range(1,s+1)))//6 # Chai Wah Wu, Oct 21 2023

[sum(k*sigma(k,1) for k in (1..n)) for n in (1..50)] # G. C. Greubel, Nov 21 2018

Sum {k=1..n} k*sigma(k), where sigma(n) = A000203: (1, 3, 4, 7, 6, 12, ...) and n*sigma(n) = A064987: (1, 6, 12, 28, ...). Equals row sums of triangle A110662. - Emeric Deutsch, Aug 12 2008
a(n) ~ n^3 * Pi^2/18. - Charles R Greathouse IV, Jun 19 2012
G.f.: x*f'(x)/(1 - x), where f(x) = Sum_{k>=1} k*x^k/(1 - x^k). - Ilya Gutkovskiy, Apr 13 2017
a(n) = Sum_{k=1..n} k^2/2 * floor(n/k) * floor(1 + n/k). - Daniel Suteu, May 29 2018
a(n) = A256533(n) - A175254(n-1), n >= 2. - Omar E. Pol, Jul 31 2018
a(n) = Sum_{k=1..s} (k*A000330(floor(n/k)) + k^2*A000217(floor(n/k))) - A000330(s)*A000217(s), where s = floor(sqrt(n)). - Daniel Suteu, Nov 26 2020
a(n) = Sum_{k=1..n} Sum_{i=1..floor(n/k)} i*k^2. - Wesley Ivan Hurt, Nov 26 2020

Corrected and extended by Emeric Deutsch, Aug 12 2008

A281372 Coefficients in q-expansion of (E_2*E_4 - E_6)/720, where E_2, E_4, E_6 are the Eisenstein series shown in A006352, A004009, A013973, respectively.

Original entry on oeis.org

0, 1, 18, 84, 292, 630, 1512, 2408, 4680, 6813, 11340, 14652, 24528, 28574, 43344, 52920, 74896, 83538, 122634, 130340, 183960, 202272, 263736, 279864, 393120, 393775, 514332, 551880, 703136, 707310, 952560, 923552, 1198368, 1230768, 1503684, 1517040, 1989396, 1874198, 2346120, 2400216, 2948400

Offset: 0

N. J. A. Sloane, Feb 04 2017

The q-expansion of the square of this expression is given in A281371.

Multiplicative because A001158 is. - Andrew Howroyd, Jul 23 2018

Seiichi Manyama, Table of n, a(n) for n = 0..1000

Cf. A001158, A003840, A006352, A004009, A013973, A064987, A145094, A281371, A328259.

[0] cat [n*DivisorSigma(3, n): n in [1..50]]; // Vincenzo Librandi, Mar 01 2018

with(gfun):
with(numtheory); M:=100;
E := proc(k) local n, t1; global M;
t1 := 1-(2*k/bernoulli(k))*add(sigma[k-1](n)*q^n, n=1..M+1);
series(t1, q, M+1); end;
e2:=E(2); e4:=E(4); e6:=E(6);
t1:=series((e2*e4-e6)/720,q,M+1);
seriestolist(t1);
# alternative program
seq(add(sigma[4](d)*phi(n/d), d in divisors(n)), n = 1..100); # Peter Bala, Jan 20 2024

Table[If[n==0, 0, n * DivisorSigma[3, n]], {n, 0, 40}] (* Indranil Ghosh, Mar 11 2017 *)
terms = 41; Ei[n_] = 1-(2n/BernoulliB[n]) Sum[k^(n-1) x^k/(1-x^k), {k, terms}]; CoefficientList[(Ei[2] Ei[4] - Ei[6])/720 + O[x]^terms, x] (* Jean-François Alcover, Mar 01 2018 *)

for(n=0, 40, print1(if(n==0, 0, n * sigma(n, 3)), ", ")) \\ Indranil Ghosh, Mar 11 2017

a(n) = A145094(n)/240 for n > 0. - Seiichi Manyama, Feb 04 2017
G.f.: phi_{4, 1}(x) where phi_{r, s}(x) = Sum_{n, m>0} m^r * n^s * x^{m*n}. - Seiichi Manyama, Feb 04 2017
a(n) = n*A001158(n) for n > 0. - Seiichi Manyama, Feb 18 2017
G.f.: x*f'(x), where f(x) = Sum_{k>=1} k^3*x^k/(1 - x^k). - Ilya Gutkovskiy, Aug 31 2017
Sum_{k=1..n} a(k) ~  Pi^4 * n^5 / 450. - Vaclav Kotesovec, May 09 2022
From Amiram Eldar, Oct 30 2023: (Start)
Multiplicative with a(p^e) = p^e * (p^(3*e+3)-1)/(p^3-1).
Dirichlet g.f.: zeta(s-1)*zeta(s-4). (End)
a(n) = Sum_{k = 1..n} sigma_4( gcd(k, n) ) = Sum_{d divides n} sigma_4(d) * phi(n/d). - Peter Bala, Jan 19 2024
a(n) = Sum_{1 <= i, j, k, l <= n} sigma_1( gcd(i, j, k, l, n) ) = Sum_{d divides n} sigma_1(d) * J_4(n/d), where the Jordan totient function J_4(n) = A059377(n). - Peter Bala, Jan 22 2024

A349745 Numbers k for which k * gcd(sigma(k), A003961(k)) is equal to sigma(k) * gcd(k, A003961(k)), where A003961 shifts the prime factorization one step towards larger primes, and sigma is the sum of divisors function.

Original entry on oeis.org

1, 120, 216, 672, 2464, 22176, 228480, 523776, 640640, 837760, 5581440, 5765760, 7539840, 12999168, 19603584, 33860736, 38342304, 71344000, 95472000, 102136320, 197308800, 220093440, 345080736, 459818240, 807009280, 975576960, 1476304896, 1510831360, 1773584640

Offset: 1

Antti Karttunen, Nov 29 2021

nonn

Numbers k for which k * A342671(k) = A000203(k) * A322361(k).
Numbers k such that gcd(A064987(k), A191002(k)) = gcd(A064987(k), A341529(k)).
Obviously, all odd terms in this sequence must be squares.
All the terms k of A005820 that satisfy A007949(k) < A007814(k) [i.e., whose 3-adic valuation is strictly less than their 2-adic valuation] are also terms of this sequence. Incidentally, the first six known terms of A005820 satisfy this condition, while on the other hand, any hypothetical odd 3-perfect number would be excluded from this sequence. Also, as a corollary, any hypothetical 3-perfect numbers of the form 4u+2 must not be multiples of 3 if they are to appear here. Similarly for any k which occurs in A349169, for 2*k to occur in this sequence, it shouldn't be a multiple of 3 and k should also be a term of A191218. See question 2 and its partial answer in A349169.
From Antti Karttunen, Feb 13-20 2022: (Start)
Question: Are all terms/2 (A351548) abundant, from n > 1 onward?
Note that of the 65 known 5-multiperfect numbers, all others except these three 1245087725796543283200, 1940351499647188992000, 4010059765937523916800 are also included in this sequence. The three exceptions are distinguished by the fact that their 3 and 5-adic valuations are equal. In 62 others the former is larger.
If k satisfying the condition were of the form 4u+2, then it should be one of the terms of A191218 doubled as only then both k and sigma(k) are of the form 4u+2, with equal 2-adic valuations for both. More precisely, one of the terms of A351538.
(End)

Amiram Eldar, Table of n, a(n) for n = 1..52 (terms below 10^11)
Antti Karttunen, Ratio A342671(n)/A322361(n) plotted with OEIS Plot2 tool
Index entries for sequences computed from indices in prime factorization
Index entries for sequences related to sigma(n)

Cf. A000203, A003961, A005101, A005820, A007814, A007949, A046060, A064987, A191002, A191218, A248150, A322361, A341529, A342671, A347870, A351534, A351536, A351538.
Cf. also A349169, A349746, A351458, A351549 for other variants.
Subsequence of A351554 and also of its subsequence A351551.
Cf. A351459 (subsequence, intersection with A351458), A351548 (terms halved).

f1[p_, e_] := (p^(e + 1) - 1)/(p - 1); f2[p_, e_] := NextPrime[p]^e; q[1] = True; q[n_] := n * GCD[(s = Times @@ f1 @@@ (f = FactorInteger[n])), (r = Times @@ f2 @@@ f)] == s*GCD[n, r]; Select[Range[10^6], q] (* Amiram Eldar, Nov 29 2021 *)

A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
isA349745(n) = { my(s=sigma(n),u=A003961(n)); (n*gcd(s,u) == (s*gcd(n,u))); };

For all n >= 1, A007814(A000203(a(n))) = A007814(a(n)). [sigma preserves the 2-adic valuation of the terms of this sequence]

A244580 Square array read by antidiagonals related to the symmetric representation of sigma.

Original entry on oeis.org

1, 2, 2, 3, 2, 3, 4, 3, 3, 4, 5, 4, 4, 4, 5, 6, 5, 4, 4, 5, 6, 7, 6, 5, 6, 5, 6, 7, 8, 7, 6, 6, 6, 6, 7, 8, 9, 8, 7, 6, 6, 6, 7, 8, 9, 10, 9, 8, 7, 8, 8, 7, 8, 9, 10, 11, 10, 9, 8, 8, 8, 8, 8, 9, 10, 11, 12, 11, 10, 9, 8, 9, 9, 8, 9, 10, 11, 12

Offset: 1

Omar E. Pol, Jul 04 2014

The number of parts k in the square array is equal to A000203(k) hence the sum of parts k is equal to A064987(k).
The structure has a three-dimensional representation using polycubes. T(n,k) is the height of a column. The total area in the horizontal level z gives A000203(z).
The main diagonal gives A244367.

			.                         _ _ _ _ _ _ _ _ _
1,2,3,4,5,6,7,8,9...     |_| | | | | | | | |
2,2,3,4,5,6,7,8,9...     |_ _|_| | | | | | |
3,3,4,4,5,6,7,8,9...     |_ _|  _|_| | | | |
4,4,4,6,6,6,7,8,9...     |_ _ _|    _|_| | |
5,5,5,6,6,8,8,8,9...     |_ _ _|  _|  _ _|_|
6,6,6,6,8,8,9...         |_ _ _ _|  _| |
7,7,7,7,8,9,9...         |_ _ _ _| |_ _|
8,8,8,8,8...             |_ _ _ _ _|
9,9,9,9,9...             |_ _ _ _ _|
.

Cf. A000203, A064987, A143128, A185283, A196020, A236104, A235791, A237270, A237271, A237591, A237593, A071562, A244367.

A351458 Numbers k for which k * gcd(sigma(k), A276086(k)) is equal to sigma(k) * gcd(k, A276086(k)), where A276086 is the primorial base exp-function, and sigma gives the sum of divisors of its argument.

Original entry on oeis.org

1, 10, 56, 9196, 9504, 56160, 121176, 239096, 354892, 411264, 555520, 716040, 804384, 904704, 1063348, 1387386, 1444352, 1454112, 1884800, 2708640, 3317248, 3548920, 4009824, 4634784, 6179712, 6795360, 7285248, 14511744, 16328466, 28377216, 29855232, 31940280, 37444736, 42711552, 49762944, 52815744

Offset: 1

Antti Karttunen, Feb 13 2022

nonn

Numbers k such that k * A324644(k) = A000203(k) * A324198(k).
Numbers k such that gcd(A064987(k), A324580(k)) = gcd(A064987(k), A351252(k)).
Numbers k such that their abundancy index [sigma(k)/k] is equal to A324644(k)/A324198(k). See A364286.
A324644 gives odd values for even numbers and for the odd squares. A324198 is odd on all arguments, therefore on odd squares the above equation reduces to odd * odd = odd * odd, and on odd nonsquares as odd * even = even * odd. It is an open question whether there are any odd terms after the initial a(1)=1.
If k is even, but not a multiple of 3, then A276086(k) is a multiple of 3, but not even (i.e., is an odd multiple of 3). If for such k also sigma(k) = 3*k, then A007949(A324644(k)) = min(A007949(sigma(k)), A007949(A276086(k))) = 1, while A007949(A324198(k)) = min(A007949(k), A007949(A276086(k))) = 0, therefore all such k's do occur in this sequence, for example, the two known terms of A005820 (3-perfect numbers) that are not multiples of three: 459818240, 51001180160, but also any hypothetical term of A005820 of the form 4u+2, where 2u+1 is not multiple of 3, and which by necessity is then also an odd perfect number.
Similarly, of the 65 known 5-multiperfect numbers (A046060), those 20 that are not multiples of five are included in this sequence. Note that all 65 are multiples of six.
It is conjectured that the intersection of this sequence with the multiperfect numbers (A007691) gives A323653, see comments in the latter.
For all even terms k of this sequence, A007814(A000203(k)) = A007814(k), sigma preserves the 2-adic valuation, and A007949(A000203(k)) >= A007949(k), i.e., does not decrease the 3-adic valuation. The condition is equivalence (=) when k is a multiple of 6. With odd terms, any hypothetical odd perfect number x would yield a one greater 2-adic valuation for sigma(x) than for x, but would satisfy the main condition of this sequence. - Corrected Feb 17 2022
If k is a nonsquare positive odd number (in A088828), then it must be a term of A191218. - Antti Karttunen, Mar 10 2024

Cf. A000203, A005820, A007691, A007814, A007949, A046060, A088828, A191218, A276086, A323653, A324198, A324644, A327938, A364286, A351459 (intersection with a similar sequence, A349745).

Cf. also A351549.

A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
isA351458(n) = { my(s=sigma(n), z=A276086(n)); (n*gcd(s,z))==(s*gcd(n,z)); };

A006530(n) = if(1==n, n, my(f=factor(n)); f[#f~, 1]); \\ Works OK with rationals also!
isA351458(n) = { my(orgn=n, s=sigma(n), abi=s/n, p=2, q=A006530(abi), d, e1, e2); while((1!=abi)&&(p<=q), d = n%p; e1 = min(d, valuation(s, p)); e2 = min(d, valuation(orgn, p)); d = e1-e2; if(valuation(abi,p)!=d, return(0), abi /= (p^d)); n = n\p; p = nextprime(1+p)); (abi==1); }; \\ (This implementation does not require the construction of largish intermediate numbers, A276086, but might still be slower and return a few false positives on the long run, so please check the results with the above program). - Antti Karttunen, Feb 19 2022

cp's OEIS Frontend

A069097 Moebius transform of A064987, n*sigma(n).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A379499 Square array A(n, k) = A064987(A246278(n, k)), read by falling antidiagonals; A064987(n) = n*sigma(n), applied to the prime shift array.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Formula

A038040 a(n) = n*d(n), where d(n) = number of divisors of n (A000005).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

MuPAD

PARI

PARI

PARI

Python

SageMath

Formula

A001001 Number of sublattices of index n in generic 3-dimensional lattice.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

Formula

A349169 Numbers k such that k * gcd(sigma(k), A003961(k)) is equal to the odd part of {sigma(k) * gcd(k, A003961(k))}, where A003961 shifts the prime factorization one step towards larger primes, and sigma is the sum of divisors function.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A143128 a(n) = Sum_{k=1..n} k*sigma(k).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs