%I #24 Feb 09 2024 11:16:10
%S 1,1,2,3,5,9,15,26,45,77,133,229,394,679,1169,2013,3467,5970,10281,
%T 17705,30489,52505,90418,155707,268141,461761,795191,1369386,2358197,
%U 4061013,6993405,12043229,20739450,35715071,61504345,105915637,182395603,314100514
%N Arises from a construction of equiangular lines in complex space of dimension 2.
%H Vincenzo Librandi, <a href="/A239909/b239909.txt">Table of n, a(n) for n = 1..1000</a>
%H G. McConnell, <a href="http://arxiv.org/abs/1402.7330">Some non-standard ways to generate SIC-POVMs in dimensions 2 and 3</a>, arXiiv preprint arXiv:1402.7330, 2014, p. 4.
%H <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (1,1,1,-1).
%F From _Vincenzo Librandi_ Apr 09 2014: (Start)
%F G.f.: x*(1-x^3)/(x^4-x^3-x^2-x+1).
%F a(n) = a(n-1) + a(n-2) + a(n-3) - a(n-4) for n>4.
%F a(n) = a(n-1) + 2*a(n-3) + A116732(n-5) for n>4. (End)
%t LinearRecurrence[{1, 1, 1, -1}, {1, 1, 2, 3}, 40] (* or *) CoefficientList[Series[(1 - x^3)/(x^4 - x^3 - x^2 - x + 1), {x, 0, 100}], x] (* _Vincenzo Librandi_, Apr 09 2014 *)
%o (Magma) I:=[1,1,2,3]; [n le 4 select I[n] else Self(n-1)+Self(n-2)+Self(n-3)-Self(n-4): n in [1..50]];
%Y Cf. A116732.
%K nonn,easy
%O 1,3
%A _N. J. A. Sloane_, Apr 09 2014