|
%I #26 Sep 08 2022 08:45:16
%S 1,1,1,1,1,1,2,3,4,8,17,50,107,239,1103,3775,14463,55283,256666,
%T 2059753,9820288,55075036,503857819,4083736906,44590046729,
%U 335845998321,3581731774609,68868876045617,782035904796497,11680434156713849,194342679446776442
%N a(n) = (a(n-1) * a(n-5) + a(n-3)^2) / a(n-6) with a(0) = ... = a(5) = 1, a(n) = a(5-n) for all n in Z.
%C Sequence defined by recursion derived from a genus 2 curve.
%C Similar to the Somos-6 and Somos-7 sequences with many bilinear identities.
%H Seiichi Manyama, <a href="/A102276/b102276.txt">Table of n, a(n) for n = 0..211</a>
%H A. J. van der Poorten, <a href="http://arXiv.org/abs/math.NT/0412372">Curves of Genus 2, Continued Fractions and Somos Sequences</a>, arXiv:math/0412372 [math.NT], 2004.
%H A. J. van der Poorten, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL8/Poorten2/vdp89.html">Curves of Genus 2, Continued Fractions and Somos Sequences</a>, J. Integer Seqs., 8 (2005), #05.3.4.
%H <a href="/index/Tu#2wis">Index entries for two-way infinite sequences</a>
%F a(n) = A256858(2*n - 5) for all n in Z. - _Michael Somos_, Apr 13 2015
%F Let b(n) = A256916(n). Then 0 = a(n) * b(n) - a(n-2) * b(n+2) + a(n-3) * b(n+3) for all n in Z. - _Michael Somos_, Apr 13 2015
%F 0 = a(n) * a(n+6) - a(n+1) * a(n+5) - a(n+3) * a(n+3) for all n in Z. - _Michael Somos_, Apr 13 2015
%F 0 = a(n) * a(n+9) + a(n+2) * a(n+7) - a(n+3) * a(n+6) - 9 * a(n+4) * a(n+5) for all n in Z. - _Michael Somos_, Apr 13 2015
%t Join[{1, 1, 1, 1, 1}, RecurrenceTable[{a[n] == (a[n-1]*a[n-5] + a[n-3]^2)/a[n-6], a[6] == 1, a[7] == 2, a[8] == 3, a[9] == 4, a[10] == 8, a[11] == 17}, a, {n, 6, 60}]] (* _G. C. Greubel_, Aug 03 2018 *)
%o (PARI) {a(n) = my(an); if( n<0, a(5-n), n++; an = vector(n,i,1); for(k=7, n, an[k] = (an[k-1]*an[k-5] + an[k-3]^2) / an[k-6]); an[n])};
%o (Magma) I:=[1, 2, 3, 4, 8, 17]; [1, 1, 1, 1, 1] cat [n le 6 select I[n] else (Self(n-1)*Self(n-5) + Self(n-3)^2)/Self(n-6): n in [1..30]]; // _G. C. Greubel_, Aug 03 2018
%Y Cf. A006720, A006722, A256858, A256916, A018896, A271341, A271835, A271831, A271837, A271838, A271839.
%K nonn
%O 0,7
%A _Michael Somos_, Jan 02 2005
|
