%I #31 Sep 08 2022 08:45:03

%S 16,123,361,778,1428,2371,3673,5406,7648,10483,14001,18298,23476,

%T 29643,36913,45406,55248,66571,79513,94218,110836,129523,150441,

%U 173758,199648,228291,259873,294586,332628,374203,419521,468798,522256,580123,642633,710026,782548,860451

%N Surround numbers of an n X 1 rectangle.

%H Harry J. Smith, <a href="/A060633/b060633.txt">Table of n, a(n) for n = 1..1000</a>

%H E. J. Friedman, <a href="https://erich-friedman.github.io/mathmagic/0599.html">Math. Magic</a>

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

%F a(n) = (n^4 + 22*n^3 + 105*n^2 - 56*n - 8)/4.

%F G.f.: x*(-16 - 43*x + 94*x^2 - 43*x^3 + 2*x^4) / (x-1)^5. - _R. J. Mathar_, Jan 30 2011

%p for n from 1 to 50 do printf(`%d,`,(n^4+22*n^3+105*n^2-56*n-8)/4) od:

%t Table[(n^4 + 22 n^3 + 105 n^2 - 56 n - 8) / 4, {n, 40}] (* _Vincenzo Librandi_, Jul 03 2018 *)

%t LinearRecurrence[{5,-10,10,-5,1},{16,123,361,778,1428},50] (* _Harvey P. Dale_, Dec 24 2019 *)

%o (PARI) a(n)={(n^4 + 22*n^3 + 105*n^2 - 56*n - 8)/4} \\ _Harry J. Smith_, Jul 08 2009

%o (Magma) [(n^4 + 22*n^3 + 105*n^2 - 56*n - 8)/4: n in [1..40]]; // _Vincenzo Librandi_, Jul 03 2018

%Y Cf. A047875.

%K easy,nonn

%O 1,1

%A _Jason Earls_, Apr 15 2001

%E More terms from _James A. Sellers_, Apr 16 2001