login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 


Boustrophedon transform of A068717 with A068717(1) = -1 instead of 0.
1

%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

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified September 22 14:04 EDT 2024. Contains 376114 sequences. (Running on oeis4.)