A257293 Numbers n such that T(n) + T(n+1) + ... + T(n+12) is a square, where T = A000217 (triangular numbers).
3, 29, 75, 432, 998, 3624, 8310, 44717, 102443, 370269, 848195, 4561352, 10448838, 37764464, 86508230, 465213837, 1065679683, 3851605709, 8822991915, 47447250672, 108688879478, 392826018504, 899858667750, 4839154355357, 11085200027723, 40064402282349
Offset: 1
Links
- Colin Barker, Table of n, a(n) for n = 1..1000
- Index entries for linear recurrences with constant coefficients, signature (1,0,0,102,-102,0,0,-1,1).
Programs
-
Magma
I:=[3,29,75,432,998,3624,8310,44717,102443]; [n le 9 select I[n] else Self(n-1)+102*Self(n-4)-102*Self(n-5)-Self(n-8)+Self(n-9): n in [1..40]]; // Vincenzo Librandi, May 05 2015
-
Mathematica
Select[Range[10^5],IntegerQ[Sqrt[(#^2+13*#+56)*13/2]]&] (* Ivan N. Ianakiev, May 04 2015 *) LinearRecurrence[{1, 0, 0, 102, -102, 0, 0, -1, 1}, {3, 29, 75, 432, 998, 3624, 8310, 44717, 102443}, 50] (* Vincenzo Librandi, May 05 2015 *)
-
PARI
for(n=0,10^8,issquare(binomial(n+14,3)-binomial(n+1,3))&&print1(n","))
-
PARI
Vec(x*(3*x^8+7*x^7+6*x^6+26*x^5-260*x^4-357*x^3-46*x^2-26*x-3) / ((x-1)*(x^4-10*x^2-1)*(x^4+10*x^2-1)) + O(x^100)) \\ Colin Barker, May 04 2015
Formula
G.f.: x*(3*x^8+7*x^7+6*x^6+26*x^5-260*x^4-357*x^3-46*x^2-26*x-3) / ((x-1)*(x^4-10*x^2-1)*(x^4+10*x^2-1)). - Colin Barker, May 04 2015
Comments