OFFSET
0,2
PROG
(Python)
from itertools import islice, count
import sympy
def A225654gen(): # generator of terms
from sympy.abc import x
f, g, blist, c = 1/x**2+1/x+1+x+x**2, 1, [1], 1
yield c
for n in count(1):
s = [int(d, 2) for d in bin(n)[2:].split('00') if d != '']
g = (g*f).expand(modulus=2)
if len(s) == 1:
blist.append(g.subs(x, 1))
else:
blist.append(prod(blist[d] for d in s))
c += blist[-1]
yield c
A225654_list = list(islice(A225654gen(), 26)) # Chai Wah Wu, Dec 23 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Omar E. Pol, Mar 01 2015
STATUS
approved