Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #9 Jun 15 2022 01:50:18
%S 1,2,2,2,8,25,90,404,2055,11792,75053,525622,4015361,33231679,
%T 296182315,2828335731,28809181418,311788442591,3572832236720,
%U 43216177809190,550245463265240,7356239983352887,103028812072639378
%N Boustrophedon transform of A068717 with A068717(1) = -1 instead of 0.
%C Signs omitted, all terms < 0. A068717(1) = -1 means soluble (1*1 -1*0*0 = +1 resp. 0*0 -1*1*1 = -1).
%o (Python)
%o from itertools import accumulate, count, islice
%o from math import isqrt
%o from sympy import continued_fraction_periodic
%o def A068718_gen(): # generator of terms
%o yield 1
%o blist = (1,)
%o for n in count(2):
%o yield (blist := tuple(accumulate(reversed(blist),initial=0 if (a:=isqrt(n)**2==n) else (1 if len(continued_fraction_periodic(0,1,n)[1]) & 1 else int(a)-1))))[-1]
%o A068718_list = list(islice(A068718_gen(),40)) # _Chai Wah Wu_, Jun 14 2022
%Y Cf. A068716, A068717, A049240, A046951.
%K nonn,easy
%O 0,2
%A _Frank Ellermann_, Feb 26 2002