|
%I #23 Dec 31 2021 14:36:10
%S 2,3,5,2,6,3,7,4,8,5,9,6,10,7,11,8,12,9,13,10,14,11,15,12,16,13,17,14,
%T 18,15,19,16,20,17,21,18,22,19,23,20,24,21,25,22,26,23,27,24,28,25,29,
%U 26,30,27,31,28,32,29,33,30,34,31,35,32,36,33,37,34,38,35
%N a(1)=2, a(2)=3, a(3)=5 and a(n) = (a(1) + a(2) + a(3) + ... + a(n-1))/a(n-1).
%H Colin Barker, <a href="/A297996/b297996.txt">Table of n, a(n) for n = 1..1000</a>
%H <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (1,1,-1).
%F a(n) = A168230(n+1) for n >= 3.
%F From _Colin Barker_, Jan 29 2018: (Start)
%F G.f.: x*(2 + x - 4*x^3 + 2*x^4) / ((1 - x)^2*(1 + x)).
%F a(n) = n/2 for n>2 and even.
%F a(n) = (n+7)/2 for n>2 and odd.
%F a(n) = a(n-1) + a(n-2) - a(n-3) for n>5.
%F (End)
%t Nest[Append[#, Total[#]/Last[#]] &, Prime@ Range@ 3, 67] (* _Michael De Vlieger_, Jan 10 2018 *)
%t LinearRecurrence[{1,1,-1},{2,3,5,2,6},70] (* _Harvey P. Dale_, Dec 31 2021 *)
%o (PARI) lista(nn) = {va = vector(nn); for (n=1, 3, va[n] = prime(n)); for (n=4, nn, va[n] = sum(k=1, n-1, va[k])/va[n-1];); va;} \\ _Michel Marcus_, Jan 10 2018
%o (PARI) Vec(x*(2 + x - 4*x^3 + 2*x^4) / ((1 - x)^2*(1 + x)) + O(x^100)) \\ _Colin Barker_, Jan 29 2018
%Y Cf. A168230.
%K nonn,easy,less
%O 1,1
%A _Mateusz Pasternak_, Jan 10 2018
%E More terms from _Michel Marcus_, Jan 10 2018
|
