login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

a(n) = A257808(n) - A257807(n).
12

%I #45 Feb 24 2023 12:10:25

%S 0,-1,0,1,0,1,2,1,2,1,2,3,2,3,4,5,6,5,4,5,6,5,6,5,4,3,4,3,4,5,4,5,6,7,

%T 6,5,6,7,6,7,8,7,6,7,8,9,10,11,12,11,12,13,12,11,10,9,10,9,10,11,10,

%U 11,12,13,12,11,12,13,12,13,12,13,14,13,12,11,10,9,10,11,12,11,10,9,10,11,12,13,14,15,14,15,16,15,16,15,14

%N a(n) = A257808(n) - A257807(n).

%C Alternative description: Start with a(0) = 0, and then to obtain each a(n), look at each successive term in the infinite trunk of inverted binary beanstalk, from A233271(1) onward, subtracting one from a(n-1) if A233271(n) is odd, and adding one to a(n-1) if A233271(n) is even.

%C In other words, starting from zero, iterate the map x -> {x + 1 + number of nonleading zeros in the binary representation of x}, and note each time whether the result is odd or even: With odd results go one step down, and even results go one step up.

%C After the zeros at a(0), a(2) and a(4) and -1 at a(1), the terms stay strictly positive for a long time, although from the terms of A257805 it can be seen that the sequence must again fall to the negative side somewhere between n = 541110611 and n = 1051158027 (i.e., A218600(33) .. A218600(34)). Indeed the fourth zero occurs at n = 671605896, and the second negative term right after that as a(671605897) = -1.

%C The maximum positive value reached prior to the slide into negative territory is 2614822 for a(278998626) and a(278998628). - _Hans Havermann_, May 23 2015

%H Antti Karttunen, <a href="/A257806/b257806.txt">Table of n, a(n) for n = 0..8727</a>

%H Hans Havermann, <a href="http://chesswanks.com/num/a257806.png">Graph of 2*10^9 terms</a>

%H Hans Havermann, <a href="http://chesswanks.com/num/a257806zeros.png">Detail graph of the eventual crossover to negative terms and a listing of its associated zeros</a>

%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Blancmange_curve">Blancmange curve</a>

%F a(n) = A257808(n) - A257807(n).

%F a(0) = 0; and for n >= 1, a(n) = a(n-1) + (-1)^A233271(n).

%F Other identities. For all n >= 0:

%F a(A218600(n+1)) = -A257805(n).

%e We consider 0 to have no nonleading zeros, so first we get to 0 -> 0+1+0 = 1, and 1 is odd, so we go one step down from the starting value a(0)=0, and thus a(1) = -1.

%e 1 has no nonleading zeros, so we get 1 -> 1+1+0 = 2, and 2 is even, so we go one step up, and thus a(2) = 0.

%e 2 has one nonleading zero in binary "10", so we get 2 -> 2+1+1 = 4, and 4 is also even, so we go one step up, and thus a(3) = 1.

%e 4 has two nonleading zeros in binary "100", so we get 4 -> 4+2+1 = 7, 7 is odd, so we go one step down, and thus a(4) = 0.

%o (PARI)

%o A070939 = n->#binary(n)+!n; \\ From _M. F. Hasler_

%o A080791 = n->if(n<1,0,(A070939(n)-hammingweight(n)));

%o A233272 = n->(n + A080791(n) + 1);

%o A257806_write_bfile(up_to_n) = { my(n,a_n=0,b_n=0); for(n=0, up_to_n, write("b257806.txt", n, " ", a_n); b_n = A233272(b_n); a_n += ((-1)^b_n)); };

%o A257806_write_bfile(8727);

%o (Python)

%o def A257806_print_upto(n):

%o a = 0

%o b = 0

%o for n in range(n):

%o print(b, end=", ")

%o ta = a

%o c0 = 0

%o while ta>0:

%o c0 += 1-(ta&1)

%o ta >>= 1

%o a += 1 + c0

%o b += ((2*(1-(a&1))) - 1)

%o # By _Antti Karttunen_ after _Alex Ratushnyak_'s Python-code for A216431.

%o (Scheme, two alternatives, the latter using memoizing definec-macro)

%o (define (A257806 n) (- (A257808 n) (A257807 n)))

%o (definec (A257806 n) (if (zero? n) n (+ (expt -1 (A233271 n)) (A257806 (- n 1)))))

%Y Cf. A218600, A233271, A233272, A257807, A257808, A257803, A257804, A257805.

%Y Cf. also A218542, A218543, A218789 and A233270 (compare the scatter plots).

%K sign,base,look

%O 0,7

%A _Antti Karttunen_, May 12 2015