load sp01;
l=size(sp01,1);

rnd01=sparse(zeros(size(sp01)));
for i=1:size(sp01,2)
  d=sum(sp01(:,i))/l;
  r=rand(l,1);
  rnd01(:,i)=(r<=d);
end

save rnd01_1 rnd01