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.

A303384 Total area of all rectangles with dimensions s and t where s | t, n = s + t and s <= t.

This page as a plain text file.
%I A303384 #35 Dec 11 2023 08:36:40
%S A303384 0,1,2,7,4,22,6,35,26,50,10,126,12,86,100,155,16,247,18,294,172,182,
%T A303384 22,590,124,242,260,518,28,860,30,651,364,386,380,1365,36,470,484,
%U A303384 1390,40,1532,42,1134,1144,662,46,2542,342,1395,772,1526,52,2380,788
%N A303384 Total area of all rectangles with dimensions s and t where s | t, n = s + t and s <= t.
%H A303384 Vincenzo Librandi, <a href="/A303384/b303384.txt">Table of n, a(n) for n = 1..1000</a>
%H A303384 <a href="/index/Par#part">Index entries for sequences related to partitions</a>
%F A303384 a(n) = Sum_{i=1..floor(n/2)} i * (n-i) * (floor((n-i)/i) - floor((n-i-1)/i)).
%F A303384 a(n) = Sum_{d|n} d*(n-d). - _Daniel Suteu_, Jun 19 2018
%F A303384 a(n) = n*sigma(n) - sigma_2(n). - _Ridouane Oudra_, Apr 15 2021
%F A303384 From _Amiram Eldar_, Dec 11 2023: (Start)
%F A303384 a(n) = A064987(n) - A001157(n).
%F A303384 Sum_{k=1..n} a(k) ~  c * n^3 / 3, where c = zeta(2) - zeta(3) = 0.442877... . (End)
%p A303384 with(numtheory): seq(n*sigma(n) - sigma[2](n), n=1..60); # _Ridouane Oudra_, Apr 15 2021
%t A303384 Table[Sum[i (n - i) (Floor[(n - i)/i] - Floor[(n - i - 1)/i]), {i, Floor[n/2]}], {n, 80}]
%t A303384 a[n_] := n * DivisorSigma[1, n] - DivisorSigma[2, n]; Array[a, 100] (* _Amiram Eldar_, Dec 11 2023 *)
%o A303384 (Magma) [0] cat [&+[k*(n-k)*((n-k) div k)-(n-k-1) div k:  k in [1..n div 2]]: n in [2..80]]; // _Vincenzo Librandi_, Jun 07 2018
%o A303384 (PARI) a(n) = sum(i=1, n\2, i*(n-i)*((n-i)\i - (n-i-1)\i)); \\ _Michel Marcus_, Jun 07 2018
%o A303384 (PARI) a(n) = sumdiv(n, d, d*(n-d)); \\ _Daniel Suteu_, Jun 19 2018
%o A303384 (PARI) a(n) = {my(f = factor(n)); n * sigma(f) - sigma(f, 2);} \\ _Amiram Eldar_, Dec 11 2023
%o A303384 (GAP) List([1..60],n->Sum([1..Int(n/2)],i->i*(n-i)*(Int((n-i)/i)-Int((n-i-1)/i)))); # _Muniru A Asiru_, Jun 07 2018
%Y A303384 Cf. A297024, A303385.
%Y A303384 Cf. A000203, A001157.
%Y A303384 Cf. A002117, A013661.
%K A303384 nonn,easy
%O A303384 1,3
%A A303384 _Wesley Ivan Hurt_, Apr 22 2018