cp's OEIS Frontend

A333260 Number of terms in polynomial sequence s(n) = (x*s(n-1)*s(n-4) + y*s(n-2)*s(n-3))/s(n-5), with s(k) = 1 for k = 0..4.

%I A333260 #7 Mar 14 2020 03:14:24
%S A333260 1,1,1,1,1,2,3,4,7,11,16,23,33,46,64,84,109,143,184,228,283,351,429,
%T A333260 515,615,734,871,1017,1181,1376,1593,1821,2077,2372,2694,3035,3409,
%U A333260 3832,4294,4777,5299,5888,6522,7180,7891,8681,9523,10400,11337,12367,13465
%N A333260 Number of terms in polynomial sequence s(n) = (x*s(n-1)*s(n-4) + y*s(n-2)*s(n-3))/s(n-5), with s(k) = 1 for k = 0..4.
%C A333260 s(n) is a generalized Somos-5 sequence (A006721) having coefficients x, y in the recurrence numerator sum of products.
%F A333260 a(n) = a(4-n) for all n in Z.
%e A333260 a(7) = 4 because s(7) = x^3 + x^2*y + 2*x*y*z + y^2*z has 4 terms.
%t A333260 a[ n_] := If[0 <= n <= 4, 1, RecurrenceTable[{s[m]*s[m - 5] == x*s[m - 1]*s[m - 4] + y*s[m - 2]*s[m - 3], s[0] == s[1] == s[2] == s[3] == s[4] == 1}, s, {m, Max[n, 4 - n]}] // Last // Factor // Expand // Length];
%Y A333260 Cf. A006721.
%K A333260 nonn
%O A333260 0,6
%A A333260 _Michael Somos_, Mar 13 2020