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.

A307642 a(n) = n!*Sum_{i=1..n} (Sum_{j=1..i} (i/j)).

This page as a plain text file.
%I A307642 #43 Sep 08 2022 08:46:21
%S A307642 1,8,57,428,3510,31644,312984,3380544,39664080,502927200,6858181440,
%T A307642 100135491840,1559197261440,25797280723200,452046655872000,
%U A307642 8364495012249600,162994310248089600,3336683369519001600,71596721810396160000,1606993396943155200000
%N A307642 a(n) = n!*Sum_{i=1..n} (Sum_{j=1..i} (i/j)).
%F A307642 a(n) = n! * Sum_{i=1..n} (Sum_{j=1..i} (i/j)).
%F A307642 a(n) = n * A182541(n+2).
%F A307642 a(n) = (1/4) * n * (n+1)! * (2*harmonic(n+1) - 1).
%e A307642 a(3) = 57 because a(3) = 3!*Sum_{i=1..3} (Sum_{j=1..i} (i/j)).
%t A307642 Array[#!*Sum[Sum[i/j, {j, i}], {i, #}] &, 25] (* _Michael De Vlieger_, Apr 21 2019 *)
%t A307642 Table[n*(n+1)!*(2*HarmonicNumber[n+1] -1)/4, {n, 25}] (* _G. C. Greubel_, Jul 15 2019 *)
%o A307642 (PARI) a(n)=n!*sum(i=1, n, sum(j=1, i, i/j)); \\ _Michel Marcus_, Apr 20 2019
%o A307642 (Magma) [n*Factorial(n+1)*(2*HarmonicNumber(n+1)-1)/4: n in [1..25]]; // _G. C. Greubel_, Jul 15 2019
%o A307642 (Sage) [n*factorial(n+1)*(2*harmonic_number(n+1)-1)/4 for n in (1..25)] # _G. C. Greubel_, Jul 15 2019
%o A307642 (GAP) List([1..25], n-> n*Factorial(n+1)*(1+2*Sum([2..n+1], j-> 1/j))/4 ); # _G. C. Greubel_, Jul 15 2019
%Y A307642 Cf. A001008/A002805 (harmonic), A182541.
%K A307642 nonn
%O A307642 1,2
%A A307642 _Pedro Caceres_, Apr 19 2019