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

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A048718 Binary expansion matches ((0)*0001)*(0*); or, Zeckendorf-like expansion of n using recurrence f(n) = f(n-1) + f(n-4). 8

%I #25 Dec 31 2022 15:20:16

%S 0,1,2,4,8,16,17,32,33,34,64,65,66,68,128,129,130,132,136,256,257,258,

%T 260,264,272,273,512,513,514,516,520,528,529,544,545,546,1024,1025,

%U 1026,1028,1032,1040,1041,1056,1057

%N Binary expansion matches ((0)*0001)*(0*); or, Zeckendorf-like expansion of n using recurrence f(n) = f(n-1) + f(n-4).

%C Max. 1 one-bit occur in each range of four bits.

%C Constructed from A003269 in the same way as A003714 is constructed from A000045.

%H Sebastian Karlsson, <a href="/A048718/a048718.txt">Walnut code that verifies the conjectures of Paul D. Hanna</a>

%H Walnut can be downloaded from <a href="https://cs.uwaterloo.ca/~shallit/walnut.html">https://cs.uwaterloo.ca/~shallit/walnut.html</a>.

%H <a href="/index/Con#CongruCrossDomain">Index entries for sequences defined by congruent products between domains N and GF(2)[X]</a>

%H <a href="/index/Con#CongruXOR">Index entries for sequences defined by congruent products under XOR</a>

%F a(0) = 0, a(n) = (2^(invfyy(n)-1))+a(n-fyy(invfyy(n))) where fyy(n) is fyy(n-1) + fyy(n-4) (A003269) and invfyy is its "integral" (floored down) inverse.

%F a(n) XOR 14*a(n) = 15*a(n); 3*a(n) XOR 9*a(n) = 10*a(n); 3*a(n) XOR 13*a(n) = 14*a(n); 5*a(n) XOR 9*a(n) = 12*a(n); 5*a(n) XOR 11*a(n) = 14*a(n); 6*a(n) XOR 11*a(n) = 13*a(n); 7*a(n) XOR 9*a(n) = 14*a(n); 7*a(n) XOR 10*a(n) = 13*a(n); 7*a(n) XOR 11*a(n) = 12*a(n); 12*a(n) XOR 21*a(n) = 25*a(n); 12*a(n) XOR 37*a(n) = 41*a(n); etc. (conjectures). - _Paul D. Hanna_, Jan 22 2006

%F The conjectures can be verified using the Walnut theorem-prover (see links). - _Sebastian Karlsson_, Dec 31 2022

%t filterQ[n_] := With[{bb = IntegerDigits[n, 2]}, !MemberQ[{{1, 1}, {1, 0, 1}, {1, 1, 0}, {1, 1, 1}}, bb] && SequencePosition[bb, {a_, b_, c_, d_} /; Count[{a, b, c, d}, 1] > 1] == {}];

%t Select[Range[0, 1057], filterQ] (* _Jean-François Alcover_, Dec 31 2020 *)

%o (PARI) is(n)=!bitand(n, 14*n) \\ _Charles R Greathouse IV_, Oct 03 2016

%Y Cf. A048715, A048719, A115422, A115423, A115424.

%K nonn,base,easy

%O 0,3

%A _Antti Karttunen_, Mar 30 1999

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | 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 April 24 20:08 EDT 2024. Contains 371963 sequences. (Running on oeis4.)