|
%I #21 Jan 02 2023 12:30:54
%S 1,3,2,3,7,4,5,11,6,7,15,8,9,19,10,11,23,12,13,27,14,15,31,16,17,35,
%T 18,19,39,20,21,43,22,23,47,24,25,51,26,27,55,28,29,59,30,31,63,32,33,
%U 67,34,35,71,36,37,75,38,39,79,40
%N Partial sums (and absolute value of first differences) of A332056: if odd (resp. even) add (resp. subtract) the partial sum to get the next term.
%C The terms show a 3-quasiperiodic pattern (2m-1, 4m-1, 2m), m = 1, 2, 3, ...
%C Or: group positive integers by pairs, then insert the sum of the pair between the two terms.
%H Antti Karttunen, <a href="/A332057/b332057.txt">Table of n, a(n) for n = 1..19683</a>
%H Eric Angelini, <a href="http://list.seqfan.eu/oldermail/seqfan/2020-February/020509.html">Re: Add or subtract my cumulative sum of terms</a>, SeqFan list, Feb 24 2020.
%H <a href="/index/Rec#order_06">Index entries for linear recurrences with constant coefficients</a>, signature (0,0,2,0,0,-1).
%F a(3k-2) = 2k - 1, a(3k-1) = 4k - 1, a(3k) = 2k, for all k >= 1.
%F From _Colin Barker_, Feb 25 2020: (Start)
%F G.f.: x*(1 + x)*(1 + 2*x + x^3) / ((1 - x)^2*(1 + x + x^2)^2).
%F a(n) = 2*a(n-3) - a(n-6) for n>6.
%F (End)
%o (PARI) apply( {A332057(n)=n<<max(n%3,1)\/3}, [1..99])
%o (PARI) Vec(x*(1 + x)*(1 + 2*x + x^3) / ((1 - x)^2*(1 + x + x^2)^2) + O(x^60)) \\ _Colin Barker_, Feb 26 2020
%Y Cf. A332056.
%K nonn,easy
%O 1,2
%A _Eric Angelini_ and _M. F. Hasler_, Feb 24 2020
|
