function s=bin(n,k)
for ii=k-1:-1:0,
    s(k-ii)=floor(n/2^ii);
    if floor(n/2^ii)~=0,
        n=mod(n,2^ii);
    end
end