cp's OEIS Frontend

A062187 a(n+1) = a(n) - a(floor(n/2)), with a(0)=0, a(1)=1.

%I A062187 #11 Jun 09 2022 02:24:38
%S A062187 0,1,1,0,-1,-2,-3,-3,-3,-2,-1,1,3,6,9,12,15,18,21,23,25,26,27,26,25,
%T A062187 22,19,13,7,-2,-11,-23,-35,-50,-65,-83,-101,-122,-143,-166,-189,-214,
%U A062187 -239,-265,-291,-318,-345,-371,-397,-422,-447,-469,-491,-510,-529,-542,-555,-562,-569,-567,-565,-554,-543,-520,-497,-462
%N A062187 a(n+1) = a(n) - a(floor(n/2)), with a(0)=0, a(1)=1.
%C A062187 Period of oscillations above and below the axis more than doubles at each cycle.
%F A062187 G.f. A(x) satisfies: A(x) = x * (1 - (1 + x)*A(x^2))/(1 - x). - _Ilya Gutkovskiy_, May 04 2019
%e A062187 a(6) = a(5) - a(2) = -2 - 1 = -3.
%e A062187 a(7) = a(6) - a(3) = -3 - 0 = -3.
%o A062187 (Python)
%o A062187 from itertools import islice
%o A062187 from collections import deque
%o A062187 def A062187_gen(): # generator of terms
%o A062187     aqueue, f, b, a = deque([1]), True, 0, 1
%o A062187     yield from (0,1)
%o A062187     while True:
%o A062187         a -= b
%o A062187         yield a
%o A062187         aqueue.append(a)
%o A062187         if f: b = aqueue.popleft()
%o A062187         f = not f
%o A062187 A061287_list = list(islice(A062187_gen(),40)) # _Chai Wah Wu_, Jun 08 2022
%Y A062187 Cf. A033485, A062186, A062188.
%K A062187 sign
%O A062187 0,6
%A A062187 _Henry Bottomley_, Jun 13 2001