A065359(n) = my (b=Vecrev(binary(n))); -sum(k=1, #b, if (b[k], (-1)^k, 0))

g(z) = z/gcd(real(z), imag(z))

{
	z=0;
	for (n=0, 2^14,
		print (n " " imag(z));
		z+=g((1+I)^A065359(n));
	);
}

quit