A175113 a(n) = ((2*n + 1)^6 + 1)/2.
1, 365, 7813, 58825, 265721, 885781, 2413405, 5695313, 12068785, 23522941, 42883061, 74017945, 122070313, 193710245, 297411661, 443751841, 645733985, 919132813, 1282863205, 1759371881, 2375052121, 3160681525, 4151882813
Offset: 0
Links
- Vincenzo Librandi, Table of n, a(n) for n = 0..1000
- R. J. Mathar, Point counts of D_k and some A_k and E_k integer lattices inside hypercubes, arXiv:1002.3844, Variable V_6^(g)(n).
- Index entries for linear recurrences with constant coefficients, signature (7,-21,35,-35,21,-7,1).
Programs
-
Magma
I:=[1, 365, 7813, 58825, 265721, 885781, 2413405]; [n le 7 select I[n] else 7*Self(n-1) - 21*Self(n-2) + 35*Self(n-3) - 35*Self(n-4) + 21*Self(n-5) - 7*Self(n-6) + Self(n-7): n in [1..40]]; // Vincenzo Librandi, Dec 20 2012
-
Mathematica
CoefficientList[Series[(1 + 358*x + 5279*x^2 + 11764*x^3 + 5279*x^4 + 358*x^5 + x^6)/(1 - x)^7, {x, 0, 40}], x] (* Vincenzo Librandi, Dec 20 2012 *)
Formula
a(n)= 7*a(n-1) -21*a(n-2) +35*a(n-3) -35*a(n-4) +21*a(n-5) -7*a(n-6) +a(n-7).
G.f.: (1+358*x+5279*x^2+11764*x^3+5279*x^4+358*x^5+x^6)/(1-x)^7.
a(n) = (2*n^2+2*n+1)*(16*n^4+32*n^3+20*n^2+4*n+1). - Bruno Berselli, Dec 27 2010
Comments