login
Expansion of g.f.: (1+x^6+x^7)/(1-2*x+x^6-x^7-x^8).
1

%I #17 Aug 31 2023 13:39:26

%S 1,2,4,8,16,32,64,128,255,508,1012,2016,4016,8000,15937,31749,63249,

%T 126002,251016,500064,996207,1984602,3953641,7876278,15690791,

%U 31258536,62271945,124055559,247138286,492338537,980816202,1953940937,3892559256,7754593434,15448376086,30775607480,61309875581,122138964964

%N Expansion of g.f.: (1+x^6+x^7)/(1-2*x+x^6-x^7-x^8).

%C a(n) is the number of binary words of length n which have no 00010100-matches.

%H G. C. Greubel, <a href="/A234589/b234589.txt">Table of n, a(n) for n = 0..1000</a>

%H B. K. Miceli, J, Remmel, <a href="http://puma.dimai.unifi.it/23_3/miceli_remmel.pdf">Minimal Overlapping Embeddings and Exact Matches in Words</a>, PU. M. A., Vol. 23 (2012), No. 3, pp. 291-315.

%H <a href="/index/Rec#order_08">Index entries for linear recurrences with constant coefficients</a>, signature (2,0,0,0,0,-1,1,1).

%p seq(coeff(series((1+x^6+x^7)/(1-2*x+x^6-x^7-x^8), x, n+1), x, n), n = 0..40); # _G. C. Greubel_, Sep 13 2019

%t CoefficientList[Series[(1+x^6+x^7)/(1-2*x+x^6-x^7-x^8), {x,0,40}], x] (* _G. C. Greubel_, Sep 13 2019 *)

%t LinearRecurrence[{2,0,0,0,0,-1,1,1},{1,2,4,8,16,32,64,128},40] (* _Harvey P. Dale_, Aug 31 2023 *)

%o (PARI) my(x='x+O('x^40)); Vec((1+x^6+x^7)/(1-2*x+x^6-x^7-x^8)) \\ _G. C. Greubel_, Sep 13 2019

%o (Magma) R<x>:=PowerSeriesRing(Integers(), 40); Coefficients(R!( (1+x^6+x^7)/(1-2*x+x^6-x^7-x^8) )); // _G. C. Greubel_, Sep 13 2019

%o (Sage)

%o def A234589_list(prec):

%o P.<x> = PowerSeriesRing(ZZ, prec)

%o return P((1+x^6+x^7)/(1-2*x+x^6-x^7-x^8)).list()

%o A234589_list(40) # _G. C. Greubel_, Sep 13 2019

%o (GAP) a:=[1,2,4,8,16,32,64,128];; for n in [9..40] do a[n]:=2*a[n-1]-a[n-6]+a[n-7]+a[n-8]; od; a; # _G. C. Greubel_, Sep 13 2019

%Y Similar to but different from A172317.

%Y Cf. A234590, A234591, A234592.

%K nonn,easy

%O 0,2

%A _N. J. A. Sloane_, Jan 01 2014