login
Somos's sequence {b(9,n)} defined in comment in A078495: a(0)=a(1)=...=a(20)=1; for n>=21, a(n)=(a(n-1)*a(n-20)+a(n-10)*a(n-11))/a(n-21).
1

%I #25 May 04 2016 08:50:51

%S 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,3,4,5,6,7,8,9,10,11,13,

%T 19,31,51,81,123,179,251,341,451,1045,1691,2451,3459,4977,7467,11679,

%U 18755,30349,48763,100474,282777,679512,1391391,2547414,4327101

%N Somos's sequence {b(9,n)} defined in comment in A078495: a(0)=a(1)=...=a(20)=1; for n>=21, a(n)=(a(n-1)*a(n-20)+a(n-10)*a(n-11))/a(n-21).

%H Seiichi Manyama, <a href="/A272038/b272038.txt">Table of n, a(n) for n = 0..578</a>

%o (Ruby)

%o def b(k, n)

%o b = Array.new(2 * k + 3, 1)

%o (2 * k + 3..n).each{|i|

%o j = (b[i - 1] * b[i - 2 * k - 2] + b[i - k - 1] * b[i - k - 2]) / b[i - 2 * k - 3].to_r

%o j = j.to_i if j.denominator == 1

%o b[i] = j

%o }

%o b[0..n]

%o end

%o p b(9, n) # _Seiichi Manyama_, May 04 2016

%Y Cf. A006721, A078495, A268199, A271954, A271948, A271949, A271950, A271952, A271955.

%K nonn

%O 0,22

%A _Vladimir Shevelev_ and _Peter J. C. Moses_, Apr 18 2016