%I #32 Apr 26 2023 06:19:05
%S 216,2744,10648,27000,54872,97336,157464,238328,343000,474552,636056,
%T 830584,1061208,1331000,1643032,2000376,2406104,2863288,3375000,
%U 3944312,4574296,5268024,6028568,6859000,7762392,8741816,9800344,10941048,12167000,13481272,14886936
%N a(n) = (8*n + 6)^3.
%C 4*n + 3 = (8*n + 6) / 2 is never a square, as 3 is not a quadratic residue modulo 4. Using this, we can show that each term has an even square part and an even squarefree part, neither part being a power of 2. (Less than 2% of integers have this property - see A339245.) - _Peter Munn_, Dec 14 2020
%H Vincenzo Librandi, <a href="/A017139/b017139.txt">Table of n, a(n) for n = 0..10000</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/QuadraticResidue.html">Quadratic Residue</a>.
%H <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (4,-6,4,-1).
%F From _R. J. Mathar_, Mar 22 2010: (Start)
%F G.f.: 8*(27 + 235*x + 121*x^2 + x^3)/(x-1)^4.
%F a(n) = 8*A016839(n). (End)
%F a(0)=216, a(1)=2744, a(2)=10648, a(3)=27000, a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - _Harvey P. Dale_, Dec 11 2012
%F a(n) = A017137(n)^3 = A000578(A017137(n)). - _Peter Munn_, Dec 20 2020
%F Sum_{n>=0} 1/a(n) = 7*zeta(3)/128 - Pi^2/512. - _Amiram Eldar_, Apr 26 2023
%t Table[(8*n+6)^3,{n,0,5!}] (* _Vladimir Joseph Stephan Orlovsky_, Mar 17 2010 *)
%t LinearRecurrence[{4,-6,4,-1},{216,2744,10648,27000},30] (* _Harvey P. Dale_, Dec 11 2012 *)
%o (Magma) [(8*n+6)^3: n in [0..35]]; // _Vincenzo Librandi_, Jul 22 2011
%Y A000578, A016839, A017137 are used in a formula defining this sequence.
%Y Subsequence of A339245.
%Y Cf. A002117, A017138, A017140.
%K nonn,easy
%O 0,1
%A _N. J. A. Sloane_
%E More terms from _Vladimir Joseph Stephan Orlovsky_, Mar 17 2010
|