login
Numbers of the form m*(4*m+1)+1, where m = 0,-1,1,-2,2,-3,3,...
5

%I #37 Sep 08 2022 08:46:15

%S 1,4,6,15,19,34,40,61,69,96,106,139,151,190,204,249,265,316,334,391,

%T 411,474,496,565,589,664,690,771,799,886,916,1009,1041,1140,1174,1279,

%U 1315,1426,1464,1581,1621,1744,1786,1915,1959,2094,2140,2281,2329,2476,2526

%N Numbers of the form m*(4*m+1)+1, where m = 0,-1,1,-2,2,-3,3,...

%C Also, numbers m such that 16*m-15 is a square. Therefore, the terms 1 and 4 are the only squares in this sequence.

%H Bruno Berselli, <a href="/A266883/b266883.txt">Table of n, a(n) for n = 0..1000</a>

%H <a href="/index/Rec#order_05">Index entries for linear recurrences with constant coefficients</a>, signature (1,2,-2,-1,1).

%F O.g.f.: (1 + 3*x + 3*x^3 + x^4)/((1 + x)^2*(1 - x)^3).

%F E.g.f.: (5 + 8*x + 4*x^2)*exp(x)/4 -(1 - 2*x)*exp(-x)/4.

%F a(n) = a(-n-1) = n*(n + 1) + 1 - ((2*n + 1)*(-1)^n - 1)/4 = (2*n + 1)*floor((n + 1)/2) + 1.

%F a(n) = A002061(n+1) + A001057(n) = A074378(n)+1.

%F a(n+1) + a(n+2) = A049486(n+3).

%t Table[n (n + 1) + 1 - ((2 n + 1) (-1)^n - 1)/4, {n, 0, 50}]

%t LinearRecurrence[{1, 2, -2, -1, 1}, {1, 4, 6, 15, 19}, 60] (* _Vincenzo Librandi_, Jan 06 2016 *)

%o (PARI) vector(50, n, n--; n*(n+1)+1-((2*n+1)*(-1)^n-1)/4)

%o (PARI) Vec((1+3*x+3*x^3+x^4)/((1+x)^2*(1-x)^3) + O(x^100)) \\ _Altug Alkan_, Jan 06 2016

%o (Sage) [n*(n+1)+1-((2*n+1)*(-1)^n-1)/4 for n in range(50)]

%o (Python) [n*(n+1)+1-((2*n+1)*(-1)**n-1)/4 for n in range(60)]

%o (Magma) [n*(n+1)+1-((2*n+1)*(-1)^n-1)/4: n in [0..50]];

%o (Magma) I:=[1,4,6,15,19]; [n le 5 select I[n] else Self(n-1) + 2*Self(n-2) -2*Self(n-3)-Self(n-4)+Self(n-5): n in [1..60]]; // _Vincenzo Librandi_, Jan 06 2016

%Y Cf. A001057, A049486, A127365, A130472.

%Y Cf. A002061: m*(4*m+2)+1 for m = 0,0,-1,1,-2,2,-3,3, ...

%Y Cf. A174114: m*(4*m+3)+1 for m = 0,-1,1,-2,2,-3,3,-4,4, ...

%Y Cf. A054556: m*(4*m+1)+1 for nonpositive m.

%Y Cf. A054567: m*(4*m+1)+1 for nonnegative m.

%Y Cf. A074378: numbers m such that 16*m+1 is a square.

%K nonn,easy

%O 0,2

%A _Bruno Berselli_, Jan 05 2016