cp's OEIS Frontend

A124072 First differences of A129819.

%I A124072 #21 Sep 17 2024 14:01:30
%S A124072 0,1,0,2,1,3,1,4,2,5,2,6,3,7,3,8,4,9,4,10,5,11,5,12,6,13,6,14,7,15,7,
%T A124072 16,8,17,8,18,9,19,9,20,10,21,10,22,11,23,11,24,12,25,12,26
%N A124072 First differences of A129819.
%C A124072 A129819 and its repeated differences are
%C A124072   0.0.1..1..3..4..7...8..12..14.19..21.27....
%C A124072   ..0.1..0..2..1..3...1...4...2..5...2..6....
%C A124072   ....1.-1..2.-1..2..-2...3..-2..3..-3..4....
%C A124072   ......-2..3.-3..3..-4...5..-5..5..-6..7....
%C A124072   ..........5.-6..6..-7...9.-10.10.-11.13...
%C A124072   ...........-11.12.-13..16.-19.20.-21.24.-27
%C A124072   ...............23.-25..29.-35.39.-41.45.-51
%C A124072 The left edge is A130668.
%C A124072 I discovered the array 1 1 -2 1 -3 2 in studying the singular points of planar polynomial differential systems (inspired by the reference).
%H A124072 G. C. Greubel, <a href="/A124072/b124072.txt">Table of n, a(n) for n = 0..5000</a>
%H A124072 Paul Curtz, <a href="http://www.numdam.org/item?id=ASENS_1980_4_13_3_293_0">Stabilite locale des systemes quadratiques</a>, Ann. sc. Ecole Norm. Sup. vol 13 no 3 (1980) pp 293-302.
%H A124072 <a href="/index/Rec#order_06">Index entries for linear recurrences with constant coefficients</a>, signature (0,1,0,1,0,-1).
%F A124072 a(2n) = A004526(n).
%F A124072 a(2n+1) = A000027(n+1) .
%F A124072 G.f.: x*(1+x^2+x^3)/((1-x)^2*(1+x)^2*(1+x^2)). - _R. J. Mathar_, Feb 25 2009
%F A124072 From _G. C. Greubel_, Sep 17 2024: (Start)
%F A124072 a(n) = (1/8)*(3*n + 1 - (-1)^n*(n + 3) + i^n*(1 + (-1)^n)).
%F A124072 E.g.f.: (1/4)*( cos(x) - (1-2*x)*cosh(x) + (2+x)*sinh(x) ). (End)
%t A124072 a[n_?OddQ] := (n+1)/2; a[n_?EvenQ] := Floor[n^2/16] - Floor[(n-2)^2/16]; Table[a[n], {n, 0, 51}] (* _Jean-François Alcover_, Aug 13 2012 *)
%t A124072 LinearRecurrence[{0,1,0,1,0,-1}, {0,1,0,2,1,3}, 61] (* _G. C. Greubel_, Sep 17 2024 *)
%o A124072 (Magma)
%o A124072 R<x>:=PowerSeriesRing(Integers(), 61);
%o A124072 [0] cat Coefficients(R!( x*(1+x^2+x^3)/((1-x^2)*(1-x^4)) )); // _G. C. Greubel_, Sep 17 2024
%o A124072 (SageMath)
%o A124072 def A124072_list(prec):
%o A124072     P.<x> = PowerSeriesRing(ZZ, prec)
%o A124072     return P( x*(1+x^2+x^3)/((1-x^2)*(1-x^4)) ).list()
%o A124072 A124072_list(60) # _G. C. Greubel_, Sep 17 2024
%Y A124072 Cf. A000027, A004526, A129819, A130668.
%K A124072 nonn,easy
%O A124072 0,4
%A A124072 _Paul Curtz_, Jun 26 2007
%E A124072 Partially edited by _R. J. Mathar_, Jul 07 2008