cp's OEIS Frontend

A271952 Somos's sequence {b(6,n)} defined in comment in A078495: a(0)=a(1)=...=a(14)=1; for n>=15, a(n)=(a(n-1)*a(n-14)+a(n-7)*a(n-8))/a(n-15).

%I A271952 #20 Sep 08 2022 08:46:16
%S A271952 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,3,4,5,6,7,8,10,16,28,48,78,120,176,
%T A271952 432,728,1120,1736,2832,4864,8576,20224,63808,162624,348224,668288,
%U A271952 1204736,2114560,6175744,13394432,34860544,104595968,304683008,807587840
%N A271952 Somos's sequence {b(6,n)} defined in comment in A078495: a(0)=a(1)=...=a(14)=1; for n>=15, a(n)=(a(n-1)*a(n-14)+a(n-7)*a(n-8))/a(n-15).
%H A271952 Seiichi Manyama, <a href="/A271952/b271952.txt">Table of n, a(n) for n = 0..440</a>
%t A271952 a[k_,n_]:=a[k,n]=If[n>2k+2,(a[k,(n-1)]*a[k,(n-2k-2)]+a[k,(n-k-1)]*a[k,(n-k-2)])/a[k,(n-2k-3)],1];
%t A271952 Map[a[6,#]&,Range[0,50]] (* _Peter J. C. Moses_, Apr 17 2016 *)
%o A271952 (PARI) {a(n) = if(n<= 15, 1, (a(n-1)*a(n-14) + a(n-7)*a(n-8))/a(n-15))}; for(n=1,50, print1(a(n), ", ")) \\ _G. C. Greubel_, Jul 30 2018
%o A271952 (Magma) [n le 15 select 1 else (Self(n-1)*Self(n-14)+Self(n-7)*Self(n-8))/Self(n-15): n in [1..60]]; // _G. C. Greubel_, Jul 30 2018
%Y A271952 Cf. A006721, A078495, A268199, A271948, A271949, A271950.
%K A271952 nonn
%O A271952 0,16
%A A271952 _Vladimir Shevelev_ and _Peter J. C. Moses_, Apr 17 2016