|
%I #5 Mar 07 2016 03:07:04
%S 8,32,33,34,35,36,64,65,128,129,130,131,132,133,134,135,136,137,138,
%T 139,140,141,142,143,144,145,146,147,148,149,150,152,160,161,256,257,
%U 258,259,260,261,262,263,264,265,266,267,268,269,270,272,273,274,275
%N Positions of negative terms in A059876.
%p negative_terms(A059876); negative_terms := proc(b) local a,i; a := []; for i from 1 to nops(b) do if(b[i] < 0) then a := [op(a),i]; fi; od; RETURN(a); end;
%t a[n_] := With[{s = Floor[Log[2, n]]}, (-1)^(n+1) + Sum[(-1)^(Floor[n/2^i] + 1)*Prime[i], {i, 1, s}] + If[1 == n, 1, Mod[s + 1, 2]*Prime[s]]];
%t A059880 = Position[Array[a, 500], _?Negative] // Flatten (* _Jean-François Alcover_, Mar 07 2016 *)
%K nonn
%O 1,1
%A _Antti Karttunen_, Feb 05 2001
|
