%I #18 Jan 05 2024 17:49:26
%S 0,4,6,98,142,2162,3128,47476,68684,1042320,1507930,22883574,33105786,
%T 502396318,726819372,11029835432,15956920408,242153983196,
%U 350325429614,5316357794890,7691202531110,116717717504394,168856130254816
%N a(n) = 22*a(n-2) - a(n-4) + 10, for n > 3, with a(0) = 0, a(1) = 4, a(2) = 6, a(3) = 98.
%C It appears this sequence gives all nonnegative m such that 120*m^2 + 120*m + 1 is a square.
%H G. C. Greubel, <a href="/A105037/b105037.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,22,-22,-1,1).
%F a(n) = 22*a(n-2) - a(n-4) + 10, for n > 3.
%F G.f.: 2*x*(2 + x + 2*x^2)/((1-x)*(1-22*x^2+x^4)). - Maksym Voznyy (voznyy(AT)mail.ru), Aug 11 2009
%t LinearRecurrence[{1,22,-22,-1,1}, {0,4,6,98,142}, 41] (* _G. C. Greubel_, Mar 14 2023 *)
%o (Magma) R<x>:=PowerSeriesRing(Integers(), 40); [0] cat Coefficients(R!( 2*x*(2+x+2*x^2)/((1-x)*(1-22*x^2+x^4)) )); // _G. C. Greubel_, Mar 14 2023
%o (SageMath)
%o @CachedFunction
%o def a(n): # a = A105037
%o if (n<5): return (0,4,6,98,142)[n]
%o else: return a(n-1) +22*a(n-2) -22*a(n-3) -a(n-4) +a(n-5)
%o [a(n) for n in range(41)] # _G. C. Greubel_, Mar 14 2023
%Y Cf. A077421.
%K nonn,easy
%O 0,2
%A _Gerald McGarvey_, Apr 03 2005