cp's OEIS Frontend

A255654 Partial sums of A247649.

%I A255654 #15 Dec 24 2021 08:17:09
%S A255654 1,6,11,18,23,40,47,66,71,96,113,132,139,170,189,214,219,244,269,304,
%T A255654 321,382,401,472,479,514,545,586,605,676,701,778,783,808,833,868,893,
%U A255654 978,1013,1108,1125,1210,1271,1342,1361,1452,1523,1600,1607,1642,1677,1726,1757,1864,1905,2026,2045,2140,2211,2296,2321,2434,2511,2614
%N A255654 Partial sums of A247649.
%H A255654 <a href="/index/Ce#cell">Index entries for sequences related to cellular automata</a>
%o A255654 (Python)
%o A255654 from itertools import islice, count
%o A255654 import sympy
%o A255654 def A225654gen(): # generator of terms
%o A255654     from sympy.abc import x
%o A255654     f, g, blist, c = 1/x**2+1/x+1+x+x**2, 1, [1], 1
%o A255654     yield c
%o A255654     for n in count(1):
%o A255654         s = [int(d,2) for d in bin(n)[2:].split('00') if d != '']
%o A255654         g = (g*f).expand(modulus=2)
%o A255654         if len(s) == 1:
%o A255654             blist.append(g.subs(x,1))
%o A255654         else:
%o A255654             blist.append(prod(blist[d] for d in s))
%o A255654         c += blist[-1]
%o A255654         yield c
%o A255654 A225654_list = list(islice(A225654gen(),26)) # _Chai Wah Wu_, Dec 23 2021
%Y A255654 Cf. A247649.
%K A255654 nonn
%O A255654 0,2
%A A255654 _Omar E. Pol_, Mar 01 2015