cp's OEIS Frontend

A158057 First differences of A051870: 16*n + 1.

%I A158057 #41 Apr 13 2025 08:19:48
%S A158057 1,17,33,49,65,81,97,113,129,145,161,177,193,209,225,241,257,273,289,
%T A158057 305,321,337,353,369,385,401,417,433,449,465,481,497,513,529,545,561,
%U A158057 577,593,609,625,641,657,673,689,705,721,737,753,769,785,801,817,833,849
%N A158057 First differences of A051870: 16*n + 1.
%C A158057 The identity (16*n+1)^2 - (16*n^2+2*n)*(4)^2 = 1 can be written as a(n+1)^2 - A158056(n)*(4)^2 = 1. - _Vincenzo Librandi_, Feb 09 2012
%C A158057 This sequence gives the 18-gonal (or octadecagonal) gnomonic numbers. Name suggested by _Todd Silvestri_, Nov 22 2014
%C A158057 All elements are odd and contains subsequence A249356. - _Todd Silvestri_, Nov 22 2014
%H A158057 Vincenzo Librandi, <a href="/A158057/b158057.txt">Table of n, a(n) for n = 0..1000</a>
%H A158057 OEIS Wiki, <a href="http://oeis.org/wiki/Gnomonic_numbers">Gnomonic numbers</a>.
%H A158057 <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (2,-1).
%F A158057 a(n) = 16*n + 1.
%F A158057 a(n) = 2*a(n-1) - a(n-2), a(0) = 1, a(1) = 17.
%F A158057 G.f.: (1+15*x)/(1-x)^2. - _Vincenzo Librandi_, Nov 23 2014
%F A158057 E.g.f.: (1 + 16*x)*exp(x). - _G. C. Greubel_, Sep 18 2019 [corrected by _Elmo R. Oliveira_, Apr 12 2025]
%F A158057 a(n) = A017077(2*n) = A016813(4*n). - _Elmo R. Oliveira_, Apr 12 2025
%p A158057 seq(16*n+1, n=0..60); # _G. C. Greubel_, Sep 18 2019
%t A158057 LinearRecurrence[{2,-1}, {1,17}, 60]
%t A158057 Table[16*n+1,{n,0,60}] (* _Vladimir Joseph Stephan Orlovsky_, Mar 10 2010 *)
%t A158057 a[n_Integer/;n>=0]:=16 n+1 (* _Todd Silvestri_, Nov 22 2014 *)
%t A158057 CoefficientList[Series[(1+15x)/(1-x)^2, {x,0,60}], x] (* _Vincenzo Librandi_, Nov 23 2014 *)
%o A158057 (Magma) I:=[1, 17]; [n le 2 select I[n] else 2*Self(n-1)-Self(n-2): n in [1..50]];
%o A158057 (PARI) a(n)=n<<4+1 \\ _Charles R Greathouse IV_, Dec 23 2011
%o A158057 (Sage) [16*n+1 for n in (0..60)] # _G. C. Greubel_, Sep 18 2019
%o A158057 (GAP) List([0..60], n-> 16*n+1); # _G. C. Greubel_, Sep 18 2019
%Y A158057 Cf. A016813, A017077, A051870, A158056, A249356.
%K A158057 nonn,easy
%O A158057 0,2
%A A158057 _Vincenzo Librandi_, Mar 12 2009
%E A158057 Name clarified and offset changed by _Todd Silvestri_, Nov 22 2014
%E A158057 Edited by _Vincenzo Librandi_ Nov 23 2014
%E A158057 Edited: Offset changed to 0 according to the
%E A158057 Todd Silvestri proposal. Name changed. - _Wolfdieter Lang_, Nov 29 2014