cp's OEIS Frontend

A046162 Reduced numerators of (n-1)^2/(n^2 + n + 1).

%I A046162 #43 Feb 16 2025 08:32:39
%S A046162 0,1,4,3,16,25,12,49,64,27,100,121,48,169,196,75,256,289,108,361,400,
%T A046162 147,484,529,192,625,676,243,784,841,300,961,1024,363,1156,1225,432,
%U A046162 1369,1444,507,1600,1681,588,1849,1936,675,2116,2209,768,2401,2500
%N A046162 Reduced numerators of (n-1)^2/(n^2 + n + 1).
%C A046162 Arises in Routh's theorem.
%C A046162 With offset 0, multiplicative with a(3^e) = 3^(2e-1), a(p^e) = p^(2e) otherwise. - _David W. Wilson_, Jun 12 2005, corrected by _Robert Israel_, Apr 28 2017
%H A046162 Robert Israel, <a href="/A046162/b046162.txt">Table of n, a(n) for n = 1..10000</a>
%H A046162 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/RouthsTheorem.html">Routh's Theorem</a>.
%F A046162 G.f.: x^2*(1 + 4*x + 3*x^2 + 13*x^3 + 13*x^4 + 3*x^5 + 4*x^6 + x^7)/(1 - x^3)^3.
%F A046162 a(n) = (n-1)^2/3 if n-1 == 0 (mod 3), (n-1)^2 otherwise. - _David W. Wilson_, Jun 12 2005, corrected by _Robert Israel_, Apr 28 2017
%F A046162 From _Amiram Eldar_, Aug 11 2022: (Start)
%F A046162 a(n) = numerator((n-1)^2/3).
%F A046162 Sum_{n>=2} 1/a(n) = 11*Pi^2/54. (End)
%F A046162 From _Amiram Eldar_, Dec 30 2022: (Start)
%F A046162 With offset 0, Dirichlet g.f.: zeta(s-2)*(1-6/3^s).
%F A046162 Sum_{k=1..n} a(k) ~ 7*n^3/27. (End)
%p A046162 seq(numer((n-1)^2/(n^2+n+1)), n=1..51) ; # _Zerinvary Lajos_, Jun 04 2008
%p A046162 seq(denom(3/n^2-2), n=0..76) ; # _Zerinvary Lajos_, Jun 04 2008
%t A046162 a[n_] := Numerator[(n - 1)^2/(n^2 + n + 1)]; Array[a, 50] (* _Amiram Eldar_, Aug 11 2022 *)
%o A046162 (Magma) [Numerator((n-1)^2/3): n in [1..70]]; // _G. C. Greubel_, Oct 27 2022
%o A046162 (SageMath) [numerator((n-1)^2/3) for n in range(1,71)] # _G. C. Greubel_, Oct 27 2022
%Y A046162 Cf. A046163 (denominators), A147560.
%K A046162 nonn,mult
%O A046162 1,3
%A A046162 _Eric W. Weisstein_