|
%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
|
