A108620 - OEIS

%I #26 Mar 15 2024 15:20:10

%S 2,4,4,2,0,-3,-6,-4,2,6,6,2,-5,-10,-6,4,10,8,0,-9,-12,-4,11,20,12,-6,

%T -16,-10,6,19,17,1,-16,-18,-2,16,21,9,-12,-22,-10,12,25,17,-5,-22,-18,

%U 4,22,20,0,-21,-24,-4,23,32,12,-21,-36,-16,23,44,24,-18,-40,-22,18,43,29,-11,-40,-30,10,40,32,-6,-39,-35,3,40,40,2,-36

%N 2*a(n) = A108618(n) + 3*A108619(n).

%C Floretion Algebra Multiplication Program, FAMP Code: 1vessum(*)seq[ + .5'i + .5'j + .5'k + .5e]

%H Creighton Dement, <a href="/A108620/b108620.txt">Table of n, a(n) for n = 0..10000</a>

%H Rémy Sigrist, <a href="/A108620/a108620.png">Colored scatterplot of a(n) for n = 0..10000</a> (where the color is function of n mod 6)

%t a[0] = b[0] = 1;

%t f[n_] := Sign[n]*Mod[n, 2];

%t a[n_] := a[n] = (1/2)*(a[n-1] - 3*b[n-1]) + 3*f[(1/2)*(a[n-1] + b[n-1])] + f[(1/2)*(a[n-1] - 3*b[n-1])] + 1;

%t b[n_] := b[n] = (1/2)*(a[n-1] + b[n-1]) + 1;

%t A108620 = Table[(a[n] + 3*b[n])/2, {n, 0, 100}] (* _Jean-François Alcover_, Feb 25 2015, after _Benoit Jubin_ *)

%Y Cf. A108618, A108619.

%K easy,sign,look

%O 0,1

%A _Creighton Dement_, Jun 22 2005