%I #35 Sep 08 2022 08:45:38
%S 0,21,78,171,300,465,666,903,1176,1485,1830,2211,2628,3081,3570,4095,
%T 4656,5253,5886,6555,7260,8001,8778,9591,10440,11325,12246,13203,
%U 14196,15225,16290,17391,18528,19701,20910,22155,23436,24753,26106,27495
%N a(n) = 3*n*(6*n + 1).
%H Harvey P. Dale, <a href="/A144314/b144314.txt">Table of n, a(n) for n = 0..1000</a>
%H <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (3,-3,1).
%F a(n) = A000217(6*n) = A014105(3*n) = A081266(2*n).
%F G.f.: 3*x*(7+5*x)/(1-x)^3. - _Vincenzo Librandi_, Dec 18 2014
%F From _Wesley Ivan Hurt_, Dec 16 2015: (Start)
%F a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n>2.
%F a(n) = 3 * A049453(n). (End)
%F E.g.f.: 3*exp(x)*x*(7 + 6*x). - _Stefano Spezia_, Jun 29 2021
%p A144314:=n->3*n*(6*n+1): seq(A144314(n), n=0..70); # _Wesley Ivan Hurt_, Dec 16 2015
%t Table[3n(6n+1),{n,0,40}] (* or *) LinearRecurrence[{3,-3,1},{0,21,78},40] (* _Harvey P. Dale_, Dec 17 2014 *)
%t CoefficientList[Series[x (21 + 15 x) / (1 - x)^3, {x, 0, 50}], x] (* _Vincenzo Librandi_, Dec 18 2014 *)
%o (Magma) [18*n^2+3*n: n in [0..50]]; // _Vincenzo Librandi_, Dec 18 2014
%o (PARI) a(n)=3*n*(6*n+1) \\ _Charles R Greathouse IV_, Oct 07 2015
%Y Cf. A000217, A014105, A033585, A049453, A081266, A144312.
%K nonn,easy
%O 0,2
%A _Reinhard Zumkeller_, Sep 17 2008