Back to Gallery
  • Remix
  • Share
  • New Entry
on 5 Oct 2021
  • 7
  • 20
  • 0
  • 0
  • 243
N = 2047;
f = zeros(N,1);
f(1) = 0;
b = @(x) dec2bin(x,ceil(log2(N)))-'0';
for i=2:N
bi = b(i);
no = [];
for j=1:i
if any(b(j).*bi>0)
no = [no f(j)];
end
end
if isempty(no)
f(i) = 0;
else
k = 0:max(no)+1;
l = k(find(~ismember(k,no),1,'first'));
f(i) = l;
end
end
plot(f,'.');

Image

Remix Tree