function s = BR_get_samp_P(P)
%

two_pow_n = length(P);
n = log2(two_pow_n);
r = rand;
p = 0;
for k = 1:two_pow_n
    p = p + P(k);
    if r<p
        s = BR_int2bits(k-1,n);
        return
    end
end
s = ones(1,n);
